|
HARM
harm and utilities
|
|
HARM
harm and utilities
|
U->P inversion definitions for all utoprim files. More...
#include "decs.h"
Go to the source code of this file.
Macros | |
| #define | GAMMA (gamideal) /* Adiabatic index used for the state equation */ |
| #define | CHANGEDTOOLDER 0 |
| #define | CRAZYDEBUG 0 |
| #define | OPTIMIZED 1 |
| #define | NEWCONVERGE 1 |
| #define | DOHISTOGRAM 0 |
| #define | GAMMASQCHECKRESID (1e2) |
| at what ^2 does the inversion change character to become very accurate? More... | |
| #define | CYCLESTOP 200 |
| #define | NUMCYCLES 200 |
| #define | ITERDAMPSTART 10 |
| Note that if using -pc64 -mp that error in inversion seems to be limited for doubles to 1E-11 instead of 1E-15. More... | |
| #define | MAX_NEWT_ITER 30 /* Max. # of Newton-Raphson iterations for find_root_2D(); */ |
| #define | NEWT_TOL (1E5*NUMEPSILON) /* Min. of tolerance allowed for Newton-Raphson iterations */ |
| #define | NEWT_TOL_ULTRAREL (1E1*NUMEPSILON) /* Min. of tolerance allowed for Newton-Raphson iterations */ |
| #define | MIN_NEWT_TOL (1E-2) |
| #define | MINERROREXPECTED (NUMEPSILON*10.0) |
| minimum relative error expected in used to allow ultrarelativistic results to have larger relative errors since limited by machine precision More... | |
| #define | EXTRA_NEWT_ITER 1 |
| always do 1 extra to avoid error estimate issues as above More... | |
| #define | EXTRA_NEWT_ITER_ULTRAREL 2 |
| #define | CYCLE_BREAK_PERIOD 1000 /* change newton step by random factor every this number of newton iterations*/ |
| #define | CHECK_FOR_STALL 0 /* check for stationary newton stepping */ |
| #define | NEWT_TOL2 (NUMEPSILON) /* TOL of new DBL gnr2 method */ |
| #define | MIN_NEWT_TOL2 (1E4*NUMEPSILON) /* TOL of new DBL gnr2 method */ |
| #define | SCALEMAX 1.0e2 /* Max. value of the factor used to scale the Newton step */ |
| #define | TOL_LINE_STEP NEWT_TOL /* Minimum value of Max(dx/x) in line search algorithm */ |
| #define | GRADMIN (1E4*NUMEPSILON) /* Magnitude of gradient below which we say we are at a local min. */ |
| #define | NEWT_FUNC_TOL 1.0e-5 /* Max. ratio of the final and initial resid magnitudes to be considered converged */ |
| #define | GAMMASQ_TOO_BIG (GAMMAFAIL*GAMMAFAIL) |
| #define | UTSQ_TOO_BIG ((GAMMAFAIL-1.0)*(GAMMAFAIL-1.0)) |
| #define | UT_TOO_BIG (GAMMAFAIL-1.0) |
| #define | VSQ_TOO_BIG ((GAMMASQ_TOO_BIG-1.0)/(GAMMASQ_TOO_BIG)) |
| #define | FAIL_VAL BIG |
| #define | VSQNEGCHECK 2 |
| seems code can recover gracefully when this happens, so don't worry about it which check 0: no check 1: original check 2: new check allowing somewhat negative result 3: just set to 0 if negative More... | |
| #define | UTSQNEGLIMIT (-0.9*UTSQ_TOO_BIG) |
| For <=-1 the equations become complex. More... | |
| #define | MAXNEGUTSQ (fabs(UTSQNEGLIMIT)) |
| positive value of max negative utsq More... | |
| #define | MAXNEGVSQ (0.9*VSQ_TOO_BIG) |
| #define | RHO 0 |
| #define | UU 1 |
| #define | UTCON1 2 |
| #define | UTCON2 3 |
| #define | UTCON3 4 |
| #define | BCON1 5 |
| #define | BCON2 6 |
| #define | BCON3 7 |
| #define | QCOV0 1 |
| #define | QCOV1 2 |
| #define | QCOV2 3 |
| #define | QCOV3 4 |
| #define | MYMAX(a, b) ( ((a) > (b)) ? (a) : (b) ) |
U->P inversion definitions for all utoprim files.
Definition in file u2p_defs.h.
| #define BCON1 5 |
Definition at line 171 of file u2p_defs.h.
| #define BCON2 6 |
Definition at line 172 of file u2p_defs.h.
| #define BCON3 7 |
Definition at line 173 of file u2p_defs.h.
| #define CHANGEDTOOLDER 0 |
Definition at line 24 of file u2p_defs.h.
| #define CHECK_FOR_STALL 0 /* check for stationary newton stepping */ |
Definition at line 104 of file u2p_defs.h.
| #define CRAZYDEBUG 0 |
Definition at line 26 of file u2p_defs.h.
| #define CYCLE_BREAK_PERIOD 1000 /* change newton step by random factor every this number of newton iterations*/ |
Definition at line 103 of file u2p_defs.h.
| #define CYCLESTOP 200 |
Definition at line 43 of file u2p_defs.h.
| #define DOHISTOGRAM 0 |
Definition at line 34 of file u2p_defs.h.
| #define EXTRA_NEWT_ITER 1 |
always do 1 extra to avoid error estimate issues as above
Definition at line 99 of file u2p_defs.h.
| #define EXTRA_NEWT_ITER_ULTRAREL 2 |
Definition at line 100 of file u2p_defs.h.
| #define FAIL_VAL BIG |
Definition at line 137 of file u2p_defs.h.
Definition at line 9 of file u2p_defs.h.
Definition at line 125 of file u2p_defs.h.
| #define GAMMASQCHECKRESID (1e2) |
at what ^2 does the inversion change character to become very accurate?
Definition at line 37 of file u2p_defs.h.
| #define GRADMIN (1E4*NUMEPSILON) /* Magnitude of gradient below which we say we are at a local min. */ |
Definition at line 113 of file u2p_defs.h.
| #define ITERDAMPSTART 10 |
Note that if using -pc64 -mp that error in inversion seems to be limited for doubles to 1E-11 instead of 1E-15.
Definition at line 64 of file u2p_defs.h.
| #define MAX_NEWT_ITER 30 /* Max. # of Newton-Raphson iterations for find_root_2D(); */ |
Definition at line 65 of file u2p_defs.h.
| #define MAXNEGUTSQ (fabs(UTSQNEGLIMIT)) |
positive value of max negative utsq
Definition at line 152 of file u2p_defs.h.
| #define MAXNEGVSQ (0.9*VSQ_TOO_BIG) |
Definition at line 154 of file u2p_defs.h.
| #define MIN_NEWT_TOL (1E-2) |
Definition at line 71 of file u2p_defs.h.
| #define MIN_NEWT_TOL2 (1E4*NUMEPSILON) /* TOL of new DBL gnr2 method */ |
Definition at line 108 of file u2p_defs.h.
| #define MINERROREXPECTED (NUMEPSILON*10.0) |
minimum relative error expected in used to allow ultrarelativistic results to have larger relative errors since limited by machine precision
Definition at line 83 of file u2p_defs.h.
Definition at line 182 of file u2p_defs.h.
| #define NEWCONVERGE 1 |
Definition at line 32 of file u2p_defs.h.
| #define NEWT_FUNC_TOL 1.0e-5 /* Max. ratio of the final and initial resid magnitudes to be considered converged */ |
Definition at line 114 of file u2p_defs.h.
| #define NEWT_TOL (1E5*NUMEPSILON) /* Min. of tolerance allowed for Newton-Raphson iterations */ |
Definition at line 66 of file u2p_defs.h.
| #define NEWT_TOL2 (NUMEPSILON) /* TOL of new DBL gnr2 method */ |
Definition at line 107 of file u2p_defs.h.
| #define NEWT_TOL_ULTRAREL (1E1*NUMEPSILON) /* Min. of tolerance allowed for Newton-Raphson iterations */ |
Definition at line 67 of file u2p_defs.h.
| #define NUMCYCLES 200 |
Definition at line 44 of file u2p_defs.h.
| #define OPTIMIZED 1 |
Definition at line 30 of file u2p_defs.h.
| #define QCOV0 1 |
Definition at line 176 of file u2p_defs.h.
| #define QCOV1 2 |
Definition at line 177 of file u2p_defs.h.
| #define QCOV2 3 |
Definition at line 178 of file u2p_defs.h.
| #define QCOV3 4 |
Definition at line 179 of file u2p_defs.h.
| #define RHO 0 |
Definition at line 161 of file u2p_defs.h.
| #define SCALEMAX 1.0e2 /* Max. value of the factor used to scale the Newton step */ |
Definition at line 111 of file u2p_defs.h.
Definition at line 112 of file u2p_defs.h.
| #define UT_TOO_BIG (GAMMAFAIL-1.0) |
Definition at line 128 of file u2p_defs.h.
| #define UTCON1 2 |
Definition at line 168 of file u2p_defs.h.
| #define UTCON2 3 |
Definition at line 169 of file u2p_defs.h.
| #define UTCON3 4 |
Definition at line 170 of file u2p_defs.h.
Definition at line 127 of file u2p_defs.h.
| #define UTSQNEGLIMIT (-0.9*UTSQ_TOO_BIG) |
For <=-1 the equations become complex.
Definition at line 148 of file u2p_defs.h.
| #define UU 1 |
Definition at line 165 of file u2p_defs.h.
| #define VSQ_TOO_BIG ((GAMMASQ_TOO_BIG-1.0)/(GAMMASQ_TOO_BIG)) |
Definition at line 134 of file u2p_defs.h.
| #define VSQNEGCHECK 2 |
seems code can recover gracefully when this happens, so don't worry about it which check 0: no check 1: original check 2: new check allowing somewhat negative result 3: just set to 0 if negative
Definition at line 145 of file u2p_defs.h.
1.8.3.1