OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fixed: wrong projection basis for tensor subgrid

Browse Source
This commit is contained in:
Arne Morten Kvarving 2022-05-06 11:17:04 +02:00
parent 3e6edab293
commit 73d6e85745
2 changed files with 6 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/ASM/ASMs2Dmx.C
View File

@ -203,8 +203,11 @@ bool ASMs2Dmx::generateFEMTopology ()
// we need to project on something that is not one of our bases
if (ASMmxBase::Type == ASMmxBase::REDUCED_CONT_RAISE_BASIS1 ||
ASMmxBase::Type == ASMmxBase::REDUCED_CONT_RAISE_BASIS2 ||
ASMmxBase::Type == ASMmxBase::DIV_COMPATIBLE)
projB = proj = ASMmxBase::raiseBasis(surf);
else if (ASMmxBase::Type == ASMmxBase::SUBGRID)
projB = proj = m_basis.front()->clone();
else
projB = proj = m_basis[2-ASMmxBase::geoBasis]->clone();
}

3
src/ASM/ASMs3Dmx.C
View File

@ -204,8 +204,11 @@ bool ASMs3Dmx::generateFEMTopology ()
// we need to project on something that is not one of our bases
if (ASMmxBase::Type == ASMmxBase::REDUCED_CONT_RAISE_BASIS1 ||
ASMmxBase::Type == ASMmxBase::REDUCED_CONT_RAISE_BASIS2 ||
ASMmxBase::Type == ASMmxBase::DIV_COMPATIBLE)
projB = proj = ASMmxBase::raiseBasis(svol);
else if (ASMmxBase::Type == ASMmxBase::SUBGRID)
projB = proj = m_basis.front()->clone();
else
projB = proj = m_basis[2-ASMmxBase::geoBasis]->clone();
}