SubStep3 3D weirdness
Something weird shows up with Windy cases.
SubStep3 proptotype is:
void
SubStep3( TheoreticalGrid3D const& theoreticalGrid,
ViscosityTensor const& viscosityTensor,
PolarGrid const& Rho, PolarGrid const& Energy, PolarGrid const& Pressure,
PolarGrid& EnergyNew,
real dt)
Energy
seems to be the input energy and EnergyNew
the output energy.
In sourceEuler we call SubStep3 like that:
SubStep3 (gas().theoreticalGrid(), viscosityTensor, gas().density(), gas().energy(), gas().pressure(), gas().energy(), dt);
That is, gas().energy() is both 'in' and 'out', which is weird.
When I try to ''fix it'' with:
PolarGrid energy(gas().energy().name(), disk);
SubStep3 (gas().theoreticalGrid(), viscosityTensor, gas().density(), gas().energy(), gas().pressure(), energy, dt);
swap(gas().energy(), energy);
integ3DWindAlpha1|4
and integ3DWind1|4
fails.
Where is the truth ?