I am trying to compute Ext for two ideals as follows:
Code: ring R = 0,(x(1),x(2),x(3),x(4),x(5),x(6),x(7)),dp;
ideal I = x(4)*x(6)-x(3)*x(7),
x(5)^2-x(3)*x(7),
x(4)*x(5)-x(2)*x(7),
x(3)*x(5)-x(2)*x(6),
x(1)*x(5)-x(6)*x(7),
x(3)*x(4)-x(2)*x(5),
x(1)*x(4)-x(5)*x(7),
x(1)*x(3)-x(5)*x(6),
x(1)*x(2)-x(3)*x(7),
x(2)*x(5)*x(6)-x(3)^2*x(7),
x(2)^2*x(6)^2-x(3)^3*x(7);
qring Q = std(I);
ideal D1 = x(7)^2, x(6)*x(7), x(5)*x(7), x(4)*x(7), x(1)*x(7), x(1)*x(6), x(1)^2, x(6)^3, x(5)*x(6)^2, x(4)^3;
ideal D2 = x(7)^2, x(4)*x(7), x(4)^2, x(2)*x(5)*x(7), x(2)^2*x(7), x(1)^2*x(7), x(2)^2*x(4), x(2)^3*x(5), x(2)^4, x(1)^4;
ideal CminD2 = x(4)^2, x(2)*x(7)^2, x(2)*x(5)*x(7), x(2)*x(4)*x(7), x(2)*x(3)*x(7), x(2)^2*x(7), x(3)*x(7)^3, x(3)^2*x(7)^2, x(3)^3*x(7), x(5)*x(6)*x(7)^3, x(2)*x(3)^2*x(6)^2, x(6)^2*x(7)^4, x(3)^4*x(6)^2;
ideal CminD1 = x(6)^3, x(3)*x(6)^2*x(7), x(3)^2*x(6)^2, x(3)^2*x(6)*x(7)^2, x(3)^3*x(6)*x(7), x(3)^4*x(6), x(3)^3*x(7)^3, x(3)^4*x(7)^2, x(3)^5*x(7), x(3)^6, x(2)^2*x(3)*x(6)*x(7)^3, x(2)^2*x(3)^2*x(7)^4, x(2)^4*x(6)*x(7)^4, x(2)^4*x(3)*x(7)^5, x(2)^5*x(5)*x(7)^5, x(2)^6*x(7)^6;
LIB "homolog.lib";
module M = Ext(5,syz(D1),syz(CminD2));
after quite some time this fails with the following error:
Code: > module M = Ext(5,syz(D1),syz(CminD2));
? internal error: creating matrix[7647][38235]
Singular : signal 11 (v: 4020):
current line:>> return(matrix(B));<<
Segment fault/Bus error occurred at 7f6dd04a25a0 because of 10246 (r:1477505230)
please inform the authors
trying to restart...
? error occurred in or before matrix.lib::concat line 106: ` return(matrix(B));`
? wrong type declaration. type 'help matrix;'
? leaving matrix.lib::concat
skipping text from `;` error at token `)`
? leaving matrix.lib::concat
Auf Wiedersehen.
I suspect this matrix that Singular is creating might be too big. Is this correct? If so, is there any possibility to rebuild Singular in a way such that it can handle bigger matrices?
|
: