// Compute expansion amplitudes of dW in right eigenvectors. void PolytropicPhysics::charAnalysis(FArrayBox& a_dW, const FArrayBox& a_W, const int& a_dir, const Box& a_box) { CH_assert(isDefined()); FORT_CHARANALYSISF(CHF_FRA(a_dW), CHF_CONST_FRA(a_W), CHF_CONST_INT(a_dir), CHF_BOX(a_box)); }
/** On input, a_dW contains the increments of the primitive variables. On output, it contains the increments in the characteristic variables. IMPORTANT NOTE: It is assumed that the characteristic analysis puts the smallest eigenvalue first, the largest eigenvalue last, and orders the characteristic variables accordingly. */ void LinElastPhysics::charAnalysis(FArrayBox& a_dW, const FArrayBox& a_W, const int& a_dir, const Box& a_box) { //JK pout() << "LinElastPhysics::charAnalysis" << endl; CH_assert(isDefined()); /* * This routine assumes that cp >= cs */ FORT_CHARANALYSISF(CHF_FRA(a_dW), CHF_CONST_FRA(a_W), CHF_CONST_INT(a_dir), CHF_BOX(a_box)); }