#!/usr/bin/perl -lw
use strict;

use Inline C => Config => BUILD_NOISY => 1;
use Inline C => <<'END_C',  NAME => 'ieeroundmode', CLEAN_AFTER_BUILD => 0;

unsigned int controlFP( unsigned int new, unsigned int mask ) {
    return _controlfp( new, mask );
}

END_C

use constant {
    RC_MASK => 0x00000300,
    RC_CHOP => 0x00000300,
    RC_UP   => 0x00000200,
    RC_DOWN => 0x00000100,
    RC_NEAR => 0x00000000,
};

my @n;
my $m;

@n = qw< 0.000005 0.000015 0.000025 0.000035 0.000045
         0.000055 0.000065 0.000075 0.000085 0.000095 >;
$m= $n[0];
for my $n ( @n ) {
    printf " %.22f   %.5f   %.22f   %.5f\n", $n, $n, $m, $m;
    $m += 0.00001;
}

print "Setting mode CHOP; controlFP returned: ", controlFP( RC_CHOP, RC_MASK );
@n = qw< 0.000005 0.000015 0.000025 0.000035 0.000045
         0.000055 0.000065 0.000075 0.000085 0.000095 >;
$m= $n[0];
for my $n ( @n ) {
    printf " %.22f   %.5f   %.22f   %.5f\n", $n, $n, $m, $m;
    $m += 0.00001;
}

print "Setting mode UP; controlFP returned: ", controlFP( RC_UP, RC_MASK );
@n = qw< 0.000005 0.000015 0.000025 0.000035 0.000045
         0.000055 0.000065 0.000075 0.000085 0.000095 >;
$m= $n[0];
for my $n ( @n ) {
    printf " %.22f   %.5f   %.22f   %.5f\n", $n, $n, $m, $m;
    $m += 0.00001;
}

print "Setting mode DOWN; controlFP returned: ", controlFP( RC_DOWN, RC_MASK );
@n = qw< 0.000005 0.000015 0.000025 0.000035 0.000045
         0.000055 0.000065 0.000075 0.000085 0.000095 >;
$m= $n[0];
for my $n ( @n ) {
    printf " %.22f   %.5f   %.22f   %.5f\n", $n, $n, $m, $m;
    $m += 0.00001;
}

print "Setting mode NEAR; controlFP returned: ", controlFP( RC_NEAR, RC_MASK );
@n = qw< 0.000005 0.000015 0.000025 0.000035 0.000045
         0.000055 0.000065 0.000075 0.000085 0.000095 >;
$m= $n[0];
for my $n ( @n ) {
    printf " %.22f   %.5f   %.22f   %.5f\n", $n, $n, $m, $m;
    $m += 0.00001;
}

##</code><code>##

Default / NEAR                       (compare computed value with constant)
 0.0000050000000000000004   0.00001  =0.0000050000000000000004= =0.00001=
 0.0000150000000000000000   0.00002  <0.0000150000000000000020> =0.00002=
 0.0000250000000000000010   0.00003  <0.0000250000000000000050> =0.00003=
 0.0000349999999999999970   0.00003  <0.0000350000000000000040> <0.00004>
 0.0000450000000000000030   0.00005  =0.0000450000000000000030= =0.00005=
 0.0000550000000000000020   0.00006  =0.0000550000000000000020= =0.00006=
 0.0000649999999999999940   0.00006  <0.0000650000000000000080> <0.00007>
 0.0000749999999999999930   0.00007  <0.0000750000000000000070> <0.00008>
 0.0000850000000000000060   0.00009  =0.0000850000000000000060= =0.00009=
 0.0000950000000000000050   0.00010  =0.0000950000000000000050= =0.00010=
CHOP / DOWN (compare value to Default / NEAR case)
>0.0000049999999999999996< >0.00000< >0.0000049999999999999996< >0.00000<
>0.0000149999999999999990< >0.00001< >0.0000150000000000000000< =0.00002=
>0.0000249999999999999980< >0.00002< >0.0000250000000000000010< =0.00003=
=0.0000349999999999999970= =0.00003= >0.0000349999999999999970< >0.00003<
>0.0000449999999999999960< >0.00004< >0.0000449999999999999960< >0.00004<
>0.0000549999999999999950< >0.00005< >0.0000549999999999999950< >0.00005<
>0.0000649999999999999940< >0.00006< >0.0000649999999999999940< >0.00006<
=0.0000749999999999999930= >0.00007< >0.0000749999999999999930< >0.00007<
>0.0000849999999999999930< >0.00008< >0.0000849999999999999930< >0.00008<
>0.0000949999999999999920< >0.00009< >0.0000949999999999999920< >0.00009<
UP (compare value to Default / NEAR case)
=0.0000050000000000000004= =0.00001= =0.0000050000000000000004= =0.00001=
=0.0000150000000000000000= =0.00002= =0.0000150000000000000020= =0.00002=
=0.0000250000000000000010= =0.00003= =0.0000250000000000000050> =0.00003=
<0.0000350000000000000040> <0.00004> <0.0000350000000000000100> =0.00004=
=0.0000450000000000000030= =0.00005= <0.0000450000000000000160> =0.00005=
=0.0000550000000000000020= =0.00006= <0.0000550000000000000220> =0.00006=
<0.0000650000000000000080> <0.00007> <0.0000650000000000000350> =0.00007=
<0.0000750000000000000070> <0.00008> <0.0000750000000000000480> =0.00008=
=0.0000850000000000000060= =0.00009= <0.0000850000000000000600> =0.00009=
=0.0000950000000000000050= =0.00010= <0.0000950000000000000730> =0.00010=