want to classify all combination of system of polynomials with hilbert series
maple code got out of memory
AllMatrices := proc (A::set, k::posint, n::posint) local B, C, E: B := [[]]: C := proc () B := [seq(seq([A[i], op(B[j])], i = 1 .. nops(A)), j = 1 .. nops(B))]: end proc: E := (C@@(k*n))(B): seq(Matrix(k, n, E[m]), m = 1 .. nops(A)^(k*n)); end proc: mm := AllMatrices({0, 1}, 3, 3);
GetRing := proc(sol) ringequation := 0; mono1 := 0; for j from 1 to 3 do mono1 := 1; for i from 1 to nops(sol[1][j]) do mono1 := mono1*op(i, sol[1][j]); od: ringequation := ringequation + mono1; od: return ringequation; end proc;
with(LinearAlgebra): polylistresult := []; for i from 1 to nops([mm]) do sol := MatrixMatrixMultiply(Matrix([[a,b,c]]), op(i,[mm])); sol := GetRing(sol); polylistresult := [op(polylistresult), sol]; od: gaga := polylistresult; nops(polylistresult); polylistresult := [op(polylistresult), op(gaga)]; nops(polylistresult); polylistresult := [op(polylistresult), op(gaga)]; nops(polylistresult); with(combinat): gg := permute(polylistresult, 3); Error, (in combinat:-permute) Maple was unable to allocate enough memory to complete this computation. Please see ?alloc
want to classify all combination of system of polynomials with hilbert series
maple code got out of memory
AllMatrices := proc (A::set, k::posint, n::posint) local B, C, E: B := [[]]: C := proc () B := [seq(seq([A[i], op(B[j])], i = 1 .. nops(A)), j = 1 .. nops(B))]: end proc: E := (C@@(k*n))(B): seq(Matrix(k, n, E[m]), m = 1 .. nops(A)^(k*n)); end proc: mm := AllMatrices({0, 1}, 3, 3);
GetRing := proc(sol) ringequation := 0; mono1 := 0; for j from 1 to 3 do mono1 := 1; for i from 1 to nops(sol[1][j]) do mono1 := mono1*op(i, sol[1][j]); od: ringequation := ringequation + mono1; od: return ringequation; end proc;
with(LinearAlgebra): polylistresult := []; for i from 1 to nops([mm]) do sol := MatrixMatrixMultiply(Matrix([[a,b,c]]), op(i,[mm])); sol := GetRing(sol); polylistresult := [op(polylistresult), sol]; od: gaga := polylistresult; nops(polylistresult); polylistresult := [op(polylistresult), op(gaga)]; nops(polylistresult); polylistresult := [op(polylistresult), op(gaga)]; nops(polylistresult); with(combinat): gg := permute(polylistresult, 3); Error, (in combinat:-permute) Maple was unable to allocate enough memory to complete this computation. Please see ?alloc
|