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

Fix another bug in ExtractCoeffsToCoeffs with variable p

parent 569da6a0
......@@ -2236,27 +2236,19 @@ namespace Nektar
int i;
int offset = 0;
// check if the same and if so just copy over coeffs
if(fromExpList->GetNcoeffs() == m_ncoeffs)
{
Vmath::Vcopy(m_ncoeffs,fromCoeffs,1,toCoeffs,1);
}
else
for(i = 0; i < (*m_exp).size(); ++i)
{
for(i = 0; i < (*m_exp).size(); ++i)
std::vector<unsigned int> nummodes;
int eid = m_offset_elmt_id[i];
for(int j= 0; j < fromExpList->GetExp(eid)->GetNumBases(); ++j)
{
std::vector<unsigned int> nummodes;
int eid = m_offset_elmt_id[i];
for(int j= 0; j < fromExpList->GetExp(eid)->GetNumBases(); ++j)
{
nummodes.push_back(fromExpList->GetExp(eid)->GetBasisNumModes(j));
}
nummodes.push_back(fromExpList->GetExp(eid)->GetBasisNumModes(j));
}
(*m_exp)[eid]->ExtractDataToCoeffs(&fromCoeffs[offset], nummodes,0,
&toCoeffs[m_coeff_offset[eid]]);
(*m_exp)[eid]->ExtractDataToCoeffs(&fromCoeffs[offset], nummodes,0,
&toCoeffs[m_coeff_offset[eid]]);
offset += fromExpList->GetExp(eid)->GetNcoeffs();
}
offset += fromExpList->GetExp(eid)->GetNcoeffs();
}
}
......
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