Core MUSCL reconstruction kernel (component-wise, left-biased stencil).
Each component k of f_hat is computed as: ΔL = s(k,2) - s(k,1), ΔR = s(k,3) - s(k,2) r = ΔR / ΔL (when |ΔL| > tiny; else f_hat(k) = s(k,2)) f_hat(k) = s(k,2) + 0.5 · lim_fn(r) · ΔL
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(in) | :: | f_stencil(:,:) | |||
| real(kind=wp), | intent(out) | :: | f_hat(:) | |||
| procedure(limiter_iface) | :: | lim_fn |
TVD limiter function phi(r) |