ASKSAGE: Sage Q&A Forum - RSS feedhttps://ask.sagemath.org/questions/Q&A Forum for SageenCopyright Sage, 2010. Some rights reserved under creative commons license.Wed, 30 Dec 2020 08:55:58 +0100Kernel died when computing large matrixhttps://ask.sagemath.org/question/55000/kernel-died-when-computing-large-matrix/ I have the following
m=5; n=3
A = matrix(QQ, (m^n)^2+1, (m^m)^2)
A[0, :] = ones_matrix(1, (m^m)^2)
for i in [1..(m^n)^2]: A[i,(m^(2*(m-n)))*(i-1):(m^(2*(m-n)))*i] = ones_matrix(1, m^(2*(m-n)))
E = A.right_kernel()
And I get the message: the kernel died...
Any suggestion to big matrices?
Wed, 30 Dec 2020 02:18:22 +0100https://ask.sagemath.org/question/55000/kernel-died-when-computing-large-matrix/Answer by Emmanuel Charpentier for <p>I have the following</p>
<pre><code>m=5; n=3
A = matrix(QQ, (m^n)^2+1, (m^m)^2)
A[0, :] = ones_matrix(1, (m^m)^2)
for i in [1..(m^n)^2]: A[i,(m^(2*(m-n)))*(i-1):(m^(2*(m-n)))*i] = ones_matrix(1, m^(2*(m-n)))
E = A.right_kernel()
</code></pre>
<p>And I get the message: the kernel died...
Any suggestion to big matrices?</p>
https://ask.sagemath.org/question/55000/kernel-died-when-computing-large-matrix/?answer=55003#post-id-55003 sage: m=5;n=3
sage: ((m^n)^2+1)*((m^m)^2)
152597656250
sage: (((m^n)^2+1)*((m^m)^2)).log(10).n()
11.1835478633337
You are trying to allocate a matrix of 152 billions rationals (and some change...). Do you have enough memory for that ?Wed, 30 Dec 2020 08:55:58 +0100https://ask.sagemath.org/question/55000/kernel-died-when-computing-large-matrix/?answer=55003#post-id-55003