walker::MixDirichletHydroTimeScale class

MixDirichlet coefficients policity: mean(rho) forced const in time

User-defined parameters b' and kappa' are functions of an externally, e.g., DNS-, provided hydrodynamics time scale ensuring decay in the evolution of <y_alpha^2>. Additionally, S_alpha is constrained to make d<rho>/dt = 0. Additionally, we pull in a hydrodynamic timescale from an external function.

Constructors, destructors, conversion operators

MixDirichletHydroTimeScale(tk::ctr::ncomp_t ncomp, ctr::NormalizationType norm, const std::vector<kw::sde_b::info::expect::type>& b_, const std::vector<kw::sde_S::info::expect::type>& S_, const std::vector<kw::sde_kappa::info::expect::type>& kprime_, const std::vector<kw::sde_rho::info::expect::type>& rho_, std::vector<kw::sde_b::info::expect::type>& b, std::vector<kw::sde_kappa::info::expect::type>& kprime, std::vector<kw::sde_S::info::expect::type>& S, std::vector<kw::sde_rho::info::expect::type>& rho, std::vector<kw::sde_r::info::expect::type>& r, std::vector<kw::sde_kappa::info::expect::type>& k)
Constructor: initialize coefficients.

Public functions

void update(char depvar, ncomp_t ncomp, ctr::NormalizationType norm, std::size_t density_offset, std::size_t volume_offset, const std::map<tk::ctr::Product, tk::real>& moments, const std::vector<kw::sde_rho::info::expect::type>& rho, const std::vector<kw::sde_r::info::expect::type>& r, const std::vector<kw::sde_kappa::info::expect::type>& kprime, const std::vector<kw::sde_b::info::expect::type>& b, std::vector<kw::sde_kappa::info::expect::type>& k, std::vector<kw::sde_kappa::info::expect::type>& S) const
Update coefficients.

Function documentation

walker::MixDirichletHydroTimeScale::MixDirichletHydroTimeScale(tk::ctr::ncomp_t ncomp, ctr::NormalizationType norm, const std::vector<kw::sde_b::info::expect::type>& b_, const std::vector<kw::sde_S::info::expect::type>& S_, const std::vector<kw::sde_kappa::info::expect::type>& kprime_, const std::vector<kw::sde_rho::info::expect::type>& rho_, std::vector<kw::sde_b::info::expect::type>& b, std::vector<kw::sde_kappa::info::expect::type>& kprime, std::vector<kw::sde_S::info::expect::type>& S, std::vector<kw::sde_rho::info::expect::type>& rho, std::vector<kw::sde_r::info::expect::type>& r, std::vector<kw::sde_kappa::info::expect::type>& k)

Constructor: initialize coefficients.

Parameters
ncomp in Number of scalar components in this SDE system
norm in Normalization type (N=heavy or N=light)
b_ in Vector used to initialize coefficient vector b
S_ in Vector used to initialize coefficient vector S
kprime_ in Vector used to initialize coefficient vector kprime and k
rho_ in Vector used to initialize coefficient vector rho and r
b in/out Coefficient vector to be initialized
kprime in/out Coefficient vector to be initialized
S in/out Coefficient vector to be initialized
rho in/out Coefficient vector to be initialized
r in/out Coefficient vector to be initialized
k in/out Coefficient vector to be initialized

void walker::MixDirichletHydroTimeScale::update(char depvar, ncomp_t ncomp, ctr::NormalizationType norm, std::size_t density_offset, std::size_t volume_offset, const std::map<tk::ctr::Product, tk::real>& moments, const std::vector<kw::sde_rho::info::expect::type>& rho, const std::vector<kw::sde_r::info::expect::type>& r, const std::vector<kw::sde_kappa::info::expect::type>& kprime, const std::vector<kw::sde_b::info::expect::type>& b, std::vector<kw::sde_kappa::info::expect::type>& k, std::vector<kw::sde_kappa::info::expect::type>& S) const

Update coefficients.

Parameters
depvar in Dependent variable
ncomp in Number of scalar components in this SDE system
norm in Normalization type (N=heavy or N=light)
density_offset in Offset of particle density in solution array relative to YN
volume_offset in Offset of particle specific volume in solution array relative to YN
moments in Map of statistical moments estimated
rho in Coefficient vector
r in Coefficient Vector
kprime in Coefficient vector
b in Coefficient vector
k in/out Coefficient vector to be updated
S in/out Coefficient vector to be updated