blob: 2635ca933c12414309aa7e7b8cbcf1701122f70d (
plain)
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
|
/*@@
@header WAVEG_guts.h
@date Jul 98
@author Gabrielle Allen
@desc
Macro to calculate wave operator acting on the metric
That is g^lm g_ij,lm
@enddesc
@@*/
#ifndef WAVEG_GUTS
#define WAVEG_GUTS
#include "UPPERMET_guts.h"
#include "DDG_guts.h"
#ifdef FCODE
WAVEG_DDGXX = UPPERMET_UXX*DDG_DXXDGXX +
& 2D0*UPPERMET_UXY*DDG_DXYDGXX +
& 2D0*UPPERMET_UXZ*DDG_DXZDGXX +
& UPPERMET_UYY*DDG_DYYDGXX +
& 2D0*UPPERMET_UYZ*DDG_DYZDGXX +
& UPPERMET_UZZ*DDG_DZZDGXX
WAVEG_DDGXY = UPPERMET_UXX*DDG_DXXDGXY +
& 2D0*UPPERMET_UXY*DDG_DXYDGXY +
& 2D0*UPPERMET_UXZ*DDG_DXZDGXY +
& UPPERMET_UYY*DDG_DYYDGXY +
& 2D0*UPPERMET_UYZ*DDG_DYZDGXY +
& UPPERMET_UZZ*DDG_DZZDGXY
WAVEG_DDGXZ = UPPERMET_UXX*DDG_DXXDGXZ +
& 2D0*UPPERMET_UXY*DDG_DXYDGXZ +
& 2D0*UPPERMET_UXZ*DDG_DXZDGXZ +
& UPPERMET_UYY*DDG_DYYDGXZ +
& 2D0*UPPERMET_UYZ*DDG_DYZDGXZ +
& UPPERMET_UZZ*DDG_DZZDGXZ
WAVEG_DDGYY = UPPERMET_UXX*DDG_DXXDGYY +
& 2D0*UPPERMET_UXY*DDG_DXYDGYY +
& 2D0*UPPERMET_UXZ*DDG_DXZDGYY +
& UPPERMET_UYY*DDG_DYYDGYY +
& 2D0*UPPERMET_UYZ*DDG_DYZDGYY +
& UPPERMET_UZZ*DDG_DZZDGYY
WAVEG_DDGYZ = UPPERMET_UXX*DDG_DXXDGYZ +
& 2D0*UPPERMET_UXY*DDG_DXYDGYZ +
& 2D0*UPPERMET_UXZ*DDG_DXZDGYZ +
& UPPERMET_UYY*DDG_DYYDGYZ +
& 2D0*UPPERMET_UYZ*DDG_DYZDGYZ +
& UPPERMET_UZZ*DDG_DZZDGYZ
WAVEG_DDGZZ = UPPERMET_UXX*DDG_DXXDGZZ +
& 2D0*UPPERMET_UXY*DDG_DXYDGZZ +
& 2D0*UPPERMET_UXZ*DDG_DXZDGZZ +
& UPPERMET_UYY*DDG_DYYDGZZ +
& 2D0*UPPERMET_UYZ*DDG_DYZDGZZ +
& UPPERMET_UZZ*DDG_DZZDGZZ
#endif
#ifdef CCODE
WAVEG_DDGXX = UPPERMET_UXX*DDG_DXXDGXX +
2*UPPERMET_UXY*DDG_DXYDGXX +
2*UPPERMET_UXZ*DDG_DXZDGXX +
UPPERMET_UYY*DDG_DYYDGXX +
2*UPPERMET_UYZ*DDG_DYZDGXX +
UPPERMET_UZZ*DDG_DYYDGXX;
WAVEG_DDGXY = UPPERMET_UXX*DDG_DXXDGXY +
2*UPPERMET_UXY*DDG_DXYDGXY +
2*UPPERMET_UXZ*DDG_DXZDGXY +
UPPERMET_UYY*DDG_DYYDGXY +
2*UPPERMET_UYZ*DDG_DYZDGXY +
UPPERMET_UZZ*DDG_DYYDGXY;
WAVEG_DDGXZ = UPPERMET_UXX*DDG_DXXDGXZ +
2*UPPERMET_UXY*DDG_DXYDGXZ +
2*UPPERMET_UXZ*DDG_DXZDGXZ +
UPPERMET_UYY*DDG_DYYDGXZ +
2*UPPERMET_UYZ*DDG_DYZDGXZ +
UPPERMET_UZZ*DDG_DYYDGXZ;
WAVEG_DDGYY = UPPERMET_UXX*DDG_DXXDGYY +
2*UPPERMET_UXY*DDG_DXYDGYY +
2*UPPERMET_UXZ*DDG_DXZDGYY +
UPPERMET_UYY*DDG_DYYDGYY +
2*UPPERMET_UYZ*DDG_DYZDGYY +
UPPERMET_UZZ*DDG_DYYDGYY;
WAVEG_DDGYZ = UPPERMET_UXX*DDG_DXXDGYZ +
2*UPPERMET_UXY*DDG_DXYDGYZ +
2*UPPERMET_UXZ*DDG_DXZDGYZ +
UPPERMET_UYY*DDG_DYYDGYZ +
2*UPPERMET_UYZ*DDG_DYZDGYZ +
UPPERMET_UZZ*DDG_DYYDGYZ;
WAVEG_DDGZZ = UPPERMET_UXX*DDG_DXXDGZZ +
2*UPPERMET_UXY*DDG_DXYDGZZ +
2*UPPERMET_UXZ*DDG_DXZDGZZ +
UPPERMET_UYY*DDG_DYYDGZZ +
2*UPPERMET_UYZ*DDG_DYZDGZZ +
UPPERMET_UZZ*DDG_DYYDGZZ;
#endif
#endif
|
