Changeset 20179


Ignore:
Timestamp:
02/16/16 15:11:27 (9 years ago)
Author:
jbondzio
Message:

BUG: cover case of 0 segment length

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/exp/exp2levelsetfunction.m

    r20178 r20179  
    4343                w=[segment.w.x segment.w.y];
    4444                verts=[x y];
    45                 t=[x-v(1) y-v(2)]*(w-v)'/norm(w-v)^2; %projection of verts on line defined by v and w
     45                dist_vw2=norm(w-v)^2;
     46                if(dist_vw2==0.),       t=zeros(size(x)); %cover case where segment has length 0
     47                else t=[x-v(1) y-v(2)]*(w-v)'/dist_vw2; end %projection of verts on line defined by v and w
    4648                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
    4749                dist=sqrt(sum(abs(dist_vec).^2,2));
Note: See TracChangeset for help on using the changeset viewer.