Commit b4a970ac authored by Douglas Serson's avatar Douglas Serson

Add checks to make sure npz is consistent with HomModesZ

parent 790e7416
......@@ -78,6 +78,19 @@ namespace Nektar
m_homogeneousBasis->GetNumPoints() /
m_StripZcomm->GetSize());
ASSERTL0( m_homogeneousBasis->GetNumPoints() %
m_StripZcomm->GetSize() == 0,
"HomModesZ should be a multiple of npz.");
if ( (m_homogeneousBasis->GetBasisType() !=
LibUtilities::eFourierHalfModeRe)
&& (m_homogeneousBasis->GetBasisType() !=
LibUtilities::eFourierHalfModeIm) )
{
ASSERTL0( m_planes.num_elements() % 2 == 0,
"HomModesZ/npz should be an even integer.");
}
if(m_useFFT)
{
m_FFT = LibUtilities::GetNektarFFTFactory().CreateInstance(
......
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