|
| From: | Dmitri A. Sergatskov |
| Subject: | [Octave-bug-tracker] [bug #56037] qp exits with an indexing error for a specific input with OpenBLAS |
| Date: | Wed, 3 Apr 2019 15:22:33 -0400 (EDT) |
| User-agent: | Mozilla/5.0 (X11; Fedora; Linux x86_64; rv:66.0) Gecko/20100101 Firefox/66.0 |
Follow-up Comment #15, bug #56037 (project octave):
I do not get error code 3 any more, but the answers are different depending on
blas version. Perhaps the problem is ill-posed...
(default, with openblas)
octave -q -f
octave:1> try1
x =
0.50000
0.00000
0.50000
0.00000
0.00000
-0.00000
0.00000
0.00000
info = 1.4050e-31
obj =
scalar structure containing the fields:
solveiter = 8
info = 0
lambda =
-1.7790e-47
0.0000e+00
0.0000e+00
4.9304e-31
0.0000e+00
0.0000e+00
0.0000e+00
4.9304e-31
0.0000e+00
4.9304e-31
0.0000e+00
4.9304e-31
0.0000e+00
4.9304e-31
0.0000e+00
4.9304e-31
0.0000e+00
With atlas:
LD_PRELOAD=/usr/lib64/atlas/libtatlas.so octave -q -f
octave:1> try1
x =
0.50000
-0.00000
0.00000
0.00000
0.35884
-0.00000
0.00000
0.14116
info = 3.0830e-33
obj =
scalar structure containing the fields:
solveiter = 14
info = 0
lambda =
4.5213e-17
0.0000e+00
0.0000e+00
4.5213e-17
0.0000e+00
4.5213e-17
0.0000e+00
0.0000e+00
0.0000e+00
0.0000e+00
0.0000e+00
4.5213e-17
0.0000e+00
4.5213e-17
0.0000e+00
0.0000e+00
0.0000e+00
with reference blas:
LD_PRELOAD=/usr/lib64/libblas.so.3 octave -q -f
octave:1> try1
x =
0.12500
0.23988
0.00000
0.00000
0.00000
0.00000
0.00000
0.63512
info = 0.56250
obj =
scalar structure containing the fields:
solveiter = 6
info = 0
lambda =
0.50000
0.00000
0.00000
0.00000
0.00000
1.00000
0.00000
1.00000
0.00000
1.00000
0.00000
1.00000
0.00000
1.00000
0.00000
0.00000
0.00000
Dmitri.
_______________________________________________________
Reply to this item at:
<https://savannah.gnu.org/bugs/?56037>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/
| [Prev in Thread] | Current Thread | [Next in Thread] |