\\ Load the arithmetic operations

\r arith.gp

\\ Now some tests

for(i=0, 10, \
n=random(90);m=random(100); \
u=pinput(random(10^n)); \
v=pinput(random(10^m)); \
print1("size(u)=");print(size(u)); \
print1("size(v)=");print(size(v)); \
time=0;w=division(u,v);print1("divtime=");print(time); \
print1("division(");print1(poutput(u));print1(",");print1(poutput(v));print1(")=");print([poutput(w[1]),poutput(w[2])]); \
print1("error=");print([poutput(v)\poutput(u),poutput(v)%poutput(u)]-[poutput(w[1]),poutput(w[2])]); \
time=0;x=multiplication(u,v);print1("multime=");print(time); \
print1("multiplication(");print1(poutput(u));print1(",");print1(poutput(v));print1(")=");print(poutput(x)); \
print1("error=");print(poutput(u)*poutput(v)-poutput(x)); \
)