Add a missing template function to TridiagonalMatrix.

diff --git a/source/lac/tridiagonal_matrix.cc b/source/lac/tridiagonal_matrix.cc
index 16937b54457d2dd0acbb95e88f8a92d832096b8c..cf9b9fe3b0f9e9d6f6584cf370d4c3d691270b68 100644 (file)
--- a/source/lac/tridiagonal_matrix.cc
+++ b/source/lac/tridiagonal_matrix.cc
@@ -35,6 +35,7 @@ TridiagonalMatrix<number>::TridiagonalMatrix(
 {}
 
 
+
 template <typename number>
 void
 TridiagonalMatrix<number>::reinit(
@@ -49,6 +50,7 @@ TridiagonalMatrix<number>::reinit(
 }
 
 
+
 template <typename number>
 bool
 TridiagonalMatrix<number>::all_zero() const
@@ -73,6 +75,7 @@ TridiagonalMatrix<number>::all_zero() const
 }
 
 
+
 template <typename number>
 void
 TridiagonalMatrix<number>::vmult (
@@ -137,6 +140,7 @@ TridiagonalMatrix<number>::vmult_add (
 }
 
 
+
 template <typename number>
 void
 TridiagonalMatrix<number>::Tvmult (
@@ -181,6 +185,7 @@ TridiagonalMatrix<number>::Tvmult (
 }
 
 
+
 template <typename number>
 void
 TridiagonalMatrix<number>::Tvmult_add (
@@ -191,6 +196,7 @@ TridiagonalMatrix<number>::Tvmult_add (
 }
 
 
+
 template <typename number>
 number
 TridiagonalMatrix<number>::matrix_scalar_product(
@@ -218,6 +224,7 @@ TridiagonalMatrix<number>::matrix_scalar_product(
 }
 
 
+
 template <typename number>
 number
 TridiagonalMatrix<number>::matrix_norm_square(
@@ -227,9 +234,10 @@ TridiagonalMatrix<number>::matrix_norm_square(
 }
 
 
-template <>
+
+template <typename number>
 void
-TridiagonalMatrix<double>::compute_eigenvalues()
+TridiagonalMatrix<number>::compute_eigenvalues()
 {
 #ifdef DEAL_II_WITH_LAPACK
   Assert(state == matrix, ExcState(state));
@@ -247,6 +255,7 @@ TridiagonalMatrix<double>::compute_eigenvalues()
 }
 
 
+
 template <typename number>
 number
 TridiagonalMatrix<number>::eigenvalue(const size_type i) const

diff --git a/tests/lapack/tridiagonal_ev1.cc b/tests/lapack/tridiagonal_ev1.cc
index ce27790bb9629f4aaf9599900d475fb509394e41..c19ee60f5593440219a548d13da26e069d0c406d 100644 (file)
--- a/tests/lapack/tridiagonal_ev1.cc
+++ b/tests/lapack/tridiagonal_ev1.cc
@@ -52,4 +52,9 @@ int main()
   test_laplacian<double>(20);
   test_laplacian<double>(40);
   test_laplacian<double>(80);
+
+  test_laplacian<float>(10);
+  test_laplacian<float>(20);
+  test_laplacian<float>(40);
+  test_laplacian<float>(80);
 }

diff --git a/tests/lapack/tridiagonal_ev1.output b/tests/lapack/tridiagonal_ev1.output
index eca9a0403264ec1f1acc7dc68b1af9142e770c09..224f45af3a002e650735976ec22db75d6dc0e2ee 100644 (file)
--- a/tests/lapack/tridiagonal_ev1.output
+++ b/tests/lapack/tridiagonal_ev1.output
@@ -3,3 +3,7 @@ DEAL::  9.78870 38.1966 82.4429 138.197 200.000  cond 39.8635
 DEAL:: 9.84933 39.1548 87.1948 152.786 234.315  cond 161.448
 DEAL:: 9.86453 39.3973 88.4163 156.619 243.585  cond 647.789
 DEAL:: 9.86834 39.4581 88.7238 157.589 245.948  cond 2593.16
+DEAL:: 9.78869 38.1966 82.4429 138.197 200.000  cond 39.8635
+DEAL:: 9.84912 39.1547 87.1948 152.786 234.315  cond 161.451
+DEAL:: 9.86457 39.3974 88.4164 156.619 243.586  cond 647.787
+DEAL:: 9.86855 39.4586 88.7249 157.592 245.949  cond 2593.10