forked from Reference-LAPACK/lapack-www
-
Notifications
You must be signed in to change notification settings - Fork 0
/
errata_from_342_to_350.txt
106 lines (93 loc) · 4.69 KB
/
errata_from_342_to_350.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
Bugs corrected from 3.4.2 to 3.5.0
==================================
== [green]*bug112* :: DSYEVR does not seem to be thread-safe
* *CORRECTED*, rev 1436
* Present in LAPACK 3.4.2
* Bug report by Daniel Strobusch on LAPACK Mailing - August 19th 2013
* Links: http://stackoverflow.com/q/18216314/1458343
* Fix: LAPACK needs to be compiled with option "-frecursive" so that local
arrays are allocated on the stack. * Corrected by Julien on Oct 4 2013
== [green]*bug110* :: possible infinite loop in ZLARTG when input vector has NaN
* *CORRECTED*, rev 1437
* Present in LAPACK 3.4.2
* Bug report by Mark Hoemmen and Jack Poulson
* See http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=5&t=3267&p=9067
* Mark Hoemmen writes: Jack Poulson found that ZLARTG has an infinite loop related
* to NaN .EQ. 0. It's the line below "ELSE IF( SCALE.LE.SAFMN2 )" between 10 and
* 20. When G is NaN, then the loop goes forever.
* Corrected by Rodney on Oct 5 2013
== [green]*bug0109* :: [xDRVRFP] Test was failing
* *CORRECTED*, rev 1425
* Present in LAPACK 3.4.2
* bug report by Hong Xu (INTEL) on July 21st
* Patch sent by Hong Xu (INTEL)
* In the input parameter checking, LDZ was not checked correctly
* LDZ was checked against M instead of N essentially
* Corrected by Julien on July 22 2013
== [green]*bug0108* :: [xDRVRFP] Test was failing
* *CORRECTED*, rev 1416
* Present in LAPACK 3.4.2
* bug report by Elena Ivanova (ORACLE)
* Patch sent by Elena Ivanova (ORACLE)
* add MAX() to ensure index is at least 1
* Corrected by Rodney on June 27 2013
== [green]*bug0107* :: [CSD routines] fixed array index problems
* *CORRECTED*, rev 1414
* Present in LAPACK 3.4.2
* add MAX() to ensure index is at least 1
* this fixes array bounds violations in the test suite when array bounds checking is enabled
* Corrected by Rodney on June 04 2013
== [green]*bug0106* :: [xGEBAL] very large backward error for certain Hessenberg matrices (Watkins)
* *CORRECTED*, rev 1413
* Present in LAPACK 3.4.2
* bug report on the forum topic 4221
* Now using 2-norm to compute vector norms of row and column for balancing algorithm.
* This seems to fix problems with balancing causing very large backward error for certain
* Hessenberg matrices, including the often cited example of Watkins.
* Corrected by Bradley Lowery, Rodney James, Julien Langou on May 26 2013
== [green]*bug0105* :: [LAPACKE] lapacke.h problem
* *CORRECTED*, rev 1388
* Present in LAPACKE 3.4.2
* bug report on the forum topic 4221
* Patch sent by bmoloney
* A few C++ library includes are inside of the 'extern C block which causes many template with C
* linkage errors when the header is included in a C++ file.
* Corrected by Julie on April 03 2013
== [green]*bug0104* :: dgeqp3 fall over when asked to do a factorization of the last two columns.
* *CORRECTED*, rev 1385
* Present in LAPACK 3.4.2
* bug report by Tyrone Rees (STFC Rutherford Appleton Laboratory) on March 22 2013
* Patch sent by Tyrone Rees
* Corrected by Rodney on March 22 2013
== [green]*bug0103* :: [dlasd4/slasd4] fixes bug where TAU2 can be uninitialized.
* *CORRECTED*, rev 1384
* Present in LAPACK 3.4.2
* bug report by Sergey Kuznetsov (INTEL) on March 15 2013
* Patch sent by Sergey Kuznetsov
* Corrected by Rodney on March 22 2013
== [green]*bug0102* :: convergence problem with some particular pencils in xGGEV
* *CORRECTED*, rev 1351
* Present in LAPACK 3.4.2
* bug report by Mathworks on October 19 2012
* Patch sent by Mathworks
* Corrected by Rodney on Thu Oct 19
== [green]*bug0101* :: [S,D, C, Z]hsein: NaN in H can cause an infinite loop
* *CORRECTED*, rev 1350
* Present in LAPACK 3.4.2
* bug report by Alexander Kobotov (INTEL) on October 3 2012
* See link:http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=3928[forum topic 3928]
* added check for NaN after the norm of each block of H is computed
* if NaN is detected, returns with INFO=-6 (H is the 6th parameter)
* Corrected by Rodney on Thu Oct 11
== [green]*bug0057* :: Non-eigenvectors from DGEEV, problem in balancing
* *CORRECTED*
* Same bug than bug0106
* bug report sent by Numpy team on April 22nd, 2010
* see link:https://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=2&t=1818[forum topic 1818]
* see follow-up on lapackers on April 22nd from Jim "bug in DGEHRD"
"The balancing is choosing to scale some rows/columns by factors like 1e18,
which rather amplifies tiny rounding errors in the (presumably correctly)
computed eigenvectors of the balanced matrix, when balancing is undone at the end.
The expert driver interface dgeevx would let the user disable balancing, until
we figure out whether this is a bug or a feature, I'm not sure which."
* Corrected by Bradley Lowery, Rodney James, Julien Langou on May 26 2013