| # | Line 82 | Line 82 | register RAY *r; | |
|---|---|---|
| 82 | objerror(m, WARNING, "compute error"); | |
| 83 | return; | |
| 84 | } | |
| 85 | + | multv3(disp, disp, mxf->fore.xfm); |
| 86 | + | multv3(disp, disp, r->rofx); |
| 87 | for (i = 0; i < 3; i++) | |
| 88 | < | r->pert[i] += ( disp[0]*mxf->fore.xfm[0][i] + |
| 87 | < | disp[1]*mxf->fore.xfm[1][i] + |
| 88 | < | disp[2]*mxf->fore.xfm[2][i] ) |
| 89 | < | / mxf->fore.sca; |
| 88 | > | r->pert[i] += disp[i] / (mxf->fore.sca * r->rofs); |
| 89 | return; | |
| 90 | memerr: | |
| 91 | error(SYSTEM, "out of memory in t_func"); | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |