para_and_paraenohybrid.h

Go to the documentation of this file.

//
// PARA settings:
//

#define PMC 100
// PMC: originally used in PPM paper -- Xiaoyue says fails on nonlinear tests since has no limiter -- but they used it, right?

#define MINM_STEEPENER 101

#define PARALINE2EXTREME 1
#define PARALINE2LIM MC

#if(WHICHPARA==PARA2)
#define PARA2LIM VANL




#elif(WHICHPARA==PARA3)

#define PARA2LIM MC
#elif(WHICHPARA==PARA4)

#define PARA2LIM MC

#endif

#define JONPARASMOOTH 0 // default 0

//
//
// HYBRID PPM SETTINGS (only used when lim=WENO5*)
//


// whether to use para with WENO
// CHANGINGMARK: Not correct with MPI -- e.g. produces monopoles for any FLUXB
//#define PARAMODWENO 1   // default: 1 // now global variable

#define FULLHYBRID 1   // default: 1

#define DOPPMCONTACTSTEEPMODWENO 0    // default: 0

#define DOPPMSTEEPVARTYPEMODWENO 0    // default: 0
#define DOPPMREDUCEMODWENO 1    // default: 0 (stay speculative) -- well, causes problems.  In general strong shock will be too oscillatory so need to keep and hope that shock indicator only activates in strong shocks

#define USEMCFORLOWERORDERWENO 0   // default: 0










//
// PPM settings for para4gen() used by *all* para methods (lim=PARA,PARAFLAT,PARALINE)
//

#define JONPARAREDUCE 0  // default: 0

#define JONPARASTEEP 0  // default: 0

#define AVGINPUT 0  // default: 0

#define WHICHLIMITERTOUSEFORLR 0  // default: 0

#define PARAGENDQALLOWEXTREMUM 0  // default: 0 // bit too speculative in general -- and with new JONPARASMOOTH==1 not necessary really



//
// PPM SETTINGS for checkparamonotonicity()
//
 
#define NONMONOLIM 0  // default: 0




//
// FOR PARAPL (i.e. lim=PARAFLAT) only
//

#define DOPPMREDUCE 1 // default: 1

#define DOPPMCONTACTSTEEP 0 // default: 0


//
// PPM SETTINGS for ftilde() (lim=PARAFLAT and PARALINE only)
//

//#define SP0 0.65 // for strong shocks
#define SP0 0.75 // FLASH

#define SPSMOOTH 0.5  // default: 0.5

#define SV0 (0.0)

//
//
// PPM SETTINGS for parasteep() (only for lim=PARAFLAT and PARALINE)
//

#define DOPPMSTEEPVARTYPE 0   // default: 0