OLD | NEW |
1 // Copyright (c) 2012, the Dart project authors. Please see the AUTHORS file | 1 // Copyright (c) 2012, the Dart project authors. Please see the AUTHORS file |
2 // for details. All rights reserved. Use of this source code is governed by a | 2 // for details. All rights reserved. Use of this source code is governed by a |
3 // BSD-style license that can be found in the LICENSE file. | 3 // BSD-style license that can be found in the LICENSE file. |
4 | 4 |
5 import "dart:typed_data"; | 5 import "dart:typed_data"; |
6 | 6 |
7 // A VM patch of the dart:math library. | 7 // A VM patch of the dart:math library. |
8 | 8 |
9 // If [x] is an [int] and [exponent] is a non-negative [int], the result is | 9 // If [x] is an [int] and [exponent] is a non-negative [int], the result is |
10 // an [int], otherwise the result is a [double]. | 10 // an [int], otherwise the result is a [double]. |
11 @patch num pow(num x, num exponent) { | 11 @patch |
| 12 num pow(num x, num exponent) { |
12 if ((x is int) && (exponent is int) && (exponent >= 0)) { | 13 if ((x is int) && (exponent is int) && (exponent >= 0)) { |
13 return _intPow(x, exponent); | 14 return _intPow(x, exponent); |
14 } | 15 } |
15 return _doublePow(x.toDouble(), exponent.toDouble()); | 16 return _doublePow(x.toDouble(), exponent.toDouble()); |
16 } | 17 } |
17 | 18 |
18 double _doublePow(double base, double exponent) { | 19 double _doublePow(double base, double exponent) { |
19 if (exponent == 0.0) { | 20 if (exponent == 0.0) { |
20 return 1.0; // ECMA-262 15.8.2.13 | 21 return 1.0; // ECMA-262 15.8.2.13 |
21 } | 22 } |
22 // Speed up simple cases. | 23 // Speed up simple cases. |
23 if (exponent == 1.0) return base; | 24 if (exponent == 1.0) return base; |
24 if (exponent == 2.0) return base * base; | 25 if (exponent == 2.0) return base * base; |
25 if (exponent == 3.0) return base * base * base; | 26 if (exponent == 3.0) return base * base * base; |
26 | 27 |
27 if (base == 1.0) return 1.0; | 28 if (base == 1.0) return 1.0; |
28 | 29 |
29 if (base.isNaN || exponent.isNaN) { | 30 if (base.isNaN || exponent.isNaN) { |
30 return double.NAN; | 31 return double.NAN; |
(...skipping 18 matching lines...) Expand all Loading... |
49 } | 50 } |
50 exponent >>= 1; | 51 exponent >>= 1; |
51 // Skip unnecessary operation (can overflow to Mint or Bigint). | 52 // Skip unnecessary operation (can overflow to Mint or Bigint). |
52 if (exponent != 0) { | 53 if (exponent != 0) { |
53 base *= base; | 54 base *= base; |
54 } | 55 } |
55 } | 56 } |
56 return result; | 57 return result; |
57 } | 58 } |
58 | 59 |
59 @patch double atan2(num a, num b) => _atan2(a.toDouble(), b.toDouble()); | 60 @patch |
60 @patch double sin(num x) => _sin(x.toDouble()); | 61 double atan2(num a, num b) => _atan2(a.toDouble(), b.toDouble()); |
61 @patch double cos(num x) => _cos(x.toDouble()); | 62 @patch |
62 @patch double tan(num x) => _tan(x.toDouble()); | 63 double sin(num x) => _sin(x.toDouble()); |
63 @patch double acos(num x) => _acos(x.toDouble()); | 64 @patch |
64 @patch double asin(num x) => _asin(x.toDouble()); | 65 double cos(num x) => _cos(x.toDouble()); |
65 @patch double atan(num x) => _atan(x.toDouble()); | 66 @patch |
66 @patch double sqrt(num x) => _sqrt(x.toDouble()); | 67 double tan(num x) => _tan(x.toDouble()); |
67 @patch double exp(num x) => _exp(x.toDouble()); | 68 @patch |
68 @patch double log(num x) => _log(x.toDouble()); | 69 double acos(num x) => _acos(x.toDouble()); |
| 70 @patch |
| 71 double asin(num x) => _asin(x.toDouble()); |
| 72 @patch |
| 73 double atan(num x) => _atan(x.toDouble()); |
| 74 @patch |
| 75 double sqrt(num x) => _sqrt(x.toDouble()); |
| 76 @patch |
| 77 double exp(num x) => _exp(x.toDouble()); |
| 78 @patch |
| 79 double log(num x) => _log(x.toDouble()); |
69 | 80 |
70 double _atan2(double a, double b) native "Math_atan2"; | 81 double _atan2(double a, double b) native "Math_atan2"; |
71 double _sin(double x) native "Math_sin"; | 82 double _sin(double x) native "Math_sin"; |
72 double _cos(double x) native "Math_cos"; | 83 double _cos(double x) native "Math_cos"; |
73 double _tan(double x) native "Math_tan"; | 84 double _tan(double x) native "Math_tan"; |
74 double _acos(double x) native "Math_acos"; | 85 double _acos(double x) native "Math_acos"; |
75 double _asin(double x) native "Math_asin"; | 86 double _asin(double x) native "Math_asin"; |
76 double _atan(double x) native "Math_atan"; | 87 double _atan(double x) native "Math_atan"; |
77 double _sqrt(double x) native "Math_sqrt"; | 88 double _sqrt(double x) native "Math_sqrt"; |
78 double _exp(double x) native "Math_exp"; | 89 double _exp(double x) native "Math_exp"; |
79 double _log(double x) native "Math_log"; | 90 double _log(double x) native "Math_log"; |
80 | 91 |
81 | |
82 // TODO(iposva): Handle patch methods within a patch class correctly. | 92 // TODO(iposva): Handle patch methods within a patch class correctly. |
83 @patch class Random { | 93 @patch |
84 | 94 class Random { |
85 @patch factory Random([int seed]) { | 95 @patch |
| 96 factory Random([int seed]) { |
86 var state = _Random._setupSeed((seed == null) ? _Random._nextSeed() : seed); | 97 var state = _Random._setupSeed((seed == null) ? _Random._nextSeed() : seed); |
87 // Crank a couple of times to distribute the seed bits a bit further. | 98 // Crank a couple of times to distribute the seed bits a bit further. |
88 return new _Random._withState(state).._nextState() | 99 return new _Random._withState(state) |
89 .._nextState() | 100 .._nextState() |
90 .._nextState() | 101 .._nextState() |
91 .._nextState(); | 102 .._nextState() |
| 103 .._nextState(); |
92 } | 104 } |
93 | 105 |
94 @patch factory Random.secure() { | 106 @patch |
| 107 factory Random.secure() { |
95 return new _SecureRandom(); | 108 return new _SecureRandom(); |
96 } | 109 } |
97 } | 110 } |
98 | 111 |
99 | |
100 class _Random implements Random { | 112 class _Random implements Random { |
101 // Internal state of the random number generator. | 113 // Internal state of the random number generator. |
102 final Uint32List _state; | 114 final Uint32List _state; |
103 static const _kSTATE_LO = 0; | 115 static const _kSTATE_LO = 0; |
104 static const _kSTATE_HI = 1; // Unused in Dart code. | 116 static const _kSTATE_HI = 1; // Unused in Dart code. |
105 | 117 |
106 _Random._withState(this._state); | 118 _Random._withState(this._state); |
107 | 119 |
108 // The algorithm used here is Multiply with Carry (MWC) with a Base b = 2^32. | 120 // The algorithm used here is Multiply with Carry (MWC) with a Base b = 2^32. |
109 // http://en.wikipedia.org/wiki/Multiply-with-carry | 121 // http://en.wikipedia.org/wiki/Multiply-with-carry |
110 // The constant A is selected from "Numerical Recipes 3rd Edition" p.348 B1. | 122 // The constant A is selected from "Numerical Recipes 3rd Edition" p.348 B1. |
111 | 123 |
112 // Implements: | 124 // Implements: |
113 // var state = | 125 // var state = |
114 // ((_A * (_state[_kSTATE_LO])) + _state[_kSTATE_HI]) & ((1 << 64) - 1); | 126 // ((_A * (_state[_kSTATE_LO])) + _state[_kSTATE_HI]) & ((1 << 64) - 1); |
115 // _state[_kSTATE_LO] = state & ((1 << 32) - 1); | 127 // _state[_kSTATE_LO] = state & ((1 << 32) - 1); |
116 // _state[_kSTATE_HI] = state >> 32; | 128 // _state[_kSTATE_HI] = state >> 32; |
117 // This is a native to prevent 64-bit operations in Dart, which | 129 // This is a native to prevent 64-bit operations in Dart, which |
118 // fail with --throw_on_javascript_int_overflow. | 130 // fail with --throw_on_javascript_int_overflow. |
119 // TODO(regis): Implement in Dart and remove Random_nextState in math.cc. | 131 // TODO(regis): Implement in Dart and remove Random_nextState in math.cc. |
120 void _nextState() native "Random_nextState"; | 132 void _nextState() native "Random_nextState"; |
121 | 133 |
122 int nextInt(int max) { | 134 int nextInt(int max) { |
123 const limit = 0x3FFFFFFF; | 135 const limit = 0x3FFFFFFF; |
124 if ((max <= 0) || ((max > limit) && (max > _POW2_32))) { | 136 if ((max <= 0) || ((max > limit) && (max > _POW2_32))) { |
125 throw new RangeError.range(max, 1, _POW2_32, "max", | 137 throw new RangeError.range( |
126 "Must be positive and <= 2^32"); | 138 max, 1, _POW2_32, "max", "Must be positive and <= 2^32"); |
127 } | 139 } |
128 if ((max & -max) == max) { | 140 if ((max & -max) == max) { |
129 // Fast case for powers of two. | 141 // Fast case for powers of two. |
130 _nextState(); | 142 _nextState(); |
131 return _state[_kSTATE_LO] & (max - 1); | 143 return _state[_kSTATE_LO] & (max - 1); |
132 } | 144 } |
133 | 145 |
134 var rnd32; | 146 var rnd32; |
135 var result; | 147 var result; |
136 do { | 148 do { |
(...skipping 29 matching lines...) Expand all Loading... |
166 // Get a seed from the VM's random number provider. | 178 // Get a seed from the VM's random number provider. |
167 static Uint32List _initialSeed() native "Random_initialSeed"; | 179 static Uint32List _initialSeed() native "Random_initialSeed"; |
168 | 180 |
169 static int _nextSeed() { | 181 static int _nextSeed() { |
170 // Trigger the PRNG once to change the internal state. | 182 // Trigger the PRNG once to change the internal state. |
171 _prng._nextState(); | 183 _prng._nextState(); |
172 return _prng._state[_kSTATE_LO]; | 184 return _prng._state[_kSTATE_LO]; |
173 } | 185 } |
174 } | 186 } |
175 | 187 |
176 | |
177 class _SecureRandom implements Random { | 188 class _SecureRandom implements Random { |
178 _SecureRandom() { | 189 _SecureRandom() { |
179 // Throw early in constructor if entropy source is not hooked up. | 190 // Throw early in constructor if entropy source is not hooked up. |
180 _getBytes(1); | 191 _getBytes(1); |
181 } | 192 } |
182 | 193 |
183 // Return count bytes of entropy as a positive integer; count <= 8. | 194 // Return count bytes of entropy as a positive integer; count <= 8. |
184 static int _getBytes(int count) native "SecureRandom_getBytes"; | 195 static int _getBytes(int count) native "SecureRandom_getBytes"; |
185 | 196 |
186 int nextInt(int max) { | 197 int nextInt(int max) { |
187 RangeError.checkValueInInterval( | 198 RangeError.checkValueInInterval( |
188 max, 1, _POW2_32, "max", "Must be positive and <= 2^32"); | 199 max, 1, _POW2_32, "max", "Must be positive and <= 2^32"); |
189 final byteCount = ((max - 1).bitLength + 7) >> 3; | 200 final byteCount = ((max - 1).bitLength + 7) >> 3; |
190 if (byteCount == 0) { | 201 if (byteCount == 0) { |
191 return 0; // Not random if max == 1. | 202 return 0; // Not random if max == 1. |
192 } | 203 } |
193 var rnd; | 204 var rnd; |
194 var result; | 205 var result; |
195 do { | 206 do { |
196 rnd = _getBytes(byteCount); | 207 rnd = _getBytes(byteCount); |
197 result = rnd % max; | 208 result = rnd % max; |
198 } while ((rnd - result + max) > (1 << (byteCount << 3))); | 209 } while ((rnd - result + max) > (1 << (byteCount << 3))); |
199 return result; | 210 return result; |
200 } | 211 } |
201 | 212 |
202 double nextDouble() { | 213 double nextDouble() { |
203 return (_getBytes(7) >> 3) / _POW2_53_D; | 214 return (_getBytes(7) >> 3) / _POW2_53_D; |
204 } | 215 } |
205 | 216 |
206 bool nextBool() { | 217 bool nextBool() { |
207 return _getBytes(1).isEven; | 218 return _getBytes(1).isEven; |
208 } | 219 } |
209 | 220 |
210 // Constants used by the algorithm. | 221 // Constants used by the algorithm. |
211 static const _POW2_32 = 1 << 32; | 222 static const _POW2_32 = 1 << 32; |
212 static const _POW2_53_D = 1.0 * (1 << 53); | 223 static const _POW2_53_D = 1.0 * (1 << 53); |
213 } | 224 } |
214 | |
OLD | NEW |