Archives/MikeMcl_decimal.js
1
0
Fork 0
mirror of https://github.com/MikeMcl/decimal.js.git synced 2024-10-27 20:34:12 +00:00
Code Issues Projects Releases Wiki Activity

use a method less prone to cancellation near x=1

Browse Source
This commit is contained in:
Alden 2023-02-13 03:50:30 -05:00
parent e634266aad
commit 2b1298cd69
2 changed files with 5 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
decimal.js
View File

@ -745,13 +745,12 @@
Ctor.precision = pr + 6; Ctor.precision = pr + 6;
Ctor.rounding = 1; Ctor.rounding = 1;
x = x.asin(); x = (new Ctor(1)).minus(x).div(x.plus(1)).sqrt().atan();
halfPi = getPi(Ctor, pr + 4, rm).times(0.5);
Ctor.precision = pr; Ctor.precision = pr;
Ctor.rounding = rm; Ctor.rounding = rm;
return halfPi.minus(x); return x.times(2);
}; };

5
decimal.mjs
View File

@ -741,13 +741,12 @@ P.inverseCosine = P.acos = function () {
Ctor.precision = pr + 6; Ctor.precision = pr + 6;
Ctor.rounding = 1; Ctor.rounding = 1;
x = x.asin(); x = (new Ctor(1)).minus(x).div(x.plus(1)).sqrt().atan();
halfPi = getPi(Ctor, pr + 4, rm).times(0.5);
Ctor.precision = pr; Ctor.precision = pr;
Ctor.rounding = rm; Ctor.rounding = rm;
return halfPi.minus(x); return x.times(2);
}; };