Singular https://www.singular.uni-kl.de/forum/ |
|
how to compute all combination of system of polynomials https://www.singular.uni-kl.de/forum/viewtopic.php?f=10&t=2457 |
Page 1 of 1 |
Author: | beginner [ Fri Sep 12, 2014 11:41 am ] |
Post subject: | how to compute all combination of system of polynomials |
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 |
Page 1 of 1 | All times are UTC + 1 hour [ DST ] |
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group http://www.phpbb.com/ |