Changeset 20179

Timestamp:

02/16/16 15:11:27 (9 years ago)

Author:

jbondzio

Message:

BUG: cover case of 0 segment length

File:

• issm/trunk-jpl/src/m/exp/exp2levelsetfunction.m (modified) (1 diff)

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

@@ -43 +43 @@
                 w=[segment.w.x segment.w.y];
                 verts=[x y];
-                t=[x-v(1) y-v(2)]*(w-v)'/norm(w-v)^2; %projection of verts on line defined by v and w
+                dist_vw2=norm(w-v)^2;
+                if(dist_vw2==0.),       t=zeros(size(x)); %cover case where segment has length 0
+                else t=[x-v(1) y-v(2)]*(w-v)'/dist_vw2; end %projection of verts on line defined by v and w
                 dist_vec=(ones(length(x),1)*v+max(0,min(1,t))*(w-v))-verts; %vector of shortest distance between verts and segment v-w
                 dist=sqrt(sum(abs(dist_vec).^2,2));