-
Notifications
You must be signed in to change notification settings - Fork 0
/
GetPI.e
56 lines (42 loc) · 3.03 KB
/
GetPI.e
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
-- Copyright (c) 2016-2023 James Cook
-- PI functions:
-- GetPI, GetTwoPI, GetHalfPI, GetQuarterPI
include ../../eunumber/minieun/common.e
include ../../eunumber/minieun/defaults.e
include ../../eunumber/minieun/eun.e
include ../../eunumber/minieun/AdjustRound.e
include ../../eunumber/eun/EunMultiply.e
include EunArcTan.e
-- constant strPI =
-- "3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679821480865132823066470938446095505822317253594081284811174502841027019385211055596446229489549303819644288109756659334461284756482337867831652712019091456485669234603486104543266482133936072602491412737245870066063155881748815209209628292540917153643678925903600113305305488204665213841469519415116094330572703657595919530921861173819326117931051185480744623799627495673518857527248912279381830119491298336733624406566430860213949463952247371907021798609437027705392171762931767523846748184676694051320005681271452635608277857713427577896091736371787214684409012249534301465495853710507922796892589235420199561121290219608640344181598136297747713099605187072113499999983729780499510597317328160963185950244594553469083026425223082533446850352619311881710100031378387528865875332083814206171776691473035982534904287554687311595628638823537875937519577818577805321712268066130019278766111959092164201989"
-- quarterPI = EunMultiply(ToEun(strPI, 10, length(strPI) - 1), NewEun({2, 5}, -1, length(strPI) - 1, 10))
global sequence quarterPI = {}
-- NOTE: To precalculate, put the largest value for targetLength first, then use the same radix for all your calculations, before switching to another radix.
-- You can also use the "SwapQuarterPI()" function below:
global function SwapQuarterPI(sequence s = quarterPI) -- SwapQuarterPI when changing radixes, then swap back when changing back.
object oldvalue = quarterPI
quarterPI = s
return oldvalue
end function
global function GetQuarterPI(TargetLength targetLength = defaultTargetLength, AtomRadix radix = defaultRadix, PositiveInteger multBy = 1)
sequence ret
-- targetLength += adjustPrecision
if not length(quarterPI) or not length(quarterPI[1]) or quarterPI[3] <= targetLength or quarterPI[4] != radix then
quarterPI = ArcTanExp({1}, 0, targetLength + 1, radix)
end if
ret = AdjustRound(quarterPI[1], quarterPI[2], targetLength, radix, NO_SUBTRACT_ADJUST)
if multBy != 1 then
object tmp = AdjustRound({multBy}, 0, targetLength, radix, 0) -- 0 makes it use Carry()
ret = EunMultiply(ret, tmp)
end if
return ret
end function
global function GetHalfPI(TargetLength targetLength = defaultTargetLength, integer radix = defaultRadix)
return GetQuarterPI(targetLength, radix, 2)
end function
global function GetPI(TargetLength targetLength = defaultTargetLength, integer radix = defaultRadix)
return GetQuarterPI(targetLength, radix, 4)
end function
global function GetTwoPI(TargetLength targetLength = defaultTargetLength, integer radix = defaultRadix)
return GetQuarterPI(targetLength, radix, 8)
end function