Changeset 18486

Timestamp:

09/09/14 02:31:03 (11 years ago)

Author:

Mathieu Morlighem

Message:

NEW: make sure to add background vertices only if they are not outside of the geometric boundary

File:

• issm/trunk-jpl/src/c/bamg/Mesh.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/bamg/Mesh.cpp

@@ -3143 +3143 @@
 
                 /*First, insert old points if requested*/
-                if (KeepVertices && (&Bh != this) && (nbv+Bh.nbv< maxnbv)){
+                if(KeepVertices && (&Bh != this) && (nbv+Bh.nbv< maxnbv)){
                         if (verbose>5) _printf_("         Inserting initial mesh points\n");
+                        bool pointsoutside = false;
                         for (i=0;i<Bh.nbv;i++){ 
                                 BamgVertex &bv=Bh[i];
+                                /*Do not insert if the point is outside*/
+                                long long det3[3];
+                                Triangle* tcvj=TriangleFindFromCoord(bv.i,det3);
+                                if(tcvj->det<0){
+                                        pointsoutside = true;
+                                        continue;
+                                }
                                 if (!bv.GeomEdgeHook){
                                         vertices[nbv].r              = bv.r;
@@ -3153 +3161 @@
                                 }
                         }
+                        if(pointsoutside) _printf_("WARNING: One or more points of the initial mesh fall outside of the geometric boundary\n");
                         Bh.CreateSingleVertexToTriangleConnectivity();     
                         InsertNewPoints(nbvold,NbTSwap,bamgopts->random);