Changeset 7922

Timestamp:

04/15/11 08:05:34 (14 years ago)

Author:

Mathieu Morlighem

Message:

Better checks and comments

File:

• issm/trunk/src/c/objects/Loads/Friction.cpp (modified) (3 diffs)

issm/trunk/src/c/objects/Loads/Friction.cpp

@@ -110 +110 @@
         else _error_("element_type %s not supported yet",element_type);
 
+        /*Checks that s-1>0 if v=0*/
+        if(vmag==0 && (s-1)<0) _error_("velocity is 0 ans (s-1)=%g<0, alpha_complement is Inf",s-1);
+
         alpha2=pow(drag_coefficient,2)*pow(Neff,r)*pow(vmag,(s-1));
         _assert_(!isnan(alpha2));
@@ -171 +174 @@
         }
         else _error_("element_type %s not supported yet",element_type);
+
+        /*Checks that s-1>0 if v=0*/
+        if(vmag==0 && (s-1)<0) _error_("velocity is 0 ans (s-1)=%g<0, alpha_complement is Inf",s-1);
 
         alpha2=pow(drag_coefficient,2)*pow(Neff,r)*pow(vmag,(s-1));
@@ -230 +236 @@
         vmag=sqrt(pow(vx,2)+pow(vy,2));
 
-        /*Checks*/
+        /*Checks that s-1>0 if v=0*/
         if(vmag==0 && (s-1)<0) _error_("velocity is 0 ans (s-1)=%g<0, alpha_complement is Inf",s-1);