From c5a2d90ed460d16b4332a1148b059be05a660de3 Mon Sep 17 00:00:00 2001 From: "Heidi K. Thornquist" Date: Tue, 16 Oct 2018 14:42:35 -0600 Subject: [PATCH] Replace PTEQR with STEQR test for #2410 Replace working tridiagonal eigensolver with the one that causes the seg fault on ride/white/waterman. --- packages/teuchos/numerics/test/LAPACK/cxx_main.cpp | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/packages/teuchos/numerics/test/LAPACK/cxx_main.cpp b/packages/teuchos/numerics/test/LAPACK/cxx_main.cpp index 80ca2e8bbbdd..7440693879a1 100644 --- a/packages/teuchos/numerics/test/LAPACK/cxx_main.cpp +++ b/packages/teuchos/numerics/test/LAPACK/cxx_main.cpp @@ -212,7 +212,7 @@ int main(int argc, char* argv[]) #endif - if (verbose) std::cout << "PTEQR test ... "; + if (verbose) std::cout << "STEQR test ... "; const int n_steqr = 10; std::vector diagonal(n_steqr); @@ -227,20 +227,20 @@ int main(int argc, char* argv[]) std::vector scalar_dummy(1,0.0); std::vector mag_dummy(4*n_steqr,0.0); - L.PTEQR (char_N, n_steqr, &diagonal[0], &subdiagonal[0], + L.STEQR (char_N, n_steqr, &diagonal[0], &subdiagonal[0], &scalar_dummy[0], n_steqr, &mag_dummy[0], &info); if (info != 0) { - if (verbose) std::cout << "PTEQR: compute symmetric tridiagonal eigenvalues: " - << "LAPACK's _PTEQR failed with info = " + if (verbose) std::cout << "STEQR: compute symmetric tridiagonal eigenvalues: " + << "LAPACK's _STEQR failed with info = " << info; numberFailedTests++; } - MagnitudeType lambda_min = diagonal[n_steqr-1]; - MagnitudeType lambda_max = diagonal[0]; + MagnitudeType lambda_min = diagonal[0]; + MagnitudeType lambda_max = diagonal[n_steqr-1]; MagnitudeType exp_lambda_min = STM::one(); MagnitudeType exp_lambda_max = STM::one()*n_steqr;