get_dt_normal_vector_kappa.m

LevelSetMethods水平集matlab程序

function [dt] = get_dt_normal_vector_kappa(alpha, dx, dy, H1_abs, H2_abs, b, dx2, dy2)
%
% Calculate the Euler time step.
%
% Author: Baris Sumengen  [e-mail]
% http://vision.ece.ucsb.edu/~sumengen/
%

if alpha <= 0 | alpha >= 1 
    error('alpha needs to be between 0 and 1!');
end

maxs = max(max( H1_abs/dx + H2_abs/dy + (2*b)/dx2 + (2*b)/dy2 ));
dt = alpha/(maxs+(maxs==0));