mod_lsm_reinitialization_hj_first_order Module Reference

Reinitialization method using the classic HJ algorithm.

Functions/Subroutines
subroutine	compute_levelset_reinitialization_hj_first_order (my_levelset, mask, n_iteration)
subroutine	integrate_hj_equation_array (phi_0, phi_n, mask, dtau, phi_np1)
	Integrate with Euler time scheme the HJ equation.
subroutine	integrate_hj_equation_array_weno5 (phi_0, phi_n, mask, dtau, phi_np1)

Function/Subroutine Documentation

integrate_hj_equation_array()

subroutine mod_lsm_reinitialization_hj_first_order::integrate_hj_equation_array	(	double precision, dimension(:,:,:), intent(in)	phi_0,
		double precision, dimension(:,:,:), intent(in)	phi_n,
		integer, dimension(:,:,:), intent(in)	mask,
		double precision, intent(in)	dtau,
		double precision, dimension(:,:,:), intent(inout)	phi_np1 )

Note

Before we were using WENO5 when the band==1 and WENO3 in the band==2 (nothing elsewhere) but this seems to be very unstable. Now we are only using WENO5 in the whole band.