Skip to content

Commit

Permalink
Replace PTEQR with STEQR test for trilinos#2410
Browse files Browse the repository at this point in the history
Replace working tridiagonal eigensolver with the one that causes the
seg fault on ride/white/waterman.
  • Loading branch information
hkthorn authored and tjfulle committed Dec 6, 2018
1 parent a7010dc commit c5a2d90
Showing 1 changed file with 6 additions and 6 deletions.
12 changes: 6 additions & 6 deletions packages/teuchos/numerics/test/LAPACK/cxx_main.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -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<MagnitudeType> diagonal(n_steqr);
Expand All @@ -227,20 +227,20 @@ int main(int argc, char* argv[])
std::vector<double> scalar_dummy(1,0.0);
std::vector<MagnitudeType> 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;

Expand Down

0 comments on commit c5a2d90

Please sign in to comment.