reconstruct_with_proj Subroutine

private subroutine reconstruct_with_proj(recon, stencil, r_mat, r_inv, proj, flux_out)

Reconstruct a face flux or state with optional characteristic projection.

If proj is .true.: 1. Transform each stencil column to characteristic space: K^{-1} * f. 2. Reconstruct the face value in characteristic space. 3. Map back to physical space: K * flux_out. If proj is .false.: Call recon directly on the physical-space stencil.

See course notes, Sec. "Generalization to Systems", Steps 3-4.

Arguments

Type IntentOptional Attributes Name
procedure(reconstructor_iface) :: recon

Active reconstruction procedure pointer.

real(kind=wp), intent(in) :: stencil(:,:)
real(kind=wp), intent(in) :: r_mat(:,:)
real(kind=wp), intent(in) :: r_inv(:,:)
logical, intent(in) :: proj

Enable characteristic-space projection.

real(kind=wp), intent(out) :: flux_out(:)