From 73d6e85745c9de9f530d189ca0df88f0fac40293 Mon Sep 17 00:00:00 2001 From: Arne Morten Kvarving Date: Fri, 6 May 2022 11:17:04 +0200 Subject: [PATCH] fixed: wrong projection basis for tensor subgrid --- src/ASM/ASMs2Dmx.C | 3 +++ src/ASM/ASMs3Dmx.C | 3 +++ 2 files changed, 6 insertions(+) diff --git a/src/ASM/ASMs2Dmx.C b/src/ASM/ASMs2Dmx.C index 51617d23..2f1770bb 100644 --- a/src/ASM/ASMs2Dmx.C +++ b/src/ASM/ASMs2Dmx.C @@ -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(); } diff --git a/src/ASM/ASMs3Dmx.C b/src/ASM/ASMs3Dmx.C index 63f6ad3f..2bcf5926 100644 --- a/src/ASM/ASMs3Dmx.C +++ b/src/ASM/ASMs3Dmx.C @@ -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(); }