vec(x) = [component(x,1), component(x,2), component(x,3)];
vec( qfbred(Qfb(6,6,-1),1) )
default(realprecision,38)
q=Qfb(7, 30, -14)^2;
qfbpowraw(q,-1)
q*q^2
q^0
q^1
q^-1
q^-(2^64+1)
q=Qfb(2, 1, 3); q2=q*q;
q3=qfbcompraw(q,q2)
qfbpowraw(q,3)
qfbred(q3,1)
q=Qfb(1009, 60, 99108027750247771)
qfbnupow(q, 8839368315)
qfbred(Qfb(22000957029, 25035917443, 7122385192))
p=2^64+13;
qfbprimeform(-4,p)
qfbprimeform(5,p)
