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
|
CCTK_REAL_VEC kzero = ToReal(0.0);
CCTK_REAL_VEC kone = ToReal(1.0);
CCTK_REAL_VEC ktwo = ToReal(2.0);
CCTK_REAL_VEC xx = vec_load(x[index]);
CCTK_REAL_VEC x2 = kmul(xx, xx);
CCTK_REAL_VEC xinv = kdiv(kone, xx);
CCTK_REAL_VEC x2inv = kdiv(kone, x2);
CCTK_REAL_VEC absx = kfabs(xx);
CCTK_REAL_VEC eps = ToReal(1e-8);
CCTK_BOOLEAN_VEC origin = kcmplt(absx, eps);
PDstandardNth2trK = kzero;
PDstandardNth2alpha = kzero;
PDstandardNth22alpha = kifthen(origin, PDstandardNth11alpha, kmul(PDstandardNth1alpha, xinv));
PDstandardNth12alpha = kzero;
PDstandardNth23alpha = kzero;
PDstandardNth2beta1 = kzero;
PDstandardNth22beta1 = kifthen(origin, kmul(PDstandardNth11beta1, khalf),
ksub(kmul(PDstandardNth1beta1, xinv), kmul(beta1L, x2inv)));
PDstandardNth12beta1 = kzero;
PDstandardNth23beta1 = kzero;
PDstandardNth2beta2 = kifthen(origin, PDstandardNth1beta1, kmul(beta1L, xinv));
PDstandardNth22beta2 = kzero;
PDstandardNth12beta2 = kifthen(origin, kmul(PDstandardNth11beta1, khalf),
ksub(kmul(PDstandardNth1beta1, xinv), kmul(beta1L, x2inv)));
PDstandardNth23beta2 = kifthen(origin, PDstandardNth13beta1, kmul(PDstandardNth3beta1, xinv));
PDstandardNth2beta3 = kzero;
PDstandardNth22beta3 = kifthen(origin, PDstandardNth11beta3, kmul(PDstandardNth1beta3, xinv));
PDstandardNth12beta3 = kzero;
PDstandardNth23beta3 = kzero;
PDstandardNth2gt11 = kzero;
PDstandardNth22gt11 = kifthen(origin, PDstandardNth11gt22,
ksub(kmul(PDstandardNth1gt11, xinv), kmul(ktwo, kmul(ksub(gt11L, gt22L), x2inv))));
PDstandardNth12gt11 = kzero;
PDstandardNth23gt11 = kzero;
PDstandardNth2gt12 = kifthen(origin, ksub(PDstandardNth1gt11, PDstandardNth1gt22),
kmul(ksub(gt11L, gt22L), xinv));
PDstandardNth22gt12 = kzero;
PDstandardNth12gt12 = kifthen(origin, kmul(ksub(PDstandardNth11gt11, PDstandardNth11gt22), khalf),
ksub(kmul(ksub(PDstandardNth1gt11, PDstandardNth1gt22), xinv), kmul(ksub(gt11L, gt22L), x2inv)));
PDstandardNth23gt12 = kifthen(origin, ksub(PDstandardNth13gt11, PDstandardNth13gt22),
kmul(ksub(PDstandardNth3gt11, PDstandardNth3gt22), xinv));
PDstandardNth2gt13 = kzero;
PDstandardNth22gt13 = kifthen(origin, kmul(PDstandardNth11gt13, khalf),
ksub(kmul(PDstandardNth1gt13, xinv), kmul(gt13L, x2inv)));
PDstandardNth12gt13 = kzero;
PDstandardNth23gt13 = kzero;
PDstandardNth2gt22 = kzero;
PDstandardNth22gt22 = kifthen(origin, PDstandardNth11gt11,
kadd(kmul(PDstandardNth1gt22, xinv), kmul(ktwo, kmul(ksub(gt11L, gt22L), x2inv))));
PDstandardNth12gt22 = kzero;
PDstandardNth23gt22 = kzero;
PDstandardNth2gt23 = kifthen(origin, PDstandardNth1gt13, kmul(gt13L, xinv));
PDstandardNth22gt23 = kzero;
PDstandardNth12gt23 = kifthen(origin, kmul(PDstandardNth11gt13, khalf),
ksub(kmul(PDstandardNth1gt13, xinv), kmul(gt13L, x2inv)));
PDstandardNth23gt23 = kifthen(origin, PDstandardNth13gt13, kmul(PDstandardNth3gt13, xinv));
PDstandardNth2gt33 = kzero;
PDstandardNth22gt33 = kifthen(origin, PDstandardNth11gt33, kmul(PDstandardNth1gt33, xinv));
PDstandardNth12gt33 = kzero;
PDstandardNth23gt33 = kzero;
PDstandardNth2phi = kzero;
PDstandardNth22phi = kifthen(origin, PDstandardNth11phi, kmul(PDstandardNth1phi, xinv));
PDstandardNth12phi = kzero;
PDstandardNth23phi = kzero;
PDstandardNth2Xt1 = kzero;
PDstandardNth2Xt2 = kifthen(origin, PDstandardNth1Xt1, kmul(Xt1L, xinv));
PDstandardNth2Xt3 = kzero;
//PDstandardNth2At11 = origin ? -2 * PDstandardNth1At12 : -2 * At12[index] / xx;
//PDstandardNth2At13 = origin ? -PDstandardNth1At23 : - At23[index] / xx;
//PDstandardNth2At22 = origin ? 2 * PDstandardNth1At12 : 2 * At12[index] / xx;
PDstandardNth2At11 = kzero;
PDstandardNth2At12 = kifthen(origin, ksub(PDstandardNth1At11, PDstandardNth1At22),
kmul(ksub(At11L, At22L), xinv));
PDstandardNth2At13 = kzero;
PDstandardNth2At22 = kzero;
PDstandardNth2At23 = kifthen(origin, PDstandardNth1At13, kmul(At13L, xinv));
PDstandardNth2At33 = kzero;
|
