24 |
|
|
25 |
|
double comf2, comf3, comf4; |
26 |
|
double dcsidx, dcsidy, dcsidz, dcskdx, dcskdy, dcskdz; |
27 |
< |
double dcsjdx, dcsjdy, dcsjdz; |
27 |
> |
// double dcsjdx, dcsjdy, dcsjdz; |
28 |
|
double dadxi, dadyi, dadzi; |
29 |
< |
double dadxk, dadyk, dadzk, dadxj, dadyj, dadzj; |
29 |
> |
double dadxk, dadyk, dadzk;//, dadxj, dadyj, dadzj; |
30 |
|
double daxi, dayi, dazi, daxk, dayk, dazk, daxj, dayj, dazj; |
31 |
|
|
32 |
|
|
80 |
|
dcskdy = dy*denom - comf3*gy; |
81 |
|
dcskdz = dz*denom - comf3*gz; |
82 |
|
|
83 |
< |
dcsjdx = -dcsidx - dcskdx; |
84 |
< |
dcsjdy = -dcsidy - dcskdy; |
85 |
< |
dcsjdz = -dcsidz - dcskdz; |
83 |
> |
// dcsjdx = -dcsidx - dcskdx; |
84 |
> |
// dcsjdy = -dcsidy - dcskdy; |
85 |
> |
// dcsjdz = -dcsidz - dcskdz; |
86 |
|
|
87 |
|
dadxi = -sinai*dcsidx; |
88 |
|
dadyi = -sinai*dcsidy; |
92 |
|
dadyk = -sinai*dcskdy; |
93 |
|
dadzk = -sinai*dcskdz; |
94 |
|
|
95 |
< |
dadxj = -dadxi - dadxk; |
96 |
< |
dadyj = -dadyi - dadyk; |
97 |
< |
dadzj = -dadzi - dadzk; |
95 |
> |
// dadxj = -dadxi - dadxk; |
96 |
> |
// dadyj = -dadyi - dadyk; |
97 |
> |
// dadzj = -dadzi - dadzk; |
98 |
|
|
99 |
|
daxi = comf4*dadxi; |
100 |
|
dayi = comf4*dadyi; |