Skip to content

Replace STEQR with PTEQR in Stokhos for issue #3542

James Willenbring requested to merge etphipp:stokhos_steqr_pteqr into develop

Created by: etphipp

The LAPACK function STEQR doesn't appear to work on IBM Power systems. In #3658 a solution was found by replacing STEQR with PTEQR, which seems to not have problems on this platform. This PR replaces usage of STEQR with PTEQR in Stokhos. However this isn't completely trivial because PTEQR is for symmetric, positive-definite systems where STEQR is for symmetric, indefinite systems (and therefore is a completely different algorithm, which is probably why it works when STEQR doesn't). So to use PTEQR I had to compute a shift to make the system positive definite. This caused a few tests to fail due to numerical differences with very tight tolerances, so I had to loosen some of those tolerances to get things to pass. After these changes, all of the Stokhos tests pass on white. For issue #3542 (closed).

Merge request reports