| # | Line 35 | Line 35 | void Bend::calcForce() { | |
|---|---|---|
| 35 | ||
| 36 | double sinTheta = sqrt(1.0 - cosTheta * cosTheta); | |
| 37 | ||
| 38 | < | if (fabs(sinTheta) < 1.0E - 12) { |
| 39 | < | sinTheta = 1.0E - 12; |
| 38 | > | if (fabs(sinTheta) < 1.0E -12) { |
| 39 | > | sinTheta = 1.0E -12; |
| 40 | } | |
| 41 | ||
| 42 | double commonFactor1 = -firstDerivative / sinTheta * d12inv; | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |