long-decimal 0.01.01 → 0.01.02
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- data/README +17 -94
- data/VERSION +1 -1
- data/doc/classes/Integer.html +239 -0
- data/doc/classes/Integer.src/M000036.html +18 -0
- data/doc/classes/Integer.src/M000039.html +145 -0
- data/doc/classes/LongDecimal.html +453 -399
- data/doc/classes/LongDecimal.src/M000040.html +4 -4
- data/doc/classes/LongDecimal.src/M000041.html +4 -4
- data/doc/classes/LongDecimal.src/M000042.html +4 -6
- data/doc/classes/LongDecimal.src/M000043.html +4 -106
- data/doc/classes/LongDecimal.src/M000044.html +4 -21
- data/doc/classes/LongDecimal.src/M000045.html +4 -11
- data/doc/classes/LongDecimal.src/M000046.html +6 -27
- data/doc/classes/LongDecimal.src/M000047.html +103 -16
- data/doc/classes/LongDecimal.src/M000048.html +21 -4
- data/doc/classes/LongDecimal.src/M000049.html +10 -14
- data/doc/classes/LongDecimal.src/M000050.html +22 -36
- data/doc/classes/LongDecimal.src/M000051.html +34 -4
- data/doc/classes/LongDecimal.src/M000052.html +17 -6
- data/doc/classes/LongDecimal.src/M000053.html +4 -6
- data/doc/classes/LongDecimal.src/M000054.html +15 -4
- data/doc/classes/LongDecimal.src/M000055.html +44 -12
- data/doc/classes/LongDecimal.src/M000056.html +4 -7
- data/doc/classes/LongDecimal.src/M000057.html +8 -4
- data/doc/classes/LongDecimal.src/M000058.html +6 -15
- data/doc/classes/LongDecimal.src/M000059.html +4 -13
- data/doc/classes/LongDecimal.src/M000060.html +15 -4
- data/doc/classes/LongDecimal.src/M000061.html +7 -5
- data/doc/classes/LongDecimal.src/M000062.html +4 -4
- data/doc/classes/LongDecimal.src/M000063.html +15 -5
- data/doc/classes/LongDecimal.src/M000064.html +13 -5
- data/doc/classes/LongDecimal.src/M000065.html +4 -4
- data/doc/classes/LongDecimal.src/M000067.html +4 -8
- data/doc/classes/LongDecimal.src/M000068.html +5 -9
- data/doc/classes/LongDecimal.src/M000069.html +5 -9
- data/doc/classes/LongDecimal.src/M000070.html +4 -9
- data/doc/classes/LongDecimal.src/M000071.html +4 -4
- data/doc/classes/LongDecimal.src/M000072.html +6 -17
- data/doc/classes/LongDecimal.src/M000073.html +7 -7
- data/doc/classes/LongDecimal.src/M000074.html +7 -7
- data/doc/classes/LongDecimal.src/M000075.html +7 -19
- data/doc/classes/LongDecimal.src/M000076.html +4 -8
- data/doc/classes/LongDecimal.src/M000077.html +19 -5
- data/doc/classes/LongDecimal.src/M000078.html +7 -7
- data/doc/classes/LongDecimal.src/M000079.html +7 -7
- data/doc/classes/LongDecimal.src/M000080.html +19 -7
- data/doc/classes/LongDecimal.src/M000081.html +8 -4
- data/doc/classes/LongDecimal.src/M000082.html +5 -7
- data/doc/classes/LongDecimal.src/M000083.html +8 -6
- data/doc/classes/LongDecimal.src/M000084.html +9 -4
- data/doc/classes/LongDecimal.src/M000085.html +9 -4
- data/doc/classes/LongDecimal.src/M000086.html +4 -9
- data/doc/classes/LongDecimal.src/M000087.html +6 -8
- data/doc/classes/LongDecimal.src/M000088.html +7 -5
- data/doc/classes/LongDecimal.src/M000089.html +4 -9
- data/doc/classes/LongDecimal.src/M000090.html +4 -4
- data/doc/classes/LongDecimal.src/M000091.html +9 -4
- data/doc/classes/LongDecimal.src/M000092.html +9 -4
- data/doc/classes/LongDecimal.src/M000093.html +5 -4
- data/doc/classes/LongDecimal.src/M000094.html +7 -57
- data/doc/classes/LongDecimal.src/M000095.html +4 -4
- data/doc/classes/LongDecimal.src/M000096.html +4 -4
- data/doc/classes/LongDecimal.src/M000097.html +4 -5
- data/doc/classes/LongDecimal.src/M000098.html +4 -5
- data/doc/classes/LongDecimal.src/M000099.html +59 -5
- data/doc/classes/LongDecimal.src/M000100.html +4 -5
- data/doc/classes/LongDecimal.src/M000101.html +4 -4
- data/doc/classes/LongDecimal.src/M000104.html +19 -0
- data/doc/classes/LongDecimal.src/M000105.html +19 -0
- data/doc/classes/LongDecimal.src/M000106.html +18 -0
- data/doc/classes/LongDecimal.src/M000107.html +18 -0
- data/doc/classes/LongDecimal.src/M000108.html +18 -0
- data/doc/classes/LongDecimal.src/M000109.html +18 -0
- data/doc/classes/LongDecimalBase.html +69 -66
- data/doc/classes/LongDecimalBase.src/M000113.html +4 -9
- data/doc/classes/LongDecimalBase.src/M000114.html +4 -5
- data/doc/classes/LongDecimalBase.src/M000115.html +4 -9
- data/doc/classes/LongDecimalBase.src/M000116.html +4 -9
- data/doc/classes/LongDecimalBase.src/M000117.html +4 -5
- data/doc/classes/LongDecimalBase.src/M000118.html +23 -0
- data/doc/classes/LongDecimalBase.src/M000121.html +23 -0
- data/doc/classes/LongDecimalBase.src/M000122.html +19 -0
- data/doc/classes/LongDecimalBase.src/M000123.html +18 -0
- data/doc/classes/LongDecimalQuot.html +33 -27
- data/doc/classes/LongDecimalQuot.src/M000003.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000004.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000005.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000006.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000007.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000008.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000009.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000010.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000011.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000012.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000013.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000014.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000015.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000016.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000017.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000018.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000019.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000020.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000021.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000022.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000023.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000024.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000025.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000026.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000027.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000028.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000029.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000030.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000031.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000034.html +1 -1
- data/doc/classes/LongDecimalQuot.src/M000035.html +1 -1
- data/doc/classes/LongDecimalRoundingMode/RoundingModeClass.html +25 -24
- data/doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000156.html +22 -0
- data/doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000157.html +18 -0
- data/doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000158.html +18 -0
- data/doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000159.html +18 -0
- data/doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000161.html +25 -0
- data/doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000162.html +18 -0
- data/doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000163.html +18 -0
- data/doc/classes/LongDecimalRoundingMode/RoundingModeClass.src/M000164.html +18 -0
- data/doc/classes/LongDecimalRoundingMode/ZeroRoundingModeClass.html +178 -0
- data/doc/classes/LongDecimalRoundingMode/ZeroRoundingModeClass.src/M000165.html +25 -0
- data/doc/classes/LongDecimalRoundingMode/ZeroRoundingModeClass.src/M000166.html +18 -0
- data/doc/classes/LongDecimalRoundingMode.html +41 -4
- data/doc/classes/LongMath/CacheKey.html +10 -10
- data/doc/classes/LongMath/CacheKey.src/M000155.html +36 -0
- data/doc/classes/LongMath/CacheKey.src/M000160.html +36 -0
- data/doc/classes/LongMath.html +246 -235
- data/doc/classes/LongMath.src/M000124.html +5 -18
- data/doc/classes/LongMath.src/M000125.html +4 -18
- data/doc/classes/LongMath.src/M000126.html +4 -5
- data/doc/classes/LongMath.src/M000127.html +5 -29
- data/doc/classes/LongMath.src/M000128.html +4 -5
- data/doc/classes/LongMath.src/M000129.html +14 -53
- data/doc/classes/LongMath.src/M000130.html +17 -20
- data/doc/classes/LongMath.src/M000131.html +5 -36
- data/doc/classes/LongMath.src/M000132.html +29 -6
- data/doc/classes/LongMath.src/M000133.html +5 -33
- data/doc/classes/LongMath.src/M000134.html +47 -41
- data/doc/classes/LongMath.src/M000135.html +21 -7
- data/doc/classes/LongMath.src/M000136.html +36 -4
- data/doc/classes/LongMath.src/M000137.html +6 -4
- data/doc/classes/LongMath.src/M000138.html +31 -70
- data/doc/classes/LongMath.src/M000139.html +45 -38
- data/doc/classes/LongMath.src/M000140.html +7 -4
- data/doc/classes/LongMath.src/M000141.html +4 -4
- data/doc/classes/LongMath.src/M000142.html +4 -6
- data/doc/classes/LongMath.src/M000143.html +67 -14
- data/doc/classes/LongMath.src/M000144.html +39 -14
- data/doc/classes/LongMath.src/M000145.html +4 -44
- data/doc/classes/LongMath.src/M000146.html +5 -339
- data/doc/classes/LongMath.src/M000147.html +6 -25
- data/doc/classes/LongMath.src/M000148.html +14 -75
- data/doc/classes/LongMath.src/M000149.html +33 -0
- data/doc/classes/LongMath.src/M000150.html +58 -0
- data/doc/classes/LongMath.src/M000151.html +406 -0
- data/doc/classes/LongMath.src/M000152.html +63 -0
- data/doc/classes/LongMath.src/M000153.html +117 -0
- data/doc/classes/LongMath.src/M000154.html +150 -0
- data/doc/classes/LongMath.src/M000155.html +63 -0
- data/doc/classes/LongMath.src/M000156.html +18 -0
- data/doc/classes/LongMath.src/M000157.html +19 -0
- data/doc/classes/LongMath.src/M000158.html +18 -0
- data/doc/classes/LongMath.src/M000159.html +19 -0
- data/doc/classes/Numeric.html +23 -21
- data/doc/classes/Numeric.src/M000106.html +18 -0
- data/doc/classes/Numeric.src/M000110.html +23 -0
- data/doc/classes/Numeric.src/M000111.html +18 -0
- data/doc/classes/Rational.html +18 -16
- data/doc/classes/Rational.src/M000107.html +23 -0
- data/doc/classes/Rational.src/M000112.html +23 -0
- data/doc/created.rid +1 -1
- data/doc/dot/f_0.dot +23 -1
- data/doc/dot/f_0.png +0 -0
- data/doc/dot/m_0_0.png +0 -0
- data/doc/dot/m_0_1.dot +13 -0
- data/doc/dot/m_0_1.png +0 -0
- data/doc/files/lib/long-decimal_rb.html +15 -13
- data/doc/files/lib/long-decimal_rb.src/M000001.html +1 -1
- data/doc/files/lib/long-decimal_rb.src/M000002.html +1 -1
- data/doc/fr_class_index.html +2 -0
- data/doc/fr_method_index.html +131 -124
- data/lib/long-decimal.rb +373 -49
- data/test/testlongdecimal.rb +689 -7
- data/test/testlongdeclib.rb +51 -11
- data/test/testrandlib.rb +23 -7
- data/test/testrandom.rb +4 -4
- data/test/testrandpower.rb +11 -9
- metadata +48 -3
|
@@ -0,0 +1,406 @@
|
|
|
1
|
+
<?xml version="1.0" encoding="iso-8859-1"?>
|
|
2
|
+
<!DOCTYPE html
|
|
3
|
+
PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
|
|
4
|
+
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
|
|
5
|
+
|
|
6
|
+
<html>
|
|
7
|
+
<head>
|
|
8
|
+
<title>log_raw (LongMath)</title>
|
|
9
|
+
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1" />
|
|
10
|
+
<link rel="stylesheet" href="../.././rdoc-style.css" type="text/css" media="screen" />
|
|
11
|
+
</head>
|
|
12
|
+
<body class="standalone-code">
|
|
13
|
+
<pre><span class="ruby-comment cmt"># File lib/long-decimal.rb, line 3023</span>
|
|
14
|
+
<span class="ruby-keyword kw">def</span> <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">log_raw</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">prec</span>, <span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">mode</span>)
|
|
15
|
+
|
|
16
|
+
<span class="ruby-comment cmt"># puts("log_raw start x=#{x}")</span>
|
|
17
|
+
|
|
18
|
+
<span class="ruby-comment cmt"># we have to rely on iprec being at least 10</span>
|
|
19
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">TypeError</span>, <span class="ruby-node">"iprec=#{iprec} out of range"</span> <span class="ruby-keyword kw">unless</span> (<span class="ruby-identifier">iprec</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Fixnum</span>) <span class="ruby-operator">&&</span> <span class="ruby-identifier">iprec</span> <span class="ruby-operator">>=</span> <span class="ruby-value">10</span>
|
|
20
|
+
|
|
21
|
+
<span class="ruby-comment cmt"># dprec = [ iprec - 1, (prec + 1) << 1 ].min</span>
|
|
22
|
+
<span class="ruby-comment cmt"># dprec >= 9</span>
|
|
23
|
+
<span class="ruby-identifier">dprec</span> = <span class="ruby-identifier">iprec</span> <span class="ruby-operator">-</span> <span class="ruby-value">1</span>
|
|
24
|
+
|
|
25
|
+
<span class="ruby-comment cmt"># result is stored in y</span>
|
|
26
|
+
<span class="ruby-identifier">y</span> = <span class="ruby-value">0</span>
|
|
27
|
+
<span class="ruby-comment cmt"># sign of result</span>
|
|
28
|
+
<span class="ruby-identifier">s</span> = <span class="ruby-value">1</span>
|
|
29
|
+
<span class="ruby-comment cmt"># make sure x is >= 1</span>
|
|
30
|
+
<span class="ruby-identifier">mode1</span> = <span class="ruby-identifier">mode</span>
|
|
31
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">x</span> <span class="ruby-operator"><</span> <span class="ruby-value">1</span>) <span class="ruby-keyword kw">then</span>
|
|
32
|
+
<span class="ruby-identifier">mode1</span> = <span class="ruby-identifier">mode1</span>.<span class="ruby-identifier">minverse</span>
|
|
33
|
+
<span class="ruby-identifier">x</span> = (<span class="ruby-value">1</span> <span class="ruby-operator">/</span> <span class="ruby-identifier">x</span>).<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">mode1</span>)
|
|
34
|
+
<span class="ruby-identifier">s</span> = <span class="ruby-value">-1</span>
|
|
35
|
+
<span class="ruby-keyword kw">end</span>
|
|
36
|
+
<span class="ruby-comment cmt"># puts("log_raw prepared t=#{Time.new-t0} x=#{x} s=#{s} mode=#{mode} mode1=#{mode1}")</span>
|
|
37
|
+
|
|
38
|
+
<span class="ruby-comment cmt"># number that are beyond the usual range of Float need to be</span>
|
|
39
|
+
<span class="ruby-comment cmt"># handled specially to reduce to something expressable as Float</span>
|
|
40
|
+
<span class="ruby-identifier">exp_keys</span> = [ <span class="ruby-constant">MAX_EXP_ABLE</span>.<span class="ruby-identifier">to_i</span>, <span class="ruby-value">100</span>, <span class="ruby-value">10</span>, <span class="ruby-value">1</span> ]
|
|
41
|
+
<span class="ruby-identifier">exp_keys</span>.<span class="ruby-identifier">each</span> <span class="ruby-keyword kw">do</span> <span class="ruby-operator">|</span><span class="ruby-identifier">exp_key</span><span class="ruby-operator">|</span>
|
|
42
|
+
<span class="ruby-identifier">exp_val</span> = <span class="ruby-identifier">exp</span>(<span class="ruby-identifier">exp_key</span>, <span class="ruby-identifier">iprec</span>)
|
|
43
|
+
<span class="ruby-keyword kw">while</span> (<span class="ruby-identifier">x</span> <span class="ruby-operator">></span> <span class="ruby-identifier">exp_val</span>) <span class="ruby-keyword kw">do</span>
|
|
44
|
+
<span class="ruby-identifier">x</span> = (<span class="ruby-identifier">x</span> <span class="ruby-operator">/</span> <span class="ruby-identifier">exp_val</span>).<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">mode1</span>)
|
|
45
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">s</span> <span class="ruby-operator"><</span> <span class="ruby-value">0</span>) <span class="ruby-keyword kw">then</span>
|
|
46
|
+
<span class="ruby-identifier">y</span> <span class="ruby-operator">-=</span> <span class="ruby-identifier">exp_key</span>
|
|
47
|
+
<span class="ruby-keyword kw">else</span>
|
|
48
|
+
<span class="ruby-identifier">y</span> <span class="ruby-operator">+=</span> <span class="ruby-identifier">exp_key</span>
|
|
49
|
+
<span class="ruby-keyword kw">end</span>
|
|
50
|
+
<span class="ruby-keyword kw">end</span>
|
|
51
|
+
<span class="ruby-keyword kw">end</span>
|
|
52
|
+
<span class="ruby-comment cmt"># puts("log_raw divided t=#{Time.new-t0} x=#{x} y=#{y} s=#{s}")</span>
|
|
53
|
+
|
|
54
|
+
<span class="ruby-identifier">factor</span> = <span class="ruby-value">1</span>
|
|
55
|
+
<span class="ruby-identifier">sprec</span> = (<span class="ruby-identifier">iprec</span> <span class="ruby-operator">*</span> <span class="ruby-value">1.5</span>).<span class="ruby-identifier">round</span>
|
|
56
|
+
<span class="ruby-comment cmt"># delta = LongDecimal(1, (iprec.to_f**(1/3)).round)</span>
|
|
57
|
+
<span class="ruby-comment cmt"># delta = LongDecimal(1, (iprec.to_f**0.3).round)</span>
|
|
58
|
+
<span class="ruby-identifier">delta</span> = <span class="ruby-constant">LongDecimal</span>(<span class="ruby-value">1</span>, (<span class="ruby-identifier">iprec</span>.<span class="ruby-identifier">to_f</span><span class="ruby-operator">**</span><span class="ruby-value">0</span><span class="ruby-value">.45</span>).<span class="ruby-identifier">round</span>)
|
|
59
|
+
<span class="ruby-comment cmt"># delta = LongDecimal(1, LongMath.sqrtw(iprec))</span>
|
|
60
|
+
<span class="ruby-comment cmt"># delta = LongDecimal(1, LongMath.sqrtw(LongMath.sqrtw(iprec+1)+1))</span>
|
|
61
|
+
<span class="ruby-keyword kw">while</span> (<span class="ruby-identifier">x</span> <span class="ruby-operator">-</span> <span class="ruby-value">1</span>).<span class="ruby-identifier">abs</span> <span class="ruby-operator">></span> <span class="ruby-identifier">delta</span> <span class="ruby-keyword kw">do</span>
|
|
62
|
+
<span class="ruby-identifier">x</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">sqrt</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">sprec</span>, <span class="ruby-identifier">mode1</span>)
|
|
63
|
+
<span class="ruby-identifier">factor</span> <span class="ruby-operator">*=</span> <span class="ruby-value">2</span>
|
|
64
|
+
<span class="ruby-keyword kw">end</span>
|
|
65
|
+
<span class="ruby-comment cmt"># puts("log_raw rooted t=#{Time.new-t0} x=#{x} y=#{y} s=#{s} f=#{factor}")</span>
|
|
66
|
+
|
|
67
|
+
<span class="ruby-identifier">ss</span> = <span class="ruby-value">1</span>
|
|
68
|
+
<span class="ruby-identifier">mode2</span> = <span class="ruby-identifier">mode1</span>.<span class="ruby-identifier">ainverse</span>
|
|
69
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">x</span> <span class="ruby-operator"><</span> <span class="ruby-value">1</span>)
|
|
70
|
+
<span class="ruby-identifier">puts</span>(<span class="ruby-value str">"x<1 inverting"</span>)
|
|
71
|
+
<span class="ruby-identifier">mode2</span> = <span class="ruby-identifier">mode2</span>.<span class="ruby-identifier">ainverse</span>
|
|
72
|
+
<span class="ruby-identifier">x</span> = (<span class="ruby-value">1</span> <span class="ruby-operator">/</span> <span class="ruby-identifier">x</span>).<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">mode2</span>)
|
|
73
|
+
<span class="ruby-identifier">ss</span> = <span class="ruby-value">-1</span>
|
|
74
|
+
<span class="ruby-keyword kw">end</span>
|
|
75
|
+
<span class="ruby-comment cmt"># puts("log_raw 2nd prep t=#{Time.new-t0} x=#{x} y=#{y} s=#{s} ss=#{ss} f=#{factor}")</span>
|
|
76
|
+
|
|
77
|
+
<span class="ruby-identifier">sum</span> = <span class="ruby-value">0</span>
|
|
78
|
+
<span class="ruby-identifier">z</span> = <span class="ruby-value">1</span> <span class="ruby-operator">-</span> <span class="ruby-identifier">x</span>
|
|
79
|
+
<span class="ruby-identifier">i</span> = <span class="ruby-value">1</span>
|
|
80
|
+
<span class="ruby-identifier">p</span> = <span class="ruby-value">1</span>.<span class="ruby-identifier">to_ld</span>
|
|
81
|
+
<span class="ruby-identifier">d</span> = <span class="ruby-value">1</span>.<span class="ruby-identifier">to_ld</span>
|
|
82
|
+
<span class="ruby-keyword kw">until</span> <span class="ruby-identifier">p</span>.<span class="ruby-identifier">abs</span>.<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-identifier">dprec</span>, <span class="ruby-constant">LongDecimal</span><span class="ruby-operator">::</span><span class="ruby-constant">ROUND_DOWN</span>).<span class="ruby-identifier">zero?</span> <span class="ruby-keyword kw">do</span>
|
|
83
|
+
<span class="ruby-identifier">p</span> = (<span class="ruby-identifier">p</span> <span class="ruby-operator">*</span> <span class="ruby-identifier">z</span>).<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">mode2</span>)
|
|
84
|
+
<span class="ruby-identifier">d</span> = (<span class="ruby-identifier">p</span> <span class="ruby-operator">/</span> <span class="ruby-identifier">i</span>).<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">mode2</span>)
|
|
85
|
+
<span class="ruby-identifier">i</span> <span class="ruby-operator">+=</span> <span class="ruby-value">1</span>
|
|
86
|
+
<span class="ruby-identifier">sum</span> <span class="ruby-operator">+=</span> <span class="ruby-identifier">d</span>
|
|
87
|
+
|
|
88
|
+
<span class="ruby-keyword kw">end</span>
|
|
89
|
+
<span class="ruby-identifier">sum</span> <span class="ruby-operator">*=</span> <span class="ruby-identifier">ss</span>
|
|
90
|
+
<span class="ruby-comment cmt"># puts("log_raw added t=#{Time.new-t0} x=#{x} y=#{y} sum=#{sum} s=#{s} ss=#{ss} f=#{factor}")</span>
|
|
91
|
+
|
|
92
|
+
<span class="ruby-identifier">y</span> <span class="ruby-operator">-=</span> ((<span class="ruby-identifier">s</span> <span class="ruby-operator">*</span> <span class="ruby-identifier">factor</span>) <span class="ruby-operator">*</span> <span class="ruby-identifier">sum</span>).<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">mode</span>.<span class="ruby-identifier">ainverse</span>)
|
|
93
|
+
<span class="ruby-comment cmt"># puts("log_raw done t=#{Time.new-t0} x=#{x} y=#{y} sum=#{sum} s=#{s} ss=#{ss} f=#{factor}")</span>
|
|
94
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-identifier">y</span>
|
|
95
|
+
<span class="ruby-keyword kw">end</span>
|
|
96
|
+
|
|
97
|
+
<span class="ruby-identifier">private</span>
|
|
98
|
+
|
|
99
|
+
<span class="ruby-comment cmt">#</span>
|
|
100
|
+
<span class="ruby-comment cmt"># internal helper method for calculating the internal precision for power</span>
|
|
101
|
+
<span class="ruby-comment cmt">#</span>
|
|
102
|
+
<span class="ruby-keyword kw">def</span> <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">calc_iprec_for_power</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">y</span>, <span class="ruby-identifier">prec</span>)
|
|
103
|
+
|
|
104
|
+
<span class="ruby-identifier">logx_f</span> = <span class="ruby-keyword kw">nil</span>
|
|
105
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">x</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator"><=</span> <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">MAX_FLOATABLE</span>)
|
|
106
|
+
<span class="ruby-identifier">x_f</span> = <span class="ruby-identifier">x</span>.<span class="ruby-identifier">to_f</span>
|
|
107
|
+
<span class="ruby-identifier">logx_f</span> = <span class="ruby-constant">Math</span>.<span class="ruby-identifier">log</span>(<span class="ruby-identifier">x_f</span>.<span class="ruby-identifier">abs</span>)
|
|
108
|
+
<span class="ruby-keyword kw">else</span>
|
|
109
|
+
<span class="ruby-identifier">logx_f</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">log</span>(<span class="ruby-identifier">x</span>, <span class="ruby-value">15</span>, <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">ROUND_UP</span>)
|
|
110
|
+
<span class="ruby-keyword kw">end</span>
|
|
111
|
+
|
|
112
|
+
<span class="ruby-identifier">y_f</span> = <span class="ruby-keyword kw">nil</span>
|
|
113
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator"><=</span> <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">MAX_FLOATABLE</span>) <span class="ruby-keyword kw">then</span>
|
|
114
|
+
<span class="ruby-identifier">y_f</span> = <span class="ruby-identifier">y</span>.<span class="ruby-identifier">to_f</span>
|
|
115
|
+
<span class="ruby-keyword kw">else</span>
|
|
116
|
+
<span class="ruby-identifier">y_f</span> = <span class="ruby-identifier">y</span>.<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-value">15</span>, <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">ROUND_UP</span>)
|
|
117
|
+
<span class="ruby-keyword kw">end</span>
|
|
118
|
+
|
|
119
|
+
<span class="ruby-identifier">logx_y_f</span> = <span class="ruby-identifier">logx_f</span> <span class="ruby-operator">*</span> <span class="ruby-identifier">y_f</span>
|
|
120
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">logx_y_f</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator">></span> <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">MAX_FLOATABLE</span>) <span class="ruby-keyword kw">then</span>
|
|
121
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">ArgumentError</span>, <span class="ruby-node">"power would be way too big: y*log(x)=#{logx_y_f}"</span>;
|
|
122
|
+
<span class="ruby-keyword kw">end</span>
|
|
123
|
+
<span class="ruby-identifier">logx_y_f</span> = <span class="ruby-identifier">logx_y_f</span>.<span class="ruby-identifier">to_f</span> <span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">logx_y_f</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Float</span>
|
|
124
|
+
|
|
125
|
+
<span class="ruby-identifier">iprec_x</span> = <span class="ruby-identifier">calc_iprec_for_exp</span>(<span class="ruby-identifier">logx_y_f</span>.<span class="ruby-identifier">abs</span>.<span class="ruby-identifier">ceil</span>, <span class="ruby-identifier">prec</span>, <span class="ruby-identifier">logx_y_f</span> <span class="ruby-operator"><</span> <span class="ruby-value">0</span>)
|
|
126
|
+
<span class="ruby-identifier">iprec_y</span> = <span class="ruby-identifier">iprec_x</span>
|
|
127
|
+
<span class="ruby-identifier">iprec</span> = <span class="ruby-identifier">iprec_x</span> <span class="ruby-operator">+</span> <span class="ruby-value">2</span>
|
|
128
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">logx_f</span> <span class="ruby-operator"><</span> <span class="ruby-value">0</span>)
|
|
129
|
+
<span class="ruby-identifier">iprec_x</span> <span class="ruby-operator">-=</span> (<span class="ruby-identifier">logx_f</span><span class="ruby-operator">/</span><span class="ruby-constant">LOG10</span>).<span class="ruby-identifier">round</span>
|
|
130
|
+
<span class="ruby-keyword kw">end</span>
|
|
131
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y_f</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator"><</span> <span class="ruby-value">1</span>)
|
|
132
|
+
<span class="ruby-identifier">logy_f</span> = <span class="ruby-keyword kw">nil</span>
|
|
133
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y_f</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Float</span>) <span class="ruby-keyword kw">then</span>
|
|
134
|
+
<span class="ruby-identifier">logy_f</span> = <span class="ruby-constant">Math</span>.<span class="ruby-identifier">log</span>(<span class="ruby-identifier">y_f</span>.<span class="ruby-identifier">abs</span>)
|
|
135
|
+
<span class="ruby-keyword kw">else</span>
|
|
136
|
+
<span class="ruby-identifier">logy_f</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">log</span>(<span class="ruby-identifier">y_f</span>.<span class="ruby-identifier">abs</span>, <span class="ruby-value">15</span>, <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">ROUND_UP</span>)
|
|
137
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">logy_f</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator">></span> <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">MAX_FLOATABLE</span>) <span class="ruby-keyword kw">then</span>
|
|
138
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">ArgumentError</span>, <span class="ruby-node">"exponent would be way too big: y=#{y} logy_f=#{logy_f}"</span>;
|
|
139
|
+
<span class="ruby-keyword kw">end</span>
|
|
140
|
+
<span class="ruby-identifier">logy_f</span> = <span class="ruby-identifier">logy_f</span>.<span class="ruby-identifier">to_f</span>
|
|
141
|
+
<span class="ruby-keyword kw">end</span>
|
|
142
|
+
<span class="ruby-comment cmt"># puts("x=#{x} y=#{y} x_f=#{x_f} y_f=#{y_f} logx_f=#{logx_f} logy_f=#{logy_f} logx_y_f=#{logx_y_f}\n")</span>
|
|
143
|
+
<span class="ruby-identifier">iprec_y</span> <span class="ruby-operator">-=</span> (<span class="ruby-identifier">logy_f</span><span class="ruby-operator">/</span><span class="ruby-constant">LOG10</span>).<span class="ruby-identifier">round</span>
|
|
144
|
+
<span class="ruby-keyword kw">end</span>
|
|
145
|
+
<span class="ruby-comment cmt"># puts("x=#{x} y=#{y} x_f=#{x_f} y_f=#{y_f} logx_f=#{logx_f} logy_f=#{logy_f} logx_y_f=#{logx_y_f}\n")</span>
|
|
146
|
+
<span class="ruby-comment cmt"># puts("\niprec: x=#{x} y=#{y} iprec=#{iprec} iprec_x=#{iprec_x} iprec_y=#{iprec_y}\n")</span>
|
|
147
|
+
[ <span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">iprec_x</span>, <span class="ruby-identifier">iprec_y</span>, <span class="ruby-identifier">logx_y_f</span> ]
|
|
148
|
+
|
|
149
|
+
<span class="ruby-keyword kw">end</span>
|
|
150
|
+
|
|
151
|
+
<span class="ruby-identifier">public</span>
|
|
152
|
+
|
|
153
|
+
<span class="ruby-comment cmt">#</span>
|
|
154
|
+
<span class="ruby-comment cmt"># calc the power of x with exponent y to the given precision as</span>
|
|
155
|
+
<span class="ruby-comment cmt"># LongDecimal. Only supports values of y such that exp(y) still</span>
|
|
156
|
+
<span class="ruby-comment cmt"># fits into a float (y <= 709)</span>
|
|
157
|
+
<span class="ruby-comment cmt">#</span>
|
|
158
|
+
<span class="ruby-keyword kw">def</span> <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">power</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">y</span>, <span class="ruby-identifier">prec</span>, <span class="ruby-identifier">mode</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">standard_mode</span>)
|
|
159
|
+
|
|
160
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">TypeError</span>, <span class="ruby-node">"x=#{x} must be numeric"</span> <span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">x</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Numeric</span>
|
|
161
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">TypeError</span>, <span class="ruby-node">"y=#{y} must be numeric"</span> <span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">y</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Numeric</span>
|
|
162
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">TypeError</span>, <span class="ruby-node">"x=#{x.inspect} must not be greater #{MAX_FLOATABLE}"</span> <span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">x</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator"><=</span> <span class="ruby-constant">MAX_FLOATABLE</span>
|
|
163
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">TypeError</span>, <span class="ruby-node">"y=#{y.inspect} must not be greater #{MAX_FLOATABLE}"</span> <span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">y</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator"><=</span> <span class="ruby-constant">MAX_FLOATABLE</span>
|
|
164
|
+
<span class="ruby-comment cmt"># raise TypeError, "y=#{y.inspect} must not be greater #{MAX_EXP_ABLE}" unless y <= MAX_EXP_ABLE</span>
|
|
165
|
+
<span class="ruby-comment cmt"># raise TypeError, "x=#{x.inspect} must not negative" unless x >= 0 || (y.kind_of? Integer) || (y.kind_of? LongDecimalBase) && y.is_int?</span>
|
|
166
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">TypeError</span>, <span class="ruby-node">"x=#{x.inspect} must not negative"</span> <span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">x</span> <span class="ruby-operator">>=</span> <span class="ruby-value">0</span>
|
|
167
|
+
<span class="ruby-identifier">check_is_prec</span>(<span class="ruby-identifier">prec</span>, <span class="ruby-value str">"prec"</span>)
|
|
168
|
+
<span class="ruby-identifier">check_is_mode</span>(<span class="ruby-identifier">mode</span>, <span class="ruby-value str">"mode"</span>)
|
|
169
|
+
|
|
170
|
+
<span class="ruby-comment cmt"># handle the special cases where base or exponent are 0 or 1 explicitely</span>
|
|
171
|
+
<span class="ruby-keyword kw">if</span> <span class="ruby-identifier">y</span>.<span class="ruby-identifier">zero?</span> <span class="ruby-keyword kw">then</span>
|
|
172
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-constant">LongDecimal</span>.<span class="ruby-identifier">one!</span>(<span class="ruby-identifier">prec</span>)
|
|
173
|
+
<span class="ruby-keyword kw">elsif</span> <span class="ruby-identifier">x</span>.<span class="ruby-identifier">zero?</span> <span class="ruby-keyword kw">then</span>
|
|
174
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-constant">LongDecimal</span>.<span class="ruby-identifier">zero!</span>(<span class="ruby-identifier">prec</span>)
|
|
175
|
+
<span class="ruby-keyword kw">elsif</span> <span class="ruby-identifier">y</span>.<span class="ruby-identifier">one?</span> <span class="ruby-keyword kw">then</span>
|
|
176
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-identifier">x</span>.<span class="ruby-identifier">to_ld</span>(<span class="ruby-identifier">prec</span>, <span class="ruby-identifier">mode</span>)
|
|
177
|
+
<span class="ruby-keyword kw">elsif</span> <span class="ruby-identifier">x</span>.<span class="ruby-identifier">one?</span> <span class="ruby-keyword kw">then</span>
|
|
178
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-constant">LongDecimal</span>.<span class="ruby-identifier">one!</span>(<span class="ruby-identifier">prec</span>)
|
|
179
|
+
<span class="ruby-keyword kw">end</span>
|
|
180
|
+
|
|
181
|
+
<span class="ruby-comment cmt"># els</span>
|
|
182
|
+
<span class="ruby-comment cmt"># could be result with our precision</span>
|
|
183
|
+
<span class="ruby-comment cmt"># x ** y <= 10**-s/2 <=> y * log(x) <= -s log(10) - log(2)</span>
|
|
184
|
+
|
|
185
|
+
<span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">iprec_x</span>, <span class="ruby-identifier">iprec_y</span>, <span class="ruby-identifier">logx_y_f</span> = <span class="ruby-identifier">calc_iprec_for_power</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">y</span>, <span class="ruby-identifier">prec</span>)
|
|
186
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">x</span> <span class="ruby-operator"><</span> <span class="ruby-value">1</span> <span class="ruby-operator">&&</span> <span class="ruby-identifier">y</span> <span class="ruby-operator">></span> <span class="ruby-value">0</span> <span class="ruby-operator">||</span> <span class="ruby-identifier">x</span> <span class="ruby-operator">></span> <span class="ruby-value">1</span> <span class="ruby-operator">&&</span> <span class="ruby-identifier">y</span> <span class="ruby-operator"><</span> <span class="ruby-value">0</span>) <span class="ruby-keyword kw">then</span>
|
|
187
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">logx_y_f</span> <span class="ruby-operator"><=</span> <span class="ruby-operator">-</span> <span class="ruby-identifier">prec</span> <span class="ruby-operator">*</span> <span class="ruby-constant">LOG10</span> <span class="ruby-operator">-</span> <span class="ruby-constant">LOG2</span>) <span class="ruby-keyword kw">then</span>
|
|
188
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-constant">LongDecimal</span>.<span class="ruby-identifier">zero!</span>(<span class="ruby-identifier">prec</span>)
|
|
189
|
+
<span class="ruby-keyword kw">end</span>
|
|
190
|
+
<span class="ruby-keyword kw">end</span>
|
|
191
|
+
<span class="ruby-comment cmt"># puts("x=#{x} y=#{y} iprec=#{iprec} iprec_x=#{iprec_x} iprec_y=#{iprec_y} prec=#{prec}")</span>
|
|
192
|
+
|
|
193
|
+
<span class="ruby-keyword kw">unless</span> (<span class="ruby-identifier">x</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimalBase</span>) <span class="ruby-operator">||</span> (<span class="ruby-identifier">x</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Integer</span>)
|
|
194
|
+
<span class="ruby-identifier">x</span> = <span class="ruby-identifier">x</span>.<span class="ruby-identifier">to_ld</span>(<span class="ruby-identifier">iprec_x</span>, <span class="ruby-identifier">mode</span>)
|
|
195
|
+
<span class="ruby-keyword kw">end</span>
|
|
196
|
+
<span class="ruby-keyword kw">unless</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimalBase</span>) <span class="ruby-operator">||</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Integer</span>)
|
|
197
|
+
<span class="ruby-identifier">y</span> = <span class="ruby-identifier">y</span>.<span class="ruby-identifier">to_ld</span>(<span class="ruby-identifier">iprec_y</span>, <span class="ruby-identifier">mode</span>)
|
|
198
|
+
<span class="ruby-keyword kw">end</span>
|
|
199
|
+
|
|
200
|
+
<span class="ruby-comment cmt"># try shortcut if exponent is an integer</span>
|
|
201
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimalBase</span>) <span class="ruby-operator">&&</span> <span class="ruby-identifier">y</span>.<span class="ruby-identifier">is_int?</span> <span class="ruby-keyword kw">then</span>
|
|
202
|
+
<span class="ruby-identifier">y</span> = <span class="ruby-identifier">y</span>.<span class="ruby-identifier">to_i</span>
|
|
203
|
+
<span class="ruby-keyword kw">end</span>
|
|
204
|
+
<span class="ruby-keyword kw">unless</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Integer</span>)
|
|
205
|
+
<span class="ruby-identifier">y2</span> = <span class="ruby-identifier">y</span><span class="ruby-operator">*</span><span class="ruby-value">2</span>
|
|
206
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y2</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimalBase</span>) <span class="ruby-operator">&&</span> <span class="ruby-identifier">y2</span>.<span class="ruby-identifier">is_int?</span> <span class="ruby-keyword kw">then</span>
|
|
207
|
+
<span class="ruby-identifier">y2</span> = <span class="ruby-identifier">y2</span>.<span class="ruby-identifier">to_i</span>
|
|
208
|
+
<span class="ruby-identifier">puts</span>(<span class="ruby-node">"y2=#{y2}"</span>)
|
|
209
|
+
<span class="ruby-keyword kw">end</span>
|
|
210
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y2</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Integer</span>)
|
|
211
|
+
<span class="ruby-identifier">x</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">sqrt</span>(<span class="ruby-identifier">x</span>, <span class="ruby-value">2</span><span class="ruby-operator">*</span><span class="ruby-identifier">iprec_x</span>, <span class="ruby-identifier">mode</span>)
|
|
212
|
+
<span class="ruby-identifier">y</span> = <span class="ruby-identifier">y2</span>
|
|
213
|
+
<span class="ruby-keyword kw">end</span>
|
|
214
|
+
<span class="ruby-keyword kw">end</span>
|
|
215
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Integer</span>)
|
|
216
|
+
<span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">x</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimal</span>
|
|
217
|
+
<span class="ruby-comment cmt"># x = x.to_ld(prec)</span>
|
|
218
|
+
<span class="ruby-identifier">x</span> = <span class="ruby-identifier">x</span>.<span class="ruby-identifier">to_ld</span>(<span class="ruby-identifier">iprec_x</span>)
|
|
219
|
+
<span class="ruby-keyword kw">end</span>
|
|
220
|
+
<span class="ruby-comment cmt"># z = x ** y</span>
|
|
221
|
+
<span class="ruby-identifier">z</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">ipower</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">y</span>, <span class="ruby-value">2</span><span class="ruby-operator">*</span><span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">mode</span>)
|
|
222
|
+
<span class="ruby-comment cmt"># puts("x=#{x} y=#{y} z=#{z} y int")</span>
|
|
223
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-identifier">z</span>.<span class="ruby-identifier">to_ld</span>(<span class="ruby-identifier">prec</span>, <span class="ruby-identifier">mode</span>)
|
|
224
|
+
<span class="ruby-keyword kw">end</span>
|
|
225
|
+
|
|
226
|
+
<span class="ruby-comment cmt"># it can be assumed that the exponent is not an integer, so it should</span>
|
|
227
|
+
<span class="ruby-comment cmt"># be converted into LongDecimal</span>
|
|
228
|
+
<span class="ruby-keyword kw">unless</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimal</span>)
|
|
229
|
+
<span class="ruby-comment cmt"># y = y.to_ld(prec, mode)</span>
|
|
230
|
+
<span class="ruby-identifier">y</span> = <span class="ruby-identifier">y</span>.<span class="ruby-identifier">to_ld</span>(<span class="ruby-identifier">iprec_y</span>, <span class="ruby-identifier">mode</span>)
|
|
231
|
+
<span class="ruby-keyword kw">end</span>
|
|
232
|
+
|
|
233
|
+
<span class="ruby-comment cmt"># if x < 1 && y < 0 then</span>
|
|
234
|
+
<span class="ruby-comment cmt"># working with x < 1 should be improved, less precision needed</span>
|
|
235
|
+
<span class="ruby-keyword kw">if</span> <span class="ruby-identifier">x</span> <span class="ruby-operator"><</span> <span class="ruby-value">1</span> <span class="ruby-keyword kw">then</span>
|
|
236
|
+
<span class="ruby-comment cmt"># since we do not allow x < 0 and we have handled x = 0 already,</span>
|
|
237
|
+
<span class="ruby-comment cmt"># we can be sure that x is no integer, so it has been converted</span>
|
|
238
|
+
<span class="ruby-comment cmt"># if necessary to LongDecimalBase</span>
|
|
239
|
+
<span class="ruby-identifier">y</span> = <span class="ruby-operator">-</span><span class="ruby-identifier">y</span>
|
|
240
|
+
<span class="ruby-identifier">x</span> = (<span class="ruby-value">1</span><span class="ruby-operator">/</span><span class="ruby-identifier">x</span>).<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-identifier">iprec_x</span><span class="ruby-operator">*</span><span class="ruby-value">2</span>, <span class="ruby-identifier">mode</span>)
|
|
241
|
+
<span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">iprec_x</span>, <span class="ruby-identifier">iprec_y</span> = <span class="ruby-identifier">calc_iprec_for_power</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">y</span>, <span class="ruby-identifier">prec</span>)
|
|
242
|
+
<span class="ruby-keyword kw">end</span>
|
|
243
|
+
|
|
244
|
+
<span class="ruby-comment cmt"># exponent is split in two parts, an integer part and a</span>
|
|
245
|
+
<span class="ruby-comment cmt"># LongDecimal with absolute value <= 0.5</span>
|
|
246
|
+
<span class="ruby-identifier">y0</span> = <span class="ruby-identifier">y</span>.<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-value">0</span>, <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">standard_imode</span>).<span class="ruby-identifier">to_i</span>
|
|
247
|
+
<span class="ruby-comment cmt"># z0 = x**y0</span>
|
|
248
|
+
<span class="ruby-identifier">z0</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">ipower</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">y0</span>, <span class="ruby-value">2</span><span class="ruby-operator">*</span><span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">mode</span>)
|
|
249
|
+
<span class="ruby-identifier">y1</span> = <span class="ruby-identifier">y</span> <span class="ruby-operator">-</span> <span class="ruby-identifier">y0</span>
|
|
250
|
+
<span class="ruby-identifier">prec_extra</span> = <span class="ruby-value">0</span>
|
|
251
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y0</span> <span class="ruby-operator">></span> <span class="ruby-value">0</span>)
|
|
252
|
+
<span class="ruby-identifier">prec_extra</span> = (<span class="ruby-identifier">y0</span><span class="ruby-operator">*</span><span class="ruby-constant">Math</span>.<span class="ruby-identifier">log10</span>(<span class="ruby-identifier">x</span>.<span class="ruby-identifier">to_f</span>).<span class="ruby-identifier">abs</span>).<span class="ruby-identifier">ceil</span>
|
|
253
|
+
<span class="ruby-keyword kw">end</span>
|
|
254
|
+
<span class="ruby-comment cmt"># z1 = LongMath.power_internal(x, y1, prec+prec_extra , mode)</span>
|
|
255
|
+
<span class="ruby-identifier">z1</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">power_internal</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">y1</span>, <span class="ruby-identifier">prec</span><span class="ruby-operator">+</span><span class="ruby-identifier">prec_extra</span> <span class="ruby-operator">+</span> <span class="ruby-value">4</span>, <span class="ruby-identifier">mode</span>)
|
|
256
|
+
<span class="ruby-identifier">z</span> = <span class="ruby-identifier">z0</span> <span class="ruby-operator">*</span> <span class="ruby-identifier">z1</span>
|
|
257
|
+
<span class="ruby-comment cmt"># puts("x=#{x} y=#{y} z=#{z} y not int")</span>
|
|
258
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-identifier">z</span>.<span class="ruby-identifier">to_ld</span>(<span class="ruby-identifier">prec</span>, <span class="ruby-identifier">mode</span>)
|
|
259
|
+
<span class="ruby-keyword kw">end</span>
|
|
260
|
+
|
|
261
|
+
<span class="ruby-comment cmt">#</span>
|
|
262
|
+
<span class="ruby-comment cmt"># internal functionality to calculate the y-th power of x assuming</span>
|
|
263
|
+
<span class="ruby-comment cmt"># that y is an integer</span>
|
|
264
|
+
<span class="ruby-comment cmt"># prec is a hint on how much internal precision is needed at most</span>
|
|
265
|
+
<span class="ruby-comment cmt"># final rounding is left to the caller</span>
|
|
266
|
+
<span class="ruby-comment cmt">#</span>
|
|
267
|
+
<span class="ruby-keyword kw">def</span> <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">ipower</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">y</span>, <span class="ruby-identifier">prec</span>, <span class="ruby-identifier">mode</span>)
|
|
268
|
+
|
|
269
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">TypeError</span>, <span class="ruby-node">"x=#{x} must be numeric"</span> <span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">x</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Numeric</span>
|
|
270
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">TypeError</span>, <span class="ruby-node">"y=#{y} must be integer"</span> <span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">y</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Integer</span>
|
|
271
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">TypeError</span>, <span class="ruby-node">"x=#{x.inspect} must not be greater #{MAX_FLOATABLE}"</span> <span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">x</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator"><=</span> <span class="ruby-constant">MAX_FLOATABLE</span>
|
|
272
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">TypeError</span>, <span class="ruby-node">"y=#{y.inspect} must not be greater #{MAX_FLOATABLE}"</span> <span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">y</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator"><=</span> <span class="ruby-constant">MAX_FLOATABLE</span>
|
|
273
|
+
<span class="ruby-identifier">check_is_prec</span>(<span class="ruby-identifier">prec</span>, <span class="ruby-value str">"prec"</span>)
|
|
274
|
+
<span class="ruby-identifier">check_is_mode</span>(<span class="ruby-identifier">mode</span>, <span class="ruby-value str">"mode"</span>)
|
|
275
|
+
|
|
276
|
+
<span class="ruby-identifier">cnt</span> = <span class="ruby-value">0</span>
|
|
277
|
+
|
|
278
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">zero?</span>)
|
|
279
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-value">1</span>
|
|
280
|
+
<span class="ruby-keyword kw">elsif</span> <span class="ruby-operator">!</span> (<span class="ruby-identifier">x</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimalBase</span>) <span class="ruby-operator">||</span> <span class="ruby-identifier">x</span>.<span class="ruby-identifier">scale</span> <span class="ruby-operator">*</span> <span class="ruby-identifier">y</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator"><=</span> <span class="ruby-identifier">prec</span>
|
|
281
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-identifier">x</span> <span class="ruby-operator">**</span> <span class="ruby-identifier">y</span>
|
|
282
|
+
<span class="ruby-keyword kw">elsif</span> (<span class="ruby-identifier">y</span> <span class="ruby-operator"><</span> <span class="ruby-value">0</span>)
|
|
283
|
+
<span class="ruby-identifier">l</span> = <span class="ruby-constant">Math</span>.<span class="ruby-identifier">log10</span>(<span class="ruby-identifier">x</span>.<span class="ruby-identifier">abs</span>.<span class="ruby-identifier">to_f</span>)
|
|
284
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">l</span> <span class="ruby-operator">></span> <span class="ruby-value">0</span>)
|
|
285
|
+
<span class="ruby-identifier">prec</span> <span class="ruby-operator">+=</span> (<span class="ruby-value">2</span><span class="ruby-operator">*</span><span class="ruby-identifier">l</span>).<span class="ruby-identifier">ceil</span>
|
|
286
|
+
<span class="ruby-keyword kw">end</span>
|
|
287
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-value">1</span><span class="ruby-operator">/</span><span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">ipower</span>(<span class="ruby-identifier">x</span>, <span class="ruby-operator">-</span><span class="ruby-identifier">y</span>, <span class="ruby-identifier">prec</span>, <span class="ruby-identifier">mode</span>)
|
|
288
|
+
<span class="ruby-keyword kw">else</span>
|
|
289
|
+
<span class="ruby-comment cmt"># y > 0</span>
|
|
290
|
+
<span class="ruby-comment cmt"># puts("ipower y>0 x=#{x} y=#{y} prec=#{prec}")</span>
|
|
291
|
+
<span class="ruby-identifier">z</span> = <span class="ruby-identifier">x</span>
|
|
292
|
+
<span class="ruby-keyword kw">while</span> <span class="ruby-keyword kw">true</span> <span class="ruby-keyword kw">do</span>
|
|
293
|
+
|
|
294
|
+
<span class="ruby-identifier">cnt</span><span class="ruby-operator">+</span><span class="ruby-operator">+</span>
|
|
295
|
+
<span class="ruby-identifier">y</span> <span class="ruby-operator">-=</span> <span class="ruby-value">1</span>
|
|
296
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">zero?</span>)
|
|
297
|
+
<span class="ruby-keyword kw">break</span>
|
|
298
|
+
<span class="ruby-keyword kw">end</span>
|
|
299
|
+
<span class="ruby-keyword kw">while</span> (<span class="ruby-identifier">y</span> <span class="ruby-operator">&</span> <span class="ruby-value">0x01</span>) <span class="ruby-operator">==</span> <span class="ruby-value">0</span> <span class="ruby-keyword kw">do</span>
|
|
300
|
+
|
|
301
|
+
<span class="ruby-identifier">cnt</span><span class="ruby-operator">+</span><span class="ruby-operator">+</span>
|
|
302
|
+
<span class="ruby-identifier">y</span> = <span class="ruby-identifier">y</span> <span class="ruby-operator">>></span> <span class="ruby-value">1</span>
|
|
303
|
+
<span class="ruby-identifier">x</span> = (<span class="ruby-identifier">x</span><span class="ruby-operator">*</span><span class="ruby-identifier">x</span>)
|
|
304
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">x</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimalBase</span>)
|
|
305
|
+
<span class="ruby-identifier">x</span> = <span class="ruby-identifier">x</span>.<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-identifier">prec</span>, <span class="ruby-identifier">mode</span>)
|
|
306
|
+
<span class="ruby-keyword kw">end</span>
|
|
307
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">cnt</span> <span class="ruby-operator">></span> <span class="ruby-value">1000</span>)
|
|
308
|
+
<span class="ruby-identifier">puts</span>(<span class="ruby-node">"ipower x=#{x} y=#{y} cnt=#{cnt} z=#{z}"</span>)
|
|
309
|
+
<span class="ruby-identifier">cnt</span> = <span class="ruby-value">0</span>
|
|
310
|
+
<span class="ruby-keyword kw">end</span>
|
|
311
|
+
|
|
312
|
+
<span class="ruby-keyword kw">end</span>
|
|
313
|
+
<span class="ruby-identifier">z</span> = (<span class="ruby-identifier">z</span><span class="ruby-operator">*</span><span class="ruby-identifier">x</span>)
|
|
314
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">z</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimalBase</span>)
|
|
315
|
+
<span class="ruby-identifier">z</span> = <span class="ruby-identifier">z</span>.<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-identifier">prec</span>, <span class="ruby-identifier">mode</span>)
|
|
316
|
+
<span class="ruby-keyword kw">end</span>
|
|
317
|
+
|
|
318
|
+
<span class="ruby-keyword kw">end</span>
|
|
319
|
+
<span class="ruby-keyword kw">end</span>
|
|
320
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-identifier">z</span>
|
|
321
|
+
<span class="ruby-keyword kw">end</span>
|
|
322
|
+
|
|
323
|
+
<span class="ruby-comment cmt">#</span>
|
|
324
|
+
<span class="ruby-comment cmt"># internal functionality of exp. exposes some more parameters, that</span>
|
|
325
|
+
<span class="ruby-comment cmt"># should usually be set to defaut values, in order to allow better testing.</span>
|
|
326
|
+
<span class="ruby-comment cmt"># do not actually call this method unless you are testing exp.</span>
|
|
327
|
+
<span class="ruby-comment cmt"># create a bug report, if the default settings for the parameters do</span>
|
|
328
|
+
<span class="ruby-comment cmt"># not work correctly</span>
|
|
329
|
+
<span class="ruby-comment cmt">#</span>
|
|
330
|
+
<span class="ruby-keyword kw">def</span> <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">power_internal</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">y</span>, <span class="ruby-identifier">prec</span> = <span class="ruby-keyword kw">nil</span>, <span class="ruby-identifier">final_mode</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">standard_mode</span>, <span class="ruby-identifier">iprec</span> = <span class="ruby-keyword kw">nil</span>, <span class="ruby-identifier">mode</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">standard_imode</span>)
|
|
331
|
+
|
|
332
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">prec</span>.<span class="ruby-identifier">nil?</span>) <span class="ruby-keyword kw">then</span>
|
|
333
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">x</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimalBase</span>) <span class="ruby-operator">&&</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimalBase</span>)
|
|
334
|
+
<span class="ruby-identifier">prec</span> = [<span class="ruby-identifier">x</span>.<span class="ruby-identifier">scale</span>, <span class="ruby-identifier">y</span>.<span class="ruby-identifier">scale</span>].<span class="ruby-identifier">max</span>
|
|
335
|
+
<span class="ruby-keyword kw">elsif</span> (<span class="ruby-identifier">x</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimalBase</span>)
|
|
336
|
+
<span class="ruby-identifier">prec</span> = <span class="ruby-identifier">x</span>.<span class="ruby-identifier">scale</span>
|
|
337
|
+
<span class="ruby-keyword kw">elsif</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimalBase</span>)
|
|
338
|
+
<span class="ruby-identifier">prec</span> = <span class="ruby-identifier">y</span>.<span class="ruby-identifier">scale</span>
|
|
339
|
+
<span class="ruby-keyword kw">else</span>
|
|
340
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">ArgumentError</span>, <span class="ruby-node">"precision must be supplied either as precision of x=#{x} or explicitely"</span>
|
|
341
|
+
<span class="ruby-keyword kw">end</span>
|
|
342
|
+
<span class="ruby-keyword kw">end</span>
|
|
343
|
+
<span class="ruby-identifier">check_is_prec</span>(<span class="ruby-identifier">prec</span>, <span class="ruby-value str">"prec"</span>)
|
|
344
|
+
|
|
345
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">final_mode</span>.<span class="ruby-identifier">nil?</span>)
|
|
346
|
+
<span class="ruby-identifier">final_mode</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">standard_mode</span>
|
|
347
|
+
<span class="ruby-keyword kw">end</span>
|
|
348
|
+
<span class="ruby-identifier">check_is_mode</span>(<span class="ruby-identifier">final_mode</span>, <span class="ruby-value str">"final_mode"</span>)
|
|
349
|
+
<span class="ruby-identifier">check_is_mode</span>(<span class="ruby-identifier">mode</span>, <span class="ruby-value str">"mode"</span>)
|
|
350
|
+
|
|
351
|
+
<span class="ruby-keyword kw">if</span> <span class="ruby-identifier">y</span>.<span class="ruby-identifier">zero?</span> <span class="ruby-keyword kw">then</span>
|
|
352
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-constant">LongDecimal</span>.<span class="ruby-identifier">one!</span>(<span class="ruby-identifier">prec</span>)
|
|
353
|
+
<span class="ruby-keyword kw">elsif</span> <span class="ruby-identifier">x</span>.<span class="ruby-identifier">zero?</span> <span class="ruby-keyword kw">then</span>
|
|
354
|
+
<span class="ruby-keyword kw">return</span> <span class="ruby-constant">LongDecimal</span>.<span class="ruby-identifier">zero!</span>(<span class="ruby-identifier">prec</span>)
|
|
355
|
+
<span class="ruby-keyword kw">end</span>
|
|
356
|
+
|
|
357
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">iprec</span>.<span class="ruby-identifier">nil?</span>) <span class="ruby-keyword kw">then</span>
|
|
358
|
+
<span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">iprec_x</span>, <span class="ruby-identifier">iprec_y</span> = <span class="ruby-identifier">calc_iprec_for_power</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">y</span>, <span class="ruby-identifier">prec</span>)
|
|
359
|
+
<span class="ruby-keyword kw">end</span>
|
|
360
|
+
<span class="ruby-keyword kw">unless</span> (<span class="ruby-identifier">x</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimal</span>)
|
|
361
|
+
<span class="ruby-comment cmt"># x = x.to_ld(iprec, mode)</span>
|
|
362
|
+
<span class="ruby-identifier">x</span> = <span class="ruby-identifier">x</span>.<span class="ruby-identifier">to_ld</span>(<span class="ruby-identifier">iprec_x</span>, <span class="ruby-identifier">mode</span>)
|
|
363
|
+
<span class="ruby-keyword kw">end</span>
|
|
364
|
+
<span class="ruby-keyword kw">unless</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">LongDecimal</span>)
|
|
365
|
+
<span class="ruby-comment cmt"># y = y.to_ld(iprec, mode)</span>
|
|
366
|
+
<span class="ruby-identifier">y</span> = <span class="ruby-identifier">y</span>.<span class="ruby-identifier">to_ld</span>(<span class="ruby-identifier">iprec_y</span>, <span class="ruby-identifier">mode</span>)
|
|
367
|
+
<span class="ruby-keyword kw">end</span>
|
|
368
|
+
|
|
369
|
+
<span class="ruby-comment cmt"># logx = log(x, iprec, mode)</span>
|
|
370
|
+
<span class="ruby-identifier">logx</span> = <span class="ruby-identifier">log</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">iprec</span> <span class="ruby-operator">+</span> <span class="ruby-value">20</span>, <span class="ruby-identifier">mode</span>)
|
|
371
|
+
<span class="ruby-identifier">logx_y</span> = <span class="ruby-identifier">logx</span><span class="ruby-operator">*</span><span class="ruby-identifier">y</span>
|
|
372
|
+
<span class="ruby-comment cmt"># xy = exp_internal(logx_y, prec + 1, mode)</span>
|
|
373
|
+
<span class="ruby-comment cmt"># xy = exp_internal(logx_y, prec + 4, mode)</span>
|
|
374
|
+
<span class="ruby-identifier">xy</span> = <span class="ruby-identifier">exp_internal</span>(<span class="ruby-identifier">logx_y</span>, <span class="ruby-identifier">prec</span> <span class="ruby-operator">+</span> <span class="ruby-value">3</span>, <span class="ruby-identifier">mode</span>)
|
|
375
|
+
<span class="ruby-identifier">xy</span>.<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-identifier">prec</span>, <span class="ruby-identifier">final_mode</span>)
|
|
376
|
+
|
|
377
|
+
<span class="ruby-keyword kw">end</span> <span class="ruby-comment cmt"># power_internal</span>
|
|
378
|
+
|
|
379
|
+
<span class="ruby-ivar">@@standard_mode</span> = <span class="ruby-constant">ROUND_HALF_UP</span>
|
|
380
|
+
|
|
381
|
+
<span class="ruby-keyword kw">def</span> <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">standard_mode</span>
|
|
382
|
+
<span class="ruby-ivar">@@standard_mode</span>
|
|
383
|
+
<span class="ruby-keyword kw">end</span>
|
|
384
|
+
|
|
385
|
+
<span class="ruby-keyword kw">def</span> <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">standard_mode=</span>(<span class="ruby-identifier">x</span>)
|
|
386
|
+
<span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">check_is_mode</span>
|
|
387
|
+
<span class="ruby-ivar">@@standard_mode</span> = <span class="ruby-identifier">x</span>
|
|
388
|
+
<span class="ruby-keyword kw">end</span>
|
|
389
|
+
|
|
390
|
+
<span class="ruby-ivar">@@standard_imode</span> = <span class="ruby-constant">ROUND_HALF_EVEN</span>
|
|
391
|
+
|
|
392
|
+
<span class="ruby-keyword kw">def</span> <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">standard_imode</span>
|
|
393
|
+
<span class="ruby-ivar">@@standard_imode</span>
|
|
394
|
+
<span class="ruby-keyword kw">end</span>
|
|
395
|
+
|
|
396
|
+
<span class="ruby-keyword kw">def</span> <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">standard_imode=</span>(<span class="ruby-identifier">x</span>)
|
|
397
|
+
<span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">check_is_mode</span>
|
|
398
|
+
<span class="ruby-ivar">@@standard_imode</span> = <span class="ruby-identifier">x</span>
|
|
399
|
+
<span class="ruby-keyword kw">end</span>
|
|
400
|
+
|
|
401
|
+
<span class="ruby-keyword kw">end</span> <span class="ruby-comment cmt"># LongMath</span>
|
|
402
|
+
|
|
403
|
+
<span class="ruby-comment cmt"># end of file long-decimal.rb</span>
|
|
404
|
+
</pre>
|
|
405
|
+
</body>
|
|
406
|
+
</html>
|
|
@@ -0,0 +1,63 @@
|
|
|
1
|
+
<?xml version="1.0" encoding="iso-8859-1"?>
|
|
2
|
+
<!DOCTYPE html
|
|
3
|
+
PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
|
|
4
|
+
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
|
|
5
|
+
|
|
6
|
+
<html>
|
|
7
|
+
<head>
|
|
8
|
+
<title>calc_iprec_for_power (LongMath)</title>
|
|
9
|
+
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1" />
|
|
10
|
+
<link rel="stylesheet" href="../.././rdoc-style.css" type="text/css" media="screen" />
|
|
11
|
+
</head>
|
|
12
|
+
<body class="standalone-code">
|
|
13
|
+
<pre><span class="ruby-comment cmt"># File lib/long-decimal.rb, line 3111</span>
|
|
14
|
+
<span class="ruby-keyword kw">def</span> <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">calc_iprec_for_power</span>(<span class="ruby-identifier">x</span>, <span class="ruby-identifier">y</span>, <span class="ruby-identifier">prec</span>)
|
|
15
|
+
|
|
16
|
+
<span class="ruby-identifier">logx_f</span> = <span class="ruby-keyword kw">nil</span>
|
|
17
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">x</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator"><=</span> <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">MAX_FLOATABLE</span>)
|
|
18
|
+
<span class="ruby-identifier">x_f</span> = <span class="ruby-identifier">x</span>.<span class="ruby-identifier">to_f</span>
|
|
19
|
+
<span class="ruby-identifier">logx_f</span> = <span class="ruby-constant">Math</span>.<span class="ruby-identifier">log</span>(<span class="ruby-identifier">x_f</span>.<span class="ruby-identifier">abs</span>)
|
|
20
|
+
<span class="ruby-keyword kw">else</span>
|
|
21
|
+
<span class="ruby-identifier">logx_f</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">log</span>(<span class="ruby-identifier">x</span>, <span class="ruby-value">15</span>, <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">ROUND_UP</span>)
|
|
22
|
+
<span class="ruby-keyword kw">end</span>
|
|
23
|
+
|
|
24
|
+
<span class="ruby-identifier">y_f</span> = <span class="ruby-keyword kw">nil</span>
|
|
25
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator"><=</span> <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">MAX_FLOATABLE</span>) <span class="ruby-keyword kw">then</span>
|
|
26
|
+
<span class="ruby-identifier">y_f</span> = <span class="ruby-identifier">y</span>.<span class="ruby-identifier">to_f</span>
|
|
27
|
+
<span class="ruby-keyword kw">else</span>
|
|
28
|
+
<span class="ruby-identifier">y_f</span> = <span class="ruby-identifier">y</span>.<span class="ruby-identifier">round_to_scale</span>(<span class="ruby-value">15</span>, <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">ROUND_UP</span>)
|
|
29
|
+
<span class="ruby-keyword kw">end</span>
|
|
30
|
+
|
|
31
|
+
<span class="ruby-identifier">logx_y_f</span> = <span class="ruby-identifier">logx_f</span> <span class="ruby-operator">*</span> <span class="ruby-identifier">y_f</span>
|
|
32
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">logx_y_f</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator">></span> <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">MAX_FLOATABLE</span>) <span class="ruby-keyword kw">then</span>
|
|
33
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">ArgumentError</span>, <span class="ruby-node">"power would be way too big: y*log(x)=#{logx_y_f}"</span>;
|
|
34
|
+
<span class="ruby-keyword kw">end</span>
|
|
35
|
+
<span class="ruby-identifier">logx_y_f</span> = <span class="ruby-identifier">logx_y_f</span>.<span class="ruby-identifier">to_f</span> <span class="ruby-keyword kw">unless</span> <span class="ruby-identifier">logx_y_f</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Float</span>
|
|
36
|
+
|
|
37
|
+
<span class="ruby-identifier">iprec_x</span> = <span class="ruby-identifier">calc_iprec_for_exp</span>(<span class="ruby-identifier">logx_y_f</span>.<span class="ruby-identifier">abs</span>.<span class="ruby-identifier">ceil</span>, <span class="ruby-identifier">prec</span>, <span class="ruby-identifier">logx_y_f</span> <span class="ruby-operator"><</span> <span class="ruby-value">0</span>)
|
|
38
|
+
<span class="ruby-identifier">iprec_y</span> = <span class="ruby-identifier">iprec_x</span>
|
|
39
|
+
<span class="ruby-identifier">iprec</span> = <span class="ruby-identifier">iprec_x</span> <span class="ruby-operator">+</span> <span class="ruby-value">2</span>
|
|
40
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">logx_f</span> <span class="ruby-operator"><</span> <span class="ruby-value">0</span>)
|
|
41
|
+
<span class="ruby-identifier">iprec_x</span> <span class="ruby-operator">-=</span> (<span class="ruby-identifier">logx_f</span><span class="ruby-operator">/</span><span class="ruby-constant">LOG10</span>).<span class="ruby-identifier">round</span>
|
|
42
|
+
<span class="ruby-keyword kw">end</span>
|
|
43
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y_f</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator"><</span> <span class="ruby-value">1</span>)
|
|
44
|
+
<span class="ruby-identifier">logy_f</span> = <span class="ruby-keyword kw">nil</span>
|
|
45
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">y_f</span>.<span class="ruby-identifier">kind_of?</span> <span class="ruby-constant">Float</span>) <span class="ruby-keyword kw">then</span>
|
|
46
|
+
<span class="ruby-identifier">logy_f</span> = <span class="ruby-constant">Math</span>.<span class="ruby-identifier">log</span>(<span class="ruby-identifier">y_f</span>.<span class="ruby-identifier">abs</span>)
|
|
47
|
+
<span class="ruby-keyword kw">else</span>
|
|
48
|
+
<span class="ruby-identifier">logy_f</span> = <span class="ruby-constant">LongMath</span>.<span class="ruby-identifier">log</span>(<span class="ruby-identifier">y_f</span>.<span class="ruby-identifier">abs</span>, <span class="ruby-value">15</span>, <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">ROUND_UP</span>)
|
|
49
|
+
<span class="ruby-keyword kw">if</span> (<span class="ruby-identifier">logy_f</span>.<span class="ruby-identifier">abs</span> <span class="ruby-operator">></span> <span class="ruby-constant">LongMath</span><span class="ruby-operator">::</span><span class="ruby-constant">MAX_FLOATABLE</span>) <span class="ruby-keyword kw">then</span>
|
|
50
|
+
<span class="ruby-identifier">raise</span> <span class="ruby-constant">ArgumentError</span>, <span class="ruby-node">"exponent would be way too big: y=#{y} logy_f=#{logy_f}"</span>;
|
|
51
|
+
<span class="ruby-keyword kw">end</span>
|
|
52
|
+
<span class="ruby-identifier">logy_f</span> = <span class="ruby-identifier">logy_f</span>.<span class="ruby-identifier">to_f</span>
|
|
53
|
+
<span class="ruby-keyword kw">end</span>
|
|
54
|
+
<span class="ruby-comment cmt"># puts("x=#{x} y=#{y} x_f=#{x_f} y_f=#{y_f} logx_f=#{logx_f} logy_f=#{logy_f} logx_y_f=#{logx_y_f}\n")</span>
|
|
55
|
+
<span class="ruby-identifier">iprec_y</span> <span class="ruby-operator">-=</span> (<span class="ruby-identifier">logy_f</span><span class="ruby-operator">/</span><span class="ruby-constant">LOG10</span>).<span class="ruby-identifier">round</span>
|
|
56
|
+
<span class="ruby-keyword kw">end</span>
|
|
57
|
+
<span class="ruby-comment cmt"># puts("x=#{x} y=#{y} x_f=#{x_f} y_f=#{y_f} logx_f=#{logx_f} logy_f=#{logy_f} logx_y_f=#{logx_y_f}\n")</span>
|
|
58
|
+
<span class="ruby-comment cmt"># puts("\niprec: x=#{x} y=#{y} iprec=#{iprec} iprec_x=#{iprec_x} iprec_y=#{iprec_y}\n")</span>
|
|
59
|
+
[ <span class="ruby-identifier">iprec</span>, <span class="ruby-identifier">iprec_x</span>, <span class="ruby-identifier">iprec_y</span>, <span class="ruby-identifier">logx_y_f</span> ]
|
|
60
|
+
|
|
61
|
+
<span class="ruby-keyword kw">end</span></pre>
|
|
62
|
+
</body>
|
|
63
|
+
</html>
|