Hi Volker, Thank you for your message. I followed your suggestion and run Sage with gdb. I agree with you that linbox and blas (ATLAS) are suspected.
Best regards ---------------------------------------------------------------------- | Sage Version 4.6, Release Date: 2010-10-30 | | Type notebook() for the GUI, and license() for information. | ---------------------------------------------------------------------- /home/iwao/sage-4.6-linux-64bit-ubuntu_10.04.1_lts-x86_64-Linux/local/bin/sage-ipython GNU gdb (GDB) 7.1-ubuntu Copyright (C) 2010 Free Software Foundation, Inc. License GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html> This is free software: you are free to change and redistribute it. There is NO WARRANTY, to the extent permitted by law. Type "show copying" and "show warranty" for details. This GDB was configured as "x86_64-linux-gnu". For bug reporting instructions, please see: <http://www.gnu.org/software/gdb/bugs/>... Reading symbols from /home/iwao/sage-4.6-linux-64bit-ubuntu_10.04.1_lts-x86_64-Linux/local/bin/python...done. [Thread debugging using libthread_db enabled] Python 2.6.4 (r264, Oct 31 2010, 17:42:01) [GCC 4.4.3] on linux2 Type "help", "copyright", "credits" or "license" for more information. sage: identity_matrix(103).det() Program received signal SIGSEGV, Segmentation fault. 0x00007fffec9f818d in ATL_dJIK40x40x40TN40x40x0_a1_b1 () from /home/iwao/sage-4.6-linux-64bit-ubuntu_10.04.1_lts-x86_64-Linux/local/lib/libatlas.so (gdb) bt #0 0x00007fffec9f818d in ATL_dJIK40x40x40TN40x40x0_a1_b1 () from /home/iwao/sage-4.6-linux-64bit-ubuntu_10.04.1_lts-x86_64-Linux/local/lib/libatlas.so #1 0x00007fffecada34a in ATL_dmmJIK2 () from /home/iwao/sage-4.6-linux-64bit-ubuntu_10.04.1_lts-x86_64-Linux/local/lib/libatlas.so #2 0x00007fffecadadea in ATL_dmmJIK () from /home/iwao/sage-4.6-linux-64bit-ubuntu_10.04.1_lts-x86_64-Linux/local/lib/libatlas.so #3 0x00007fffecad2366 in ATL_dgemm () from /home/iwao/sage-4.6-linux-64bit-ubuntu_10.04.1_lts-x86_64-Linux/local/lib/libatlas.so #4 0x00007fffd8467d76 in ClassicMatmul<LinBox::UnparametricField<double> > ( F=..., ta=<value optimized out>, tb=<value optimized out>, m=<value optimized out>, n=52, k=<value optimized out>, alpha=46336, A=0x4404888, lda=103, B=0x43fa5f8, ldb=103, beta=<value optimized out>, C=0x4404a20, ldc=103, kmax=51, base=LinBox::FFLAS::FflasDouble) at ../../linbox/fflas/fflas_fgemm.inl:244 #5 LinBox::FFLAS::ClassicMatmul<LinBox::Modular<double> > (F=..., ta=<value optimized out>, tb=<value optimized out>, m=<value optimized out>, n=52, k=<value optimized out>, alpha=46336, A=0x4404888, lda=103, B=0x43fa5f8, ldb=103, beta=<value optimized out>, C=0x4404a20, ldc=103, kmax=51, base=LinBox::FFLAS::FflasDouble) at ../../linbox/fflas/fflas_fgemm.inl:422 ---Type <return> to continue, or q <return> to quit--- #6 0x00007fffd84aec64 in LinBox::FFLAS::WinoMain<LinBox::Modular<double> > ( F=..., ta=40, tb=<value optimized out>, m=52, n=52, k=51, alpha=46336, A=0x4404888, lda=103, B=0x43fa5f8, ldb=103, beta=1, C=0x4404a20, ldc=103, kmax=51, w=0, base=LinBox::FFLAS::FflasDouble) at ../../linbox/fflas/fflas_fgemm.inl:1276 #7 0x00007fffd84b1ecd in LinBox::FFLAS::fgemm<LinBox::Modular<double> > ( F=..., ta=<value optimized out>, tb=<value optimized out>, m=52, n=51, k=51, alpha=1, A=0x4404888, lda=103, B=0x43fa5f8, ldb=103, beta=1, C=0x4404a20, ldc=103) at ../../linbox/fflas/fflas.h:294 #8 0x00007fffd842fce8 in LUdivine<LinBox::Modular<double> > (F=..., Diag=<value optimized out>, trans=<value optimized out>, M=<value optimized out>, N=<value optimized out>, A=<value optimized out>, lda=103, P=0x43c5020, Q=0x43c5360, LuTag=LinBox::FFPACK::FfpackLQUP, cutoff=<value optimized out>) at ../../linbox/ffpack/ffpack_ludivine.inl:511 #9 0x00007fffd8430428 in Rank<LinBox::Modular<double> > ( modulus=<value optimized out>, matrix=<value optimized out>, nrows=103, ncols=103) at ../../linbox/ffpack/ffpack.h:84 #10 linbox_modn_dense_rank (modulus=<value optimized out>, matrix=<value optimized out>, nrows=103, ncols=103) at linbox-sage.C:119 #11 0x00007fffd8754fd9 in __pyx_f_4sage_4libs_6linbox_6linbox_17Linbox_modn_dense_rank (__pyx_v_self=<value optimized out>) at sage/libs/linbox/linbox.cpp:2569 #12 0x00007fffda6547ba in __pyx_pf_4sage_6matrix_17matrix_modn_dense_17Matrix_m---Type <return> to continue, or q <return> to quit--- odn_dense_rank (__pyx_v_self=<value optimized out>, unused=<value optimized out>) at sage/matrix/matrix_modn_dense.c:11420 #13 0x00007ffff7a6a163 in PyObject_Call (func=0x4174248, arg=0x28, kw=0x28) at Objects/abstract.c:2492 #14 0x00007fffd9f8ce87 in __pyx_pf_4sage_6matrix_20matrix_integer_dense_20Matrix_integer_dense__rank_modp (__pyx_v_self=<value optimized out>, __pyx_args=<value optimized out>, __pyx_kwds=<value optimized out>) at sage/matrix/matrix_integer_dense.c:24877 #15 0x00007ffff7a6a163 in PyObject_Call (func=0x4174200, arg=0x28, kw=0x28) at Objects/abstract.c:2492 #16 0x00007fffd9fa2365 in __pyx_pf_4sage_6matrix_20matrix_integer_dense_20Matrix_integer_dense_rank (__pyx_v_self=0x3d70290, unused=<value optimized out>) at sage/matrix/matrix_integer_dense.c:24590 #17 0x00007ffff7b0e2ed in call_function (f=0x43bf090, throwflag=<value optimized out>) at Python/ceval.c:3690 #18 PyEval_EvalFrameEx (f=0x43bf090, throwflag=<value optimized out>) at Python/ceval.c:2389 #19 0x00007ffff7b1026d in PyEval_EvalCodeEx (co=0x419d918, globals=<value optimized out>, locals=<value optimized out>, args=0x2, argcount=<value optimized out>, kws=0x414f3d8, kwcount=2, defs=0x414f3c8, defcount=2, closure=0x0) at Python/ceval.c:2968 #20 0x00007ffff7a967cb in function_call (func=0x7ffff652bd70, arg=0x7ffff7f5f6d0, kw=0xb93c80) at Objects/funcobject.c:524 #21 0x00007ffff7a6a163 in PyObject_Call (func=0x7ffff652bd70, arg=0x28, ---Type <return> to continue, or q <return> to quit--- kw=0x28) at Objects/abstract.c:2492 #22 0x00007ffff7b086c3 in PyEval_CallObjectWithKeywords (func=0x7ffff652bd70, arg=0x7ffff7f5f6d0, kw=0x28) at Python/ceval.c:3575 #23 0x00007fffd9f84c32 in __pyx_pf_4sage_6matrix_20matrix_integer_dense_20Matrix_integer_dense_determinant (__pyx_v_self=0x3d70290, __pyx_args=<value optimized out>, __pyx_kwds=<value optimized out>) at sage/matrix/matrix_integer_dense.c:25287 #24 0x00007ffff7a6a163 in PyObject_Call (func=0x4174290, arg=0x28, kw=0x28) at Objects/abstract.c:2492 #25 0x00007ffff7b086c3 in PyEval_CallObjectWithKeywords (func=0x4174290, arg=0x7ffff7f96050, kw=0x28) at Python/ceval.c:3575 #26 0x00007fffdb5866e1 in __pyx_pf_4sage_6matrix_7matrix2_6Matrix_det ( __pyx_v_self=<value optimized out>, __pyx_args=0x7ffff7f96050, __pyx_kwds=<value optimized out>) at sage/matrix/matrix2.c:7757 #27 0x00007ffff7b0e592 in call_function (f=0x4284610, throwflag=<value optimized out>) at Python/ceval.c:3706 #28 PyEval_EvalFrameEx (f=0x4284610, throwflag=<value optimized out>) at Python/ceval.c:2389 #29 0x00007ffff7b1026d in PyEval_EvalCodeEx (co=0x41a1b70, globals=<value optimized out>, locals=<value optimized out>, args=0x0, argcount=<value optimized out>, kws=0x0, kwcount=0, defs=0x0, defcount=0, closure=0x0) at Python/ceval.c:2968 #30 0x00007ffff7b10342 in PyEval_EvalCode (co=0x28, globals=0x28, locals=0x28) at Python/ceval.c:522 ---Type <return> to continue, or q <return> to quit--- #31 0x00007ffff7b0f619 in exec_statement (f=0x4284220, throwflag=<value optimized out>) at Python/ceval.c:4401 #32 PyEval_EvalFrameEx (f=0x4284220, throwflag=<value optimized out>) at Python/ceval.c:1717 #33 0x00007ffff7b1026d in PyEval_EvalCodeEx (co=0xac6918, globals=<value optimized out>, locals=<value optimized out>, args=0x4283b20, argcount=<value optimized out>, kws=0x2, kwcount=0, defs=0x0, defcount=0, closure=0x0) at Python/ceval.c:2968 #34 0x00007ffff7b0e60d in fast_function (f=0x4283980, throwflag=<value optimized out>) at Python/ceval.c:3802 #35 call_function (f=0x4283980, throwflag=<value optimized out>) at Python/ceval.c:3727 #36 PyEval_EvalFrameEx (f=0x4283980, throwflag=<value optimized out>) at Python/ceval.c:2389 #37 0x00007ffff7b1026d in PyEval_EvalCodeEx (co=0xac6828, globals=<value optimized out>, locals=<value optimized out>, args=0x2, argcount=<value optimized out>, kws=0xb29d20, kwcount=0, defs=0xb29d10, defcount=2, closure=0x0) at Python/ceval.c:2968 #38 0x00007ffff7b0e60d in fast_function (f=0x6d83f0, throwflag=<value optimized out>) at Python/ceval.c:3802 #39 call_function (f=0x6d83f0, throwflag=<value optimized out>) at Python/ceval.c:3727 #40 PyEval_EvalFrameEx (f=0x6d83f0, throwflag=<value optimized out>) at Python/ceval.c:2389 ---Type <return> to continue, or q <return> to quit--- #41 0x00007ffff7b0ee88 in fast_function (f=0x6da190, throwflag=<value optimized out>) at Python/ceval.c:3792 #42 call_function (f=0x6da190, throwflag=<value optimized out>) at Python/ceval.c:3727 #43 PyEval_EvalFrameEx (f=0x6da190, throwflag=<value optimized out>) at Python/ceval.c:2389 #44 0x00007ffff7b1026d in PyEval_EvalCodeEx (co=0xac63f0, globals=<value optimized out>, locals=<value optimized out>, args=0xb93c28, argcount=<value optimized out>, kws=0x2, kwcount=0, defs=0xb27c68, defcount=1, closure=0x0) at Python/ceval.c:2968 #45 0x00007ffff7b0e60d in fast_function (f=0xb93aa0, throwflag=<value optimized out>) at Python/ceval.c:3802 #46 call_function (f=0xb93aa0, throwflag=<value optimized out>) at Python/ceval.c:3727 #47 PyEval_EvalFrameEx (f=0xb93aa0, throwflag=<value optimized out>) at Python/ceval.c:2389 #48 0x00007ffff7b1026d in PyEval_EvalCodeEx (co=0xac6120, globals=<value optimized out>, locals=<value optimized out>, args=0x681e60, argcount=<value optimized out>, kws=0x2, kwcount=0, defs=0xb27c28, defcount=1, closure=0x0) at Python/ceval.c:2968 #49 0x00007ffff7b0e60d in fast_function (f=0x681cd0, throwflag=<value optimized out>) at Python/ceval.c:3802 #50 call_function (f=0x681cd0, throwflag=<value optimized out>) at Python/ceval.c:3727 ---Type <return> to continue, or q <return> to quit--- #51 PyEval_EvalFrameEx (f=0x681cd0, throwflag=<value optimized out>) at Python/ceval.c:2389 #52 0x00007ffff7b1026d in PyEval_EvalCodeEx (co=0x7ffff7e2ba08, globals=<value optimized out>, locals=<value optimized out>, args=0x2, argcount=<value optimized out>, kws=0x8c1780, kwcount=2, defs=0x8c1770, defcount=2, closure=0x0) at Python/ceval.c:2968 #53 0x00007ffff7b0e60d in fast_function (f=0x6bfb10, throwflag=<value optimized out>) at Python/ceval.c:3802 #54 call_function (f=0x6bfb10, throwflag=<value optimized out>) at Python/ceval.c:3727 #55 PyEval_EvalFrameEx (f=0x6bfb10, throwflag=<value optimized out>) at Python/ceval.c:2389 #56 0x00007ffff7b1026d in PyEval_EvalCodeEx (co=0x7ffff7f04558, globals=<value optimized out>, locals=<value optimized out>, args=0x0, argcount=<value optimized out>, kws=0x0, kwcount=0, defs=0x0, defcount=0, closure=0x0) at Python/ceval.c:2968 #57 0x00007ffff7b10342 in PyEval_EvalCode (co=0x28, globals=0x28, locals=0x28) at Python/ceval.c:522 #58 0x00007ffff7b304e0 in run_mod (fp=0x661b30, filename=<value optimized out>, start=<value optimized out>, globals=<value optimized out>, locals=0x63d270, closeit=0, flags=0x7fffffffd9e0) at Python/pythonrun.c:1335 #59 PyRun_FileExFlags (fp=0x661b30, filename=<value optimized out>, start=<value optimized out>, globals=<value optimized out>, ---Type <return> to continue, or q <return> to quit--- locals=0x63d270, closeit=0, flags=0x7fffffffd9e0) at Python/pythonrun.c:1321 #60 0x00007ffff7b306ac in PyRun_SimpleFileExFlags (fp=<value optimized out>, filename=0x7fffffffea5d "/home/iwao/sage-4.6-linux-64bit-ubuntu_10.04.1_lts-x86_64-Linux/local/bin/sage-ipython", closeit=0, flags=0x7fffffffd9e0) at Python/pythonrun.c:931 #61 0x00007ffff7b3cb77 in RunStartupFile (argc=<value optimized out>, argv=<value optimized out>) at Modules/main.c:142 #62 Py_Main (argc=<value optimized out>, argv=<value optimized out>) at Modules/main.c:558 #63 0x00007ffff6e17c4d in __libc_start_main () from /lib/libc.so.6 #64 0x00000000004006d9 in _start () (gdb) (gdb) quit A debugging session is active. Inferior 1 [process 14077] will be killed. Quit anyway? (y or n) y i...@octa:~$ exit On Tue, Dec 28, 2010 at 7:43 PM, Volker Braun <vbraun.n...@gmail.com> wrote: > You are not running out of memory in that computation if you have 3GB of > ram. It seems like something in linbox (which is used for matrices of that > size) dies during the computation. I also recently ran into some issues with > linbox and blas, and was wondering if they are related. Can you run sage in > the debugger (start with "sage -gdb"), run your example, and then when it > crashes post the output of "bt" (stack backtrace) in gdb? > Volker > > -- > To post to this group, send email to sage-support@googlegroups.com > To unsubscribe from this group, send email to > sage-support+unsubscr...@googlegroups.com > For more options, visit this group at > http://groups.google.com/group/sage-support > URL: http://www.sagemath.org > -- --- Iwao KIMURA Dept. Math., University of Toyama, Japan. -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org
