From 600656b8e1deb89e8f89d5902c5bea1f4572396b Mon Sep 17 00:00:00 2001 From: guido Date: Wed, 25 Jun 2003 12:25:22 +0000 Subject: [PATCH] test conjugate_add git-svn-id: https://svn.dealii.org/trunk@7811 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/lac/sparse_matrices.cc | 42 +++++++++++++++++++ .../lac/sparse_matrices.output | 15 +++++++ 2 files changed, 57 insertions(+) diff --git a/tests/lac/sparse_matrices.cc b/tests/lac/sparse_matrices.cc index 981fe3ff14..8f8ffc379e 100644 --- a/tests/lac/sparse_matrices.cc +++ b/tests/lac/sparse_matrices.cc @@ -168,6 +168,47 @@ void check_ez_iterator() } +void check_conjugate(std::ostream& out) +{ + SparseMatrixEZ B(3,2); + SparseMatrixEZ A1(2,2); + SparseMatrixEZ A2(3,3); + SparseMatrixEZ C1(3,3); + SparseMatrixEZ C2(2,2); + + B.set(0, 0, 0.); + B.set(0, 1, 1.); + B.set(1, 0, 2.); + B.set(1, 1, 3.); + B.set(2, 0, 4.); + B.set(2, 1, 5.); + + A1.set(0, 0, 1.); + A1.set(0, 1, 2.); + A1.set(1, 0, 3.); + A1.set(1, 1, 4.); + + A2.set(0, 0, 1.); + A2.set(0, 1, 2.); + A2.set(0, 2, 3.); + A2.set(1, 0, 4.); + A2.set(1, 1, 5.); + A2.set(1, 2, 6.); + A2.set(2, 0, 7.); + A2.set(2, 1, 8.); + A2.set(2, 2, 9.); + + C1.conjugate_add(A1,B); + C2.conjugate_add(A2,B,true); + + out << "First conjugate" << std::endl; + C1.print(out); + + out << "Second conjugate" << std::endl; + C2.print(out); +} + + int main() { std::ofstream logfile("sparse_matrices.output"); @@ -190,6 +231,7 @@ int main() #endif check_ez_iterator(); + check_conjugate(logfile); FDMatrix testproblem (size, size); unsigned int dim = (size-1)*(size-1); diff --git a/tests/results/i686-pc-linux-gnu+gcc3.2/lac/sparse_matrices.output b/tests/results/i686-pc-linux-gnu+gcc3.2/lac/sparse_matrices.output index d0fcd083aa..0e2be7d4dd 100644 --- a/tests/results/i686-pc-linux-gnu+gcc3.2/lac/sparse_matrices.output +++ b/tests/results/i686-pc-linux-gnu+gcc3.2/lac/sparse_matrices.output @@ -47,6 +47,21 @@ DEAL:: 2 0 0 3.0 DEAL:: 4 0 3 17. DEAL::Repeat row 2 DEAL:: 2 0 0 3.0 +First conjugate +0 0 4.0 +0 1 18. +0 2 32. +1 0 16. +1 1 70. +1 2 1.2e+02 +2 0 28. +2 1 1.2e+02 +2 2 2.2e+02 +Second conjugate +0 0 2.8e+02 +0 1 4.0e+02 +1 0 3.8e+02 +1 1 5.5e+02 DEAL::Structure DEAL::Assemble DEAL:5-SparseMatrix:Richardson::Starting value 4.0 -- 2.39.5