Impure wrapper function for f_dst_chi2_cdf_core
.
Handles optional arguments and invalid values for arguments.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=wp), | intent(in) | :: | x |
sample position |
||
real(kind=wp), | intent(in) | :: | df |
degrees of freedom |
||
real(kind=wp), | intent(in), | optional | :: | loc |
location parameter |
|
real(kind=wp), | intent(in), | optional | :: | scale |
scale parameter |
|
character(len=*), | intent(in), | optional | :: | tail |
tail options |
resulting CDF value
impure function f_dst_chi2_cdf(x, df, loc, scale, tail) result(p) ! ==== Description !! Impure wrapper function for `f_dst_chi2_cdf_core`. !! Handles optional arguments and invalid values for arguments. ! ==== Declarations real(wp) , intent(in) :: x !! sample position real(wp) , intent(in) :: df !! degrees of freedom real(wp) , intent(in), optional :: loc !! location parameter real(wp) , intent(in), optional :: scale !! scale parameter character(len=*), intent(in), optional :: tail !! tail options real(wp) :: loc_w !! final value for loc real(wp) :: scale_w !! final value for scale character(len=16) :: tail_w !! final tail option real(wp) :: p !! resulting CDF value ! ==== Instructions ! ---- handle input ! assume loc = 0, overwrite if specified loc_w = 0.0_wp if (present(loc)) loc_w = loc ! assume scale = 1, overwrite if specified scale_w = 1.0_wp if (present(scale)) scale_w = scale ! assume left-tailed, overwrite if specified tail_w = "left" if (present(tail)) tail_w = tail ! check if scale value is valid if (scale_w .le. 0.0_wp) then ! write error message and assign sentinel value if invalid call s_err_print(fsml_error(1)) p = c_sentinel_r return endif ! check if degrees of freedom value is valid if (df .le. 0.0_wp) then ! write error message and assign sentinel value if invalid call s_err_print(fsml_error(1)) p = c_sentinel_r return endif ! check if tail options are valid if (tail_w .ne. "left" .and. tail_w .ne. "right" .and. & &tail_w .ne. "two" .and. tail_w .ne. "confidence") then ! write error message and assign sentinel value if invalid call s_err_print(fsml_error(2)) p = c_sentinel_r return endif ! ----compute CDF ! call pure function to calculate probability integral p = f_dst_chi2_cdf_core(x, df, loc_w, scale_w, tail_w) end function f_dst_chi2_cdf