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
|
/* Add the contributions to T_munu that we gathered via the AddToTmunu
mechanism, using the CalcTmunu.inc mechanism. */
/* Only add if we are not included from this thorn; otherwise we would
be adding to ourselves. */
#ifndef TMUNUBASE_SETTMUNU
#ifdef CCODE
if (* stress_energy_2_state)
{
#ifndef NDEBUG
if ((ijk) < 0 || (ijk) >= cctk_lsh[0] * cctk_lsh[1] * cctk_lsh[2])
{
CCTK_WARN (0, "Illegal array index");
}
#endif
Ttt += eT2tt[ijk];
Ttx += eT2tx[ijk];
Tty += eT2ty[ijk];
Ttz += eT2tz[ijk];
Txx += eT2xx[ijk];
Txy += eT2xy[ijk];
Txz += eT2xz[ijk];
Tyy += eT2yy[ijk];
Tyz += eT2yz[ijk];
Tzz += eT2zz[ijk];
}
#endif
#ifdef FCODE
if (stress_energy_2_state .ne. 0) then
Ttt = Ttt + eT2tt(i,j,k)
Ttx = Ttx + eT2tx(i,j,k)
Tty = Tty + eT2ty(i,j,k)
Ttz = Ttz + eT2tz(i,j,k)
Txx = Txx + eT2xx(i,j,k)
Txy = Txy + eT2xy(i,j,k)
Txz = Txz + eT2xz(i,j,k)
Tyy = Tyy + eT2yy(i,j,k)
Tyz = Tyz + eT2yz(i,j,k)
Tzz = Tzz + eT2zz(i,j,k)
end if
#endif
#endif /* ifndef TMUNUBASE_SETTMUNU */
|