Beispiel #1
0
void Problem::checkGCEResult(const LifeV::Real& time)
{
    LifeV::Real dispNorm=M_fsi->displacement().norm2();
	if (time==0.000 && (dispNorm-110316)/dispNorm*(dispNorm-110316)/dispNorm<1e-5) Problem::resultCorrect(time);
	else if (time==0.001 && (dispNorm-99468.8)/dispNorm*(dispNorm-99468.8)/dispNorm<1e-5)  Problem::resultCorrect(time);
	else if (time==0.002 && (dispNorm-90957)/dispNorm*(dispNorm-90957)/dispNorm<1e-5)  Problem::resultCorrect(time);
	else if (time==0.003 && (dispNorm-90070.5)/dispNorm*(dispNorm-90070.5)/dispNorm<1e-5)  Problem::resultCorrect(time);
	else if (time==0.004 && (dispNorm-88162.4)/dispNorm*(dispNorm-88162.4)/dispNorm<1e-5)  Problem::resultCorrect(time);
}
Beispiel #2
0
void Problem::checkCEResult(const LifeV::Real& time)
{
	LifeV::Real dispNorm=M_fsi->displacement().norm2();
    if (time==0.000 && (dispNorm-106344)/dispNorm*(dispNorm-106344)/dispNorm<1e-3) Problem::resultCorrect(time);
    else if (time==0.001 && (dispNorm-147017)/dispNorm*(dispNorm-147017)/dispNorm<1e-3) Problem::resultCorrect(time);
    else if (time==0.002 && (dispNorm-108341)/dispNorm*(dispNorm-108341)/dispNorm<1e-3) Problem::resultCorrect(time);
    else if (time==0.003 && (dispNorm-106092)/dispNorm*(dispNorm-106092)/dispNorm<1e-3) Problem::resultCorrect(time);
    else if (time==0.004 && (dispNorm-105614)/dispNorm*(dispNorm-105614)/dispNorm<1e-3) Problem::resultCorrect(time);

}
Beispiel #3
0
void Problem::checkResultGCE (const LifeV::Real& time)
{

    //Extract the previous solution
    LifeV::Real dispNorm = M_fsi->displacement().norm2();
    if (time == 0.006 &&      (dispNorm - 89122.6) / dispNorm * (dispNorm - 89122.6) / dispNorm < 1e-5)
    {
        resultCorrect (time);
    }
    else if (time == 0.007 && (dispNorm - 83260.4) / dispNorm * (dispNorm - 83260.4) / dispNorm < 1e-5)
    {
        resultCorrect (time);
    }
    else if (time == 0.008 && (dispNorm - 79342.5) / dispNorm * (dispNorm - 79342.5) / dispNorm < 1e-5)
    {
        resultCorrect (time);
    }
}
Beispiel #4
0
void Problem::checkResultGI (const LifeV::Real& time)
{

    //Extract the previous solution
    LifeV::Real dispNorm = M_fsi->displacement().norm2();
    if (time == 0.006 &&      ( (dispNorm - 104264) / dispNorm * (dispNorm - 104264) / dispNorm < 1e-5 ) )
    {
        resultCorrect (time);
    }
    else if (time == 0.007 && ( (dispNorm - 101014) / dispNorm * (dispNorm - 101014) / dispNorm < 1e-5 ) )
    {
        resultCorrect (time);
    }
    else if (time == 0.008 && ( (dispNorm - 98183.7) / dispNorm * (dispNorm - 98183.7) / dispNorm < 1e-5 ) )
    {
        resultCorrect (time);
    }
}