Dear Gorzel:
The following is a toy model.
LIB "nctools.lib"; LIB "central.lib"; LIB "ring.lib"; LIB "linalg.lib";
// Define the ring ring r0 = complex,(x1,x2),dp; ring r1 = (0,t), (y1,y2), dp; matrix C[2][2]; matrix D[2][2]; C[1,2]=1; def S=nc_algebra(C,D); setring S;S; ideal Basis = PBW_eqDeg(1);
// Define a matrix for a linear transformation. matrix A[2][2]; A[1,1]=1; A[1,2]=t^2; A[2,1]=t^2; A[2,2]=t^4;
int rankA = mat_rk(A); write (":w rankA.m", rankA);
// Now generate the ideal of images.
ring r2 = 0, (z1,z2), dp;
fetchall(S);
poly h(1)= A[1,1]*z1 + A[1,2]*z2; poly h(2)= A[2,1]*z1 + A[2,2]*z2;
map F = r0, h(1..2);
ideal Image = h(1..2); module M = linearMapKernel(Image); M; ideal K = linearCombinations(Basis, M); K; write (".w K.m", K);
exit;
Now, we want to do the same, but the matrix we are working with is 2100 X 1001, so it seems to require a lot of memory. The system administrator suggest us to either try to run Singular using the 8 Cores of the machine (i.e. a parallel process) and /or run Singular in silent mode. To achieve the last goal we wrote a file, called toy-model.m, with the same instructions as above, and try to run it as Singular toy-model.m & hopping for it to run at the background, without showing anything, just writing the answer in a file, but we fail...
Thanks in advance for your help.
PS. You can also see a previous note of october 2, 2009, posted also by Nieto, called vectorspace and relations, answered by Oleksandr
Dear Gorzel:
The following is a toy model.
LIB "nctools.lib"; LIB "central.lib"; LIB "ring.lib"; LIB "linalg.lib";
// Define the ring ring r0 = complex,(x1,x2),dp; ring r1 = (0,t), (y1,y2), dp; matrix C[2][2]; matrix D[2][2]; C[1,2]=1; def S=nc_algebra(C,D); setring S;S; ideal Basis = PBW_eqDeg(1);
// Define a matrix for a linear transformation. matrix A[2][2]; A[1,1]=1; A[1,2]=t^2; A[2,1]=t^2; A[2,2]=t^4;
int rankA = mat_rk(A); write (":w rankA.m", rankA);
// Now generate the ideal of images.
ring r2 = 0, (z1,z2), dp;
fetchall(S);
poly h(1)= A[1,1]*z1 + A[1,2]*z2; poly h(2)= A[2,1]*z1 + A[2,2]*z2;
map F = r0, h(1..2);
ideal Image = h(1..2); module M = linearMapKernel(Image); M; ideal K = linearCombinations(Basis, M); K; write (".w K.m", K);
exit;
Now, we want to do the same, but the matrix we are working with is 2100 X 1001, so it seems to require a lot of memory. The system administrator suggest us to either try to run Singular using the 8 Cores of the machine (i.e. a parallel process) and /or run Singular in silent mode. To achieve the last goal we wrote a file, called toy-model.m, with the same instructions as above, and try to run it as Singular toy-model.m & hopping for it to run at the background, without showing anything, just writing the answer in a file, but we fail...
Thanks in advance for your help.
PS. You can also see a previous note of october 2, 2009, posted also by Nieto, called vectorspace and relations, answered by Oleksandr
|