fields.h

Go to the documentation of this file.

/*
 *
 *    Copyright (c) 2021-2022
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 *
 */
#ifndef SRC_INCLUDE_SMASH_FIELDS_H_
#define SRC_INCLUDE_SMASH_FIELDS_H_
 
#include <iostream>
#include <typeinfo>
#include <utility>
#include <vector>
 
#include "density.h"
#include "experimentparameters.h"
#include "fourvector.h"
#include "potentials.h"
#include "threevector.h"
 
namespace smash {
 
class FieldsOnLattice {
 public:
  FieldsOnLattice()
      : A_mu_(FourVector()),
        dAmu_dxnu_({FourVector(), FourVector(), FourVector(), FourVector()}) {}
 
  FourVector A_mu() const { return A_mu_; }
 
  std::array<FourVector, 4> dAmu_dxnu() const { return dAmu_dxnu_; }
 
  ThreeVector dvecA_dt() { return dAmu_dxnu_[0].threevec(); }
 
  ThreeVector grad_A0() {
    ThreeVector A_0_grad = ThreeVector();
    for (int i = 1; i < 4; i++) {
      A_0_grad[i - 1] = dAmu_dxnu_[i].x0();
    }
    return A_0_grad;
  }
 
  ThreeVector curl_vecA() {
    ThreeVector curl_vec_A = ThreeVector();
    curl_vec_A.set_x1(dAmu_dxnu_[2].x3() - dAmu_dxnu_[3].x2());
    curl_vec_A.set_x2(dAmu_dxnu_[3].x1() - dAmu_dxnu_[1].x3());
    curl_vec_A.set_x3(dAmu_dxnu_[1].x2() - dAmu_dxnu_[2].x1());
    return curl_vec_A;
  }
 
  void overwrite_A_mu(FourVector new_A_mu) { A_mu_ = new_A_mu; }
 
  void overwrite_dAmu_dt_to_zero() {
    FourVector tmp0(0.0, 0.0, 0.0, 0.0);
 
    dAmu_dxnu_[0] = tmp0;
  }
 
  void overwrite_dAmu_dxnu(FourVector dAmu_dt, FourVector dAmu_dx,
                           FourVector dAmu_dy, FourVector dAmu_dz) {
    dAmu_dxnu_[0] = dAmu_dt;
    dAmu_dxnu_[1] = dAmu_dx;
    dAmu_dxnu_[2] = dAmu_dy;
    dAmu_dxnu_[3] = dAmu_dz;
  }
 
 private:
  FourVector A_mu_;
  std::array<FourVector, 4> dAmu_dxnu_;
};
 
typedef RectangularLattice<FieldsOnLattice> FieldsLattice;
 
void update_fields_lattice(
    RectangularLattice<FieldsOnLattice> *fields_lat,
    RectangularLattice<FourVector> *old_fields,
    RectangularLattice<FourVector> *new_fields,
    RectangularLattice<std::array<FourVector, 4>> *fields_four_grad_lattice,
    DensityLattice *jmu_B_lat, const LatticeUpdate fields_lat_update,
    const Potentials &potentials, const double time_step);
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_FIELDS_H_