From 95df5a1e10cac0f5dc829fd8fca57d063bf2bfa5 Mon Sep 17 00:00:00 2001 From: David Wells Date: Sat, 9 Sep 2017 21:17:15 -0400 Subject: [PATCH] Add a missing template function to TridiagonalMatrix. --- source/lac/tridiagonal_matrix.cc | 13 +++++++++++-- tests/lapack/tridiagonal_ev1.cc | 5 +++++ tests/lapack/tridiagonal_ev1.output | 4 ++++ 3 files changed, 20 insertions(+), 2 deletions(-) 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::TridiagonalMatrix( {} + template void TridiagonalMatrix::reinit( @@ -49,6 +50,7 @@ TridiagonalMatrix::reinit( } + template bool TridiagonalMatrix::all_zero() const @@ -73,6 +75,7 @@ TridiagonalMatrix::all_zero() const } + template void TridiagonalMatrix::vmult ( @@ -137,6 +140,7 @@ TridiagonalMatrix::vmult_add ( } + template void TridiagonalMatrix::Tvmult ( @@ -181,6 +185,7 @@ TridiagonalMatrix::Tvmult ( } + template void TridiagonalMatrix::Tvmult_add ( @@ -191,6 +196,7 @@ TridiagonalMatrix::Tvmult_add ( } + template number TridiagonalMatrix::matrix_scalar_product( @@ -218,6 +224,7 @@ TridiagonalMatrix::matrix_scalar_product( } + template number TridiagonalMatrix::matrix_norm_square( @@ -227,9 +234,10 @@ TridiagonalMatrix::matrix_norm_square( } -template <> + +template void -TridiagonalMatrix::compute_eigenvalues() +TridiagonalMatrix::compute_eigenvalues() { #ifdef DEAL_II_WITH_LAPACK Assert(state == matrix, ExcState(state)); @@ -247,6 +255,7 @@ TridiagonalMatrix::compute_eigenvalues() } + template number TridiagonalMatrix::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(20); test_laplacian(40); test_laplacian(80); + + test_laplacian(10); + test_laplacian(20); + test_laplacian(40); + test_laplacian(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 -- 2.39.5