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
|
#include <math.h>
#include <cctk.h>
#include <cctk_Arguments.h>
#include <cctk_Parameters.h>
static void
newrad (cGH const * restrict cctkGH,
CCTK_REAL const * restrict var,
CCTK_REAL * restrict rhs,
CCTK_REAL var0,
CCTK_REAL v0);
void
ML_BSSN_NewRad (CCTK_ARGUMENTS)
{
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
CCTK_REAL const v0 = sqrt (harmonicF);
newrad (cctkGH, phi , phirhs , 0.0, v0 );
newrad (cctkGH, gt11 , gt11rhs , 1.0, 1.0);
newrad (cctkGH, gt12 , gt12rhs , 0.0, 1.0);
newrad (cctkGH, gt13 , gt13rhs , 0.0, 1.0);
newrad (cctkGH, gt22 , gt22rhs , 1.0, 1.0);
newrad (cctkGH, gt23 , gt23rhs , 0.0, 1.0);
newrad (cctkGH, gt33 , gt33rhs , 1.0, 1.0);
newrad (cctkGH, Xt1 , Xt1rhs , 0.0, 1.0);
newrad (cctkGH, Xt2 , Xt2rhs , 0.0, 1.0);
newrad (cctkGH, Xt3 , Xt3rhs , 0.0, 1.0);
newrad (cctkGH, trK , trKrhs , 0.0, v0 );
newrad (cctkGH, At11 , At11rhs , 0.0, 1.0);
newrad (cctkGH, At12 , At12rhs , 0.0, 1.0);
newrad (cctkGH, At13 , At13rhs , 0.0, 1.0);
newrad (cctkGH, At22 , At22rhs , 0.0, 1.0);
newrad (cctkGH, At23 , At23rhs , 0.0, 1.0);
newrad (cctkGH, At33 , At33rhs , 0.0, 1.0);
newrad (cctkGH, alpha, alpharhs, 1.0, v0 );
newrad (cctkGH, A , Arhs , 0.0, v0 );
newrad (cctkGH, beta1, beta1rhs, 0.0, 1.0);
newrad (cctkGH, beta2, beta2rhs, 0.0, 1.0);
newrad (cctkGH, beta3, beta3rhs, 0.0, 1.0);
newrad (cctkGH, B1 , B1rhs , 0.0, 1.0);
newrad (cctkGH, B2 , B2rhs , 0.0, 1.0);
newrad (cctkGH, B3 , B3rhs , 0.0, 1.0);
}
static void
newrad (cGH const * restrict const cctkGH,
CCTK_REAL const * restrict const var,
CCTK_REAL * restrict const rhs,
CCTK_REAL const var0,
CCTK_REAL const v0)
{
DECLARE_CCTK_PARAMETERS;
NewRad_Apply (cctkGH, var, rhs, var0, v0, radpower);
}
|