Commit 8f16a3fb authored by Chris Cantwell's avatar Chris Cantwell
Browse files

Fixed misplaced ASSERTL2 checks.

parent a6b85658
......@@ -116,14 +116,19 @@ namespace Nektar
NekDouble StdExpansion2D::v_PhysEvaluate(const Array<OneD, const NekDouble>& coords, const Array<OneD, const NekDouble> & physvals)
{
Array<OneD, NekDouble> coll(2);
Array<OneD, DNekMatSharedPtr> I(2);
Array<OneD, DNekMatSharedPtr> I(2);
ASSERTL2(coords[0] > -1 - NekConstants::kNekZeroTol, "coord[0] < -1");
ASSERTL2(coords[0] < 1 + NekConstants::kNekZeroTol, "coord[0] > 1");
ASSERTL2(coords[1] > -1 - NekConstants::kNekZeroTol, "coord[1] < -1");
ASSERTL2(coords[1] < 1 + NekConstants::kNekZeroTol, "coord[1] > 1");
LocCoordToLocCollapsed(coords,coll);
I[0] = m_base[0]->GetI(coll);
I[1] = m_base[1]->GetI(coll+1);
I[0] = m_base[0]->GetI(coll);
I[1] = m_base[1]->GetI(coll+1);
return v_PhysEvaluate(I,physvals);
return v_PhysEvaluate(I,physvals);
}
NekDouble StdExpansion2D::v_PhysEvaluate(
......@@ -136,11 +141,6 @@ namespace Nektar
int nq1 = m_base[1]->GetNumPoints();
Array<OneD, NekDouble> wsp1(nq1);
ASSERTL2(coords[0] > -1 - NekConstants::kNekZeroTol, "coord[0] < -1");
ASSERTL2(coords[0] < 1 + NekConstants::kNekZeroTol, "coord[0] > 1");
ASSERTL2(coords[1] > -1 - NekConstants::kNekZeroTol, "coord[1] < -1");
ASSERTL2(coords[1] < 1 + NekConstants::kNekZeroTol, "coord[1] > 1");
// interpolate first coordinate direction
for (i = 0; i < nq1;++i)
{
......
Supports Markdown
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