79 |
|
v->vvinc[1] *= dt; |
80 |
|
v->vvinc[2] *= dt; |
81 |
|
|
82 |
< |
v->vhs2 = 1.0/DOT(v->vhinc,v->vhinc); |
83 |
< |
v->vvs2 = 1.0/DOT(v->vvinc,v->vvinc); |
82 |
> |
v->vhn2 = DOT(v->vhinc,v->vhinc); |
83 |
> |
v->vvn2 = DOT(v->vvinc,v->vvinc); |
84 |
|
|
85 |
|
return(NULL); |
86 |
|
} |
136 |
|
disp[1] *= d; |
137 |
|
disp[2] *= d; |
138 |
|
} |
139 |
< |
*xp = DOT(disp,v->vhinc)*v->vhs2 + 0.5*v->hresolu; |
140 |
< |
*yp = DOT(disp,v->vvinc)*v->vvs2 + 0.5*v->vresolu; |
139 |
> |
*xp = DOT(disp,v->vhinc)/v->vhn2 + 0.5*v->hresolu; |
140 |
> |
*yp = DOT(disp,v->vvinc)/v->vvn2 + 0.5*v->vresolu; |
141 |
|
} |
142 |
|
|
143 |
|
|