Changeset 877 for XIOS/trunk/extern/remap/src
- Timestamp:
- 07/01/16 11:50:37 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
XIOS/trunk/extern/remap/src/meshutil.cpp
r844 r877 79 79 double area=0 ; 80 80 int k ; 81 81 int count=0 ; 82 82 83 for (int j = 0; j < elt.n; j++) 83 84 { … … 108 109 109 110 area+=polygonarea(neighBaryc, 3) ; 111 count++ ; 110 112 111 113 } 112 grad=grad*(1./area) ; 113 return grad - elt.x * scalarprod(elt.x, grad); 114 if (count>0) 115 { 116 grad=grad*(1./area) ; 117 return grad - elt.x * scalarprod(elt.x, grad); 118 } 119 else return grad ; 114 120 } 115 121
Note: See TracChangeset
for help on using the changeset viewer.