Perform MP5 reconstruction on the supplied stencil.
Applies MP5 reconstruction component-wise to a 5-point stencil. The unlimited 5th-order approximation is identical to the optimal (linear) WENO5 reconstruction. A two-stage monotonicity check and curvature-based limiter enforce non-oscillatory behaviour near discontinuities while preserving 5th-order accuracy in smooth regions.
Stencil column mapping (half_stencil=3): col 1 = v_{i-2}, col 2 = v_{i-1}, col 3 = v_i (centre), col 4 = v_{i+1}, col 5 = v_{i+2}.
See Suresh & Huynh (1997), Sec. 2.3, Eqs. (2.27)-(2.29).
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(in) | :: | f_stencil(:,:) | |||
| real(kind=wp), | intent(out) | :: | f_hat(:) |