Commit 3c60283d authored by Andrew Comerford's avatar Andrew Comerford
Browse files

Minor change to low energy routine


git-svn-id: https://gforge.sci.utah.edu/svn/nektar/trunk@4044 305cdda6-5ce1-45b3-a98d-dfc68c8b3305
parent 33677192
......@@ -1671,21 +1671,16 @@ namespace Nektar
int nLocal = m_locToGloMap->GetNumLocalBndCoeffs();
int nLocalDir = m_locToGloMap->GetNumLocalDirBndCoeffs();
int nNonDir = nGlobal-nDir;
DNekScalBlkMat &M = (*GloBlkMat);
int nblks = expList->GetExpSize();
Array<OneD,unsigned int> exp_size(nblks);
for(i=0; i<expList->GetExpSize(); ++i)
{
exp_size[i]=m_transformationmatrix->GetRows();
}
const Array<OneD,const unsigned int>& exp_size = m_locToGloMap->GetNumLocalBndCoeffsPerPatch();
MatrixStorage blkmatStorage = eDIAGONAL;
DNekScalBlkMatSharedPtr RBlkMat = MemoryManager<DNekScalBlkMat>::AllocateSharedPtr(exp_size,exp_size,blkmatStorage);
DNekScalBlkMatSharedPtr TRBlkMat = MemoryManager<DNekScalBlkMat>::AllocateSharedPtr(exp_size,exp_size,blkmatStorage);
for(i=0; i<expList->GetExpSize(); ++i)
for(i=0; i<nblks; ++i)
{
RBlkMat->SetBlock(i,i,Blktmp = MemoryManager<DNekScalMat>::AllocateSharedPtr(one,m_transformationmatrix));
TRBlkMat->SetBlock(i,i,Blktmp = MemoryManager<DNekScalMat>::AllocateSharedPtr(one,m_transposedtransformationmatrix));
......@@ -1693,6 +1688,7 @@ namespace Nektar
DNekScalBlkMat &R = (*RBlkMat);
DNekScalBlkMat &RT = (*TRBlkMat);
DNekScalBlkMat &M = (*GloBlkMat);
NekVector<NekDouble> r1(nLocal,0.0);
NekVector<NekDouble> z1(nLocal,0.0);
......
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