Changeset 21931 for issm/trunk-jpl/src/c/bamg/Mesh.cpp
- Timestamp:
- 08/10/17 10:23:14 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/bamg/Mesh.cpp
r21915 r21931 1594 1594 long i1 = GetId(triangles[it][VerticesOfTriangularEdge[j][1]]); 1595 1595 k = edge4->SortAndFind(i0,i1); 1596 _assert_(k>=0); 1596 if(k<0){ 1597 delete [] colorV; 1598 _error_("This should not happen"); 1599 } 1597 1600 subdomains[i].direction = (vertices + i0 == edges[k].v[0]) ? 1 : -1; 1598 1601 subdomains[i].edge = edges+k; … … 2054 2057 else if (xy==1) dd=dxdy; 2055 2058 else if (xy==2) dd=dydy; 2056 else _error_("not supported yet"); 2059 else{ 2060 delete [] detT; 2061 delete [] Mmass; 2062 delete [] workT; 2063 delete [] workV; 2064 delete [] Mmassxx; 2065 delete [] OnBoundary; 2066 _error_("not supported yet"); 2067 } 2057 2068 // do leat 2 iteration for boundary problem 2058 2069 for (int ijacobi=0;ijacobi<Max(NbJacobi,2);ijacobi++){ … … 2103 2114 delete [] dxdy; 2104 2115 delete [] dydy; 2105 delete [] 2116 delete [] workT; 2106 2117 delete [] workV; 2107 2118 delete [] Mmassxx; 2108 delete [] 2119 delete [] OnBoundary; 2109 2120 2110 2121 } … … 2394 2405 if (nbt == nbtout || !NbSubDomTot) { 2395 2406 delete [] HeapArete; 2407 delete [] HeapTriangle; 2396 2408 _error_("The boundary is not close: all triangles are outside"); 2397 2409 } … … 2883 2895 } 2884 2896 } 2885 if(kk) _error_("See above"); 2897 if(kk){ 2898 delete [] e; 2899 _error_("See above"); 2900 } 2886 2901 2887 2902 return e; … … 4398 4413 printf("\n"); 4399 4414 if(NbVerticesOnGeomVertex >= maxnbv){ 4415 delete [] bcurve; 4400 4416 _error_("too many vertices on geometry: " << NbVerticesOnGeomVertex << " >= " << maxnbv); 4401 4417 }
Note:
See TracChangeset
for help on using the changeset viewer.