From: David Wells <wellsd2@rpi.edu>
Date: Sun, 10 Sep 2017 01:17:15 +0000 (-0400)
Subject: Add a missing template function to TridiagonalMatrix.
X-Git-Tag: v9.0.0-rc1~854^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=95df5a1e10cac0f5dc829fd8fca57d063bf2bfa5;p=dealii.git

Add a missing template function to TridiagonalMatrix.
---

diff --git a/source/lac/tridiagonal_matrix.cc b/source/lac/tridiagonal_matrix.cc
index 16937b5445..cf9b9fe3b0 100644
--- 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 ce27790bb9..c19ee60f55 100644
--- 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 eca9a04032..224f45af3a 100644
--- 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