blob: badc89c87215d8a55c791a607d59ee4680a358ad (
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
|
double xx = x[index];
int origin = fabs(xx) < 1e-8;
PDstandardNth2gt11 = 0;
PDstandardNth22gt11 = origin ? PDstandardNth11gt22 : PDstandardNth1gt11 / xx - 2 * (gt11[index] - gt22[index]) / (xx * xx);
PDstandardNth12gt11 = 0;
PDstandardNth23gt11 = 0;
PDstandardNth2gt12 = origin ? PDstandardNth1gt11 - PDstandardNth1gt22 : (gt11[index] - gt22[index]) / xx;
PDstandardNth22gt12 = 0;
PDstandardNth12gt12 = origin ? (PDstandardNth11gt11 - PDstandardNth11gt22) / 2 : (PDstandardNth1gt11 - PDstandardNth1gt22) / xx - (gt11[index] - gt22[index]) / (xx * xx);
PDstandardNth23gt12 = origin ? PDstandardNth13gt11 - PDstandardNth13gt22 : (PDstandardNth3gt11 - PDstandardNth3gt22) / xx;
PDstandardNth2gt13 = 0;
PDstandardNth22gt13 = origin ? PDstandardNth11gt13 / 2 : PDstandardNth1gt13 / xx - gt13[index] / (xx * xx);
PDstandardNth12gt13 = 0;
PDstandardNth23gt13 = 0;
PDstandardNth2gt22 = 0;
PDstandardNth22gt22 = origin ? PDstandardNth11gt11 : PDstandardNth1gt22 / xx + 2 * (gt11[index] - gt22[index]) / (xx * xx);
PDstandardNth12gt22 = 0;
PDstandardNth23gt22 = 0;
PDstandardNth2gt23 = origin ? PDstandardNth1gt13 : gt13[index] / xx;
PDstandardNth22gt23 = 0;
PDstandardNth12gt23 = origin ? PDstandardNth11gt13 / 2 : PDstandardNth1gt13 / xx - gt13[index] / (xx * xx);
PDstandardNth23gt23 = origin ? PDstandardNth13gt13 : PDstandardNth3gt13 / xx;
PDstandardNth2gt33 = 0;
PDstandardNth22gt33 = origin ? PDstandardNth11gt33 : PDstandardNth1gt33 / xx;
PDstandardNth12gt33 = 0;
PDstandardNth23gt33 = 0;
PDstandardNth2phi = 0;
PDstandardNth22phi = origin ? PDstandardNth11phi : PDstandardNth1phi / xx;
PDstandardNth12phi = 0;
PDstandardNth23phi = 0;
PDstandardNth2trK = 0;
PDstandardNth2Xt1 = 0;
PDstandardNth2Xt2 = origin ? PDstandardNth1Xt1 : Xt1[index] / xx;
PDstandardNth2Xt3 = 0;
PDstandardNth2At11 = origin ? -2 * PDstandardNth1At12 : -2 * At12[index] / xx;
PDstandardNth2At12 = origin ? PDstandardNth1At11 - PDstandardNth1At22 : (At11[index] - At22[index]) / xx;
PDstandardNth2At13 = origin ? -PDstandardNth1At23 : - At23[index] / xx;
PDstandardNth2At22 = origin ? 2 * PDstandardNth1At12 : 2 * At12[index] / xx;
PDstandardNth2At23 = origin ? PDstandardNth1At13 : At13[index] / xx;
PDstandardNth2At33 = 0;
|
