particles_data.h

#ifndef COSMO_PARTICLES_DATA
#define COSMO_PARTICLES_DATA

namespace cosmo
{

template<typename RT>
struct Particle {
  RT X[3];     // particle position
  RT U[3];     // Particle velocity (covariant / lowered)
  RT M;        // Particle Mass
};

template<typename RT>
struct ParticleRegister {
  // 4 standard "registers" for each particle
  Particle<RT> p_p;
  Particle<RT> p_a;
  Particle<RT> p_c;
  Particle<RT> p_f;
};

typedef std::vector<ParticleRegister<real_t>> particle_vec;

template<typename RT>
struct ParticleMetricPrimitives {
  // metric components
  RT rootdetg;      // rooted metric determinant
  RT alpha;         // lapse
  RT beta[3];       // shift
  RT dalpha[3];     // derivative of lapse
  RT dbeta[3][3];   // derivative of shift
  RT gi[6];         // inverse spatial metric
  RT dgi[3][6];     // Inverse metric derivative
};

}

#endif