|\^/| Maple 2022 (X86 64 LINUX) ._|\| |/|_. Copyright (c) Maplesoft, a division of Waterloo Maple Inc. 2022 \ MAPLE / All rights reserved. Maple is a trademark of <____ ____> Waterloo Maple Inc. | Type ? for help. #Author := Ayoola Jinadu # If you find a bug in the code please email me at ajinadu@sfu.ca > Sys := [ > aa * t2^2 * t1^2 * b2^2 * b1^2 + 2 * l1 * t1^2 * b2^2 * b1^2 + aa * > t1^2 * b2^2 * b1^2 - 2 * l1 * t2^2 * b2^2 * b1^2 + aa * t2^2 * b2^2 > * b1^2 + aa * b2^2 * b1^2 > + 2 * l2 * t2^2 * t1^2 * b1^2 + aa * t2^2 * t1^2 * b1^2 + 2 * l2 * > t1^2 * b1^2 + 2 * l1 * t1^2 * b1^2 + aa * t1^2 * b1^2 + 2 * l2 * > t2^2 * b1^2 - 2 * l1 * t2^2 * b1^2 > + aa * t2^2 * b1^2 + 2 * l2 * b1^2 + aa * b1^2 - 2 * l2 * t2^2 * > t1^2 * b2^2 + aa * t2^2 * t1^2 * b2^2 - 2 * l2 * t1^2 * b2^2 + 2 * > l1 * t1^2 * b2^2 + aa * t1^2 * b2^2 > - 2 * l2 * t2^2 * b2^2 - 2 * l1 * t2^2 * b2^2 + aa * t2^2 * b2^2 - > 2 * l2 * b2^2 + aa * b2^2 + aa * t2^2 * t1^2 + 2 * l1 * t1^2 + aa * > t1^2 - 2 * l1 * t2^2 + aa * t2^2 + aa, > 2 * l1 * t2 * t1^2 * b2^2 * b1^2 - 2 * l1 * t2^2 * t1 * b2^2 * b1^ > 2 - 2 * l1 * t1 * b2^2 * b1^2 + 2 * l1 * t2 * b2^2 * b1^2 + 2 * l2 > * t2^2 * t1^2 * b2 * b1^2 > + 2 * l2 * t1^2 * b2 * b1^2 + 2 * l2 * t2^2 * b2 * b1^2 + 2 * l2 > * b2 * b1^2 + 2 * l1 * t2 * t1^2 * b1^2 - 2 * l1 * t2^2 * t1 * b1^ > 2 - 2 * l1 * t1 * b1^2 + 2 * l1 * t2 * b1^2 > - 2 * l2 * t2^2 * t1^2 * b2^2 * b1 - 2 * l2 * t1^2 * b2^2 * b1 - 2 > * l2 * t2^2 * b2^2 * b1 - 2 * l2 * b2^2 * b1 - 2 * l2 * t2^2 * t1^ > 2 * b1 - 2 * l2 * t1^2 * b1 - 2 * l2 * t2^2 * > b1 - 2 * l2 * b1 + 2 * l1 * t2 * t1^2 * b2^2 - 2 * l1 * t2^2 * t1 > * b2^2 - 2 * l1 * t1 * b2^2 + 2 * l1 * t2 * b2^2 + 2 * l2 * t2^2 * > t1^2 * b2 + 2 * l2 * t1^2 * b2 > + 2 * l2 * t2^2 * b2 + 2 * l2 * b2 + 2 * l1 * t2 * t1^2 - 2 * l1 > * t2^2 * t1 - 2 * l1 * t1 + 2 * l1 * t2, > - al^2 * x * t1^2 * b1^2 - x * t1^2 * b1^2 - l3 * al^2 * t1^2 * b1^2 > - l2 * al^2 * t1^2 * b1^2 - l1 * al^2 * t1^2 * b1^2 + l3 * t1^2 * > b1^2 - l2 * t1^2 * b1^2 > - l1 * t1^2 * b1^2 - al^2 * x * b1^2 - x * b1^2 - l3 * al^2 * b1^2 > - l2 * al^2 * b1^2 + l1 * al^2 * b1^2 + l3 * b1^2 - l2 * b1^2 + l1 > * b1^2 - 4 * l3 * al * t1^2 * > b1 - 4 * l3 * al * b1 - al^2 * x * t1^2 - x * t1^2 + l3 * al^2 * t1^2 > + l2 * al^2 * t1^2 - l1 * al^2 * t1^2 - l3 * t1^2 + l2 * t1^2 - l1 > * t1^2 - al^2 * x - x > + l3 * al^2 + l2 * al^2 + l1 * al^2 - l3 + l2 + l1, > - al^2 * y * t1^2 * b1^2 - y * t1^2 * b1^2 - 2 * l3 * al * t1^2 * > b1^2 + 2 * l1 * al^2 * t1 * b1^2 + 2 * l1 * t1 * b1^2 - al^2 * y * > b1^2 - y * b1^2 - 2 * l3 * al * b1^2 + 2 * l3 * al^2 * t1^2 * b1 + 2 > * l2 * al^2 * t1^2 * b1 - 2 * l3 * t1^2 * b1 + 2 * l2 * t1^2 * b1 > + 2 * l3 * al^2 * b1 + 2 * l2 * al^2 * b1 - 2 * l3 * b1 + 2 * l2 * > b1 - al^2 * y * t1^2 - y * t1^2 + 2 * l3 * al * t1^2 + 2 * l1 * al^2 * > t1 + 2 * l1 * t1 - al^2 * y - y + 2 * l3 * al]: > elim := [t2,b1,b2]; elim := [t2, b1, b2] > read det; read newdeg; > read dixon;read minor; read dixres; read nextprime; read Bmbot; memory used=2.7MB, alloc=136.3MB, time=0.07 > with(LinearAlgebra): > M := dixonmatrix( Sys, elim ): Dixon: n=3 Dixon: matrix done memory used=4.8MB, alloc=137.9MB, time=0.10 Dixon: #delta=104256 Dixon: 32 x 32 > rank,rows,cols := minor(M): minor: 32 x 32 minor: #nonzero=580 memory used=75.6MB, alloc=144.7MB, time=0.36 minor: starting elimination > B := M[rows,cols]: # extract minor > L := StronglyConnectedBlocks(B): > BlockStructure = map(RowDimension,L); BlockStructure = [8, 8] > X := [t1, aa, al, l1, l2, l3, x, y]; X := [t1, aa, al, l1, l2, l3, x, y] > R := convert(L[1],listlist): > n := RowDimension(R); LL := R: n := 8 > deg_M := [seq( max(seq(seq(degree(LL[i][j],var),i=1..n),j=1..n)), var in X )]; deg_M := [8, 1, 4, 2, 4, 2, 1, 1] > rt := [R,X,n ]: > GlobalCArray := Array(0..n-1,0..n-1,order=C_order,datatype=integer[8]): > E := DixonRes(BB,X); Number of Monic Square Factors = 3 memory used=216.9MB, alloc=152.7MB, time=3.79 deg x1 =65,emax = 10,#probes for the first prime = 16641 2 2 2 2 2 2 2 2 2 2 E := (al l1 t1 + 2 al l1 t1 x - al l2 t1 - 2 al l2 l3 t1 2 2 2 2 2 2 2 2 2 2 2 2 - al l3 t1 + al t1 x + al t1 y - 4 al l1 t1 y + al l1 2 2 2 2 2 2 2 2 2 2 2 2 - 2 al l1 x - al l2 - 2 al l2 l3 - al l3 + al x + al y + l1 t1 2 2 2 2 2 2 2 2 2 2 + 2 l1 t1 x - l2 t1 + 2 l2 l3 t1 - l3 t1 + t1 x + t1 y 2 2 2 2 2 4 4 - 4 l1 t1 y + l1 - 2 l1 x - l2 + 2 l2 l3 - l3 + x + y ) (aa t1 2 2 4 4 2 3 2 2 2 2 2 2 2 - 4 aa l2 t1 + 2 aa t1 - 8 aa l1 t1 + 8 aa l1 t1 - 8 aa l2 t1 2 2 4 3 2 2 2 2 3 + 16 aa l1 l2 t1 + aa - 8 aa l1 + 24 aa l1 - 4 aa l2 - 32 aa l1 2 4 2 2 2 + 16 aa l1 l2 + 16 l1 - 16 l1 l2 ) (t1 + 1) > nops( expand(E) ); 450 > quit memory used=238.6MB, alloc=152.7MB, time=4.24