Commit 723c0310 authored by Spencer Sherwin's avatar Spencer Sherwin Committed by Chris Cantwell
Browse files

Small check in case elmtid goes over the number on one plane

parent b5532354
......@@ -956,6 +956,10 @@ namespace Nektar
elmtid = Vmath::Imax(n_element,cfl,1);
NekDouble CFL = cfl[elmtid];
if(m_HomogeneousType == eHomogeneous1D) // express element id with respect to plane
elmtid = elmtid%m_fields[0]->GetPlane(0)->GetExpSize();
return CFL;
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