From 044aba47289bf3923a28b4fd4ab9edeb42e1b54c Mon Sep 17 00:00:00 2001 From: abn Date: Wed, 16 Jan 2019 15:09:33 +0100 Subject: [PATCH] restoring intersect test Tmp10 with a little fix to help on non-conform mesh --- .../Geometric2D/InterpKernelGeo2DQuadraticPolygon.cxx | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/INTERP_KERNEL/Geometric2D/InterpKernelGeo2DQuadraticPolygon.cxx b/src/INTERP_KERNEL/Geometric2D/InterpKernelGeo2DQuadraticPolygon.cxx index 298a704ac..4007f2fbe 100644 --- a/src/INTERP_KERNEL/Geometric2D/InterpKernelGeo2DQuadraticPolygon.cxx +++ b/src/INTERP_KERNEL/Geometric2D/InterpKernelGeo2DQuadraticPolygon.cxx @@ -1346,8 +1346,9 @@ void QuadraticPolygon::ComputeResidual(const QuadraticPolygon& pol1, const std:: // Convert to integer connectivity: for(std::list::iterator itConstr=retPolsUnderContruction.begin();itConstr!=retPolsUnderContruction.end();itConstr++) { - if((*itConstr)->getStartNode()==(*itConstr)->getEndNode()) // take only fully closed reconstructed polygon (?? might there be others??) + if((*itConstr)->getStartNode()==(*itConstr)->getEndNode()) // take only fully closed reconstructed polygon { + (*itConstr)->cleanDegeneratedConsecutiveEdges(); (*itConstr)->appendCrudeData(mapp,0.,0.,1.,offset,addCoordsQuadratic,conn,connI); nb1.push_back(idThis); nb2.push_back(-1); for(std::list::iterator it6=pol1ZipConsumed[*itConstr].begin();it6!=pol1ZipConsumed[*itConstr].end();it6++) delete *it6; -- 2.39.2