Post a reply
Username:
Note:If not registered, provide any username. For more comfort, register here.
Subject:
Message body:
Enter your message here, it may contain no more than 60000 characters. 

Smilies
:D :) :( :o :shock: :? 8) :lol: :x :P :oops: :cry: :evil: :twisted: :roll: :wink: :!: :?: :idea: :arrow: :| :mrgreen:
Font size:
Font colour
#000000 #000040 #000080 #0000BF #0000FF
#004000 #004040 #004080 #0040BF #0040FF
#008000 #008040 #008080 #0080BF #0080FF
#00BF00 #00BF40 #00BF80 #00BFBF #00BFFF
#00FF00 #00FF40 #00FF80 #00FFBF #00FFFF
#400000 #400040 #400080 #4000BF #4000FF
#404000 #404040 #404080 #4040BF #4040FF
#408000 #408040 #408080 #4080BF #4080FF
#40BF00 #40BF40 #40BF80 #40BFBF #40BFFF
#40FF00 #40FF40 #40FF80 #40FFBF #40FFFF
#800000 #800040 #800080 #8000BF #8000FF
#804000 #804040 #804080 #8040BF #8040FF
#808000 #808040 #808080 #8080BF #8080FF
#80BF00 #80BF40 #80BF80 #80BFBF #80BFFF
#80FF00 #80FF40 #80FF80 #80FFBF #80FFFF
#BF0000 #BF0040 #BF0080 #BF00BF #BF00FF
#BF4000 #BF4040 #BF4080 #BF40BF #BF40FF
#BF8000 #BF8040 #BF8080 #BF80BF #BF80FF
#BFBF00 #BFBF40 #BFBF80 #BFBFBF #BFBFFF
#BFFF00 #BFFF40 #BFFF80 #BFFFBF #BFFFFF
#FF0000 #FF0040 #FF0080 #FF00BF #FF00FF
#FF4000 #FF4040 #FF4080 #FF40BF #FF40FF
#FF8000 #FF8040 #FF8080 #FF80BF #FF80FF
#FFBF00 #FFBF40 #FFBF80 #FFBFBF #FFBFFF
#FFFF00 #FFFF40 #FFFF80 #FFFFBF #FFFFFF
Options:
BBCode is ON
[img] is ON
[flash] is OFF
[url] is ON
Smilies are ON
Disable BBCode
Disable smilies
Do not automatically parse URLs
Confirmation of post
To prevent automated posts the board requires you to enter a confirmation code. The code is displayed in the image you should see below. If you are visually impaired or cannot otherwise read this code please contact the %sBoard Administrator%s.
Confirmation code:
Enter the code exactly as it appears. All letters are case insensitive, there is no zero.
   

Topic review - Dealing with big matrices
Author Message
  Post subject:  Re: Dealing with big matrices  Reply with quote
Thank you. Unfortunately it seems that this was just covering up another bottleneck, which I (and probably you as well) already feared :( :
Code:
> module M = Ext(5,syz(D1),syz(CminD2));

Singular error: no more memory
System 11919379k:17115977k Appl 11905349k/14029k Malloc 11880467k/0k Valloc 38912k/14029k Pages 9453/275 Regions 19:103 Internal: 9808

halt 14

I will try to find a smaller example.
Post Posted: Thu Nov 03, 2016 8:05 pm
  Post subject:  Re: Dealing with big matrices  Reply with quote
See http://www.singular.uni-kl.de:8002/trac ... 40e01a/git
Post Posted: Mon Oct 31, 2016 5:44 pm
  Post subject:  Dealing with big matrices  Reply with quote
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?
Post Posted: Thu Oct 27, 2016 4:24 pm


It is currently Fri May 13, 2022 10:58 am
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group