1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
|
/*@@
@file GRHydro_Minima.cc
@date Tue Aug 29 18:52:10 2006
@author
@desc
Sets up the scalars used for the atmosphere, after initial data.
@enddesc
@@*/
#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctk_Arguments.h"
//#include "Carpet/Carpet/src/carpet.hh"
#include "carpet.hh"
#ifdef HAVE_CARPET
using namespace Carpet;
#endif
#ifdef __cplusplus
extern "C" {
#endif
/* Scheduled functions */
void GRHydro_Rho_Minima_Setup_Final(CCTK_ARGUMENTS);
void GRHydro_Rho_Minima_Setup_Final_PUGH(CCTK_ARGUMENTS);
#ifdef __cplusplus
} /* extern "C" */
#endif
/*@@
@routine GRHydro_Rho_Minima_Setup_Final
@date Tue Aug 29 18:54:05 2006
@author Luca Baiotti
@desc
After initial data, set up the scalar GRHydro_rho_min used for checking the atmosphere.
This is computed taking into account the actual maximum rest-mass density on the grid.
Version for Carpet.
@enddesc
@calls
@calledby
@history
@endhistory
@@*/
void GRHydro_Rho_Minima_Setup_Final(CCTK_ARGUMENTS)
{
#ifdef HAVE_CARPET
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
CCTK_REAL max_rho;
static int flag = true;
if (flag) // Actually run this reduction only the first time the routine is called.
{
if (rho_abs_min > 0.0)
{
*GRHydro_rho_min = rho_abs_min;
}
else
{
// Go to global mode
BEGIN_GLOBAL_MODE (cctkGH) {
// Find the global maximum of rho
const int Reduction_Handle = CCTK_ReductionHandle ("maximum");
const CCTK_INT input_array_variable_indices= {CCTK_VarIndex("HydroBase::rho")};
const int ierr = CCTK_Reduce(cctkGH,
-1, // target processors; -1 -> all
Reduction_Handle,
1, // number of output variables
CCTK_VARIABLE_REAL,
&max_rho,
1, // number of variables to be reduced
input_array_variable_indices);
if (ierr != 0)
{
CCTK_WARN(0, "Failed to compute the global maximum of rho");
}
*GRHydro_rho_min = max_rho * rho_rel_min;
// Go back to local mode
} END_GLOBAL_MODE;
}
// After this has run once, set the flag so that this does not run again
flag = false;
} // end if (flag)
else
{
return;
}
#endif
//Debug stuff
// char warnline;
//CCTK_VInfo (CCTK_THORNSTRING, "STEP 2: compute rho max; rho min: %13.12e \n",*GRHydro_rho_min);
// printf(warnline,"STEP 2: compute rho max; rho min: %13.12e \n",*GRHydro_rho_min);
//CCTK_WARN(1,warnline);
// printf(warnline,"STEP 3: recompute ID with new atmosphere; rho min: ', GRHydro_rho_min, 'reflev: ',%i), GRHydro_rho_min, GRHydro_reflevel);
}
/*@@
@routine GRHydro_Rho_Minima_Setup_Final_PUGH
@date Tue Aug 29 18:57:42 2006
@author Luca Baiotti
@desc
As above, but for PUGH.
@enddesc
@calls
@calledby
@history
@endhistory
@@*/
void GRHydro_Rho_Minima_Setup_Final_PUGH(CCTK_ARGUMENTS)
{
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
CCTK_REAL max_rho;
if (rho_abs_min > 0.0)
{
*GRHydro_rho_min = rho_abs_min;
}
else
{
// Find the global maximum of rho
const int Reduction_Handle = CCTK_ReductionHandle("maximum");
const CCTK_INT input_array_variable_indices={CCTK_VarIndex("HydroBase::rho")};
const int ierr = CCTK_Reduce(cctkGH,
-1, // target processors; -1 -> all
Reduction_Handle,
1, // number of output variables
CCTK_VARIABLE_REAL,
&max_rho,
1, // number of variables to be reduced
input_array_variable_indices);
if (ierr != 0)
{
CCTK_WARN(0, "Failed to compute the global maximum of rho");
}
*GRHydro_rho_min = max_rho * rho_rel_min;
}
}
|
