relspeed_QBigInt.jl

using FastRationals, LinearAlgebra, BenchmarkTools

BenchmarkTools.DEFAULT_PARAMETERS.evals = 1;
BenchmarkTools.DEFAULT_PARAMETERS.samples = 100;
BenchmarkTools.DEFAULT_PARAMETERS.time_tolerance = 1.0e-8;
BenchmarkTools.DEFAULT_PARAMETERS.overhead = BenchmarkTools.estimate_overhead();

function Base.rand(::Type{Rational{I}}, n::Int=1) where {I<:Base.BitInteger}
    n = max(1, n)
    numers = rand(I, n)
    denoms = map(x->max(1,x), abs.(rand(I, n)))
    return numers .// denoms
end;

function Base.rand(::Type{Rational{BigInt}}, n::Int=1)
    n = max(1, n)
    qs = rand(Rational{Int128}, n)
    return Rational{BigInt}.(qs)
end;

q200 = rand(Rational{BigInt}, 200);
fq200 = FastRational{BigInt}.(q200);

q500 = rand(Rational{BigInt}, 500);
fq500 = FastRational{BigInt}.(q500);

q200sum = @belapsed sum($q200);
fq200sum = @belapsed sum($fq200);
q200prod = @belapsed prod($q200);
fq200prod = @belapsed prod($fq200);

q500sum = @belapsed sum($q500);
fq500sum = @belapsed sum($fq500);
q500prod = @belapsed prod($q500);
fq500prod = @belapsed prod($fq500);

relspeed_200sum =
  round( (q200sum / fq200sum), digits=1);
relspeed_200prod =
  round( (q200prod / fq200prod), digits=1);

relspeed_500sum =
  round( (q500sum / fq500sum), digits=1);
relspeed_500prod =
  round( (q500prod / fq500prod), digits=1);


q8x8 = reshape(q200[1:64], (8,8));
fq8x8 = reshape(fq200[1:64], (8,8));
q15x15 = reshape(q500[1:225], (15,15));
fq15x15 = reshape(fq500[1:225], (15,15));

q8x8mul = @belapsed $q8x8 * $q8x8;
fq8x8mul = @belapsed $fq8x8 * $fq8x8;
q15x15mul = @belapsed $q15x15 * $q15x15;
fq15x15mul = @belapsed $fq15x15 * $fq15x15;

q8x8tr = @belapsed tr($q8x8);
fq8x8tr = @belapsed tr($fq8x8);
q15x15tr = @belapsed tr($q15x15);
fq15x15tr = @belapsed tr($fq15x15);

relspeed_8x8mul =
  round( (q8x8mul / fq8x8mul), digits=1);
relspeed_8x8tr =
  round( (q8x8tr / fq8x8tr), digits=1);

relspeed_15x15mul =
  round( (q15x15mul / fq15x15mul), digits=1);
relspeed_15x15tr =
  round( (q15x15tr / fq15x15tr), digits=1);


relspeeds = string(
"\n\n\trelative speeds",
"\n\t (200)\t (500)\n\n",
"sum:   \t $relspeed_200sum \t $relspeed_500sum \n",
"prod:  \t $relspeed_200prod \t $relspeed_500prod \n");

print(relspeeds);

relspeeds = string(
"\n\n\trelative speeds",
"\n\t (8x8)\t (15x15)\n\n",
"matmul:\t $relspeed_8x8mul \t $relspeed_15x15mul\n",
"mat tr:\t $relspeed_8x8tr \t $relspeed_15x15tr \n");

print(relspeeds);
;