bugGNU Octave - Bugs: bug #45577, Solving systems using complex NaN...

 
 

bug #45577: Solving systems using complex NaN causes crash

Submitter:  Lachlan Andrew <lachlan>
Submitted:  Sun 19 Jul 2015 10:31:18 AM UTC
   
 
Category:  Octave Function Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Segfault, Bus Error, etc.
Status:  Fixed Assigned to:  None
Originator Name:  Lachlan Open/Closed:  * Closed
Release:  * dev Operating System:  * GNU/Linux
Fixed Release:  None Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Jump to the original submission

Fri 15 Jul 2016 11:56:24 PM UTC, comment #6: 

It works for me too.  Closing as fixed.

Lachlan Andrew <lachlan>
Fri 15 Jul 2016 10:58:25 PM UTC, comment #5: 

@Lachlan: Can you re-verify that


inv (single (coeff))


still causes a segfault?  Use anything after


changeset:   22085:bf5fbf347aaf
user:        Kai T. Ohlhus <k.ohlhus@gmail.com>
date:        Mon Jul 11 09:27:10 2016 +0200
summary:     Avoid segfault with complex matrices using LAPACK (bug #46330).


It works for me.

Rik <rik5>
Group administrator
Fri 08 Jul 2016 01:19:09 PM UTC, comment #4: 

I'm re-opening this because calculating


inv (single (coeff))


still crashes, with coeff given in the original post.

I'll copy Rik's changes into fCMatrix.cc soon.

Lachlan Andrew <lachlan>
Sat 10 Oct 2015 11:48:52 PM UTC, comment #3: 

The problem seems not to be in Octave, but rather in the ZGETRF function from lapack.  The code which calls ZGETRF is in liboctave/array/CMatrix.cc.  In Octave we use the exact same code sequence for double matrices, single matrices, single complex matrices, and double complex matrices.  There are only slight changes to reflect the different inputs.  Thus, for single complex matrices the routine is CGETRF and for regular double matrices it is DGETRF.  None of the other routines produce segfaults if the input matrix has a row or column of NaN values--it is only double complex matrices which call ZGETRF.

Because it is only one routine, I have worked around the issue in Octave.  The norm of the matrix is calculated before the LU decomposition and if the norm is NaN that is an indication not to call lapack and, instead, issue a warning.

Interestingly, in searching through CMatrix.cc for calls to ZGETRF I found three other instances: finverse, determinant, rcond.  Sure enough, trying to calculate the determinant or rcond of the problem matrix in this bug report produces a segfault.  So I patched those too.

Since this is a serious crash of the interpreter I fixed the problem on the stable branch.  See this changeset (http://hg.savannah.gnu.org/hgweb/octave/rev/ffc6cdcd02c5).

Closing report.

Rik <rik5>
Group administrator
Sat 12 Sep 2015 09:59:57 PM UTC, comment #2: 

This is not quite the same as bug #45484, that bug is about BLAS errors that are caught correctly on Linux. This error causes a crash even on Linux, probably a different error. It looks like a memory double-free error in a destructor:


#0  0x00007ffff4640107 in __GI_raise (sig=sig@entry=6) at ../nptl/sysdeps/unix/sysv/linux/raise.c:56
#1  0x00007ffff46414e8 in __GI_abort () at abort.c:89
#2  0x00007ffff467e214 in __libc_message (do_abort=do_abort@entry=1, fmt=fmt@entry=0x7ffff4771000 "*** Error in `%s': %s: 0x%s ***\n") at ../sysdeps/posix/libc_fatal.c:175
#3  0x00007ffff46839ee in malloc_printerr (action=1, str=0x7ffff4771108 "double free or corruption (!prev)", ptr=<optimized out>) at malloc.c:4996
#4  0x00007ffff46846f6 in _int_free (av=<optimized out>, p=<optimized out>, have_lock=0) at malloc.c:3840
#5  0x00007ffff6de13a9 in Array<std::complex<double> >::ArrayRep::~ArrayRep (this=0xc0abc0, __in_chrg=<optimized out>) at ../liboctave/array/Array.h:89
#6  Array<std::complex<double> >::~Array (this=0x7fffffffad10, __in_chrg=<optimized out>) at ../liboctave/array/Array.h:223
#7  0x00007ffff5ced2e4 in MArray<std::complex<double> >::~MArray (this=0x7fffffffad10, __in_chrg=<optimized out>) at ../liboctave/array/MArray.h:63
#8  ComplexNDArray::~ComplexNDArray (this=0x7fffffffad10, __in_chrg=<optimized out>) at ../liboctave/array/CNDArray.h:34
#9  ComplexMatrix::~ComplexMatrix (this=0x7fffffffad10, __in_chrg=<optimized out>) at ../liboctave/array/CMatrix.h:38
#10 ComplexMatrix::fsolve (this=this@entry=0x7fffffffb610, mattype=..., b=..., info=@0x7fffffffb594: -2, rcon=@0x7fffffffb598: 0, sing_handler=0x7ffff785d730 <solve_singularity_warning(double)>, calc_cond=true)
    at ../liboctave/array/CMatrix.cc:2215
#11 0x00007ffff5cf2752 in ComplexMatrix::solve (this=this@entry=0x7fffffffb610, mattype=..., b=..., info=@0x7fffffffb594: -2, rcon=@0x7fffffffb598: 0, sing_handler=0x7ffff785d730 <solve_singularity_warning(double)>,
    singular_fallback=true, transt=blas_no_trans) at ../liboctave/array/CMatrix.cc:2372
#12 0x00007ffff7862411 in xleftdiv (a=..., b=..., typ=..., transt=transt@entry=blas_no_trans) at ../libinterp/corefcn/xdiv.cc:402
#13 0x00007ffff6df4afb in oct_binop_ldiv (a1=..., a2=...) at ../libinterp/operators/op-cm-cm.cc:107
#14 0x00007ffff728fb49 in do_binary_op (op=octave_value::op_ldiv, v1=..., v2=...) at ../libinterp/octave-value/ov.cc:2021
#15 0x00007ffff72f24d6 in tree_binary_expression::rvalue1 (this=0x8a2550) at ../libinterp/parse-tree/pt-binop.cc:146
#16 0x00007ffff7305c7e in tree_evaluator::visit_statement (this=0x7ffff7dd4440 <std_evaluator>, stmt=...) at ../libinterp/parse-tree/pt-eval.cc:744
#17 0x00007ffff7304f49 in tree_evaluator::visit_statement_list (this=0x7ffff7dd4440 <std_evaluator>, lst=...) at ../libinterp/parse-tree/pt-eval.cc:794
#18 0x00007ffff7818310 in main_loop () at ../libinterp/corefcn/toplev.cc:599
#19 0x00007ffff6dc90bf in octave_execute_interpreter () at ../libinterp/octave.cc:919
#20 0x00007ffff462cb45 in __libc_start_main (main=0x400920 <main(int, char**)>, argc=9, argv=0x7fffffffc048, init=<optimized out>, fini=<optimized out>, rtld_fini=<optimized out>, stack_end=0x7fffffffc038) at libc-start.c:287
#21 0x00000000004009af in _start ()


Mike Miller <mtmiller>
Group Member
Wed 22 Jul 2015 04:57:31 PM UTC, comment #1: 

Marking as a duplicate of bug #45484.

Rik <rik5>
Group administrator
Sun 19 Jul 2015 10:31:18 AM UTC, original submission:  

The following code causes Octave 4.0.0 to crash:


coeffs = [
   NaN - NaN*i   NaN - NaN*i   NaN - NaN*i   NaN - NaN*i
   NaN - NaN*i   NaN - NaN*i   NaN - NaN*i   NaN - NaN*i
   NaN + NaN*i     0 +   0i     0 +   0i     0 +   0i
     0 +   0i   NaN + NaN*i     0 +   0i     0 +   0i
];

rhs = [ NaN - NaN*i
   NaN - NaN*i
   NaN - NaN*i
   NaN - NaN*i];

coeffs \ rhs


The diagnostics are:


octave: no graphical display found
octave: disabling GUI features
warning: matrix singular to machine precision
warning: called from
    crash_small.m at line 15 column 8
*** Error in `/home/lachlan/src/octave-4.0.0-play/src/.libs/lt-octave-cli': double free or corruption (out): 0x000000000178abf0 ***
======= Backtrace: =========
/lib/x86_64-linux-gnu/libc.so.6(+0x7aa16)[0x7f97093d8a16]
/lib/x86_64-linux-gnu/libc.so.6(+0x7b793)[0x7f97093d9793]
/home/lachlan/src/octave-4.0.0-play/libinterp/.libs/liboctinterp.so.3(_ZN5ArrayISt7complexIdEED1Ev+0x39)[0x7f970b8e0809]
/home/lachlan/src/octave-4.0.0-play/liboctave/.libs/liboctave.so.3(_ZNK13ComplexMatrix6fsolveER10MatrixTypeRKS_RiRdPFvdEb+0x36b)[0x7f970a7c046b]
/home/lachlan/src/octave-4.0.0-play/liboctave/.libs/liboctave.so.3(_ZNK13ComplexMatrix5solveER10MatrixTypeRKS_RiRdPFvdEb15blas_trans_type+0x4d5)[0x7f970a7c5bd5]
/home/lachlan/src/octave-4.0.0-play/libinterp/.libs/liboctinterp.so.3(_Z8xleftdivRK13ComplexMatrixS1_R10MatrixType15blas_trans_type+0x69)[0x7f970c12d569]
/home/lachlan/src/octave-4.0.0-play/libinterp/.libs/liboctinterp.so.3(+0x30a23d)[0x7f970b8ed23d]
/home/lachlan/src/octave-4.0.0-play/libinterp/.libs/liboctinterp.so.3(_Z12do_binary_opN12octave_value9binary_opERKS_S2_+0xa8)[0x7f970bc959e8]
/home/lachlan/src/octave-4.0.0-play/libinterp/.libs/liboctinterp.so.3(_ZN22tree_binary_expression7rvalue1Ei+0x263)[0x7f970bcefe03]
/home/lachlan/src/octave-4.0.0-play/libinterp/.libs/liboctinterp.so.3(_ZN14tree_evaluator15visit_statementER14tree_statement+0x108)[0x7f970bcf97b8]
/home/lachlan/src/octave-4.0.0-play/libinterp/.libs/liboctinterp.so.3(_ZN14tree_evaluator20visit_statement_listER19tree_statement_list+0x51)[0x7f970bcf8e01]
/home/lachlan/src/octave-4.0.0-play/libinterp/.libs/liboctinterp.so.3(_ZN18octave_user_script17do_multi_index_opEiRK17octave_value_list+0x26a)[0x7f970bc8c07a]
/home/lachlan/src/octave-4.0.0-play/libinterp/.libs/liboctinterp.so.3(_Z11source_fileRKSsS0_bbS0_+0x2ab)[0x7f970bd37f5b]
/home/lachlan/src/octave-4.0.0-play/libinterp/.libs/liboctinterp.so.3(+0x2e8eb3)[0x7f970b8cbeb3]
/home/lachlan/src/octave-4.0.0-play/libinterp/.libs/liboctinterp.so.3(octave_execute_interpreter+0x515)[0x7f970b8cf8b5]
/lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf5)[0x7f970937f995]
/home/lachlan/src/octave-4.0.0-play/src/.libs/lt-octave-cli[0x40098f]
======= Memory map: ========
00400000-00401000 r-xp 00000000 08:06 21712327                           /home/lachlan/src/octave-4.0.0-play/src/.libs/lt-octave-cli
00600000-00601000 rw-p 00000000 08:06 21712327                           /home/lachlan/src/octave-4.0.0-play/src/.libs/lt-octave-cli
011df000-017b1000 rw-p 00000000 00:00 0                                  [heap]
7f96f7a39000-7f96f9a39000 rw-p 00000000 00:00 0
7f96f9a39000-7f96f9a45000 r-xp 00000000 08:06 21712310                   /home/lachlan/src/octave-4.0.0-play/libinterp/dldfcn/__init_gnuplot__.oct
7f96f9a45000-7f96f9c45000 ---p 0000c000 08:06 21712310                   /home/lachlan/src/octave-4.0.0-play/libinterp/dldfcn/__init_gnuplot__.oct
7f96f9c45000-7f96f9c46000 rw-p 0000c000 08:06 21712310                   /home/lachlan/src/octave-4.0.0-play/libinterp/dldfcn/__init_gnuplot__.oct
7f96f9c46000-7f96f9c4b000 r-xp 00000000 08:06 21712309                   /home/lachlan/src/octave-4.0.0-play/libinterp/dldfcn/__init_fltk__.oct
7f96f9c4b000-7f96f9e4a000 ---p 00005000 08:06 21712309                   /home/lachlan/src/octave-4.0.0-play/libinterp/dldfcn/__init_fltk__.oct
7f96f9e4a000-7f96f9e4b000 rw-p 00004000 08:06 21712309                   /home/lachlan/src/octave-4.0.0-play/libinterp/dldfcn/__init_fltk__.oct
7f96f9e4b000-7f96fbe4c000 rw-p 00000000 00:00 0
7f96fbe4c000-7f96fbe58000 r-xp 00000000 08:05 2625033                    /lib/x86_64-linux-gnu/libnss_files-2.17.so
7f96fbe58000-7f96fc057000 ---p 0000c000 08:05 2625033                    /lib/x86_64-linux-gnu/libnss_files-2.17.so
7f96fc057000-7f96fc058000 r--p 0000b000 08:05 2625033                    /lib/x86_64-linux-gnu/libnss_files-2.17.so
7f96fc058000-7f96fc059000 rw-p 0000c000 08:05 2625033                    /lib/x86_64-linux-gnu/libnss_files-2.17.so
7f96fc059000-7f96fc063000 r-xp 00000000 08:05 2625037                    /lib/x86_64-linux-gnu/libnss_nis-2.17.so
7f96fc063000-7f96fc262000 ---p 0000a000 08:05 2625037                    /lib/x86_64-linux-gnu/libnss_nis-2.17.so
7f96fc262000-7f96fc263000 r--p 00009000 08:05 2625037                    /lib/x86_64-linux-gnu/libnss_nis-2.17.so
7f96fc263000-7f96fc264000 rw-p 0000a000 08:05 2625037                    /lib/x86_64-linux-gnu/libnss_nis-2.17.so
7f96fc264000-7f96fc279000 r-xp 00000000 08:05 2625027                    /lib/x86_64-linux-gnu/libnsl-2.17.so
7f96fc279000-7f96fc478000 ---p 00015000 08:05 2625027                    /lib/x86_64-linux-gnu/libnsl-2.17.so
7f96fc478000-7f96fc479000 r--p 00014000 08:05 2625027                    /lib/x86_64-linux-gnu/libnsl-2.17.so
7f96fc479000-7f96fc47a000 rw-p 00015000 08:05 2625027                    /lib/x86_64-linux-gnu/libnsl-2.17.so
7f96fc47a000-7f96fc47c000 rw-p 00000000 00:00 0
7f96fc47c000-7f96fc483000 r-xp 00000000 08:05 2625029                    /lib/x86_64-linux-gnu/libnss_compat-2.17.so
7f96fc483000-7f96fc682000 ---p 00007000 08:05 2625029                    /lib/x86_64-linux-gnu/libnss_compat-2.17.so
7f96fc682000-7f96fc683000 r--p 00006000 08:05 2625029                    /lib/x86_64-linux-gnu/libnss_compat-2.17.so
7f96fc683000-7f96fc684000 rw-p 00007000 08:05 2625029                    /lib/x86_64-linux-gnu/libnss_compat-2.17.so
7f96fc684000-7f9702684000 rw-p 00000000 00:00 0
7f9702684000-7f9702685000 ---p 00000000 00:00 0
7f9702685000-7f9702e85000 rw-p 00000000 00:00 0                          [stack:22247]
7f9702e85000-7f9702e86000 ---p 00000000 00:00 0
7f9702e86000-7f9703686000 rw-p 00000000 00:00 0                          [stack:22246]
7f9703686000-7f9703687000 ---p 00000000 00:00 0
7f9703687000-7f9703e87000 rw-p 00000000 00:00 0                          [stack:22245]
7f9703e87000-7f9703ec2000 r-xp 00000000 08:05 4204827                    /usr/lib/x86_64-linux-gnu/libquadmath.so.0.0.0
7f9703ec2000-7f97040c1000 ---p 0003b000 08:05 4204827                    /usr/lib/x86_64-linux-gnu/libquadmath.so.0.0.0
7f97040c1000-7f97040c2000 rw-p 0003a000 08:05 4204827                    /usr/lib/x86_64-linux-gnu/libquadmath.so.0.0.0
7f97040c2000-7f97041d7000 r-xp 00000000 08:05 4204318                    /usr/lib/x86_64-linux-gnu/libgfortran.so.3.0.0
7f97041d7000-7f97043d7000 ---p 00115000 08:05 4204318                    /usr/lib/x86_64-linux-gnu/libgfortran.so.3.0.0
7f97043d7000-7f97043d9000 rw-p 00115000 08:05 4204318                    /usr/lib/x86_64-linux-gnu/libgfortran.so.3.0.0
7f97043d9000-7f97043dc000 r-xp 00000000 08:05 2624979                    /lib/x86_64-linux-gnu/libdl-2.17.so
7f97043dc000-7f97045db000 ---p 00003000 08:05 2624979                    /lib/x86_64-linux-gnu/libdl-2.17.so
7f97045db000-7f97045dc000 r--p 00002000 08:05 2624979                    /lib/x86_64-linux-gnu/libdl-2.17.so
7f97045dc000-7f97045dd000 rw-p 00003000 08:05 2624979                    /lib/x86_64-linux-gnu/libdl-2.17.so
7f97045dd000-7f970461b000 r-xp 00000000 08:05 2625055                    /lib/x86_64-linux-gnu/libpcre.so.3.13.1
7f970461b000-7f970481a000 ---p 0003e000 08:05 2625055                    /lib/x86_64-linux-gnu/libpcre.so.3.13.1
7f970481a000-7f970481b000 rw-p 0003d000 08:05 2625055                    /lib/x86_64-linux-gnu/libpcre.so.3.13.1
7f970481b000-7f9704840000 r-xp 00000000 08:05 2625095                    /lib/x86_64-linux-gnu/libtinfo.so.5.9
7f9704840000-7f9704a3f000 ---p 00025000 08:05 2625095                    /lib/x86_64-linux-gnu/libtinfo.so.5.9
7f9704a3f000-7f9704a43000 r--p 00024000 08:05 2625095                    /lib/x86_64-linux-gnu/libtinfo.so.5.9
7f9704a43000-7f9704a44000 rw-p 00028000 08:05 2625095                    /lib/x86_64-linux-gnu/libtinfo.so.5.9
7f9704a44000-7f9704a81000 r-xp 00000000 08:05 2625075                    /lib/x86_64-linux-gnu/libreadline.so.6.2
7f9704a81000-7f9704c81000 ---p 0003d000 08:05 2625075                    /lib/x86_64-linux-gnu/libreadline.so.6.2
7f9704c81000-7f9704c83000 r--p 0003d000 08:05 2625075                    /lib/x86_64-linux-gnu/libreadline.so.6.2
7f9704c83000-7f9704c89000 rw-p 0003f000 08:05 2625075                    /lib/x86_64-linux-gnu/libreadline.so.6.2
7f9704c89000-7f9704c8b000 rw-p 00000000 00:00 0
7f9704c8b000-7f9706269000 r-xp 00000000 08:05 4206214                    /usr/lib/openblas-base/libopenblas.so.0
7f9706269000-7f9706469000 ---p 015de000 08:05 4206214                    /usr/lib/openblas-base/libopenblas.so.0
7f9706469000-7f9706480000 rw-p 015de000 08:05 4206214                    /usr/lib/openblas-base/libopenblas.so.0
7f9706480000-7f970648b000 rw-p 00000000 00:00 0
7f970648b000-7f97064b7000 r-xp 00000000 08:05 4205398                    /usr/lib/libqrupdate.so.1.1
7f97064b7000-7f97066b6000 ---p 0002c000 08:05 4205398                    /usr/lib/libqrupdate.so.1.1
7f97066b6000-7f97066b7000 rw-p 0002b000 08:05 4205398                    /usr/lib/libqrupdate.so.1.1
7f97066b7000-7f9706704000 r-xp 00000000 08:05 4205387                    /usr/lib/libarpack.so.2.0.0
7f9706704000-7f9706903000 ---p 0004d000 08:05 4205387                    /usr/lib/libarpack.so.2.0.0
7f9706903000-7f9706905000 rw-p 0004c000 08:05 4205387                    /usr/lib/libarpack.so.2.0.0
7f9706905000-7f9706937000 r-xp 00000000 08:05 4205392                    /usr/lib/libcxsparse.so.2.2.3
7f9706937000-7f9706b36000 ---p 00032000 08:05 4205392                    /usr/lib/libcxsparse.so.2.2.3
7f9706b36000-7f9706b37000 rw-p 00031000 08:05 4205392                    /usr/lib/libcxsparse.so.2.2.3
7f9706b37000-7f9706b42000 r-xp 00000000 08:05 4205390                    /usr/lib/libccolamd.so.2.7.1
7f9706b42000-7f9706d41000 ---p 0000b000 08:05 4205390                    /usr/lib/libccolamd.so.2.7.1
7f9706d41000-7f9706d42000 rw-p 0000a000 08:05 4205390                    /usr/lib/libccolamd.so.2.7.1
7f9706d42000-7f9706d48000 r-xp 00000000 08:05 4198398                    /usr/lib/libcolamd.so.2.7.1
7f9706d48000-7f9706f48000 ---p 00006000 08:05 4198398                    /usr/lib/libcolamd.so.2.7.1
7f9706f48000-7f9706f49000 rw-p 00006000 08:05 4198398                    /usr/lib/libcolamd.so.2.7.1
7f9706f49000-7f9706f53000 r-xp 00000000 08:05 4205389                    /usr/lib/libcamd.so.2.2.0
7f9706f53000-7f9707152000 ---p 0000a000 08:05 4205389                    /usr/lib/libcamd.so.2.2.0
7f9707152000-7f9707153000 rw-p 00009000 08:05 4205389                    /usr/lib/libcamd.so.2.2.0
7f9707153000-7f970715c000 r-xp 00000000 08:05 4205386                    /usr/lib/libamd.so.2.2.0
7f970715c000-7f970735b000 ---p 00009000 08:05 4205386                    /usr/lib/libamd.so.2.2.0
7f970735b000-7f970735c000 rw-p 00008000 08:05 4205386                    /usr/lib/libamd.so.2.2.0
7f970735c000-7f970740f000 r-xp 00000000 08:05 4205400                    /usr/lib/libumfpack.so.5.4.0
7f970740f000-7f970760f000 ---p 000b3000 08:05 4205400                    /usr/lib/libumfpack.so.5.4.0
7f970760f000-7f9707610000 rw-p 000b3000 08:05 4205400                    /usr/lib/libumfpack.so.5.4.0
7f9707610000-7f97076f8000 r-xp 00000000 08:05 4205391                    /usr/lib/libcholmod.so.1.7.1
7f97076f8000-7f97078f7000 ---p 000e8000 08:05 4205391                    /usr/lib/libcholmod.so.1.7.1
7f97078f7000-7f97078f8000 rw-p 000e7000 08:05 4205391                    /usr/lib/libcholmod.so.1.7.1
7f97078f8000-7f9707e8e000 r-xp 00000000 08:05 4201263                    /usr/lib/lapack/liblapack.so.3.0
7f9707e8e000-7f970808d000 ---p 00596000 08:05 4201263                    /usr/lib/lapack/liblapack.so.3.0
7f970808d000-7f970808e000 r--p 00595000 08:05 4201263                    /usr/lib/lapack/liblapack.so.3.0
7f970808e000-7f9708092000 rw-p 00596000 08:05 4201263                    /usr/lib/lapack/liblapack.so.3.0
7f9708092000-7f9708288000 r-xp 00000000 08:05 4204213                    /usr/lib/x86_64-linux-gnu/libfftw3f.so.3.3.2
7f9708288000-7f9708487000 ---p 001f6000 08:05 4204213                    /usr/lib/x86_64-linux-gnu/libfftw3f.so.3.3.2
7f9708487000-7f9708498000 r--p 001f5000 08:05 4204213                    /usr/lib/x86_64-linux-gnu/libfftw3f.so.3.3.2
7f9708498000-7f9708499000 rw-p 00206000 08:05 4204213                    /usr/lib/x86_64-linux-gnu/libfftw3f.so.3.3.2
7f9708499000-7f970849f000 r-xp 00000000 08:05 4204217                    /usr/lib/x86_64-linux-gnu/libfftw3f_threads.so.3.3.2
7f970849f000-7f970869e000 ---p 00006000 08:05 4204217                    /usr/lib/x86_64-linux-gnu/libfftw3f_threads.so.3.3.2
7f970869e000-7f970869f000 r--p 00005000 08:05 4204217                    /usr/lib/x86_64-linux-gnu/libfftw3f_threads.so.3.3.2
7f970869f000-7f97086a0000 rw-p 00006000 08:05 4204217                    /usr/lib/x86_64-linux-gnu/libfftw3f_threads.so.3.3.2
7f97086a0000-7f970888a000 r-xp 00000000 08:05 4204207                    /usr/lib/x86_64-linux-gnu/libfftw3.so.3.3.2
7f970888a000-7f9708a89000 ---p 001ea000 08:05 4204207                    /usr/lib/x86_64-linux-gnu/libfftw3.so.3.3.2
7f9708a89000-7f9708a9a000 r--p 001e9000 08:05 4204207                    /usr/lib/x86_64-linux-gnu/libfftw3.so.3.3.2
7f9708a9a000-7f9708a9b000 rw-p 001fa000 08:05 4204207                    /usr/lib/x86_64-linux-gnu/libfftw3.so.3.3.2
7f9708a9b000-7f9708aa1000 r-xp 00000000 08:05 4204211                    /usr/lib/x86_64-linux-gnu/libfftw3_threads.so.3.3.2
7f9708aa1000-7f9708ca0000 ---p 00006000 08:05 4204211                    /usr/lib/x86_64-linux-gnu/libfftw3_threads.so.3.3.2
7f9708ca0000-7f9708ca1000 r--p 00005000 08:05 4204211                    /usr/lib/x86_64-linux-gnu/libfftw3_threads.so.3.3.2
7f9708ca1000-7f9708ca2000 rw-p 00006000 08:05 4204211                    /usr/lib/x86_64-linux-gnu/libfftw3_threads.so.3.3.2
7f9708ca2000-7f9708cb9000 r-xp 00000000 08:05 2625116                    /lib/x86_64-linux-gnu/libz.so.1.2.8
7f9708cb9000-7f9708eb8000 ---p 00017000 08:05 2625116                    /lib/x86_64-linux-gnu/libz.so.1.2.8
7f9708eb8000-7f9708eb9000 r--p 00016000 08:05 2625116                    /lib/x86_64-linux-gnu/libz.so.1.2.8
7f9708eb9000-7f9708eba000 rw-p 00017000 08:05 2625116                    /lib/x86_64-linux-gnu/libz.so.1.2.8
7f9708eba000-7f9709151000 r-xp 00000000 08:05 4205311                    /usr/lib/x86_64-linux-gnu/libhdf5.so.7.0.0
7f9709151000-7f9709351000 ---p 00297000 08:05 4205311                    /usr/lib/x86_64-linux-gnu/libhdf5.so.7.0.0
7f9709351000-7f9709356000 r--p 00297000 08:05 4205311                    /usr/lib/x86_64-linux-gnu/libhdf5.so.7.0.0
7f9709356000-7f970935d000 rw-p 0029c000 08:05 4205311                    /usr/lib/x86_64-linux-gnu/libhdf5.so.7.0.0
7f970935d000-7f970935e000 rw-p 00000000 00:00 0
7f970935e000-7f9709501000 r-xp 00000000 08:05 2624965                    /lib/x86_64-linux-gnu/libc-2.17.sopanic: Aborted -- stopping myself...
attempting to save variables to 'octave-workspace'...
save to 'octave-workspace' complete
Aborted



I'm not sure if the issue is related to the following error.  The following code


coeffs = [
   NaN - NaN*i   NaN - NaN*i   NaN - NaN*i
   NaN + NaN*i     0 +   0i     0 +   0i
     0 +   0i   NaN + NaN*i     0 +   0i
];


rhs = [ NaN - NaN*i
   NaN - NaN*i
   NaN - NaN*i];

coeffs \ rhs

gives the error

warning: matrix singular to machine precision
 ** On entry to DLASCL parameter number  4 had an illegal value
error: exception encountered in Fortran subroutine zgelsd_


Lachlan Andrew <lachlan>

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach Files:
   
   
Comment:
   

Attached Files
file #34464:  octave-workspace added by lachlan (422B - application/octet-stream)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by mtmiller (Posted a comment)
  • -email is unavailable- added by rik5 (Posted a comment)
  • -email is unavailable- added by lachlan (Submitted the item)
  •  

    There are 0 votes so far. Votes easily highlight which items people would like to see resolved in priority, independently of the priority of the item set by tracker managers.

    Only group members can vote.

     

    Follow 12 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2016-07-15 lachlan StatusReady For Test Fixed
        Open/ClosedOpen Closed
    2016-07-15 rik5 StatusFixed Ready For Test
    2016-07-08 lachlan Open/ClosedClosed Open
        Release4.0.0 dev
    2015-10-10 rik5 StatusConfirmed Fixed
        Open/ClosedOpen Closed
    2015-09-12 mtmiller DependenciesRemoved dependency to bugs #45484 -
    2015-09-12 mtmiller StatusDuplicate Confirmed
    2015-07-22 rik5 StatusNone Duplicate
    2015-07-22 rik5 Dependencies- Depends on bugs #45484
    2015-07-19 lachlan Attached File- Added octave-workspace, #34464

    Back to the top

    Powered by Savane 3.13-4448.
    Corresponding source code