Changeset 1178

Timestamp:

06/30/09 10:35:49 (16 years ago)

Author:

Mathieu Morlighem

Message:

Initialization with default value

Location:

issm/trunk/src/mex

Files:

• InterpFromMesh2d/InterpFromMesh2d.cpp (modified) (2 diffs)
• InterpFromMesh3d/InterpFromMesh3d.cpp (modified) (2 diffs)

issm/trunk/src/mex/InterpFromMesh2d/InterpFromMesh2d.cpp

@@ -116 +116 @@
         /*Initialize output*/
         data_prime=NewVec(nods_prime);
+        for (i=0;i<nods_prime;i++) VecSetValue(data_prime,i,default_value,INSERT_VALUES);
 
+        /*Loop over the elements*/
         printf("\n   interpolation progress:   %5.2lf %%",0.0);
         for (i=0;i<nels_data;i++){
@@ -146 +148 @@
 
                         /*is the current point in the current element?*/
-                        if ((area_1>0 || area_1==0) && (area_2>0 || area_2==0) && (area_3>0 || area_3==0)){
+                        if (area_1>=0 && area_2>=0 && area_3>=0){
 
                                 /*Yes ! compute the value on the point*/

issm/trunk/src/mex/InterpFromMesh3d/InterpFromMesh3d.cpp

@@ -121 +121 @@
         /*Initialize output*/
         data_prime=NewVec(nods_prime);
-
+        for (i=0;i<nods_prime;i++) VecSetValue(data_prime,i,default_value,INSERT_VALUES);
+
+        /*Loop over the elements*/
         printf("\n   interpolation progress:   %5.2lf %%",0.0);
         for (i=0;i<nels_data;i++){
@@ -151 +153 @@
 
                         /*is the current point in the current 2d element?*/
-                        if ((area_1>0 || area_1==0) && (area_2>0 || area_2==0) && (area_3>0 || area_3==0)){
+                        if (area_1>=0 && area_2>=0 && area_3>=0){
 
                                 /*compute bottom and top height of the element at this 2d position*/