get_dt_normal_kappa.m

来自「LevelSetMethods水平集matlab程序」· M 代码 · 共 16 行

function [dt] = get_dt_normal_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));