HARM
harm and utilities
 All Data Structures Files Functions Variables Typedefs Macros Pages
Functions
wavespeeds.funcdeclare.h File Reference

function declrations for global use of functions in wavespeeds.c More...

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions

int get_global_wavespeeds (int dir, struct of_geom *ptrgeom, FTYPE *pr, FTYPE *wspeedtemp, FTYPE *wspeedradtemp, FTYPE *wspeedrad2temp)
 store wavespeeds somewhere More...
 
int get_global_wavespeeds_full (int dir, int is, int ie, int js, int je, int ks, int ke, int idel, int jdel, int kdel, FTYPE(*prim)[NSTORE2][NSTORE3][NPR], FTYPE(*wspeed)[COMPDIM][NUMCS][NSTORE1][NSTORE2][NSTORE3])
 complete storage of wave speeds per dimension only called for STOREWAVESPEEDS==1 More...
 
int global_vchar (FTYPE(*pointspeed)[NSTORE1][NSTORE2][NSTORE3][NUMCS], int dir, int is, int ie, int js, int je, int ks, int ke, int idel, int jdel, int kdel, FTYPE(*wspeed)[COMPDIM][NUMCS][NSTORE1][NSTORE2][NSTORE3])
 GODMARK: something wrong with comparing multiple velocities since grid/metric changes in space (i.e. More...
 
int get_wavespeeds (int dir, struct of_geom *ptrgeom, FTYPE *p_l, FTYPE *p_r, FTYPE *U_l, FTYPE *U_r, FTYPE *F_l, FTYPE *F_r, struct of_state *state_l, struct of_state *state_r, FTYPE *cminmax_l, FTYPE *cminmax_r, FTYPE *cminmax, FTYPE *ctop, FTYPE *cminmaxrad_l, FTYPE *cminmaxrad_r, FTYPE *cminmaxrad, FTYPE *ctoprad, FTYPE *cminmaxrad2_l, FTYPE *cminmaxrad2_r, FTYPE *cminmaxrad2, FTYPE *ctoprad2)
 get wave speeds for flux calculation More...
 
int vchar (FTYPE *pr, struct of_state *q, int dir, struct of_geom *geom, FTYPE *cmax, FTYPE *cmin, int *ignorecourant)
 Charles believed that near polar axis one should set vmax=vmin=0, but not right since makes flux completely non-diffusive (leads to numerical instabilities) More...
 
FTYPE chk_disp (FTYPE v)
 
void make_co_to_comov (FTYPE *ucon, FTYPE(*ecov)[NDIM], FTYPE(*econ)[NDIM])
 
void transform (FTYPE *vec, FTYPE(*t)[NDIM])
 
void coeff_set (FTYPE rho, FTYPE u)
 
int vchar_all (FTYPE *pr, struct of_state *q, int dir, struct of_geom *geom, FTYPE *cmax, FTYPE *cmin, int *ignorecourant)
 get wave speeds for flux calculation and for dt calculation More...
 
int vchar_each (FTYPE *pr, struct of_state *q, int dir, struct of_geom *geom, FTYPE *cmaxmhd, FTYPE *cminmhd, FTYPE *cmaxrad, FTYPE *cminrad, FTYPE *cmaxrad2, FTYPE *cminrad2, int *ignorecourant)
 

Detailed Description

function declrations for global use of functions in wavespeeds.c

Definition in file wavespeeds.funcdeclare.h.

Function Documentation

FTYPE chk_disp ( FTYPE  v)
void coeff_set ( FTYPE  rho,
FTYPE  u 
)
int get_global_wavespeeds ( int  dir,
struct of_geom *  ptrgeom,
FTYPE *  pr,
FTYPE *  wspeedtemp,
FTYPE *  wspeedradtemp,
FTYPE *  wspeedrad2temp 
)

store wavespeeds somewhere

Definition at line 369 of file wavespeeds.c.

Here is the call graph for this function:

Here is the caller graph for this function:

int get_global_wavespeeds_full ( int  dir,
int  is,
int  ie,
int  js,
int  je,
int  ks,
int  ke,
int  idel,
int  jdel,
int  kdel,
FTYPE(*)  prim[NSTORE2][NSTORE3][NPR],
FTYPE(*)  wspeed[COMPDIM][NUMCS][NSTORE1][NSTORE2][NSTORE3] 
)

complete storage of wave speeds per dimension only called for STOREWAVESPEEDS==1

Definition at line 322 of file wavespeeds.c.

Here is the call graph for this function:

Here is the caller graph for this function:

int get_wavespeeds ( int  dir,
struct of_geom *  ptrgeom,
FTYPE *  p_l,
FTYPE *  p_r,
FTYPE *  U_l,
FTYPE *  U_r,
FTYPE *  F_l,
FTYPE *  F_r,
struct of_state *  state_l,
struct of_state *  state_r,
FTYPE *  cminmax_l,
FTYPE *  cminmax_r,
FTYPE *  cminmax,
FTYPE *  ctop,
FTYPE *  cminmaxrad_l,
FTYPE *  cminmaxrad_r,
FTYPE *  cminmaxrad,
FTYPE *  ctoprad,
FTYPE *  cminmaxrad2_l,
FTYPE *  cminmaxrad2_r,
FTYPE *  cminmaxrad2,
FTYPE *  ctoprad2 
)

get wave speeds for flux calculation

Definition at line 50 of file wavespeeds.c.

Here is the call graph for this function:

Here is the caller graph for this function:

int global_vchar ( FTYPE(*)  pointspeed[NSTORE1][NSTORE2][NSTORE3][NUMCS],
int  dir,
int  is,
int  ie,
int  js,
int  je,
int  ks,
int  ke,
int  idel,
int  jdel,
int  kdel,
FTYPE(*)  wspeed[COMPDIM][NUMCS][NSTORE1][NSTORE2][NSTORE3] 
)

GODMARK: something wrong with comparing multiple velocities since grid/metric changes in space (i.e.

v=dx/dt means something different at each grid point) GODMARK: assumes boundary zones exist (flux method of bounding won't work) – have to apply extra limits on values (i,j,k) used here

defines an effective maximum wave speed centered on the cell interface (FACE)

might choose wavespeeds that correspond to interpolation stencil, which to first approximation is a symmetric stencil of size interporder[reallim]

Definition at line 394 of file wavespeeds.c.

Here is the caller graph for this function:

void make_co_to_comov ( FTYPE *  ucon,
FTYPE(*)  ecov[NDIM],
FTYPE(*)  econ[NDIM] 
)
void transform ( FTYPE *  vec,
FTYPE(*)  t[NDIM] 
)
int vchar ( FTYPE *  pr,
struct of_state *  q,
int  dir,
struct of_geom *  geom,
FTYPE *  cmax,
FTYPE *  cmin,
int *  ignorecourant 
)

Charles believed that near polar axis one should set vmax=vmin=0, but not right since makes flux completely non-diffusive (leads to numerical instabilities)

Definition at line 34 of file vchar.c.

Here is the call graph for this function:

Here is the caller graph for this function:

int vchar_all ( FTYPE *  pr,
struct of_state *  q,
int  dir,
struct of_geom *  geom,
FTYPE *  vmaxall,
FTYPE *  vminall,
int *  ignorecourant 
)

get wave speeds for flux calculation and for dt calculation

Definition at line 16 of file wavespeeds.c.

Here is the call graph for this function:

Here is the caller graph for this function:

int vchar_each ( FTYPE *  pr,
struct of_state *  q,
int  dir,
struct of_geom *  geom,
FTYPE *  cmaxmhd,
FTYPE *  cminmhd,
FTYPE *  cmaxrad,
FTYPE *  cminrad,
FTYPE *  cmaxrad2,
FTYPE *  cminrad2,
int *  ignorecourant 
)

Definition at line 31 of file wavespeeds.c.

Here is the call graph for this function:

Here is the caller graph for this function: