aboutsummaryrefslogtreecommitdiff
path: root/cartoon_advect.c
blob: 91da4da226a0d0b18266d997c5cdcfb5e14ac71b (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
    double xx = x[index];
    int origin = fabs(xx) < 1e-8;

    PDupwindNthAnti2At11 = origin ? -2 * PDupwindNthAnti1At12 : -2 * At12[index] / xx;
    PDupwindNthSymm2At11 = origin ? -2 * PDupwindNthSymm1At12 : -2 * At12[index] / xx;
    PDupwindNthAnti2At12 = origin ? PDupwindNthAnti1At11 - PDupwindNthAntiAt22 : (At11[index] - At22[index]) / xx;
    PDupwindNthSymm2At12 = origin ? PDupwindNthSymm1At11 - PDupwindNthSymmAt22 : (At11[index] - At22[index]) / xx;
    PDupwindNthAnti2At13 = origin ? -PDupwindNthAnti1At23 : - At23[index] / xx;
    PDupwindNthSymm2At13 = origin ? -PDupwindNthSymm1At23 : - At23[index] / xx;
    PDupwindNthAnti2At22 = origin ? 2 * PDupwindNthAnti1At12 : 2 * At12[index] / xx;
    PDupwindNthSymm2At22 = origin ? 2 * PDupwindNthSymm1At12 : 2 * At12[index] / xx;
    PDupwindNthAnti2At23 = origin ? PDupwindNth1AntiAt13 : At13[index] / xx;
    PDupwindNthSymm2At23 = origin ? PDupwindNth1SymmAt13 : At13[index] / xx;
    PDupwindNthAnti2At33 = 0;
    PDupwindNthSymm2At33 = 0;

    PDupwindNthAnti2B1 = 0;
    PDupwindNthSymm2B1 = 0;
    PDupwindNthAnti2B2 = origin ? PDstandardNth1B1 : B1[index] / xx;
    PDupwindNthSymm2B2 = origin ? PDstandardNth1B1 : B1[index] / xx;
    PDupwindNthAnti2B3 = 0;
    PDupwindNthSymm2B3 = 0;

    PDupwindNthAnti2beta1 = 0;
    PDupwindNthSymm2beta1 = 0;
    PDupwindNthAnti2beta2 = PDupwindNthAntifdOrder82(&beta2[index]);
    PDupwindNthSymm2beta2 = PDupwindNthSymmfdOrder82(&beta2[index]);
    PDupwindNthAnti2beta3 = 0;
    PDupwindNthSymm2beta3 = 0;

    PDupwindNthAnti2gt11 = 0;
    PDupwindNthSymm2gt11 = 0;
    PDupwindNthAnti2gt12 = origin ? PDstandardNthAnti1gt11 - PDstandardNthAnti1gt22 : (gt11[index] - gt22[index]) / xx;
    PDupwindNthSymm2gt12 = origin ? PDstandardNthSymm1gt11 - PDstandardNthSymm1gt22 : (gt11[index] - gt22[index]) / xx;
    PDupwindNthAnti2gt13 = 0;
    PDupwindNthSymm2gt13 = 0;
    PDupwindNthAnti2gt22 = 0;
    PDupwindNthSymm2gt22 = 0;
    PDupwindNthAnti2gt23 = origin ? PDstandardNthAnti1gt13 : gt13[index] / xx;
    PDupwindNthSymm2gt23 = origin ? PDstandardNthSymm1gt13 : gt13[index] / xx;
    PDupwindNthAnti2gt33 = 0;
    PDupwindNthSymm2gt33 = 0;

    PDupwindNthAnti2phi = 0;
    PDupwindNthSymm2phi = 0;
    PDupwindNthAnti2trK = 0;
    PDupwindNthSymm2trK = 0;
    PDupwindNthAnti2Xt2 = origin ? PDstandardNthAnti1Xt1 : Xt1[index] / xx;
    PDupwindNthSymm2Xt2 = origin ? PDstandardNthSymm1Xt1 : Xt1[index] / xx;