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
|
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);
PDupwindNthAnti2At11 = kzero;
PDupwindNthSymm2At11 = kzero;
PDupwindNthAnti2At12 = kifthen(origin, ksub(PDupwindNthAnti1At11, PDupwindNthAntiAt22), kmul(ksub(At11L, At22L), xinv));
PDupwindNthSymm2At12 = kifthen(origin, ksub(PDupwindNthSymm1At11, PDupwindNthSymmAt22), kmul(ksub(At11L, At22L), xinv));
PDupwindNthAnti2At13 = kzero;
PDupwindNthSymm2At13 = kzero;
PDupwindNthAnti2At22 = kzero;
PDupwindNthSymm2At22 = kzero;
PDupwindNthAnti2At23 = kifthen(origin, PDupwindNth1AntiAt13, kmul(At13L, xinv));
PDupwindNthSymm2At23 = kifthen(origin, PDupwindNth1SymmAt13, kmul(At13L, xinv));
PDupwindNthAnti2At33 = kzero;
PDupwindNthSymm2At33 = kzero;
PDupwindNthAnti2B1 = kzero;
PDupwindNthSymm2B1 = kzero;
PDupwindNthAnti2B2 = kifthen(origin, PDstandardNthAnti1B1, kmul(B1L, xinv));
PDupwindNthSymm2B2 = kifthen(origin, PDstandardNthSymm1B1, kmul(B1L, xinv));
PDupwindNthAnti2B3 = kzero;
PDupwindNthSymm2B3 = kzero;
PDupwindNthAnti2beta1 = kzero;
PDupwindNthSymm2beta1 = kzero;
PDupwindNthAnti2beta2 = kifthen(origin, PDstandardNthAnti1beta1, kmul(beta1L, xinv));
PDupwindNthSymm2beta2 = kifthen(origin, PDstandardNthSymm1beta1, kmul(beta1L, xinv));
PDupwindNthAnti2beta3 = kzero;
PDupwindNthSymm2beta3 = kzero;
PDupwindNthAnti2gt11 = kzero;
PDupwindNthSymm2gt11 = kzero;
PDupwindNthAnti2gt12 = kifthen(origin, ksub(PDstandardNthAnti1gt11, PDstandardNthAnti1gt22), kmul(ksub(gt11L, gt22L), xinv));
PDupwindNthSymm2gt12 = kifthen(origin, ksub(PDstandardNthSymm1gt11, PDstandardNthSymm1gt22), kmul(ksub(gt11L, gt22L), xinv));
PDupwindNthAnti2gt13 = kzero;
PDupwindNthSymm2gt13 = kzero;
PDupwindNthAnti2gt22 = kzero;
PDupwindNthSymm2gt22 = kzero;
PDupwindNthAnti2gt23 = kifthen(origin, PDstandardNthAnti1gt13, kmul(gt13L, xinv));
PDupwindNthSymm2gt23 = kifthen(origin, PDstandardNthSymm1gt13, kmul(gt13L, xinv));
PDupwindNthAnti2gt33 = kzero;
PDupwindNthSymm2gt33 = kzero;
PDupwindNthAnti2phi = kzero;
PDupwindNthSymm2phi = kzero;
PDupwindNthAnti2trK = kzero;
PDupwindNthSymm2trK = kzero;
PDupwindNthAnti2Xt2 = kifthen(origin, PDstandardNthAnti1Xt1, kmul(Xt1L, xinv));
PDupwindNthSymm2Xt2 = kifthen(origin, PDstandardNthSymm1Xt1, kmul(Xt1L, xinv));
|
