Closed jakobkroeker closed 8 years ago
def assumeLevel = 3; system("random",1752752706); LIB("primdecint.lib"); proc quotientZWrapper(ideal A,ideal B) { return(Primdecint::quotientZ(A,B)); } intvec op = 100663426,10321; ring rng = (integer),(x,y,z),(dp(3),C); short = 0; ideal I = 17*x*z+10*x-43,41*x*y-16*x,47*z^2-2*z; string(I); ideal J = 34*x-16,-2*y*z,-20; string(J); bigint trial = 4; ideal quotientZ_JI = quotientZWrapper(J,I); ideal quotient_JI = quotient(J,I); quotientZ_JI; quotient_JI; string(quotientZ_JI); // 20,4*y*z,2*x*z+12*z string(quotient_JI); // 20,2*x+12,2*y*z
looks better now
bugfix was merged