ply_lambda Function

public elemental function ply_lambda(val) result(funcVal)

\todo: as we use a relation of gamma, it might be better to use the gammln function provided by the numerical recipes, and just use the difference in an exponential function.

Arguments

TypeIntentOptionalAttributesName
real(kind=rk), intent(in) :: val

Return Value real(kind=rk)


Called by

proc~~ply_lambda~~CalledByGraph proc~ply_lambda ply_lambda proc~ply_l_int ply_l_int proc~ply_l_int->proc~ply_lambda proc~ply_l ply_l proc~ply_l_int->proc~ply_l proc~ply_l->proc~ply_lambda proc~ply_m ply_m proc~ply_m->proc~ply_lambda proc~ply_m_int ply_m_int proc~ply_m_int->proc~ply_lambda proc~ply_fpt_init ply_fpt_init proc~ply_fpt_init->proc~ply_l_int proc~ply_fpt_init->proc~ply_l proc~ply_fpt_init->proc~ply_m proc~ply_fpt_init->proc~ply_m_int proc~ply_init_legfpt ply_init_legFpt proc~ply_init_legfpt->proc~ply_fpt_init proc~ply_poly_project_fillbody ply_poly_project_fillbody proc~ply_poly_project_fillbody->proc~ply_init_legfpt

Contents


Variables

TypeVisibilityAttributesNameInitial
real(kind=rk), private, parameter:: lb_poly =1.0_rk/0.067_rk
real(kind=rk), private, parameter:: inter1 =0.058_rk
real(kind=rk), private, parameter:: inter2 =0.04_rk
real(kind=rk), private, parameter:: invbound =0.02_rk
real(kind=rk), private :: invVal
real(kind=rk), private :: a0
real(kind=rk), private :: a1
real(kind=rk), private :: a2
real(kind=rk), private :: a3
real(kind=rk), private :: a4
real(kind=rk), private :: a5