Right-eigenvector matrix K and its inverse for the flux Jacobian in
direction dir (1 = x: eigenvalues u-c,u,u,u+c; 2 = y: v-c,v,v,v+c).
Columns of K, conserved variables Q=(rho,rho u,rho v,E):
x: [1, u-c, v, H-u c] [1, u, v, q2] [0, 0, 1, v] [1, u+c, v, H+u c]
y: [1, u, v-c, H-v c] [1, u, v, q2] [0, 1, 0, u] [1, u, v+c, H+v c]
where H=(E+p)/rho, q2=0.5*(u^2+v^2). K^{-1} via LAPACK getrf/getri.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(in) | :: | q(neq2d) | |||
| integer, | intent(in) | :: | dir | |||
| real(kind=wp), | intent(out) | :: | r_mat(neq2d,neq2d) | |||
| real(kind=wp), | intent(out) | :: | r_inv(neq2d,neq2d) | |||
| real(kind=wp), | intent(in) | :: | gam |