Beispiel #1
0
// Compute the maximum wave speed
Real PolytropicPhysics::getMaxWaveSpeed(const FArrayBox& a_U,
                                        const Box&       a_box)
{
    CH_assert(isDefined());
    CH_assert(a_U.contains(a_box));

    Real speed = 0.0;

    FORT_MAXWAVESPEEDF(CHF_REAL(speed),
                       CHF_CONST_FRA(a_U),
                       CHF_BOX(a_box));

    return speed;
}
Beispiel #2
0
/// Compute the maximum wave speed
Real LinElastPhysics::getMaxWaveSpeed(const FArrayBox& a_U,
    const Box&       a_box)
{
        CH_assert(isDefined());
    CH_assert(a_U.contains(a_box));

    Real speed = 0.0;

    FORT_MAXWAVESPEEDF(CHF_REAL(speed),
                       CHF_CONST_FRA(a_U),
                       CHF_BOX(a_box));

    //JK pout() << "LinElastPhysics::getMaxWaveSpeed :: SPEED = " << speed << endl;

    return speed;
}