Test PETSc FullMatrix reinit function.

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

diff --git a/tests/petsc/full_matrix_00.cc b/tests/petsc/full_matrix_00.cc
new file mode 100644 (file)
index 0000000..cbdde4b
--- /dev/null
+++ b/tests/petsc/full_matrix_00.cc
@@ -0,0 +1,122 @@
+// ---------------------------------------------------------------------
+// $Id: full_matrix_01.cc 31602 2013-11-09 15:38:36Z young $
+//
+// Copyright (C) 2014 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// The deal.II library is free software; you can use it, redistribute
+// it, and/or modify it under the terms of the GNU Lesser General
+// Public License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+// The full text of the license can be found in the file LICENSE at
+// the top level of the deal.II distribution.
+//
+// ---------------------------------------------------------------------
+
+
+// check initialiser (aka do_reinit). .
+//
+// PETScWrappers::FullMatrix::reinit ()
+
+#include "../tests.h"
+#include <deal.II/lac/petsc_full_matrix.h>
+#include <fstream>
+#include <iostream>
+
+// Simply take an already intialised full matrix, fill in some of the
+// elements, reinitialise it to a different size and fill in some
+// elements again. To be sure this generates the correct results the
+// test is verbose (and uses tiny matrices).
+void test (PETScWrappers::FullMatrix &m)
+{
+  // asign some matrix elements
+  for (unsigned int i=0; i<m.m(); ++i)
+    for (unsigned int j=0; j<m.n(); ++j)
+      m.set (i, j, i+2*j);
+
+  m.compress (VectorOperation::add);
+
+  // things we know
+  Assert (m.m() == 3, ExcInternalError());
+  Assert (m.n() == 3, ExcInternalError());
+  
+  // Generate some output
+  deallog << "initial matrix: " << std::endl;
+  for (unsigned int i=0; i<m.m(); ++i)
+    {
+      for (unsigned int j=0; j<m.n(); ++j)
+       deallog << m(i,j) << " ";
+      deallog << std::endl;
+    }
+  deallog << std::endl;
+
+  // test reinit, this time a different size matrix
+  m.reinit (5, 5);
+
+  // set some entries
+  for (unsigned int i=0; i<m.m(); ++i)
+    for (unsigned int j=0; j<m.n(); ++j)
+      m.set (i, j, j+2*i);
+  m.compress (VectorOperation::add);
+  
+  // things we know
+  Assert (m.m() == 5, ExcInternalError());
+  Assert (m.n() == 5, ExcInternalError());
+
+  // Generate some output
+  deallog << "after reinit: " << std::endl;
+  for (unsigned int i=0; i<m.m(); ++i)
+    {
+      for (unsigned int j=0; j<m.n(); ++j)
+       deallog << m(i,j) << " ";
+      deallog << std::endl;
+    }
+  deallog << std::endl;
+  
+  // done
+  deallog << "OK" << std::endl;
+}
+
+int main (int argc,char **argv)
+{
+  std::ofstream logfile("output");
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+
+  try
+    {
+      Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
+      {
+       // Standard intialiser...
+        PETScWrappers::FullMatrix m (3,3);
+        test (m);
+      }
+
+    }
+  catch (std::exception &exc)
+    {
+      std::cerr << std::endl << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+      std::cerr << "Exception on processing: " << std::endl
+                << exc.what() << std::endl
+                << "Aborting!" << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+
+      return 1;
+    }
+  catch (...)
+    {
+      std::cerr << std::endl << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+      std::cerr << "Unknown exception!" << std::endl
+                << "Aborting!" << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+      return 1;
+    };
+}

diff --git a/tests/petsc/full_matrix_00.output b/tests/petsc/full_matrix_00.output
new file mode 100644 (file)
index 0000000..c7aa78d
--- /dev/null
+++ b/tests/petsc/full_matrix_00.output
@@ -0,0 +1,14 @@
+
+DEAL::initial matrix: 
+DEAL::0 2.00000 4.00000 
+DEAL::1.00000 3.00000 5.00000 
+DEAL::2.00000 4.00000 6.00000 
+DEAL::
+DEAL::after reinit: 
+DEAL::0 1.00000 2.00000 3.00000 4.00000 
+DEAL::2.00000 3.00000 4.00000 5.00000 6.00000 
+DEAL::4.00000 5.00000 6.00000 7.00000 8.00000 
+DEAL::6.00000 7.00000 8.00000 9.00000 10.0000 
+DEAL::8.00000 9.00000 10.0000 11.0000 12.0000 
+DEAL::
+DEAL::OK