146 lines
2.7 KiB
C
146 lines
2.7 KiB
C
|
|
#ifndef MICRODEF_H
|
|
#define MICRODEF_H
|
|
|
|
#include "macrodef.fh"
|
|
|
|
// application parameters
|
|
|
|
#define SommerType 0
|
|
|
|
#define GaussInt
|
|
|
|
#define ABEtype 0
|
|
|
|
//#define With_AHF
|
|
#define Psi4type 0
|
|
|
|
//#define Point_Psi4
|
|
|
|
#define RPS 1
|
|
|
|
#define AGM 0
|
|
|
|
#define RPB 0
|
|
|
|
#define MAPBH 1
|
|
|
|
#define PSTR 0
|
|
|
|
#define REGLEV 0
|
|
|
|
//#define USE_GPU
|
|
|
|
//#define CHECKDETAIL
|
|
|
|
//#define FAKECHECK
|
|
|
|
//
|
|
// define SommerType
|
|
// sommerfeld boundary type
|
|
// 0: bam
|
|
// 1: shibata
|
|
//
|
|
// define GaussInt
|
|
// for Using Gauss-Legendre quadrature in theta direction
|
|
//
|
|
// define ABEtype
|
|
// 0: BSSN vacuum
|
|
// 1: coupled to scalar field
|
|
// 2: Z4c vacuum
|
|
// 3: coupled to Maxwell field
|
|
//
|
|
// define With_AHF
|
|
// using Apparent Horizon Finder
|
|
//
|
|
// define Psi4type
|
|
// Psi4 calculation method
|
|
// 0: EB method
|
|
// 1: 4-D method
|
|
//
|
|
// define Point_Psi4
|
|
// for Using point psi4 or not
|
|
//
|
|
// define RPS
|
|
// RestrictProlong in Step (0) or after Step (1)
|
|
//
|
|
// define AGM
|
|
// Enforce algebra constraint
|
|
// for every RK4 sub step: 0
|
|
// only when iter_count == 3: 1
|
|
// after routine Step: 2
|
|
//
|
|
// define RPB
|
|
// Restrict Prolong using BAM style 1 or old style 0
|
|
//
|
|
// define MAPBH
|
|
// 1: move Analysis out ot 4 sub steps and treat PBH with Euler method
|
|
//
|
|
// define PSTR
|
|
// parallel structure
|
|
// 0: level by level
|
|
// 1: considering all levels
|
|
// 2: as 1 but reverse the CPU order
|
|
// 3: Frank's scheme
|
|
//
|
|
// define REGLEV
|
|
// regrid for every level or for all levels at a time
|
|
// 0: for every level;
|
|
// 1: for all
|
|
//
|
|
// define USE_GPU
|
|
// use gpu or not
|
|
//
|
|
// define CHECKDETAIL
|
|
// use checkpoint for every process
|
|
//
|
|
// define FAKECHECK
|
|
// use FakeCheckPrepare to write CheckPoint
|
|
//
|
|
|
|
////================================================================
|
|
// some basic parameters for numerical calculation
|
|
////================================================================
|
|
|
|
#define dim 3
|
|
|
|
//#define Cell or Vertex in "macrodef.fh"
|
|
|
|
#define buffer_width 6
|
|
|
|
#define SC_width buffer_width
|
|
|
|
#define CS_width (2*buffer_width)
|
|
|
|
//
|
|
// define Cell or Vertex in "macrodef.fh"
|
|
//
|
|
// define buffer_width
|
|
// buffer point number for mesh refinement interface
|
|
//
|
|
// define SC_width buffer_width
|
|
// buffer point number shell-box interface, on shell
|
|
//
|
|
// define CS_width
|
|
// buffer point number shell-box interface, on box
|
|
//
|
|
|
|
#if(buffer_width < ghost_width)
|
|
# error we always assume buffer_width>ghost_width
|
|
#endif
|
|
|
|
#define PACK 1
|
|
#define UNPACK 2
|
|
|
|
#define Mymax(a,b) (((a) > (b)) ? (a) : (b))
|
|
#define Mymin(a,b) (((a) < (b)) ? (a) : (b))
|
|
|
|
#define feq(a,b,d) (fabs(a-b)<d)
|
|
#define flt(a,b,d) ((a-b)<d)
|
|
#define fgt(a,b,d) ((a-b)>d)
|
|
|
|
#define TINY 1e-10
|
|
|
|
#endif /* MICRODEF_H */
|
|
|