Go to the source code of this file.
double** GetBackgroundField |
( |
int |
beg, |
|
|
int |
end, |
|
|
int |
where, |
|
|
Grid * |
grid |
|
) |
| |
Definition at line 5 of file bckgrnd_field.c.
16 static double **bck_fld;
22 #if (TIME_STEPPING != RK2) && (TIME_STEPPING != RK3)
23 print1 (
"! Background field splitting works with RK integrators ONLY \n");
26 #if DIVB_CONTROL == EIGHT_WAVES
27 print1 (
"! Background field splitting works with CT or GLM ONLY \n");
29 #elif DIVB_CONTROL == CONSTRAINED_TRANSPORT
30 #if (CT_EMF_AVERAGE != ARITHMETIC) && (CT_EMF_AVERAGE != UCT_HLL)
31 print1 (
"! Background field splitting works with ARITHMETIC or");
32 print1 (
" UCT_HLL averages only\n");
37 if (bck_fld == NULL) {
51 for (i = beg; i <= end; i++){
56 for (i = beg; i <= end; i++){
61 for (i = beg; i <= end; i++){
void print1(const char *fmt,...)
void BackgroundField(double x1, double x2, double x3, double *B0)
int g_i
x1 grid index when sweeping along the x2 or x3 direction.
int g_dir
Specifies the current sweep or direction of integration.
int g_j
x2 grid index when sweeping along the x1 or x3 direction.
int g_k
x3 grid index when sweeping along the x1 or x2 direction.
long int NMAX_POINT
Maximum number of points among the three directions, boundaries excluded.
#define ARRAY_2D(nx, ny, type)
#define QUIT_PLUTO(e_code)