]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Add a couple of tests for TrilinosWrappers::SparseMatrix::n_nonzero_entries(). 5680/head
authorWolfgang Bangerth <bangerth@colostate.edu>
Thu, 28 Dec 2017 19:56:39 +0000 (12:56 -0700)
committerWolfgang Bangerth <bangerth@colostate.edu>
Thu, 28 Dec 2017 19:56:39 +0000 (12:56 -0700)
tests/trilinos/n_nonzeros_01.cc [new file with mode: 0644]
tests/trilinos/n_nonzeros_01.output [new file with mode: 0644]
tests/trilinos/n_nonzeros_02.cc [new file with mode: 0644]
tests/trilinos/n_nonzeros_02.output [new file with mode: 0644]

diff --git a/tests/trilinos/n_nonzeros_01.cc b/tests/trilinos/n_nonzeros_01.cc
new file mode 100644 (file)
index 0000000..93d2c80
--- /dev/null
@@ -0,0 +1,75 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2017 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 n_nonzero_elements() for an empty matrix
+
+#include "../tests.h"
+#include <deal.II/lac/trilinos_sparsity_pattern.h>
+#include <deal.II/lac/sparsity_pattern.h>
+#include <deal.II/lac/trilinos_sparse_matrix.h>
+
+
+void test ()
+{
+  // create an empty sparsity pattern
+  TrilinosWrappers::SparsityPattern sparsity(4,5,1);
+  sparsity.compress ();
+
+  // attach a sparse matrix to it
+  TrilinosWrappers::SparseMatrix A(sparsity);
+
+  // see how many nonzero elements it reports
+  deallog << A.n_nonzero_elements() << std::endl;
+}
+
+
+
+int main (int argc, char **argv)
+{
+  initlog();
+
+  Utilities::MPI::MPI_InitFinalize mpi_initialization (argc, argv, testing_max_num_threads());
+
+  try
+    {
+      test ();
+    }
+  catch (std::exception &exc)
+    {
+      deallog << std::endl << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      deallog << "Exception on processing: " << std::endl
+              << exc.what() << std::endl
+              << "Aborting!" << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+
+      return 1;
+    }
+  catch (...)
+    {
+      deallog << std::endl << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      deallog << "Unknown exception!" << std::endl
+              << "Aborting!" << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      return 1;
+    };
+}
diff --git a/tests/trilinos/n_nonzeros_01.output b/tests/trilinos/n_nonzeros_01.output
new file mode 100644 (file)
index 0000000..199228a
--- /dev/null
@@ -0,0 +1,2 @@
+
+DEAL::0
diff --git a/tests/trilinos/n_nonzeros_02.cc b/tests/trilinos/n_nonzeros_02.cc
new file mode 100644 (file)
index 0000000..4a5c780
--- /dev/null
@@ -0,0 +1,96 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2017 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 n_nonzero_elements() for an empty matrix
+
+#include "../tests.h"
+#include <deal.II/lac/trilinos_sparsity_pattern.h>
+#include <deal.II/lac/sparsity_pattern.h>
+#include <deal.II/lac/trilinos_sparse_matrix.h>
+#include <deal.II/dofs/dof_tools.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/grid/grid_generator.h>
+
+
+void test ()
+{
+  Triangulation<2> tria;
+  GridGenerator::hyper_cube (tria);
+
+  FE_Q<2> fe(1);
+  DoFHandler<2> dof_handler (tria);
+  dof_handler.distribute_dofs(fe);
+
+  Table<2,DoFTools::Coupling> coupling (1,1);
+  coupling.fill (DoFTools::none);
+
+  // create an empty sparsity pattern
+  TrilinosWrappers::SparsityPattern sparsity;
+  sparsity.reinit (dof_handler.n_dofs(), dof_handler.n_dofs());
+  DoFTools::make_sparsity_pattern (dof_handler,
+                                   coupling, sparsity,
+                                   ConstraintMatrix(), false,
+                                   Utilities::MPI::
+                                   this_mpi_process(MPI_COMM_WORLD));
+  sparsity.compress ();
+
+  // attach a sparse matrix to it
+  TrilinosWrappers::SparseMatrix A;
+  A.reinit(sparsity);
+
+  // see how many nonzero elements it reports
+  deallog << A.n_nonzero_elements() << std::endl;
+}
+
+
+
+int main (int argc, char **argv)
+{
+  initlog();
+
+  Utilities::MPI::MPI_InitFinalize mpi_initialization (argc, argv, testing_max_num_threads());
+
+  try
+    {
+      test ();
+    }
+  catch (std::exception &exc)
+    {
+      deallog << std::endl << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      deallog << "Exception on processing: " << std::endl
+              << exc.what() << std::endl
+              << "Aborting!" << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+
+      return 1;
+    }
+  catch (...)
+    {
+      deallog << std::endl << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      deallog << "Unknown exception!" << std::endl
+              << "Aborting!" << std::endl
+              << "----------------------------------------------------"
+              << std::endl;
+      return 1;
+    };
+}
diff --git a/tests/trilinos/n_nonzeros_02.output b/tests/trilinos/n_nonzeros_02.output
new file mode 100644 (file)
index 0000000..199228a
--- /dev/null
@@ -0,0 +1,2 @@
+
+DEAL::0

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.