mirror of
https://gitlab.freedesktop.org/mesa/mesa.git
synced 2026-05-09 02:28:10 +02:00
check for tiny values in dir[] vector in calcposobjs(), prevents FP underflows later
This commit is contained in:
parent
61cb081a35
commit
a41edc31ce
6 changed files with 42 additions and 0 deletions
|
|
@ -257,6 +257,13 @@ calcposobs(void)
|
|||
dir[2] = cos(alpha * M_PI / 180.0) * sin(beta * M_PI / 180.0);
|
||||
dir[1] = cos(beta * M_PI / 180.0);
|
||||
|
||||
if (dir[0] < 1.0e-5 && dir[0] > -1.0e-5)
|
||||
dir[0] = 0;
|
||||
if (dir[1] < 1.0e-5 && dir[1] > -1.0e-5)
|
||||
dir[1] = 0;
|
||||
if (dir[2] < 1.0e-5 && dir[2] > -1.0e-5)
|
||||
dir[2] = 0;
|
||||
|
||||
obs[0] += v * dir[0];
|
||||
obs[1] += v * dir[1];
|
||||
obs[2] += v * dir[2];
|
||||
|
|
|
|||
|
|
@ -205,6 +205,13 @@ calcposobs(void)
|
|||
dir[1] = cos(alpha * M_PI / 180.0) * sin(beta * M_PI / 180.0);
|
||||
dir[2] = cos(beta * M_PI / 180.0);
|
||||
|
||||
if (dir[0] < 1.0e-5 && dir[0] > -1.0e-5)
|
||||
dir[0] = 0;
|
||||
if (dir[1] < 1.0e-5 && dir[1] > -1.0e-5)
|
||||
dir[1] = 0;
|
||||
if (dir[2] < 1.0e-5 && dir[2] > -1.0e-5)
|
||||
dir[2] = 0;
|
||||
|
||||
obs[0] += v * dir[0];
|
||||
obs[1] += v * dir[1];
|
||||
obs[2] += v * dir[2];
|
||||
|
|
|
|||
|
|
@ -102,6 +102,13 @@ calcposobs(void)
|
|||
dir[1] = cos(alpha * M_PI / 180.0) * sin(beta * M_PI / 180.0);
|
||||
dir[2] = cos(beta * M_PI / 180.0);
|
||||
|
||||
if (dir[0] < 1.0e-5 && dir[0] > -1.0e-5)
|
||||
dir[0] = 0;
|
||||
if (dir[1] < 1.0e-5 && dir[1] > -1.0e-5)
|
||||
dir[1] = 0;
|
||||
if (dir[2] < 1.0e-5 && dir[2] > -1.0e-5)
|
||||
dir[2] = 0;
|
||||
|
||||
obs[0] += v * dir[0];
|
||||
obs[1] += v * dir[1];
|
||||
obs[2] += v * dir[2];
|
||||
|
|
|
|||
|
|
@ -81,6 +81,13 @@ calcposobs(void)
|
|||
dir[2] = cos(alpha * M_PI / 180.0) * sin(beta * M_PI / 180.0);
|
||||
dir[1] = cos(beta * M_PI / 180.0);
|
||||
|
||||
if (dir[0] < 1.0e-5 && dir[0] > -1.0e-5)
|
||||
dir[0] = 0;
|
||||
if (dir[1] < 1.0e-5 && dir[1] > -1.0e-5)
|
||||
dir[1] = 0;
|
||||
if (dir[2] < 1.0e-5 && dir[2] > -1.0e-5)
|
||||
dir[2] = 0;
|
||||
|
||||
alpha1 = alpha + FOV / 2.0;
|
||||
v1[0] = sin(alpha1 * M_PI / 180.0);
|
||||
v1[1] = cos(alpha1 * M_PI / 180.0);
|
||||
|
|
|
|||
|
|
@ -164,6 +164,13 @@ calcposobs(void)
|
|||
dir[1] = cos(alpha * M_PI / 180.0) * sin(beta * M_PI / 180.0);
|
||||
dir[2] = cos(beta * M_PI / 180.0);
|
||||
|
||||
if (dir[0] < 1.0e-5 && dir[0] > -1.0e-5)
|
||||
dir[0] = 0;
|
||||
if (dir[1] < 1.0e-5 && dir[1] > -1.0e-5)
|
||||
dir[1] = 0;
|
||||
if (dir[2] < 1.0e-5 && dir[2] > -1.0e-5)
|
||||
dir[2] = 0;
|
||||
|
||||
obs[0] += v * dir[0];
|
||||
obs[1] += v * dir[1];
|
||||
obs[2] += v * dir[2];
|
||||
|
|
|
|||
|
|
@ -162,6 +162,13 @@ calcposobs(void)
|
|||
dir[1] = cos(alpha * M_PI / 180.0) * sin(beta * M_PI / 180.0);
|
||||
dir[2] = cos(beta * M_PI / 180.0);
|
||||
|
||||
if (dir[0] < 1.0e-5 && dir[0] > -1.0e-5)
|
||||
dir[0] = 0;
|
||||
if (dir[1] < 1.0e-5 && dir[1] > -1.0e-5)
|
||||
dir[1] = 0;
|
||||
if (dir[2] < 1.0e-5 && dir[2] > -1.0e-5)
|
||||
dir[2] = 0;
|
||||
|
||||
obs[0] += v * dir[0];
|
||||
obs[1] += v * dir[1];
|
||||
obs[2] += v * dir[2];
|
||||
|
|
|
|||
Loading…
Add table
Reference in a new issue