implement acos(), asin(), atan()

This commit is contained in:
Brian 2007-05-01 14:02:49 -06:00
parent bfd5cf72c4
commit 594b5ad87d
1 changed files with 14 additions and 9 deletions

View File

@ -313,9 +313,16 @@ vec4 tan(const vec4 angle)
float asin(const float x)
{
// XXX FIX ME!
// __asm float_arcsine y, x;
__retVal = 0.5;
const float a0 = 1.5707288; // PI/2?
const float a1 = -0.2121144;
const float a2 = 0.0742610;
//const float a3 = -0.0187293;
const float halfPi = 3.1415926 * 0.5;
const float y = abs(x);
// three terms seem to be enough:
__retVal = (halfPi - sqrt(1.0 - y) * (a0 + y * (a1 + a2 * y))) * sign(x);
// otherwise, try four:
//__retVal = (halfPi - sqrt(1.0 - y) * (a0 + y * (a1 + y * (a2 + y * a3)))) * sign(x);
}
vec2 asin(const vec2 v)
@ -340,8 +347,8 @@ vec4 asin(const vec4 v)
float acos(const float x)
{
// XXX FIX ME!
__retVal = 0.5;
const float halfPi = 3.1415926 * 0.5;
__retVal = halfPi - asin(x);
}
vec2 acos(const vec2 v)
@ -365,11 +372,9 @@ vec4 acos(const vec4 v)
__retVal.w = acos(v.w);
}
float atan(const float y_over_x)
float atan(const float x)
{
// XXX FIX ME
//__asm float_arctan z, y_over_x;
__retVal = 0.5;
__retVal = asin(x * inversesqrt(x * x + 1.0));
}
vec2 atan(const vec2 y_over_x)