Curvature

Curvature computation. More...

Functions/Subroutines
subroutine	mod_lsm_curvature::compute_levelset_curvature_object (my_levelset)
	Compute the curvature within the level set object.
subroutine	mod_lsm_curvature::compute_levelset_curvature_by_method (my_levelset, curvature)
	Compute the curvature with various methods.
subroutine	mod_lsm_curvature::crop_curvature (curvature)
	Crop the curvature to the maximum local possible value.
subroutine	mod_lsm_curvature::compute_levelset_curvature_div_n (my_levelset, curvature)
	Compute the curvature with the standard \( \nabla \cdot \vec{n} \) method.
subroutine	mod_lsm_curvature::compute_levelset_curvature_with_cp (my_levelset, init_curvature, closest_points, proc_id, curvature)
	Compute the curvature with the CP method.

Detailed Description

Curvature computation

.

A set of routines to compute the curvature associated to the level set field. We provide 2 methods:

• the standard method where \( \kappa_{ls} = \nabla \cdot \vec{n} \);
• the closest point method where \( \kappa = \kappa_{ls}( \mathrm{CP} (\vec{x}) ) \) by interpolation.

Function/Subroutine Documentation

compute_levelset_curvature_by_method()

subroutine mod_lsm_curvature::compute_levelset_curvature_by_method	(	class(t_levelset), intent(in), target	my_levelset,
		double precision, dimension(:,:,:), intent(inout)	curvature )

Compute the curvature with various methods.

The method is chosen within the my_levelset object.

compute_levelset_curvature_div_n()

subroutine mod_lsm_curvature::compute_levelset_curvature_div_n	(	class(t_levelset), intent(in)	my_levelset,
		double precision, dimension(:,:,:), intent(out)	curvature )

Compute the curvature with the standard \( \nabla \cdot \vec{n} \) method.

Todo

MCO: find another solution for bndc

Parameters

[in,out]	my_levelset	The input levelset
[out]	curvature	The resulting curvature field

compute_levelset_curvature_with_cp()

subroutine mod_lsm_curvature::compute_levelset_curvature_with_cp	(	class(t_levelset_cp), intent(in)	my_levelset,
		double precision, dimension(:,:,:), intent(in)	init_curvature,
		double precision, dimension(:,:,:,:), intent(in)	closest_points,
		integer, dimension(:,:,:), intent(in)	proc_id,
		double precision, dimension(:,:,:), intent(inout)	curvature )

Compute the curvature with the CP method.

Parameters

[in,out]	my_levelset	The levelset
[in]	init_curvature	The input smooth curvature (usually obtained with normal div method)
[in]	closest_points	The (pre computed) closest points
[in]	proc_id	The (pre computed) closest points associated processor ID (for MPI queries)
[in,out]	curvature	The output curvature