| # | Line 289 | Line 289 | FVECT sd, vd, vu; | |
|---|---|---|
| 289 | d = DOT(sd,vu)/sqrt(1.0-d*d); | |
| 290 | if (d >= 1.0) | |
| 291 | tau = 0.0; | |
| 292 | + | else if (d <= -1.0) |
| 293 | + | tau = 180.0; |
| 294 | else | |
| 295 | tau = acos(d) * (180./PI); | |
| 296 | return( exp( sigma*( (35.2 - tau*.31889 - 1.22*exp(-.22222*tau))*1e-3 | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines (old) |
| > | Changed lines (new) |