|
| From: | Markus Bergholz |
| Subject: | [Octave-bug-tracker] [bug #45651] Crash when inverting matrices sized 10x10 or bigger |
| Date: | Wed, 29 Jul 2015 19:36:08 +0000 |
| User-agent: | Mozilla/5.0 (X11; Linux x86_64; rv:39.0) Gecko/20100101 Firefox/39.0 |
Follow-up Comment #1, bug #45651 (project octave):
I can not confirm this, neither with 3.8.2 nor with 4.0.0
21:46 $ /usr/bin/octave
GNU Octave, version 3.8.2
Copyright (C) 2014 John W. Eaton and others.
This is free software; see the source code for copying conditions.
There is ABSOLUTELY NO WARRANTY; not even for MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. For details, type 'warranty'.
Octave was configured for "x86_64-unknown-linux-gnu".
Additional information about Octave is available at http://www.octave.org.
Please contribute if you find this software useful.
For more information, visit http://www.octave.org/get-involved.html
Read http://www.octave.org/bugs.html to learn how to submit bug reports.
For information about changes from previous versions, type 'news'.
octave:1> N=12;
octave:2> A=zeros(N);
octave:3> for n=2:N-1
> A(n,n-1)= 1;
> A(n,n )=-2;
> A(n,n+1)= 1;
> end
octave:4> A=A(2:end-1,2:end-1);
octave:5> disp(A);
-2 1 0
0 0 0 0
0 0 0
1 -2 1
0 0 0 0
0 0 0
0 1 -2
1 0 0 0
0 0 0
0 0 1
-2 1 0 0
0 0 0
0 0 0
1 -2 1 0
0 0 0
0 0 0
0 1 -2 1
0 0 0
0 0 0
0 0 1 -2
1 0 0
0 0 0
0 0 0 1
-2 1 0
0 0 0
0 0 0 0
1 -2 1
0 0 0
0 0 0 0
0 1 -2
octave:6> inv(A);
octave:7>
21:49 $ octave-cli
GNU Octave, version 4.0.0
Copyright (C) 2015 John W. Eaton and others.
This is free software; see the source code for copying conditions.
There is ABSOLUTELY NO WARRANTY; not even for MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. For details, type 'warranty'.
Octave was configured for "x86_64-unknown-linux-gnu".
Additional information about Octave is available at http://www.octave.org.
Please contribute if you find this software useful.
For more information, visit http://www.octave.org/get-involved.html
Read http://www.octave.org/bugs.html to learn how to submit bug reports.
For information about changes from previous versions, type 'news'.
octave:1> N=12;
octave:2> A=zeros(N);
octave:3> for n=2:N-1
> A(n,n-1)= 1;
> A(n,n )=-2;
> A(n,n+1)= 1;
> end
octave:4> A=A(2:end-1,2:end-1);
octave:5> disp(A);
-2 1 0
0 0 0 0
0 0 0
1 -2 1
0 0 0 0
0 0 0
0 1 -2
1 0 0 0
0 0 0
0 0 1
-2 1 0 0
0 0 0
0 0 0
1 -2 1 0
0 0 0
0 0 0
0 1 -2 1
0 0 0
0 0 0
0 0 1 -2
1 0 0
0 0 0
0 0 0 1
-2 1 0
0 0 0
0 0 0 0
1 -2 1
0 0 0
0 0 0 0
0 1 -2
octave:6> inv(A);
octave:7>
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?45651>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
| [Prev in Thread] | Current Thread | [Next in Thread] |