function [vertices] = findSaddles(dt,k1s,k2s,levelset,epsilon)

sz = size(k1s);

saddle=1; 
for slice = 1:sz(3)  
 for row = 1:sz(1)
    for col = 1:sz(2)   
      if(((k1s(row,col,slice)*k2s(row,col,slice))<epsilon) & (abs(dt(row,col,slice))<levelset))
        vertices(saddle,1)=row;
	vertices(saddle,2)=col;
	vertices(saddle,3)=slice;
        saddle = saddle+1;
      end
    end
  end
end