mod_lsm_reinitialization_hj_first_order Module Reference

Reinitialization method using the classic HJ algorithm. More...

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)

Detailed Description

Reinitialization method using the classic HJ algorithm.

Function/Subroutine Documentation

compute_levelset_reinitialization_hj_first_order()

subroutine mod_lsm_reinitialization_hj_first_order::compute_levelset_reinitialization_hj_first_order	(	class(t_levelset), intent(inout)	my_levelset,
		integer, dimension(:,:,:), intent(in)	mask,
		integer, intent(in)	n_iteration )

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 )

Integrate with Euler time scheme the HJ equation.

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.

integrate_hj_equation_array_weno5()

subroutine mod_lsm_reinitialization_hj_first_order::integrate_hj_equation_array_weno5	(	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 )