Commit f8ca5c6d authored by Chris Cantwell's avatar Chris Cantwell
Browse files

Fixed bugs in solvers due to geom factors updates.

Updated values in some particularly sensitive compressible flow solver tests.
parent 99719e31
......@@ -250,7 +250,7 @@ namespace Nektar
phys_offset = pFields[0]->GetPhys_Offset(n);
jac = LocalRegions::Expansion2D::FromStdExp(pFields[0]->GetExp(n))->GetGeom2D()->GetMetricInfo()->GetJac();
gmat = LocalRegions::Expansion2D::FromStdExp(pFields[0]->GetExp(n))->GetGeom2D()->GetMetricInfo()->GetGmat();
gmat = LocalRegions::Expansion2D::FromStdExp(pFields[0]->GetExp(n))->GetGeom2D()->GetMetricInfo()->GetDerivFactors();
if (LocalRegions::Expansion2D::FromStdExp(pFields[0]->GetExp(n))->GetGeom2D()->GetMetricInfo()->GetGtype()
== SpatialDomains::eDeformed)
......
......@@ -293,7 +293,7 @@ namespace Nektar
phys_offset = pFields[0]->GetPhys_Offset(n);
jac = LocalRegions::Expansion2D::FromStdExp(pFields[0]->GetExp(n))->GetGeom2D()->GetMetricInfo()->GetJac();
gmat = LocalRegions::Expansion2D::FromStdExp(pFields[0]->GetExp(n))->GetGeom2D()->GetMetricInfo()->GetGmat();
gmat = LocalRegions::Expansion2D::FromStdExp(pFields[0]->GetExp(n))->GetGeom2D()->GetMetricInfo()->GetDerivFactors();
if (LocalRegions::Expansion2D::FromStdExp(pFields[0]->GetExp(n))->GetGeom2D()->GetMetricInfo()->GetGtype()
== SpatialDomains::eDeformed)
......
......@@ -359,7 +359,7 @@ namespace Nektar
Array<OneD, const NekDouble> jac =
el2D->GetGeom2D()->GetMetricInfo()->GetJac();
Array<TwoD, const NekDouble> gmat =
el2D->GetGeom2D()->GetMetricInfo()->GetGmat();
el2D->GetGeom2D()->GetMetricInfo()->GetDerivFactors();
if (el2D->GetGeom2D()->GetMetricInfo()->GetGtype()
== SpatialDomains::eDeformed)
......@@ -411,7 +411,7 @@ namespace Nektar
Array<OneD, const NekDouble> jac =
el3D->GetGeom3D()->GetMetricInfo()->GetJac();
Array<TwoD, const NekDouble> gmat =
el3D->GetGeom3D()->GetMetricInfo()->GetGmat();
el3D->GetGeom3D()->GetMetricInfo()->GetDerivFactors();
if (el3D->GetGeom3D()->GetMetricInfo()->GetGtype()
== SpatialDomains::eDeformed)
......@@ -436,7 +436,7 @@ namespace Nektar
Array<OneD, const NekDouble> jac =
el3D->GetGeom3D()->GetMetricInfo()->GetJac();
Array<TwoD, const NekDouble> gmat =
el3D->GetGeom3D()->GetMetricInfo()->GetGmat();
el3D->GetGeom3D()->GetMetricInfo()->GetDerivFactors();
for (int i = 0; i < n_points; i++)
{
......
......@@ -8,16 +8,16 @@
</files>
<metrics>
<metric type="L2" id="1">
<value variable="rho" tolerance="1e-12">0.087961</value>
<value variable="rho" tolerance="1e-12">0.0879609</value>
<value variable="rhou" tolerance="1e-12">60.336</value>
<value variable="rhov" tolerance="1e-8">0.227323</value>
<value variable="rhov" tolerance="1e-8">0.227321</value>
<value variable="E" tolerance="1e-12">4924.14</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.0736791</value>
<value variable="rhou" tolerance="1e-12">61.0599</value>
<value variable="rhov" tolerance="2e-6">0.262734</value>
<value variable="E" tolerance="1e-12">4423.8</value>
<value variable="rho" tolerance="1e-12">0.0736795</value>
<value variable="rhou" tolerance="1e-12">61.0601</value>
<value variable="rhov" tolerance="2e-6">0.262767</value>
<value variable="E" tolerance="1e-12">4423.81</value>
</metric>
</metrics>
</test>
......
......@@ -8,16 +8,16 @@
</files>
<metrics>
<metric type="L2" id="1">
<value variable="rho" tolerance="1e-12">0.0878866</value>
<value variable="rho" tolerance="1e-12">0.0878864</value>
<value variable="rhou" tolerance="1e-12">60.3318</value>
<value variable="rhov" tolerance="1e-8">0.229391</value>
<value variable="E" tolerance="1e-12">4925.54</value>
<value variable="rhov" tolerance="1e-8">0.229384</value>
<value variable="E" tolerance="1e-12">4925.53</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.0743004</value>
<value variable="rhou" tolerance="1e-12">61.7703</value>
<value variable="rhov" tolerance="2e-6">0.278359</value>
<value variable="E" tolerance="1e-12">4406.13</value>
<value variable="rho" tolerance="1e-12">0.0743007</value>
<value variable="rhou" tolerance="1e-12">61.7707</value>
<value variable="rhov" tolerance="2e-6">0.278381</value>
<value variable="E" tolerance="1e-12">4406.17</value>
</metric>
</metrics>
</test>
......
......@@ -8,16 +8,16 @@
</files>
<metrics>
<metric type="L2" id="1">
<value variable="rho" tolerance="1e-12">0.0878866</value>
<value variable="rho" tolerance="1e-12">0.0878864</value>
<value variable="rhou" tolerance="1e-12">60.3318</value>
<value variable="rhov" tolerance="1e-8">0.229391</value>
<value variable="E" tolerance="1e-12">4925.54</value>
<value variable="rhov" tolerance="1e-8">0.229384</value>
<value variable="E" tolerance="1e-12">4925.53</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.0743004</value>
<value variable="rhou" tolerance="1e-12">61.7703</value>
<value variable="rhov" tolerance="2e-6">0.278359</value>
<value variable="E" tolerance="1e-12">4406.13</value>
<value variable="rho" tolerance="1e-12">0.0743007</value>
<value variable="rhou" tolerance="1e-12">61.7707</value>
<value variable="rhov" tolerance="2e-6">0.278381</value>
<value variable="E" tolerance="1e-12">4406.17</value>
</metric>
</metrics>
</test>
......
......@@ -8,16 +8,16 @@
</files>
<metrics>
<metric type="L2" id="1">
<value variable="rho" tolerance="1e-12">0.0878646</value>
<value variable="rhou" tolerance="1e-12">60.3308</value>
<value variable="rhov" tolerance="1e-8">0.230215</value>
<value variable="rho" tolerance="1e-12">0.0878644</value>
<value variable="rhou" tolerance="1e-12">60.3309</value>
<value variable="rhov" tolerance="1e-8">0.230213</value>
<value variable="E" tolerance="1e-12">4926.56</value>
</metric>
<metric type="Linf" id="2">
<value variable="rho" tolerance="1e-12">0.0743362</value>
<value variable="rhou" tolerance="1e-12">60.8039</value>
<value variable="rhov" tolerance="1e-8">0.302484</value>
<value variable="E" tolerance="1e-12">4448.58</value>
<value variable="rho" tolerance="1e-12">0.0743376</value>
<value variable="rhou" tolerance="1e-12">60.8075</value>
<value variable="rhov" tolerance="1e-8">0.302439</value>
<value variable="E" tolerance="1e-12">4448.92</value>
</metric>
</metrics>
</test>
......
......@@ -1239,7 +1239,7 @@ namespace Nektar
}
Array<TwoD, const NekDouble> gmat =
el3D->GetGeom3D()->GetMetricInfo()->GetGmat();
el3D->GetGeom3D()->GetMetricInfo()->GetDerivFactors();
if (el3D->GetGeom3D()->GetMetricInfo()->GetGtype()
== SpatialDomains::eDeformed)
......
......@@ -8,16 +8,16 @@
</files>
<metrics>
<metric type="L2" id="1">
<value variable="u" tolerance="1e-12">0</value>
<value variable="v" tolerance="1e-12">0</value>
<value variable="w" tolerance="1e-12">4.44507e-16</value>
<value variable="p" tolerance="1e-12">1.52206e-14</value>
<value variable="u" tolerance="1e-08">0</value>
<value variable="v" tolerance="1e-08">0</value>
<value variable="w" tolerance="1e-08">4.44507e-16</value>
<value variable="p" tolerance="1e-08">1.52206e-14</value>
</metric>
<metric type="Linf" id="2">
<value variable="u" tolerance="1e-12">5.28658e-16</value>
<value variable="v" tolerance="1e-12">3.37939e-16</value>
<value variable="w" tolerance="1e-12">1.39888e-14</value>
<value variable="p" tolerance="1e-12">9.61453e-14</value>
<value variable="u" tolerance="1e-08">5.28658e-16</value>
<value variable="v" tolerance="1e-08">3.37939e-16</value>
<value variable="w" tolerance="1e-08">1.39888e-14</value>
<value variable="p" tolerance="1e-08">9.61453e-14</value>
</metric>
</metrics>
</test>
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment