From: maier <maier@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Wed, 3 Jul 2013 09:01:46 +0000 (+0000)
Subject: Adapt the first three tests in tests/umfpack to test for the transpose as well.
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=0e61d878b88c45876a3794fceaaa6e97edd99aee;p=dealii-svn.git

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
---

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.