PLUTO
4.4-patch2
|
Inversion scheme for RMHD using total energy density. More...
#include "pluto.h"
Functions | |
int | RMHD_EnergySolve (Map_param *par) |
Try to recover gas pressure from conserved variables {D, m, E, B} using the algorithm outlined in Section A3 of Mignone & McKinney (2007). Specifically, we solve Eq. (A4) or (A6) (depending on the value of RMHD_REDUCED_ENERGY
) using a Newton-Raphson scheme. Here W = rho*h*lorentz^2, E, D, etc... have the same meaning as in the mentioned paper.
References
int RMHD_EnergySolve | ( | Map_param * | par | ) |
Solve f(W) = 0, where f(W) is Equation (A4) or (A6).