From aac571cf014f7b808d2c28e5b7afac72875ebcd2 Mon Sep 17 00:00:00 2001
From: Matthias Maier <tamiko@kyomu.43-1.org>
Date: Wed, 3 Jul 2013 09:01:46 +0000
Subject: [PATCH] Adapt the first three tests in tests/umfpack to test for the
 transpose as well.

git-svn-id: https://svn.dealii.org/trunk@29924 0785d39b-7218-0410-832d-ea1e28bc413d
---
 tests/umfpack/umfpack_01.cc          | 15 ++++++++++++---
 tests/umfpack/umfpack_01/cmp/generic | 24 ++++++++++++++++++++++++
 tests/umfpack/umfpack_02.cc          | 14 +++++++++++---
 tests/umfpack/umfpack_02/cmp/generic | 24 ++++++++++++++++++++++++
 tests/umfpack/umfpack_03.cc          | 15 ++++++++++++---
 tests/umfpack/umfpack_03/cmp/generic | 24 ++++++++++++++++++++++++
 6 files changed, 107 insertions(+), 9 deletions(-)

diff --git a/tests/umfpack/umfpack_01.cc b/tests/umfpack/umfpack_01.cc
index 7aa5b088b7..246ff8044e 100644
--- a/tests/umfpack/umfpack_01.cc
+++ b/tests/umfpack/umfpack_01.cc
@@ -12,6 +12,7 @@
 //----------------------------------------------------------------------
 
 // test the umfpack sparse direct solver on a mass matrix
+// test of the transpose as well
 
 #include "../tests.h"
 #include <iostream>
@@ -39,7 +40,7 @@
 
 
 template <int dim>
-void test ()
+void test (bool transpose = false)
 {
   deallog << dim << 'd' << std::endl;
   
@@ -84,9 +85,13 @@ void test ()
       for (unsigned int j=0; j<dof_handler.n_dofs(); ++j)
         solution(j) = j+j*(i+1)*(i+1);
 
-      B.vmult (b, solution);
+      if (transpose)
+        B.Tvmult (b, solution);
+      else
+        B.vmult (b, solution);
+
       x = b;
-      SparseDirectUMFPACK().solve (B, x);
+      SparseDirectUMFPACK().solve (B, x, transpose);
 
       x -= solution;
       deallog << "relative norm distance = "
@@ -111,4 +116,8 @@ int main ()
   test<1> ();
   test<2> ();
   test<3> ();
+
+  test<1> (/*transpose =*/ true);
+  test<2> (/*transpose =*/ true);
+  test<3> (/*transpose =*/ true);
 }
diff --git a/tests/umfpack/umfpack_01/cmp/generic b/tests/umfpack/umfpack_01/cmp/generic
index 8df3f884c1..5cbb5f8764 100644
--- a/tests/umfpack/umfpack_01/cmp/generic
+++ b/tests/umfpack/umfpack_01/cmp/generic
@@ -23,3 +23,27 @@ DEAL::relative norm distance = 0
 DEAL::absolute norms = 0 140414.
 DEAL::relative norm distance = 0
 DEAL::absolute norms = 0 280828.
+DEAL::1d
+DEAL::Number of dofs = 193
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 3084.00
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 7709.99
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 15420.0
+DEAL::2d
+DEAL::Number of dofs = 1889
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 94764.3
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 236911.
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 473822.
+DEAL::3d
+DEAL::Number of dofs = 1333
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 56165.6
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 140414.
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 280828.
diff --git a/tests/umfpack/umfpack_02.cc b/tests/umfpack/umfpack_02.cc
index 2554f6cbca..4eae46c954 100644
--- a/tests/umfpack/umfpack_02.cc
+++ b/tests/umfpack/umfpack_02.cc
@@ -41,7 +41,7 @@
 
 
 template <int dim>
-void test ()
+void test (bool transpose = false)
 {
   deallog << dim << 'd' << std::endl;
   
@@ -90,9 +90,13 @@ void test ()
       for (unsigned int j=0; j<dof_handler.n_dofs(); ++j)
         solution(j) = j+j*(i+1)*(i+1);
 
-      B.vmult (b, solution);
+      if (transpose)
+        B.Tvmult (b, solution);
+      else
+        B.vmult (b, solution);
+
       x = b;
-      Binv.solve (x);
+      Binv.solve (x, transpose);
 
       x -= solution;
       deallog << "relative norm distance = "
@@ -117,4 +121,8 @@ int main ()
   test<1> ();
   test<2> ();
   test<3> ();
+
+  test<1> (/*transpose =*/ true);
+  test<2> (/*transpose =*/ true);
+  test<3> (/*transpose =*/ true);
 }
diff --git a/tests/umfpack/umfpack_02/cmp/generic b/tests/umfpack/umfpack_02/cmp/generic
index 8df3f884c1..5cbb5f8764 100644
--- a/tests/umfpack/umfpack_02/cmp/generic
+++ b/tests/umfpack/umfpack_02/cmp/generic
@@ -23,3 +23,27 @@ DEAL::relative norm distance = 0
 DEAL::absolute norms = 0 140414.
 DEAL::relative norm distance = 0
 DEAL::absolute norms = 0 280828.
+DEAL::1d
+DEAL::Number of dofs = 193
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 3084.00
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 7709.99
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 15420.0
+DEAL::2d
+DEAL::Number of dofs = 1889
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 94764.3
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 236911.
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 473822.
+DEAL::3d
+DEAL::Number of dofs = 1333
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 56165.6
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 140414.
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 280828.
diff --git a/tests/umfpack/umfpack_03.cc b/tests/umfpack/umfpack_03.cc
index 4249df449c..753672cbbd 100644
--- a/tests/umfpack/umfpack_03.cc
+++ b/tests/umfpack/umfpack_03.cc
@@ -12,6 +12,7 @@
 //----------------------------------------------------------------------
 
 // test the umfpack sparse direct solver on a mass matrix that is slightly modified to make it nonsymmetric
+// test for the transpose as well
 
 #include "../tests.h"
 #include <iostream>
@@ -39,7 +40,7 @@
 
 
 template <int dim>
-void test ()
+void test (bool transpose = false)
 {
   deallog << dim << 'd' << std::endl;
 
@@ -101,9 +102,13 @@ void test ()
       for (unsigned int j=0; j<dof_handler.n_dofs(); ++j)
         solution(j) = j+j*(i+1)*(i+1);
 
-      B.vmult (b, solution);
+      if (transpose)
+        B.Tvmult (b, solution);
+      else
+        B.vmult (b, solution);
+
       x = b;
-      SparseDirectUMFPACK().solve (B, x);
+      SparseDirectUMFPACK().solve (B, x, transpose);
 
       x -= solution;
       deallog << "relative norm distance = "
@@ -128,4 +133,8 @@ int main ()
   test<1> ();
   test<2> ();
   test<3> ();
+
+  test<1> (/*transpose =*/ true);
+  test<2> (/*transpose =*/ true);
+  test<3> (/*transpose =*/ true);
 }
diff --git a/tests/umfpack/umfpack_03/cmp/generic b/tests/umfpack/umfpack_03/cmp/generic
index 8df3f884c1..5cbb5f8764 100644
--- a/tests/umfpack/umfpack_03/cmp/generic
+++ b/tests/umfpack/umfpack_03/cmp/generic
@@ -23,3 +23,27 @@ DEAL::relative norm distance = 0
 DEAL::absolute norms = 0 140414.
 DEAL::relative norm distance = 0
 DEAL::absolute norms = 0 280828.
+DEAL::1d
+DEAL::Number of dofs = 193
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 3084.00
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 7709.99
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 15420.0
+DEAL::2d
+DEAL::Number of dofs = 1889
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 94764.3
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 236911.
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 473822.
+DEAL::3d
+DEAL::Number of dofs = 1333
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 56165.6
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 140414.
+DEAL::relative norm distance = 0
+DEAL::absolute norms = 0 280828.
-- 
2.39.5