8 double BAROTROPIC_PR(
double);
10 double BodyForcePotential(
double,
double,
double);
11 void BodyForceVector(
double *,
double *,
double,
double,
double);
13 void ChangeDumpVar ();
14 void CheckConsStates (
double **,
double **,
double **,
int,
int);
15 void CheckPrimStates (
double **,
double **,
double **,
int,
int);
16 int CheckNaN (
double **,
int,
int,
int);
18 void ComputeUserVar (
const Data *,
Grid *);
23 void CreateImage (
char *);
27 #ifdef FINITE_DIFFERENCE
29 Reconstruct MP5_Reconstruct, PPM_Reconstruct, LIMO3_Reconstruct,
30 WENOZ_Reconstruct, WENO3_Reconstruct;
34 void FindShock (
const Data *,
Grid *);
39 Image *GetImage (
char *);
40 double *GetInverse_dl (
const Grid *);
44 void Init (
double *,
double,
double,
double);
57 double Median (
double a,
double b,
double c);
61 void ParseCmdLineArgs (
int,
char *argv[],
char *, Cmd_Line *);
63 char *ParGet (
const char *,
int );
64 int ParQuery (
const char *);
65 void PrimToChar (
double **,
double *,
double *);
71 void SetColorMap (
unsigned char *,
unsigned char *,
unsigned char *,
char *);
72 void SetDefaultVarNames(
Output *);
74 Riemann_Solver *SetSolver (
const char *);
78 int Setup (Input *, Cmd_Line *,
char *);
81 void SoundSpeed2 (
double **,
double *,
double *,
int,
int,
int,
Grid *);
89 void VISC_FLUX (
const State_1D *,
int beg,
int end,
Grid *);
96 char **
Array2D (
int,
int,
size_t);
97 char ***
Array3D (
int,
int,
int,
size_t);
98 char ****
Array4D (
int,
int,
int,
int,
size_t);
100 double ***
ArrayBox(
long int,
long int,
long int,
long int,
long int,
long int);
102 double ***
ArrayBoxMap (
int,
int,
int,
int,
int,
int,
double *);
103 double ***
ArrayMap (
int,
int,
int,
double *);
105 unsigned char ***ArrayCharMap (
int,
int,
int,
unsigned char *);
117 void FreeArrayCharMap(
unsigned char ***);
119 #define ARRAY_1D(nx,type) (type *)Array1D(nx,sizeof(type))
120 #define ARRAY_2D(nx,ny,type) (type **)Array2D(nx,ny,sizeof(type))
121 #define ARRAY_3D(nx,ny,nz,type) (type ***)Array3D(nx,ny,nz,sizeof(type))
122 #define ARRAY_4D(nx,ny,nz,nv,type) (type ****)Array4D(nx,ny,nz,nv,sizeof(type))
138 void print (
const char *fmt, ...);
139 void print1 (
const char *fmt, ...);
141 void Show (
double **,
int);
164 void Async_EndWriteData (Input *ini);
165 void Async_BegWriteData (
const Data *d,
Output *output,
Grid *grid);
185 void WritePPM (
double ***,
char *,
char *,
Grid *);
186 void WritePNG (
double ***,
char *,
char *,
Grid *);
194 void Jacobian (
double *v,
double *rhs,
double **dfdy);
196 #if COOLING == POWER_LAW
197 void PowerLawCooling (Data_Arr,
double,
Time_Step *,
Grid *);
200 double SolveODE_CK45 (
double *,
double *,
double *,
double,
double);
201 double SolveODE_RKF23 (
double *,
double *,
double *,
double,
double);
202 double SolveODE_RKF12 (
double *,
double *,
double *,
double,
double);
203 double SolveODE_ROS34 (
double *,
double *,
double *,
double,
double);
204 double SolveODE_RK4 (
double *,
double *,
double *,
double);
205 double SolveODE_RK2 (
double *,
double *,
double *,
double);
206 double SolveODE_Euler (
double *,
double *,
double *,
double);
215 int LUDecomp (
double **a,
int n,
int *indx,
double *d);
216 void LUBackSubst (
double **a,
int n,
int *indx,
double b[]);
218 double Length_1 (
int i,
int j,
int k,
Grid *);
219 double Length_2 (
int i,
int j,
int k,
Grid *);
220 double Length_3 (
int i,
int j,
int k,
Grid *);
222 #if UPDATE_VECTOR_POTENTIAL == YES
223 void VectorPotentialUpdate (
const Data *d,
const void *vp,
228 void CYLSOLVE (
const Data *d,
double dt,
Grid *grid);
229 void CONS_CYLSOLVE (Data_Arr,
double dt,
Grid *grid);
233 void TransformGrid (
Grid *grid);
234 double GridMap(
double,
int);