IEEE math for the masses

[p5sagit/p5-mst-13.2.git] / t / op / exp.t

#!./perl

BEGIN {
    chdir 't' if -d 't';
    @INC = '../lib';
    require './test.pl';
}

plan tests => 22;

# compile time evaluation

$s = sqrt(2);
is(substr($s,0,5), '1.414');

$s = exp(1);
is(substr($s,0,7), '2.71828');

ok(exp(log(1)) == 1);

# run time evaluation

$x1 = 1;
$x2 = 2;
$s = sqrt($x2);
is(substr($s,0,5), '1.414');

$s = exp($x1);
is(substr($s,0,7), '2.71828');

ok(exp(log($x1)) == 1);

# tests for transcendental functions

my $pi = 3.1415926535897931160;
my $pi_2 = 1.5707963267948965580;

sub round {
   my $result = shift;
   return sprintf("%.9f", $result);
}

# sin() tests
ok(sin(0) == 0.0);
ok(round(sin($pi)) == 0.0);
ok(round(sin(-1 * $pi)) == 0.0);
ok(round(sin($pi_2)) == 1.0);
ok(round(sin(-1 * $pi_2)) == -1.0);

# cos() tests
ok(cos(0) == 1.0);
ok(round(cos($pi)) == -1.0);
ok(round(cos(-1 * $pi)) == -1.0);
ok(round(cos($pi_2)) == 0.0);
ok(round(cos(-1 * $pi_2)) == 0.0);

# atan2() tests
ok(round(atan2(-0.0, 0.0)) == 0);
ok(round(atan2(0.0, 0.0)) == 0);
ok(round(atan2(-0.0, -0.0)) == round(-1 * $pi));
ok(round(atan2(0.0, -0.0)) == round($pi));
ok(round(atan2(-1.0, 0.0)) == round(-1 * $pi_2));
ok(round(atan2(1.0, 0.0)) == round($pi_2));