Follow up with ad5bc21cb50535d66d628a309d60128db96c8851

Contribution from Zlatko Drmac

1) LWORK query added;

2) few modifications in pure one sided Jacobi (XGESVJ) to remove possible error
in the really extreme cases (sigma_max close to overflow and sigma_min close to
underflow)  - note that XGESVJ is deigned to compute the singular values in the
full range; I used it (double complex) to compute SVD of certain factored
Hankel matrices with the condition number 1.0e616;

3) in the preconditioned Jacobi SVD (XGEJSV), the code I sent before to Julie
had one experimental modification that I had forgotten to remove before sending
- now this is done (the idea was to extend the computational range, but that
brings to much too risky dependence on how other lapack routines behave under
those extreme conditions).

0 files changed, 0 insertions, 0 deletions