PLUTO
4.4-patch2
|
Compute the right hand side of the primitive HD/MHD equations. More...
#include "pluto.h"
Functions | |
void | PrimRHS (double *v, double *dv, double cs2, double h, double *Adv) |
void | PrimSource (const State *gas, double **src, int beg, int end, Grid *grid) |
Implements the right hand side of the quasi-linear form of the MHD 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 * | v, |
double * | dv, | ||
double | cs2, | ||
double | h, | ||
double * | Adv | ||
) |
Compute the matrix-vector multiplication where A is the matrix of the quasi-linear form of the MHD equations.
References
[in] | v | vector of primitive variables |
[in] | dv | limited (linear) slopes |
[in] | cs2 | local sound speed |
[in] | h | local enthalpy |
[out] | AdV | matrix-vector product |
Compute source terms of the MHD equations in primitive variables. These include:
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.
[in] | gas | pointer to a Sweep structure |
[out] | src | array of source terms |
[in] | beg | initial index of computation |
[in] | end | final index of computation |
[in] | grid | pointer to a Grid structure |
where