Commit caaa38c2 authored by Douglas Serson's avatar Douglas Serson Committed by Spencer Sherwin
Browse files

Fix bug in ContField3DHomogeneous2D::v_HelmSolve

parent 726e4844
......@@ -182,10 +182,11 @@ namespace Nektar
// Fourier transform forcing function
HomogeneousFwdTrans(inarray,fce,(flags.isSet(eUseGlobal))?eGlobal:eLocal);
}
int l = 0;
for(n = 0; n < nhom_modes_z; ++n)
{
for(m = 0; m < nhom_modes_y; ++m)
for(m = 0; m < nhom_modes_y; ++m, l++)
{
beta_z = 2*M_PI*(n/2)/m_lhom_z;
beta_y = 2*M_PI*(m/2)/m_lhom_y;
......@@ -193,13 +194,13 @@ namespace Nektar
new_factors = factors;
new_factors[StdRegions::eFactorLambda] += beta;
m_lines[n]->HelmSolve(fce + cnt,
m_lines[l]->HelmSolve(fce + cnt,
e_out = outarray + cnt1,
flags, new_factors, varcoeff, dirForcing);
cnt += m_lines[n]->GetTotPoints();
cnt += m_lines[l]->GetTotPoints();
cnt1 += m_lines[n]->GetNcoeffs();
cnt1 += m_lines[l]->GetNcoeffs();
}
}
......
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