diff options
author | Nicholas Noll <nbnoll@eml.cc> | 2021-12-05 16:16:21 -0800 |
---|---|---|
committer | Nicholas Noll <nbnoll@eml.cc> | 2021-12-05 16:16:21 -0800 |
commit | 07e77936d535e58b0aeb4f2a11400c1050556739 (patch) | |
tree | fb50fad6436ecbc159505bee73a92cd289b99a07 /src/base/math/asin.c | |
parent | b48327d357e0818d1a6ae2a064cfa7d1567e1242 (diff) |
Feat: added math library
Used Plan9's libc as starting point. This cleans up dangling references
due to loss of libc.
Diffstat (limited to 'src/base/math/asin.c')
-rw-r--r-- | src/base/math/asin.c | 40 |
1 files changed, 40 insertions, 0 deletions
diff --git a/src/base/math/asin.c b/src/base/math/asin.c new file mode 100644 index 0000000..4687715 --- /dev/null +++ b/src/base/math/asin.c @@ -0,0 +1,40 @@ +/* + * asin(arg) and acos(arg) return the arcsin, arccos, + * respectively of their arguments. + * + * Arctan is called after appropriate range reduction. + */ + +#include <u.h> +#include <base.h> + +double +math·asin(double arg) +{ + double temp; + int sign; + + sign = 0; + if(arg < 0) { + arg = -arg; + sign++; + } + if(arg > 1) + return math·NaN(); + temp = math·sqrt(1 - arg*arg); + if(arg > 0.7) + temp = math·PIO2 - math·atan(temp/arg); + else + temp = math·atan(arg/temp); + if(sign) + temp = -temp; + return temp; +} + +double +math·acos(double arg) +{ + if(arg > 1 || arg < -1) + return math·NaN(); + return math·PIO2 - math·asin(arg); +} |