From 5d9d8e302d62232d7a55d50f96e7311e8cfd82ff Mon Sep 17 00:00:00 2001 From: Anthony Geay Date: Wed, 28 Dec 2022 10:20:40 +0100 Subject: [PATCH] WIP --- src/MEDLoader/MEDFileBlowStrEltUp.cxx | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/src/MEDLoader/MEDFileBlowStrEltUp.cxx b/src/MEDLoader/MEDFileBlowStrEltUp.cxx index db2b47512..aba6b59f6 100644 --- a/src/MEDLoader/MEDFileBlowStrEltUp.cxx +++ b/src/MEDLoader/MEDFileBlowStrEltUp.cxx @@ -410,6 +410,9 @@ MCAuto LocInfo::BuildMeshCommon(INTERP_KERNEL::NormalizedCellT fakeF->setGaussLocalizationOnType(gt,loc.getRefCoords(),loc.getGaussCoords(),loc.getGaussWeights()); ptsForLoc=fakeF->getLocalizationOfDiscr(); // + if(ptsForLoc->getNumberOfComponents()!=3) + ptsForLoc = ptsForLoc->changeNbOfComponents(3,0.0); + // return geoMesh; } @@ -533,7 +536,11 @@ MCAuto LocInfo::BuildMeshPipeSEG3(const DataArrayDouble *angle, geoMesh2->convertQuadraticCellsToLinear(); dir=geoMesh2->buildDirectionVectorField(); } - MCAuto rot(dir->getArray()->fromCartToSpher()); + MCAuto dirArr; + dirArr.takeRef(dir->getArray()); + if(dir->getNumberOfComponents()!=3) + dirArr = dirArr->changeNbOfComponents(3,0.0); + MCAuto rot(dirArr->fromCartToSpher()); std::size_t nbCompo(ptsForLoc->getNumberOfComponents()); MCAuto secPts(section->getCoords()->changeNbOfComponents(nbCompo,0.)); { -- 2.39.2