PLUTO
4.4-patch2
|
Compute the right hand side of the relativistic hydro (RHD) equations in primitive form. More...
#include "pluto.h"
Functions | |
void | PrimRHS (double *w, double *dw, double cs2, double h, double *Adw) |
void | PrimSource (const State *state, double **src, int beg, int end, Grid *grid) |
Implements the right hand side of the quasi-linear form of the relativistic hydro equations. In 1D this may be written as
where is the matrix of the primitive form of the equations,
is the source term.
Reference:
The function PrimRHS() implements the first term while PrimSource() implements the source term part.
void PrimRHS | ( | double * | w, |
double * | dw, | ||
double | cs2, | ||
double | h, | ||
double * | Adw | ||
) |
Compute the matrix-vector multiplication where
A
is the matrix of the quasi-linear form of the RHD equations.
[in] | w | vector of primitive variables; |
[in] | dw | limited (linear) slopes; |
[in] | cs2 | local sound speed; |
[in] | h | local enthalpy; |
[out] | Adw | matrix-vector product. |
Compute source terms of the RHD equations in primitive variables.
The rationale for choosing during which sweep a particular source term has to be incorporated should match the same criterion used during the conservative update. For instance, in polar or cylindrical coordinates, curvilinear source terms are included during the radial sweep only.