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. More...

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.