From: bangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Fri, 29 Feb 2008 22:04:35 +0000 (+0000)
Subject: Add a test where we use the ILU to compute actual inverses that we can check.
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=38e0d4816a43d5147cf981eefbd5788e79b2b1cd;p=dealii-svn.git

Add a test where we use the ILU to compute actual inverses that we can check.

git-svn-id: https://svn.dealii.org/trunk@15818 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/tests/lac/Makefile b/tests/lac/Makefile
index 7ae78a04de..b30389e76f 100644
--- a/tests/lac/Makefile
+++ b/tests/lac/Makefile
@@ -1,6 +1,6 @@
 ############################################################
 # $Id$
-# Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007 by the deal.II authors
+# Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008 by the deal.II authors
 ############################################################
 
 ############################################################
@@ -39,8 +39,7 @@ tests_x = vector-vector vector_memory \
 	solver \
 	eigen \
 	filtered_matrix \
-	sparse_ilu \
-	sparse_ilu_t \
+	sparse_ilu* \
 	lapack_fill \
 	block_minres \
 	identity_matrix* \
diff --git a/tests/lac/sparse_ilu_inverse.cc b/tests/lac/sparse_ilu_inverse.cc
new file mode 100644
index 0000000000..24f02cf90e
--- /dev/null
+++ b/tests/lac/sparse_ilu_inverse.cc
@@ -0,0 +1,87 @@
+//----------------------------  sparse_ilu_inverse.cc  ---------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2008 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  sparse_ilu_inverse.cc  ---------------------------
+
+
+// compute the inverse of a small matrix using the SparseILU with
+// infinite fill-in
+
+#include "../tests.h"
+#include <cmath>
+#include <fstream>
+#include <iomanip>
+#include <iomanip>
+#include <cstdlib>
+#include "testmatrix.h"
+#include <base/logstream.h>
+#include <lac/sparse_matrix.h>
+#include <lac/full_matrix.h>
+#include <lac/sparse_ilu.h>
+#include <lac/vector.h>
+
+
+int main()
+{
+  std::ofstream logfile("sparse_ilu_inverse/output");
+  deallog << std::fixed;
+  deallog << std::setprecision(3);
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+  
+
+  for (unsigned int N=1; N<5; ++N)
+    {
+      deallog << "N=" << N << std::endl;
+      
+      SparsityPattern structure(N, N, N);
+      for (unsigned int i=0; i<N; ++i)
+	for (unsigned int j=0; j<N; ++j)
+	  structure.add(i,j);
+      structure.compress();
+      SparseMatrix<double>  A(structure);
+      for (unsigned int i=0; i<N; ++i)
+	{
+	  A.set(i,i,2);
+	  if (i>=1)
+	    A.set(i,i-1, -1);
+	  if (i<N-1)
+	    A.set(i,i+1, -1);
+	}
+
+      SparseILU<double> ilu;
+      ilu.initialize (A, SparseILU<double>::AdditionalData());
+
+				       // now get an explicit
+				       // representation of the
+				       // inverse
+      FullMatrix<double> inverse (N,N);
+      Vector<double> tmp1(N), tmp2(N);
+      for (unsigned int i=0; i<N; ++i)
+	{
+	  tmp1 = 0;
+	  tmp1(i) = 1;
+	  ilu.vmult (tmp2, tmp1);
+	  for (unsigned int j=0; j<N; ++j)
+	    inverse(i,j) = tmp2(j);
+	}
+
+      deallog << "Matrix A:" << std::endl;
+      A.print_formatted (deallog.get_file_stream());
+  
+      deallog << "Matrix A^{-1}:" << std::endl;
+      inverse.print_formatted (deallog.get_file_stream());
+
+      deallog << std::endl;
+    }
+}
+
diff --git a/tests/lac/sparse_ilu_inverse/cmp/generic b/tests/lac/sparse_ilu_inverse/cmp/generic
new file mode 100644
index 0000000000..d156ab1dca
--- /dev/null
+++ b/tests/lac/sparse_ilu_inverse/cmp/generic
@@ -0,0 +1,28 @@
+
+DEAL::Size 4 Unknowns 9
+DEAL::Test 0
+DEAL::Residual with test vector 0:   left=0, right=0
+DEAL::Residual with test vector 1:   left=0, right=0
+DEAL::Residual with test vector 2:   left=0, right=0
+DEAL::Test 1
+DEAL::Residual with test vector 0:   left=0.308, right=0.350
+DEAL::Residual with test vector 1:   left=0.337, right=0.371
+DEAL::Residual with test vector 2:   left=0.290, right=0.327
+DEAL::Size 8 Unknowns 49
+DEAL::Test 0
+DEAL::Residual with test vector 0:   left=0, right=0
+DEAL::Residual with test vector 1:   left=0, right=0
+DEAL::Residual with test vector 2:   left=0, right=0
+DEAL::Test 1
+DEAL::Residual with test vector 0:   left=1.886, right=1.910
+DEAL::Residual with test vector 1:   left=2.123, right=2.138
+DEAL::Residual with test vector 2:   left=1.999, right=2.004
+DEAL::Size 16 Unknowns 225
+DEAL::Test 0
+DEAL::Residual with test vector 0:   left=0, right=0
+DEAL::Residual with test vector 1:   left=0, right=0
+DEAL::Residual with test vector 2:   left=0, right=0
+DEAL::Test 1
+DEAL::Residual with test vector 0:   left=5.944, right=5.953
+DEAL::Residual with test vector 1:   left=5.619, right=5.640
+DEAL::Residual with test vector 2:   left=6.079, right=6.094
diff --git a/tests/lac/sparse_ilu_inverse/cmp/mips-sgi-irix6.5+MIPSpro7.4 b/tests/lac/sparse_ilu_inverse/cmp/mips-sgi-irix6.5+MIPSpro7.4
new file mode 100644
index 0000000000..5a918f46fa
--- /dev/null
+++ b/tests/lac/sparse_ilu_inverse/cmp/mips-sgi-irix6.5+MIPSpro7.4
@@ -0,0 +1,28 @@
+
+DEAL::Size 4 Unknowns 9
+DEAL::Test 0
+DEAL::Residual with test vector 0:   left=0.000, right=0.000
+DEAL::Residual with test vector 1:   left=0.000, right=0.000
+DEAL::Residual with test vector 2:   left=0.000, right=0.000
+DEAL::Test 1
+DEAL::Residual with test vector 0:   left=0.425, right=0.447
+DEAL::Residual with test vector 1:   left=0.266, right=0.267
+DEAL::Residual with test vector 2:   left=0.297, right=0.377
+DEAL::Size 8 Unknowns 49
+DEAL::Test 0
+DEAL::Residual with test vector 0:   left=0.000, right=0.000
+DEAL::Residual with test vector 1:   left=0.000, right=0.000
+DEAL::Residual with test vector 2:   left=0.000, right=0.000
+DEAL::Test 1
+DEAL::Residual with test vector 0:   left=1.895, right=1.900
+DEAL::Residual with test vector 1:   left=1.705, right=1.711
+DEAL::Residual with test vector 2:   left=2.061, right=2.057
+DEAL::Size 16 Unknowns 225
+DEAL::Test 0
+DEAL::Residual with test vector 0:   left=0.000, right=0.000
+DEAL::Residual with test vector 1:   left=0.000, right=0.000
+DEAL::Residual with test vector 2:   left=0.000, right=0.000
+DEAL::Test 1
+DEAL::Residual with test vector 0:   left=5.900, right=5.911
+DEAL::Residual with test vector 1:   left=5.659, right=5.662
+DEAL::Residual with test vector 2:   left=5.914, right=5.921
diff --git a/tests/lac/sparse_ilu_inverse/cmp/powerpc-apple-darwin8.8.0+gcc4.0 b/tests/lac/sparse_ilu_inverse/cmp/powerpc-apple-darwin8.8.0+gcc4.0
new file mode 100644
index 0000000000..cb5e27caaf
--- /dev/null
+++ b/tests/lac/sparse_ilu_inverse/cmp/powerpc-apple-darwin8.8.0+gcc4.0
@@ -0,0 +1,28 @@
+
+DEAL::Size 4 Unknowns 9
+DEAL::Test 0
+DEAL::Residual with test vector 0:   left=0, right=0
+DEAL::Residual with test vector 1:   left=0, right=0
+DEAL::Residual with test vector 2:   left=0, right=0
+DEAL::Test 1
+DEAL::Residual with test vector 0:   left=0.295, right=0.333
+DEAL::Residual with test vector 1:   left=0.383, right=0.385
+DEAL::Residual with test vector 2:   left=0.299, right=0.285
+DEAL::Size 8 Unknowns 49
+DEAL::Test 0
+DEAL::Residual with test vector 0:   left=0, right=0
+DEAL::Residual with test vector 1:   left=0, right=0
+DEAL::Residual with test vector 2:   left=0, right=0
+DEAL::Test 1
+DEAL::Residual with test vector 0:   left=2.243, right=2.263
+DEAL::Residual with test vector 1:   left=1.707, right=1.707
+DEAL::Residual with test vector 2:   left=1.768, right=1.788
+DEAL::Size 16 Unknowns 225
+DEAL::Test 0
+DEAL::Residual with test vector 0:   left=0, right=0
+DEAL::Residual with test vector 1:   left=0, right=0
+DEAL::Residual with test vector 2:   left=0, right=0
+DEAL::Test 1
+DEAL::Residual with test vector 0:   left=5.672, right=5.683
+DEAL::Residual with test vector 1:   left=5.931, right=5.944
+DEAL::Residual with test vector 2:   left=5.988, right=5.996