Commit ce830d4e authored by Dave Moxey's avatar Dave Moxey

Merge branch 'feature/ExtractSurfaceCFS' of /opt/gitlab/repositories/nektar

parents 57a900d5 a5af0100
......@@ -1771,19 +1771,21 @@ namespace Nektar
Array<OneD, NekDouble > &soundspeed,
Array<OneD, NekDouble > &mach)
{
const int nBCEdgePts = m_fields[0]->GetTotPoints();
Vmath::Vmul(nBCEdgePts, physfield[1], 1, physfield[1], 1, mach, 1);
const int nq = m_fields[0]->GetTotPoints();
Vmath::Vmul(nq, physfield[1], 1, physfield[1], 1, mach, 1);
for (int i = 1; i < m_spacedim; ++i)
{
Vmath::Vvtvp(nBCEdgePts, physfield[1+i], 1, physfield[1+i], 1,
mach, 1, mach, 1);
Vmath::Vvtvp(nq, physfield[1+i], 1, physfield[1+i], 1,
mach, 1, mach, 1);
}
Vmath::Vdiv(nBCEdgePts, mach, 1, physfield[0], 1, mach, 1);
Vmath::Vdiv(nBCEdgePts, mach, 1, physfield[0], 1, mach, 1);
Vmath::Vdiv(nBCEdgePts, mach, 1, soundspeed, 1, mach, 1);
Vmath::Vdiv(nq, mach, 1, physfield[0], 1, mach, 1);
Vmath::Vdiv(nq, mach, 1, physfield[0], 1, mach, 1);
Vmath::Vsqrt(nq, mach, 1, mach, 1);
Vmath::Vdiv(nq, mach, 1, soundspeed, 1, mach, 1);
}
/**
......
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