PLUTO
|
Inversion scheme for RMHD using a pressure fix. More...
#include "pluto.h"
Go to the source code of this file.
Macros | |
#define | MAX_ITER 50 |
Functions | |
static double | VelocitySquareFunc (double, void *) |
int | PressureFix (Map_param *par) |
Inversion scheme for RMHD using a pressure fix.
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.
Definition in file pressure_fix.c.
#define MAX_ITER 50 |
Definition at line 20 of file pressure_fix.c.
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.
Definition at line 23 of file pressure_fix.c.
|
static |
Implement Eq (A3) of Mignone & McKinney (2007).
Definition at line 131 of file pressure_fix.c.