ppm install http://www.sisyphusion.tk/ppm/Math-GMPz.ppd
ppm install http://www.sisyphusion.tk/ppm/Math-Int128.ppd

##</code><code>##

use warnings;
use Math::Int128 qw(int128);
use Math::GMPz qw(:mpz);
use Benchmark qw(:all);

$count = 40000;

$mpz1  = Math::GMPz->new('676469752303423489');
$mpz2  = Math::GMPz->new('776469752999423489');
$i_1   = int128("$mpz1");
$i_2   = int128("$mpz2");

$mpz_sub = Math::GMPz->new('976469752313423489');
$i_sub   = int128("$mpz_sub");

$mpz_div = Math::GMPz->new('76469752313423489');
$i_div   = int128("$mpz_div");

$mpz_ret = Rmpz_init2(128);

print "
******************
**MULTIPLICATION**
******************\n\n";

cmpthese($count * 9, {
    'mul_M::I' => '$ri = Math::Int128::_mul($i_1, $i_2, 0)',
    'mul_M::G1' =>'$mpz_ret = $mpz1 * $mpz2',
    'mul_M::G2'=> 'Rmpz_mul($mpz_ret, $mpz1, $mpz2)',
});


die "Error 1:\n$ri\n$mpz_ret\n" if $ri != int128("$mpz_ret")
 || $ri != int128('525258301482620425304858018020933121');


$i_1 *= $i_1;
$i_2 *= $i_2;
$mpz1 *= $mpz1;
$mpz2 *= $mpz2;


print "
******************
*****DIVISION*****
******************\n\n";

cmpthese($count * 9, {
    'div_M::I' => '$ri = Math::Int128::_div($i_1, $i_div, 0)',
    'div_M::G1' =>'$mpz_ret = $mpz1 / $mpz_div',
    'div_M::G2'=> 'Rmpz_tdiv_q($mpz_ret, $mpz1, $mpz_div)',
});


die "Error 2:\n$ri\n$mpz_ret\n" if $ri != int128("$mpz_ret")
 || $ri != int128('5984213521522366751');

print"
******************
*****ADDITION*****
******************\n\n";

cmpthese($count * 10, {
    'add_M::I' => '$ri = Math::Int128::_add($i_1, $i_2, 0)',
    'add_M::G1' => '$mpz_ret = $mpz1  + $mpz2',
    'add_M::G2' => 'Rmpz_add($mpz_ret, $mpz1, $mpz2)',
});


die "Error 3:\n$ri\n$mpz_ret\n" if $ri != int128("$mpz_ret") 
 || $ri != int128('1060516603104440851094132036041866242');

print "
******************
****SUBTRACTION***
******************\n\n";

cmpthese($count * 10, {
    'add_M::I' => '$ri = Math::Int128::_sub($i_1, $i_sub, 0)',
    'add_M::G1' => '$mpz_ret = $mpz1 - $mpz_sub',
    'add_M::G2' => 'Rmpz_sub($mpz_ret, $mpz1, $mpz_sub)',
});


die "Error 4:\n$ri\n$mpz_ret\n" if $ri != int128("$mpz_ret") 
 || $ri != int128('457611325781455127825205517363509632');

##</code><code>##

C:\_64\pscrpt>perl bench.pl
Name "main::i_div" used only once: possible typo at bench.pl line 17.
Name "main::i_sub" used only once: possible typo at bench.pl line 14.

******************
**MULTIPLICATION**
******************

              Rate mul_M::G1  mul_M::I mul_M::G2
mul_M::G1 156726/s        --      -58%      -82%
mul_M::I  371517/s      137%        --      -58%
mul_M::G2 886700/s      466%      139%        --

******************
*****DIVISION*****
******************

              Rate div_M::G1  div_M::I div_M::G2
div_M::G1 148637/s        --      -57%      -81%
div_M::I  344168/s      132%        --      -57%
div_M::G2 794702/s      435%      131%        --

******************
*****ADDITION*****
******************

              Rate add_M::G1  add_M::I add_M::G2
add_M::G1 148810/s        --      -60%      -83%
add_M::I  376648/s      153%        --      -56%
add_M::G2 852878/s      473%      126%        --

******************
****SUBTRACTION***
******************

              Rate add_M::G1  add_M::I add_M::G2
add_M::G1 147113/s        --      -60%      -83%
add_M::I  365631/s      149%        --      -59%
add_M::G2 883002/s      500%      142%        --

C:\_64\pscrpt>

##</code><code>##

            Rate    gmpz mul_128
gmpz    900000/s      --     -3%
mul_128 927835/s      3%      --