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

macros/definitions for interppoint.para.c and interpline.para.c More...

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

Go to the source code of this file.

Macros

#define PMC   100
 
#define MINM_STEEPENER   101
 
#define PARALINE2EXTREME   1
 whether to allow extremum in para that has monotonicity check at end More...
 
#define PARALINE2LIM   MC
 
#define PARA2LIM   VANL
 
#define JONPARASMOOTH   0
 jon tests show for PARA2: PARA2LIM = MC completely fails PARA2LIM = VANL best ok PARA2LIM = PMC kinda ok More...
 
#define FULLHYBRID   1
 whether to use para for only contacts (0) or full hybrid method based upon indicators (e.g. stiffindicator) More...
 
#define DOPPMCONTACTSTEEPMODWENO   0
 see interppoint.para.c on the meaning of these contact steepener is invalid procedure for general problems with gravity and such things (e.g. More...
 
#define DOPPMSTEEPVARTYPEMODWENO   0
 not sure if field method works well enough to use generally // not used? GODMARK More...
 
#define DOPPMREDUCEMODWENO   1
 
#define USEMCFORLOWERORDERWENO   0
 whether to use MC when WENO reduces (only replaces WENO part, not PARA part) Point is that WENO3 is much more diffusion than MC may introduce slight switching problem when MC goes sharp More...
 
#define JONPARAREDUCE   0
 
#define JONPARASTEEP   0
 actually helps keep para flat, but also stepeends in perfect limit of infinite drop in density Helps #2 (double rarefaction) reach low u/rho! (not significantly better) More...
 
#define AVGINPUT   0
 whether to assume input y's are cell average or point values More...
 
#define WHICHLIMITERTOUSEFORLR   0
 don't have to use PARA to get original l,r 0 = para 1 = PARA2LIM limiter directly (then acts like MCSTEEP) More...
 
#define PARAGENDQALLOWEXTREMUM   0
 Noticed with MCSTEEP that allowing extremum causes catastrophic problems with torus problem with field assume monocheck on, or else results in crazy behavior Even with PARA this (monotonic condition check) causes problems (e.g. More...
 
#define NONMONOLIM   0
 whether to reduce to DONOR or MONO method when parabola is not monotonic DONOR: 0 PARA2LIM: 1 (not robust) More...
 
#define DOPPMREDUCE   1
 DOPPMREDUCE probably needed if using contact steepener and other such things, but does result in more diffusion result that makes benefit of general PPM not good. More...
 
#define DOPPMCONTACTSTEEP   0
 whether to use PPM contact steepener (often results in bad behvior in general even if works very well for 1D LW2003 tests) More...
 
#define SP0   0.75
 
#define SPSMOOTH   0.5
 
#define SV0   (0.0)
 
#define DOPPMSTEEPVARTYPE   0
 0 == steepend only density 1 == steepend density and cross-fields (i.e. current sheets) More...
 

Detailed Description

macros/definitions for interppoint.para.c and interpline.para.c

Definition in file para_and_paraenohybrid.h.

Macro Definition Documentation

#define AVGINPUT   0

whether to assume input y's are cell average or point values

Definition at line 88 of file para_and_paraenohybrid.h.

#define DOPPMCONTACTSTEEP   0

whether to use PPM contact steepener (often results in bad behvior in general even if works very well for 1D LW2003 tests)

Definition at line 125 of file para_and_paraenohybrid.h.

#define DOPPMCONTACTSTEEPMODWENO   0

see interppoint.para.c on the meaning of these contact steepener is invalid procedure for general problems with gravity and such things (e.g.

can steepend density profiles due to force balance! – like disks!!!)

Definition at line 56 of file para_and_paraenohybrid.h.

#define DOPPMREDUCE   1

DOPPMREDUCE probably needed if using contact steepener and other such things, but does result in more diffusion result that makes benefit of general PPM not good.

Definition at line 123 of file para_and_paraenohybrid.h.

#define DOPPMREDUCEMODWENO   1

Definition at line 59 of file para_and_paraenohybrid.h.

#define DOPPMSTEEPVARTYPE   0

0 == steepend only density 1 == steepend density and cross-fields (i.e. current sheets)

Definition at line 149 of file para_and_paraenohybrid.h.

#define DOPPMSTEEPVARTYPEMODWENO   0

not sure if field method works well enough to use generally // not used? GODMARK

Definition at line 58 of file para_and_paraenohybrid.h.

#define FULLHYBRID   1

whether to use para for only contacts (0) or full hybrid method based upon indicators (e.g. stiffindicator)

Definition at line 52 of file para_and_paraenohybrid.h.

#define JONPARAREDUCE   0

Definition at line 81 of file para_and_paraenohybrid.h.

#define JONPARASMOOTH   0

jon tests show for PARA2: PARA2LIM = MC completely fails PARA2LIM = VANL best ok PARA2LIM = PMC kinda ok

whether to use internal parabola when -1,0,+1 dq's are central for density field, look at -2..+2 if possible (PARAFLAT or PARALINE) to ensure moving contact preserved in few points This leads to more stable and accurate evolution for relatively smooth flows, but leads to much more dissipation in turbulent flows like accretion disks.

Definition at line 37 of file para_and_paraenohybrid.h.

#define JONPARASTEEP   0

actually helps keep para flat, but also stepeends in perfect limit of infinite drop in density Helps #2 (double rarefaction) reach low u/rho! (not significantly better)

Definition at line 85 of file para_and_paraenohybrid.h.

#define MINM_STEEPENER   101

Definition at line 13 of file para_and_paraenohybrid.h.

#define NONMONOLIM   0

whether to reduce to DONOR or MONO method when parabola is not monotonic DONOR: 0 PARA2LIM: 1 (not robust)

Definition at line 111 of file para_and_paraenohybrid.h.

#define PARA2LIM   VANL

Definition at line 20 of file para_and_paraenohybrid.h.

#define PARAGENDQALLOWEXTREMUM   0

Noticed with MCSTEEP that allowing extremum causes catastrophic problems with torus problem with field assume monocheck on, or else results in crazy behavior Even with PARA this (monotonic condition check) causes problems (e.g.

horizon failures and polar blow-ups) With PARAFLAT or PARALINE, enough ddq's (second derivatives) to ensure really a parabola at extrema

Definition at line 98 of file para_and_paraenohybrid.h.

#define PARALINE2EXTREME   1

whether to allow extremum in para that has monotonicity check at end

Definition at line 16 of file para_and_paraenohybrid.h.

#define PARALINE2LIM   MC

Definition at line 17 of file para_and_paraenohybrid.h.

#define PMC   100

Definition at line 10 of file para_and_paraenohybrid.h.

#define SP0   0.75

Definition at line 135 of file para_and_paraenohybrid.h.

#define SPSMOOTH   0.5

Definition at line 137 of file para_and_paraenohybrid.h.

#define SV0   (0.0)

Definition at line 139 of file para_and_paraenohybrid.h.

#define USEMCFORLOWERORDERWENO   0

whether to use MC when WENO reduces (only replaces WENO part, not PARA part) Point is that WENO3 is much more diffusion than MC may introduce slight switching problem when MC goes sharp

Definition at line 64 of file para_and_paraenohybrid.h.

#define WHICHLIMITERTOUSEFORLR   0

don't have to use PARA to get original l,r 0 = para 1 = PARA2LIM limiter directly (then acts like MCSTEEP)

Definition at line 92 of file para_and_paraenohybrid.h.