// 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));
}
