PLUTO  4.0
 All Data Structures Files Functions Variables Enumerations Macros Pages
Functions
pressure_fix.c File Reference

Inversion scheme for RMHD using a pressure fix. More...

#include "pluto.h"

Functions

static double VelocitySquareFunc (double, Map_param *)
 
int PressureFix (Map_param *par)
 

Detailed Description

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.

Authors
A. Mignone
C. Zanni
Date
Oct 5, 2012

Function Documentation

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.

double VelocitySquareFunc ( double  v2,
Map_param par 
)
static

Implement Eq (A3) of Mignone & McKinney (2007).