Skip to content

Commit 60d5dbc

Browse files
committed
Fix #52
Added in comment of xLARRV Note: VU is currently not used by this implmentation of DLARRV, VU is passed to DLARRV because it could be used compute gaps on the right end of the extremal eigenvalues. However, with not much initial accuracy in LAMBDA and VU, the formula can lead to an overestimation of the right gap and thus to inadequately early RQI 'convergence'. This is currently prevented this by forcing a small right gap. And so it turns out that VU is currently not used by this implementation of DLARRV.
1 parent 4445018 commit 60d5dbc

2 files changed

Lines changed: 16 additions & 4 deletions

File tree

SRC/dlarrv.f

Lines changed: 8 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -68,8 +68,14 @@
6868
*> \verbatim
6969
*> VU is DOUBLE PRECISION
7070
*> Upper bound of the interval that contains the desired
71-
*> eigenvalues. VL < VU. Needed to compute gaps on the left or right
72-
*> end of the extremal eigenvalues in the desired RANGE.
71+
*> eigenvalues. VL < VU.
72+
*> Note: VU is currently not used by this implmentation of DLARRV, VU is
73+
*> passed to DLARRV because it could be used compute gaps on the right end
74+
*> of the extremal eigenvalues. However, with not much initial accuracy in
75+
*> LAMBDA and VU, the formula can lead to an overestimation of the right gap
76+
*> and thus to inadequately early RQI 'convergence'. This is currently
77+
*> prevented this by forcing a small right gap. And so it turns out that VU
78+
*> is currently not used by this implementation of DLARRV.
7379
*> \endverbatim
7480
*>
7581
*> \param[in,out] D

SRC/slarrv.f

Lines changed: 8 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -68,8 +68,14 @@
6868
*> \verbatim
6969
*> VU is REAL
7070
*> Upper bound of the interval that contains the desired
71-
*> eigenvalues. VL < VU. Needed to compute gaps on the left or right
72-
*> end of the extremal eigenvalues in the desired RANGE.
71+
*> eigenvalues. VL < VU.
72+
*> Note: VU is currently not used by this implmentation of SLARRV, VU is
73+
*> passed to SLARRV because it could be used compute gaps on the right end
74+
*> of the extremal eigenvalues. However, with not much initial accuracy in
75+
*> LAMBDA and VU, the formula can lead to an overestimation of the right gap
76+
*> and thus to inadequately early RQI 'convergence'. This is currently
77+
*> prevented this by forcing a small right gap. And so it turns out that VU
78+
*> is currently not used by this implementation of SLARRV.
7379
*> \endverbatim
7480
*>
7581
*> \param[in,out] D

0 commit comments

Comments
 (0)