blob: de4c849028f82dc5a6340a628b5d5809d3bcd458 (
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
|
/*@@
@header DDA_guts.h
@date Jul 98
@author Gabrielle Allen
@desc
Macro to calculate all second spatial derivative of lapse
@enddesc
@@*/
#ifndef DDA_GUTS
#define DDA_GUTS
#ifdef FCODE
DDA_OODX2 = 1D0/(DDA_DX*DDA_DX)
DDA_OODY2 = 1D0/(DDA_DY*DDA_DY)
DDA_OODZ2 = 1D0/(DDA_DZ*DDA_DZ)
DDA_OO4DXDY = 1D0/(4D0*DDA_DX*DDA_DY)
DDA_OO4DXDZ = 1D0/(4D0*DDA_DX*DDA_DZ)
DDA_OO4DYDZ = 1D0/(4D0*DDA_DY*DDA_DZ)
DDA_DXXDA = DDA_OODX2*(DDA_A_IP - 2D0*DDA_A + DDA_A_IM)
DDA_DYYDA = DDA_OODY2*(DDA_A_JP - 2D0*DDA_A + DDA_A_JM)
DDA_DZZDA = DDA_OODZ2*(DDA_A_KP - 2D0*DDA_A + DDA_A_KM)
DDA_DXYDA = DDA_OO4DXDY*(DDA_A_IPJP-DDA_A_IPJM-DDA_A_IMJP+DDA_A_IMJM)
DDA_DXZDA = DDA_OO4DXDZ*(DDA_A_IPKP-DDA_A_IPKM-DDA_A_IMKP+DDA_A_IMKM)
DDA_DYZDA = DDA_OO4DYDZ*(DDA_A_JPKP-DDA_A_JPKM-DDA_A_JMKP+DDA_A_JMKM)
#endif
#ifdef CCODE
DDA_OODX2 = 1/(DDA_DX*DDA_DX);
DDA_OODY2 = 1/(DDA_DY*DDA_DY);
DDA_OODZ2 = 1/(DDA_DZ*DDA_DZ);
DDA_OO4DXDY = 1/(4*DDA_DX*DDA_DY);
DDA_OO4DXDZ = 1/(4*DDA_DX*DDA_DZ);
DDA_OO4DYDZ = 1/(4*DDA_DY*DDA_DZ);
DDA_DXXDA = DDA_OODX2*(DDA_A_IP - 2*DDA_A + DDA_A_IM);
DDA_DYYDA = DDA_OODY2*(DDA_A_JP - 2*DDA_A + DDA_A_JM);
DDA_DZZDA = DDA_OODZ2*(DDA_A_KP - 2*DDA_A + DDA_A_KM);
DDA_DXYDA = DDA_OO4DXDY*(DDA_A_IPJP-DDA_A_IPJM-DDA_A_IMJP+DDA_A_IMJM);
DDA_DXZDA = DDA_OO4DXDZ*(DDA_A_IPKP-DDA_A_IPKM-DDA_A_IMKP+DDA_A_IMKM);
DDA_DYZDA = DDA_OO4DYDZ*(DDA_A_JPKP-DDA_A_JPKM-DDA_A_JMKP+DDA_A_JMKM);
#endif
#endif
|