Changeset 26585

Timestamp:

11/10/21 09:19:42 (3 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixing a division by 0 in CoDiPack

File:

• issm/trunk-jpl/src/c/classes/Loads/Friction.cpp (modified) (4 diffs)

issm/trunk-jpl/src/c/classes/Loads/Friction.cpp

@@ -198 +198 @@
 
         /*Check to prevent dividing by zero if vmag==0*/
-        if(vmag==0. && (s-1.)<0.) alpha_complement=0.;
+        if(vmag==0. && (s-1.)<=0.) alpha_complement=0.;
         else alpha_complement=pow(Neff,r)*pow(vmag,(s-1));
 
@@ -345 +345 @@
 
         /*Check to prevent dividing by zero if vmag==0*/
-        if(vmag==0. && (s-1.)<0.){
+        if(vmag==0. && (s-1.)<=0.){
                 alpha2=0.;
         }
@@ -541 +541 @@
 
         /*Check to prevent dividing by zero if vmag==0*/
-        if(vmag==0. && (s-1.)<0.) alpha2=0.;
+        if(vmag==0. && (s-1.)<=0.) alpha2=0.;
         else alpha2=drag_coefficient*drag_coefficient*pow(Neff,r)*pow(vmag,(s-1.));
 
@@ -587 +587 @@
 
         /*Check to prevent dividing by zero if vmag==0*/
-        if(vmag==0. && (s-1.)<0.) alpha2=0.;
+        if(vmag==0. && (s-1.)<=0.) alpha2=0.;
         else alpha2=drag_coefficient*drag_coefficient*pow(Neff,r)*pow(vmag,(s-1.));