Singularhttp://www.singular.uni-kl.de/forum/ how to compute all combination of system of polynomialshttp://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 seriesmaple code got out of memoryAllMatrices := 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 domono1 := 1;for i from 1 to nops(sol[1][j]) domono1 := 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]) dosol := 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 Grouphttp://www.phpbb.com/