# | Line 24 | Line 24 | void Bend::calc_forces(){ | |
---|---|---|
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 | ||
# | Line 80 | Line 80 | void Bend::calc_forces(){ | |
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; | |
# | Line 92 | Line 92 | void Bend::calc_forces(){ | |
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; |
– | Removed lines |
+ | Added lines |
< | Changed lines |
> | Changed lines |