File Coverage

lib/Math/Random/NormalDistribution.pm

Criterion	Covered	Total	%
statement	23	23	100.0
branch	2	4	50.0
condition			n/a
subroutine	8	8	100.0
pod	1	1	100.0
total	34	36	94.4

package Math::Random::NormalDistribution;
# coding: UTF-8

use utf8;
use strict;
use warnings;

# ------------------------------------------------------------------------------
our $VERSION = '0.01';

use Exporter;
use base qw(Exporter);
our @EXPORT = qw(
    rand_nd_generator
);
# ------------------------------------------------------------------------------
use constant TWOPI => 2.0 * 4.0 * atan2(1.0, 1.0);

sub rand_nd_generator(;@)
{
    my ($mean, $stddev) = @_;
    $mean = 0.0 if ! defined $mean;
    $stddev = 1.0 if ! defined $stddev;

    return sub {
        return $mean + $stddev * cos(TWOPI * (1.0 - rand)) * sqrt(-2.0 * log(1.0 - rand));
    }
}
# ------------------------------------------------------------------------------
1;
__END__