Bugfix: Fix lapack stev call in TridiagonalMatrix<number>::compute_eigenvalues()

diff --git a/source/lac/tridiagonal_matrix.cc b/source/lac/tridiagonal_matrix.cc
index 986b7c7af17c2f4fd3dfcc9deea564b7f2d538f1..9c9ba64ef37dc51f4593eb6cde47fc9d53ff542c 100644 (file)
--- a/source/lac/tridiagonal_matrix.cc
+++ b/source/lac/tridiagonal_matrix.cc
@@ -241,8 +241,14 @@ TridiagonalMatrix<number>::compute_eigenvalues()
 
   const types::blas_int nn = n();
   types::blas_int       info;
-  //   stev(&N, &nn, diagonal.data(), right.data(), nullptr, &one, nullptr,
-  //   &info);
+  stev(&N,
+       &nn,
+       diagonal.data(),
+       right.data(),
+       static_cast<number *>(nullptr),
+       &one,
+       static_cast<number *>(nullptr),
+       &info);
   Assert(info == 0, ExcInternalError());
 
   state = LAPACKSupport::eigenvalues;