evolve_normal_vector_weno.m

LevelSetMethods水平集matlab程序

function [delta, H1_abs, H2_abs] = evolve_normal_vector_WENO(phi, dx, dy, Vn_ext, u_ext, v_ext)
%
% Finds the amount of evolution under a force in
% normal direction and a force based on a vector field,
% and using 5th order accurate WENO scheme.
% Does not assume that phi is approximately a signed
% distance function and uses SLLF scheme.
%
% Author: Baris Sumengen  [e-mail]
% http://vision.ece.ucsb.edu/~sumengen/
%

delta = zeros(size(phi)+6);
data_ext = zeros(size(phi)+6);
data_ext(4:end-3,4:end-3) = phi;

% Calculate the derivatives (both + and -)
phi_x_minus = zeros(size(phi)+6);
phi_x_plus = zeros(size(phi)+6);
phi_y_minus = zeros(size(phi)+6);
phi_y_plus = zeros(size(phi)+6);
% first scan the rows
for i=1:size(phi,1)
	phi_x_minus(i+3,:) = der_WENO_minus(data_ext(i+3,:), dx);
	phi_x_plus(i+3,:) = der_WENO_plus(data_ext(i+3,:), dx);	
end

% then scan the columns
for j=1:size(phi,2)
	phi_y_minus(:,j+3) = der_WENO_minus(data_ext(:,j+3), dy);	
	phi_y_plus(:,j+3) = der_WENO_plus(data_ext(:,j+3), dy);	
end


[delta, H1_abs, H2_abs] = LLF_normal_vector(dx, dy, Vn_ext, u_ext, v_ext, phi_x_minus, phi_x_plus, phi_y_minus, phi_y_plus);
H1_abs = H1_abs(4:end-3,4:end-3);
H2_abs = H2_abs(4:end-3,4:end-3);
delta = delta(4:end-3,4:end-3);