[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #54399] function LuAminusSigmaB in eigs-base.c
|
From: |
Dmitri A. Sergatskov |
|
Subject: |
[Octave-bug-tracker] [bug #54399] function LuAminusSigmaB in eigs-base.cc can be made faster for special case of sigma == 0 |
|
Date: |
Wed, 1 Aug 2018 11:23:21 -0400 (EDT) |
|
User-agent: |
Mozilla/5.0 (X11; Fedora; Linux x86_64; rv:61.0) Gecko/20100101 Firefox/61.0 |
Follow-up Comment #23, bug #54399 (project octave):
octave:1> m = ones(11, 11) * 0.15 / 11;
octave:2> m(1, 2:11) += 0.85;
octave:3> m(2:11, 1) = (1 - m(1, 1)) / 10;
octave:4> [v, lambda] = eigs(m, 1)
v =
-0.94054
-0.10742
-0.10742
-0.10742
-0.10742
-0.10742
-0.10742
-0.10742
-0.10742
-0.10742
-0.10742
lambda = 1.00000
So I do not have that bug.
Dmitri.
--
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?54399>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/