Changeset 20512

Timestamp:

04/14/16 19:41:00 (9 years ago)

Author:

Mathieu Morlighem

Message:

BUG: take care of case vx=vy=0

File:

• issm/trunk-jpl/src/m/exp/flowlines.m (modified) (2 diffs)

issm/trunk-jpl/src/m/exp/flowlines.m

@@ -88 +88 @@
 
         %velocity of the current triangle and norm it
-        ut=u(tria); vt=v(tria); normv=sqrt(ut.^2+vt.^2);
+        ut=u(tria); vt=v(tria); normv=max(eps,sqrt(ut.^2+vt.^2));
         ut=ut./normv;vt=vt./normv;
 
@@ -138 +138 @@
 
         %velocity of the current triangle and norm it
-        ut=-u(tria); vt=-v(tria); normv=sqrt(ut.^2+vt.^2);
+        ut=-u(tria); vt=-v(tria); normv=max(eps,sqrt(ut.^2+vt.^2));
         ut=ut./normv;vt=vt./normv;