Changeset 19478

Timestamp:

08/03/15 15:31:34 (10 years ago)

Author:

Mathieu Morlighem

Message:

BUG: double should be IssmDouble

File:

• issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/classes/Elements/Tria.cpp

@@ -468 +468 @@
 
         /*Intermediaries*/
-        double d,xn,yn;
+        IssmDouble d,xn,yn;
 
         /*Get current levelset and vertex coordinates*/
@@ -479 +479 @@
         /*Get distance from list of segments and reset ls*/
         for(int j=0;j<NUMVERTICES;j++){
-                double dmin = 1.e+50;
+                IssmDouble dmin = 1.e+50;
                 for(int i=0;i<numsegments;i++){
                         IssmDouble x = xyz_list[j][0];
@@ -495 +495 @@
                          *We find projection of point p onto the line.
                          *It falls where t = [(p-v) . (w-v)] / |w-v|^2*/
-                        double t = ((x-segments[4*i+0])*(segments[4*i+2]-segments[4*i+0]) + (y-segments[4*i+1])*(segments[4*i+3]-segments[4*i+1]))/l2;
+                        IssmDouble t = ((x-segments[4*i+0])*(segments[4*i+2]-segments[4*i+0]) + (y-segments[4*i+1])*(segments[4*i+3]-segments[4*i+1]))/l2;
                         if(t < 0.0){
                                 // Beyond the 'v' end of the segment