compute_eigensystem Subroutine

public pure subroutine compute_eigensystem(q_avg, r_mat, r_inv, gam)

Compute the right eigenvector matrix K and its inverse K^{-1} of the Euler Jacobian at a given state.

The eigenvectors are used for the characteristic projection in the WENO reconstruction of systems. See course notes, Eq. (10) for eigenvalues, Eq. (11) for K, and Eq. (14) for K^{-1}.

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: q_avg(neq)
real(kind=wp), intent(out) :: r_mat(neq,neq)
real(kind=wp), intent(out) :: r_inv(neq,neq)
real(kind=wp), intent(in) :: gam

Ratio of specific heats