Inversion scheme for RMHD using a pressure fix. More...
#include "pluto.h"
Functions | |
static double | VelocitySquareFunc (double, Map_param *) |
int | PressureFix (Map_param *par) |
Fix p to a small value, solve for the square of velocity by using secant 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 PressureFix | ( | Map_param * | par | ) |
Fix p to a small value, solve for the square of velocity by using secant algorithm applied to Eq (9) of Mignone, Plewa & Bodo (2005). 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.
|
static |
Implement Eq (A3) of Mignone & McKinney (2007).