| // Copyright (c) 2012, the Dart project authors. Please see the AUTHORS file |
| // for details. All rights reserved. Use of this source code is governed by a |
| // BSD-style license that can be found in the LICENSE file. |
| |
| // Test that the default PRNG does converge towards Pi when doing a Monte Carlo |
| // simulation. |
| |
| import "package:expect/expect.dart"; |
| import 'dart:math'; |
| |
| var known_bad_seeds = const [50051, 55597, 59208]; |
| |
| void main([args]) { |
| // Select a seed either from the argument passed in or |
| // otherwise a random seed. |
| var seed = -1; |
| if ((args != null) && (args.length > 0)) { |
| seed = int.parse(args[0]); |
| } else { |
| var seed_prng = new Random(); |
| while (seed == -1) { |
| seed = seed_prng.nextInt(1 << 16); |
| if (known_bad_seeds.contains(seed)) { |
| // Reset seed and try again. |
| seed = -1; |
| } |
| } |
| } |
| |
| // Setup the PRNG for the Monte Carlo simulation. |
| print("pi_test seed: $seed"); |
| var prng = new Random(seed); |
| |
| var outside = 0; |
| var inside = 0; |
| for (var i = 0; i < 600000; i++) { |
| var x = prng.nextDouble(); |
| var y = prng.nextDouble(); |
| if ((x * x) + (y * y) < 1.0) { |
| inside++; |
| } else { |
| outside++; |
| } |
| } |
| // Mmmmh, Pie! |
| var pie = 4.0 * (inside / (inside + outside)); |
| print("$pie"); |
| Expect.isTrue(((pi - 0.009) < pie) && (pie < (pi + 0.009))); |
| } |