Commit 1b4b0101 authored by Douglas Serson's avatar Douglas Serson

Fix number of points in normals calculation

parent e1411f8a
......@@ -1443,7 +1443,7 @@ namespace Nektar
// interpolate Jacobian and invert
LibUtilities::Interp1D(
from_key,jac, m_base[0]->GetPointsKey(), work);
Vmath::Sdiv(nq,1.0,&work[0],1,&work[0],1);
Vmath::Sdiv(nqe,1.0,&work[0],1,&work[0],1);
// interpolate
for (i = 0; i < GetCoordim(); ++i)
......
......@@ -961,7 +961,7 @@ namespace Nektar
// interpolate Jacobian and invert
LibUtilities::Interp1D(from_key,jac,m_base[0]->GetPointsKey(),work);
Vmath::Sdiv(nq,1.0,&work[0],1,&work[0],1);
Vmath::Sdiv(nqe,1.0,&work[0],1,&work[0],1);
// interpolate
for(i = 0; i < GetCoordim(); ++i)
......
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