PLUTO
4.4-patch2
|
Inversion scheme for RMHD using a pressure fix. More...
#include "pluto.h"
Functions | |
static double | VelocitySquareFunc (double, void *) |
int | RMHD_PressureFix (Map_param *par) |
Fix p to a small value, solve for the square of velocity by using secant or bisection algorithm applied to Eq (A3). This step involved re-computing W at each step of the iteration. Once the root has been found, we recompute total energy E. Return 0 if succesful, 1 otherwise.
int RMHD_PressureFix | ( | Map_param * | par | ) |
|
static |
Implement Eq (A3) of Mignone & McKinney (2007).