// Compute increment of primitive variables. void PolytropicPhysics::quasilinearUpdate(FArrayBox& a_dWdx, const FArrayBox& a_WHalf, const FArrayBox& a_W, const Real& a_scale, const int& a_dir, const Box& a_box) { CH_assert(isDefined()); CH_assert(a_dWdx.box().contains(a_box)); FORT_GETADWDXF(CHF_FRA(a_dWdx), CHF_CONST_FRA(a_WHalf), CHF_CONST_FRA(a_W), CHF_CONST_REAL(a_scale), CHF_CONST_INT(a_dir), CHF_BOX(a_box)); }
void LinElastPhysics::quasilinearUpdate(FArrayBox& a_dWdx, const FArrayBox& a_wHalf, const FArrayBox& a_W, const Real& a_scale, const int& a_dir, const Box& a_box) { //JK pout() << "LinElastPhysics::quasilinear" << endl; CH_assert(isDefined()); CH_assert(a_dWdx.box().contains(a_box)); FORT_GETADWDXF(CHF_FRA(a_dWdx), CHF_CONST_FRA(a_wHalf), CHF_CONST_FRA(a_W), CHF_CONST_REAL(a_scale), CHF_CONST_INT(a_dir), CHF_BOX(a_box)); }