linear_weno5_reconstruct Subroutine

public pure subroutine linear_weno5_reconstruct(f_stencil, f_hat)

5th-order linear WENO reconstruction on a 5-point stencil.

Equivalent to WENO5-JS in the smooth limit: applies the optimal weights d0=1/10, d1=3/5, d2=3/10 to the three 3rd-order substencil approximations without computing beta-based nonlinear corrections. The resulting combined coefficients are:

f_hat = (1/30)v1 - (13/60)v2 + (47/60)v3 + (27/60)v4 - (3/60)*v5

See Jiang & Shu (1996), Eq. (19) combined with substencil Eq. (16a-16c).

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: f_stencil(:,:)
real(kind=wp), intent(out) :: f_hat(:)