AtmosphericData Atmosphere::GetAtmosphericData(F32 altitude)
{
AtmosphericData data;
data.temperature = GetTemperature(altitude);
data.soundSpeed = GetSoundSpeed(altitude);
return data;
}
void XFunctions::StoreConfigurationValues()
{
wxString path;
path.Printf(wxT("/Aurora/XFunctions/FftSize"));
m_pAuroraCfg->Write(path, GetFFTLength());
path.Printf(wxT("/Aurora/XFunctions/WindowType"));
m_pAuroraCfg->Write(path, AFSampleCount(GetWindowType()));
path.Printf(wxT("/Aurora/XFunctions/XFunctionType"));
m_pAuroraCfg->Write(path, AFSampleCount(GetXFunctionType()));
path.Printf(wxT("/Aurora/XFunctions/ProbeFreeField"));
m_pAuroraCfg->Write(path, IsProbeFreeField());
path.Printf(wxT("/Aurora/XFunctions/ProbeRigidTerminated"));
m_pAuroraCfg->Write(path, IsProbeRigidTerminated());
path.Printf(wxT("/Aurora/XFunctions/SoundSpeed"));
m_pAuroraCfg->Write(path, GetSoundSpeed());
path.Printf(wxT("/Aurora/XFunctions/ProbeMicDistance"));
m_pAuroraCfg->Write(path, GetProbeDistance());
path.Printf(wxT("/Aurora/XFunctions/ProbeMaxFreq"));
m_pAuroraCfg->Write(path, GetProbeMaxFreq());
path.Printf(wxT("/Aurora/XFunctions/FollowingFilter"));
m_pAuroraCfg->Write(path, IsFollowingFilterEnabled());
path.Printf(wxT("/Aurora/XFunctions/FollowingBandwidth"));
m_pAuroraCfg->Write(path, GetFollowingBandwidth());
path.Printf(wxT("/Aurora/XFunctions/Normalize"));
m_pAuroraCfg->Write(path, IsNormalizeSet());
path.Printf(wxT("/Aurora/XFunctions/ShiftToHalfWindow"));
m_pAuroraCfg->Write(path, IsShiftToHalfWindowSet());
path.Printf(wxT("/Aurora/XFunctions/CoherenceWeighting"));
m_pAuroraCfg->Write(path, IsCoherenceWeightingSet());
path.Printf(wxT("/Aurora/XFunctions/HilbertTransform"));
m_pAuroraCfg->Write(path, IsHilbertTransformSet());
path.Printf(wxT("/Aurora/XFunctions/TimeReversal"));
m_pAuroraCfg->Write(path, IsTimeReversalSet());
path.Printf(wxT("/Aurora/XFunctions/SaveMultispectrum"));
m_pAuroraCfg->Write(path, IsSaveMultispectrumSet());
path.Printf(wxT("/Aurora/XFunctions/DiracPulse"));
m_pAuroraCfg->Write(path, IsDiracPulseSet());
path.Printf(wxT("/Aurora/XFunctions/TriggerLevel"));
m_pAuroraCfg->Write(path, GetTriggerLevel());
// Writes immediately changes
m_pAuroraCfg->Flush();
}
void EulerADCFE::DoOdeRhs(
const Array<OneD, const Array<OneD, NekDouble> > &inarray,
Array<OneD, Array<OneD, NekDouble> > &outarray,
const NekDouble time)
{
int i;
int nvariables = inarray.num_elements();
int npoints = GetNpoints();
Array<OneD, Array<OneD, NekDouble> > advVel;
Array<OneD, Array<OneD, NekDouble> > outarrayAdv(nvariables);
Array<OneD, Array<OneD, NekDouble> > outarrayDiff(nvariables);
for (i = 0; i < nvariables; ++i)
{
outarrayAdv[i] = Array<OneD, NekDouble>(npoints, 0.0);
outarrayDiff[i] = Array<OneD, NekDouble>(npoints, 0.0);
}
m_advection->Advect(nvariables, m_fields, advVel, inarray,
outarrayAdv, m_time);
for (i = 0; i < nvariables; ++i)
{
Vmath::Neg(npoints, outarrayAdv[i], 1);
}
m_diffusion->Diffuse(nvariables, m_fields, inarray, outarrayDiff);
if (m_shockCaptureType == "NonSmooth")
{
for (i = 0; i < nvariables; ++i)
{
Vmath::Vadd(npoints,
outarrayAdv[i], 1,
outarrayDiff[i], 1,
outarray[i], 1);
}
}
if(m_shockCaptureType == "Smooth")
{
const Array<OneD, int> ExpOrder = GetNumExpModesPerExp();
NekDouble pOrder = Vmath::Vmax(ExpOrder.num_elements(), ExpOrder, 1);
Array <OneD, NekDouble > a_vel (npoints, 0.0);
Array <OneD, NekDouble > u_abs (npoints, 0.0);
Array <OneD, NekDouble > pres (npoints, 0.0);
Array <OneD, NekDouble > wave_sp(npoints, 0.0);
GetPressure(inarray, pres);
GetSoundSpeed(inarray, pres, a_vel);
GetAbsoluteVelocity(inarray, u_abs);
Vmath::Vadd(npoints, a_vel, 1, u_abs, 1, wave_sp, 1);
NekDouble max_wave_sp = Vmath::Vmax(npoints, wave_sp, 1);
Vmath::Smul(npoints,
m_C2,
outarrayDiff[nvariables-1], 1,
outarrayDiff[nvariables-1], 1);
Vmath::Smul(npoints,
max_wave_sp,
outarrayDiff[nvariables-1], 1,
outarrayDiff[nvariables-1], 1);
Vmath::Smul(npoints,
pOrder,
outarrayDiff[nvariables-1], 1,
outarrayDiff[nvariables-1], 1);
for (i = 0; i < nvariables; ++i)
{
Vmath::Vadd(npoints,
outarrayAdv[i], 1,
outarrayDiff[i], 1,
outarray[i], 1);
}
Array<OneD, Array<OneD, NekDouble> > outarrayForcing(nvariables);
for (i = 0; i < nvariables; ++i)
{
outarrayForcing[i] = Array<OneD, NekDouble>(npoints, 0.0);
}
GetForcingTerm(inarray, outarrayForcing);
for (i = 0; i < nvariables; ++i)
{
// Add Forcing Term
Vmath::Vadd(npoints,
outarray[i], 1,
outarrayForcing[i], 1,
outarray[i], 1);
}
}
// Add sponge layer if defined in the session file
std::vector<SolverUtils::ForcingSharedPtr>::const_iterator x;
for (x = m_forcing.begin(); x != m_forcing.end(); ++x)
{
(*x)->Apply(m_fields, inarray, outarray, time);
}
}