CD_ItoParticle.H

Go to the documentation of this file.

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
#ifndef CD_ITOPARTICLE_H
#define CD_ITOPARTICLE_H
 
// Our includes
#include <CD_ParticleSoA.H>
#include <CD_NamespaceHeader.H>
 
struct ItoParticle
{
  ParticleReal mobility  = 0.0; 
  ParticleReal diffusion = 0.0; 
  ParticleReal energy    = 0.0; 
  ParticleReal scratch   = 0.0; 
 
  double old_x = 0.0; 
  double old_y = 0.0; 
#if CH_SPACEDIM == 3
  double old_z = 0.0; 
#endif
 
  ParticleReal vx = 0.0; 
  ParticleReal vy = 0.0; 
#if CH_SPACEDIM == 3
  ParticleReal vz = 0.0; 
#endif
 
  ParticleReal scratch_x = 0.0; 
  ParticleReal scratch_y = 0.0; 
#if CH_SPACEDIM == 3
  ParticleReal scratch_z = 0.0; 
#endif
};
 
template <>
struct ParticleTraits<ItoParticle>
{
  static constexpr auto columns = std::make_tuple(
    &ItoParticle::mobility,
    &ItoParticle::diffusion,
    &ItoParticle::energy,
    &ItoParticle::scratch,
    D_DECL(&ItoParticle::old_x, &ItoParticle::old_y, &ItoParticle::old_z),
    D_DECL(&ItoParticle::vx, &ItoParticle::vy, &ItoParticle::vz),
    D_DECL(&ItoParticle::scratch_x, &ItoParticle::scratch_y, &ItoParticle::scratch_z));
 
  static constexpr auto h5Columns = std::make_tuple(D_DECL(&ItoParticle::vx, &ItoParticle::vy, &ItoParticle::vz));
 
  static constexpr auto h5PartColumns = std::make_tuple(
    H5Part::Scalar{"mobility", &ItoParticle::mobility},
    H5Part::Scalar{"diffusion", &ItoParticle::diffusion},
    H5Part::Scalar{"energy", &ItoParticle::energy},
    H5Part::Vector{"v", D_DECL(&ItoParticle::vx, &ItoParticle::vy, &ItoParticle::vz)});
};
 
#include <CD_NamespaceFooter.H>
 
#endif