Convert between primitive and conservative variables. More...
#include "pluto.h"
Functions | |
void | PrimToCons (double **uprim, double **ucons, int ibeg, int iend) |
int | ConsToPrim (double **ucons, double **uprim, int ibeg, int iend, unsigned char *flag) |
int | EnergySolve (double *u, double *v) |
int | PressureFix (double *u, double *v) |
int | EntropySolve (double *u, double *v) |
The PrimToCons() converts an array of primitive quantities to an array of conservative variables for the RMHD equations.
The ConsToPrim() converts an array of conservative quantities to an array of primitive quantities. During the conversion, pressure is normally recovered from total energy unless zone has been tagged with FLAG_ENTROPY. In this case we recover pressure from conserved entropy:
if (FLAG_ENTROPY is TRUE) --> p = p(S) else --> p = p(E)
int ConsToPrim | ( | double ** | ucons, |
double ** | uprim, | ||
int | ibeg, | ||
int | iend, | ||
unsigned char * | flag | ||
) |
Convert from conservative to primitive variables.
[in] | ucons | array of conservative variables |
[out] | uprim | array of primitive variables |
[in] | beg | starting index of computation |
[in] | end | final index of computation |
[out] | flag | array of flags tagging zones where conversion went wrong. |
int EnergySolve | ( | double * | u, |
double * | v | ||
) |
Try to recover gas pressure from conserved variables {D, m, E, B}. Return (0) is successful, (1) otherwise.
int EntropySolve | ( | double * | u, |
double * | v | ||
) |
int PressureFix | ( | double * | u, |
double * | v | ||
) |
Fix p to a small value, solve for the square of velocity. Then recompute W by solving the velocity-momentum relation and correct energy. Return 0 if succesful, 1 otherwise.
void PrimToCons | ( | double ** | uprim, |
double ** | ucons, | ||
int | ibeg, | ||
int | iend | ||
) |
Convert primitive variables to conservative variables.
[in] | uprim | array of primitive variables |
[out] | ucons | array of conservative variables |
[in] | beg | starting index of computation |
[in] | end | final index of computation |