blob: 523d0d53776de216346c96009a0d2d634133569d (
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
|
/*@@
@header HAMADM_guts.h
@date Aug 98
@author Gabrielle Allen
@desc
Macro to calculate the spacetime part of the
Hamiltonian Constraint. That is:
R - K^i_j K^j_i + trK^2
@enddesc
@@*/
#ifndef HAMADM_GUTS
#define HAMADM_GUTS
#include "TRRICCI_guts.h"
#include "TRKK_guts.h"
#include "TRK_guts.h"
#ifdef FCODE
HAMADM_HAMADM = TRRICCI_TRRICCI - TRKK_TRKK + TRK_TRK**2
HAMADM_HAMADMABS = abs(TRRICCI_TRRICCI) + abs(TRKK_TRKK) + TRK_TRK**2
#endif
#ifdef CCODE
HAMADM_HAMADM = TRRICCI_TRRICCI - TRKK_TRKK + TRK_TRK**2
HAMADM_HAMADMABS = fabs(TRRICCI_TRRICCI) + fabs(TRKK_TRKK) + TRK_TRK**2
#endif
#endif
|