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 // Test that the default PRNG does converge towards Pi when doing a Monte Carlo | 5 // Test that the default PRNG does converge towards Pi when doing a Monte Carlo |
6 // simulation. | 6 // simulation. |
7 | 7 |
8 // Library tag to allow Dartium to run the test. | 8 // Library tag to allow Dartium to run the test. |
9 library pi_test; | 9 library pi_test; |
10 | 10 |
11 import "package:expect/expect.dart"; | |
12 import 'dart:math'; | 11 import 'dart:math'; |
13 | 12 |
14 void main() { | 13 void main() { |
15 var seed = new Random().nextInt(1<<16); | 14 var seed = new Random().nextInt(1<<16); |
16 print("pi_test seed: $seed"); | 15 print("pi_test seed: $seed"); |
17 var prng = new Random(seed); | 16 var prng = new Random(seed); |
18 var outside = 0; | 17 var outside = 0; |
19 var inside = 0; | 18 var inside = 0; |
20 for (var i = 0; i < 600000; i++) { | 19 for (var i = 0; i < 600000; i++) { |
21 var x = prng.nextDouble(); | 20 var x = prng.nextDouble(); |
22 var y = prng.nextDouble(); | 21 var y = prng.nextDouble(); |
23 if ((x*x) + (y*y) < 1.0) { | 22 if ((x*x) + (y*y) < 1.0) { |
24 inside++; | 23 inside++; |
25 } else { | 24 } else { |
26 outside++; | 25 outside++; |
27 } | 26 } |
28 } | 27 } |
29 // Mmmmh, Pie! | 28 // Mmmmh, Pie! |
30 var pie = 4.0 * (inside/(inside + outside)); | 29 var pie = 4.0 * (inside/(inside + outside)); |
31 print("$pie"); | 30 print("$pie"); |
32 Expect.isTrue(((PI - 0.009) < pie) && (pie < (PI + 0.009))); | 31 Expect.isTrue(((PI - 0.009) < pie) && (pie < (PI + 0.009))); |
33 } | 32 } |
OLD | NEW |