Commit 23443e34 authored by Dave Moxey's avatar Dave Moxey
Browse files

Fixing memory overflow bug in TetExp (occurs when using direct/iterative full

solver => matrix-free Helmholtz operation).


git-svn-id: https://gforge.sci.utah.edu/svn/nektar/trunk@4019 305cdda6-5ce1-45b3-a98d-dfc68c8b3305
parent 14913b44
......@@ -1212,8 +1212,8 @@ namespace Nektar
int nmodes1 = m_base[1]->GetNumModes();
int nmodes2 = m_base[2]->GetNumModes();
int wspsize = max(nquad0*nmodes2*(nmodes1+nquad1),
nquad0*nquad1*(nquad2+nmodes0)+
nmodes0*nmodes1*nquad2);
nquad2*nmodes0*nmodes1*(nmodes1+1)/2+
nquad2*nquad1*nmodes0);
NekDouble lambda = mkey.GetConstFactor(StdRegions::eFactorLambda);
......@@ -1298,8 +1298,8 @@ namespace Nektar
int nmodes1 = m_base[1]->GetNumModes();
int nmodes2 = m_base[2]->GetNumModes();
int wspsize = max(nquad0*nmodes2*(nmodes1+nquad1),
nquad0*nquad1*(nquad2+nmodes0)+
nmodes0*nmodes1*nquad2);
nquad2*nmodes0*nmodes1*(nmodes1+1)/2+
nquad2*nquad1*nmodes0);
const Array<OneD, const NekDouble>& base0 = m_base[0]->GetBdata ();
const Array<OneD, const NekDouble>& base1 = m_base[1]->GetBdata ();
......
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