walker::MixDirichletHomCoeffConst class

MixDirichlet coefficients policity: constants in time + mean(rho) = const

User-defined parameters b and kappaprime are constant vectors in time and, S is constrained to make $\mathrm{d}<rho>/\mathrm{d}t = 0$ .

Constructors, destructors, conversion operators

MixDirichletHomCoeffConst(tk::ctr::ncomp_type ncomp, 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, 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, std::vector<kw::sde_kappa::info::expect::type>& k, std::vector<kw::sde_kappa::info::expect::type>& S) const
Update coefficients vector S so avg(rho) = const.