Changeset 2269 for XIOS/dev/dev_ym/XIOS_COUPLING/extern/remap/src/intersection_ym.cpp

Timestamp:

12/07/21 12:00:23 (3 years ago)

Author:

ymipsl

Message:

• Solve memory leak from remapper.
• shared_ptr add add for manage nodes.

YM

File:

• XIOS/dev/dev_ym/XIOS_COUPLING/extern/remap/src/intersection_ym.cpp (modified) (2 diffs)

XIOS/dev/dev_ym/XIOS_COUPLING/extern/remap/src/intersection_ym.cpp

@@ -186 +186 @@
        if (dx*dx+dy*dy>fusion_vertex*fusion_vertex)
        {
-          intersectPolygon[nv]=intersectPolygon[n] ;
-
-          vect_points.push_back( array<double, 2>() );
-          vect_points[nv][0] = intersectPolygon[n].x;
-          vect_points[nv][1] = intersectPolygon[n].y;
-
-          nv++ ;
+         intersectPolygon[nv]=intersectPolygon[n] ;
+         vect_points.push_back( array<double, 2>() );
+         vect_points[nv][0] = intersectPolygon[n].x;
+         vect_points[nv][1] = intersectPolygon[n].y;
+         nv++ ;
        }
-      
-
     }
 
@@ -209 +205 @@
       for(int i=0;i<indices.size()/3;++i)
       {
-          Coord x0 = Ox * polyline[0][indices[3*i]][0] + Oy* polyline[0][indices[3*i]][1] + Oz ;
-          Coord x1 = Ox * polyline[0][indices[3*i+1]][0] + Oy* polyline[0][indices[3*i+1]][1] + Oz ;
-          Coord x2 = Ox * polyline[0][indices[3*i+2]][0] + Oy* polyline[0][indices[3*i+2]][1] + Oz ;
-          area2+=triarea(x0 * (1./norm(x0)),x1* (1./norm(x1)), x2* (1./norm(x2))) ;
+        Coord x0 = Ox * polyline[0][indices[3*i]][0] + Oy* polyline[0][indices[3*i]][1] + Oz ;
+        Coord x1 = Ox * polyline[0][indices[3*i+1]][0] + Oy* polyline[0][indices[3*i+1]][1] + Oz ;
+        Coord x2 = Ox * polyline[0][indices[3*i+2]][0] + Oy* polyline[0][indices[3*i+2]][1] + Oz ;
+        area2+=triarea(x0 * (1./norm(x0)),x1* (1./norm(x1)), x2* (1./norm(x2))) ;
       }