]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Add a test that currently fails.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 10 May 2011 16:27:44 +0000 (16:27 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 10 May 2011 16:27:44 +0000 (16:27 +0000)
git-svn-id: https://svn.dealii.org/trunk@23688 0785d39b-7218-0410-832d-ea1e28bc413d

tests/mpi/trilinos_sparse_matrix_mmult_01.cc [new file with mode: 0644]
tests/mpi/trilinos_sparse_matrix_mmult_01/ncpu_1/cmp/generic [new file with mode: 0644]
tests/mpi/trilinos_sparse_matrix_mmult_01/ncpu_2/cmp/generic [new file with mode: 0644]

diff --git a/tests/mpi/trilinos_sparse_matrix_mmult_01.cc b/tests/mpi/trilinos_sparse_matrix_mmult_01.cc
new file mode 100644 (file)
index 0000000..354eda9
--- /dev/null
@@ -0,0 +1,167 @@
+//----------------------------  trilinos_vector_equality_4.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 2004, 2005, 2008, 2010, 2011 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.
+//
+//----------------------------  trilinos_vector_equality_4.cc  ---------------------------
+
+
+// TrilinosWrappers::SparseMatrix::print got column indices wrong
+
+#include "../tests.h"
+#include <base/utilities.h>
+#include <base/index_set.h>
+#include <lac/trilinos_sparse_matrix.h>
+#include <lac/trilinos_sparsity_pattern.h>
+#include <fstream>
+#include <iostream>
+#include <vector>
+
+
+void test ()
+{
+  const unsigned int n_procs = Utilities::System::get_n_mpi_processes(MPI_COMM_WORLD);
+  const unsigned int my_id = Utilities::System::get_this_mpi_process(MPI_COMM_WORLD);
+
+  const unsigned int n_rows = 2;
+  const unsigned int n_cols = 2;
+
+  IndexSet row_partitioning (n_rows);
+  IndexSet col_partitioning (n_cols);
+
+  if (n_procs == 1)
+    {
+      row_partitioning.add_range(0, n_rows);
+      col_partitioning.add_range(0, n_cols);
+    }
+  else if (n_procs == 2)
+    {
+      if (my_id == 0)
+       {
+         row_partitioning.add_range(0, 1);
+         col_partitioning.add_range(0, 1);
+       }
+      else if (my_id == 1)
+       {
+         row_partitioning.add_range(1, n_rows);
+         col_partitioning.add_range(1, n_cols);
+       }
+    }
+  else
+    Assert (false, ExcNotImplemented());
+
+                                  /* A is
+
+                                     0     1
+                                     0     1
+                                  */
+  const unsigned int n_entries = 2;
+  const unsigned int line [n_entries] = {0, 1};
+  const unsigned int local_index [n_entries] = {1, 1};
+  const double local_value [n_entries] = {1.0, 1.0};
+
+  TrilinosWrappers::SparsityPattern sp (row_partitioning, col_partitioning, MPI_COMM_WORLD);
+  for (unsigned int i = 0; i < n_entries; ++i)
+    if (row_partitioning.is_element(line[i]))
+      sp.add(line[i], local_index[i]);
+  sp.compress();
+
+  TrilinosWrappers::SparseMatrix A;
+  A.clear ();
+  A.reinit (sp);
+  for (unsigned int i = 0; i<n_entries; ++i)
+    if (row_partitioning.is_element(line[i]))
+      A.add(line[i], local_index[i], local_value[i]);
+  A.compress();
+
+  if (my_id == 0)
+    {
+      Assert(A.el(0, 0) == 0, ExcMessage("Wrong element in A!"));
+      Assert(A.el(0, 1) == 1, ExcMessage("Wrong element in A!"));
+    }
+  if ((n_procs == 1) || (my_id == 1))
+    {
+      Assert(A.el(1, 0) == 0, ExcMessage("Wrong element in A!"));
+      Assert(A.el(1, 1) == 1, ExcMessage("Wrong element in A!"));
+    }
+
+  TrilinosWrappers::SparseMatrix AtA;
+  A.Tmmult (AtA, A);
+  AtA.compress ();
+
+                                  /* AtA should be
+
+                                     0     0
+                                     0     2
+                                  */
+
+                                  // checking AtA row partitioning
+  if (n_procs == 2)
+    {
+      if (my_id == 0)
+       {
+         Assert(AtA.local_range().first == 0,
+                ExcMessage("AtA Local Range is not as expected."));
+         Assert(AtA.local_range().second == 1,
+                ExcMessage("AtA Local Range is not as expected."));
+       }
+      if (my_id == 1)
+       {
+         Assert(AtA.local_range().first == 1,
+                ExcMessage("AtA Local Range is not as expected."));
+         Assert(AtA.local_range().second == 2,
+                ExcMessage("AtA Local Range is not as expected."));
+       }
+    }
+
+                                  // checking AtA elements. note that
+                                  // the el() function either returns
+                                  // the correct value, or zero in
+                                  // case the element is stored on
+                                  // other processors. consequently,
+                                  // in the following we only test
+                                  // that the values that *must* be
+                                  // zero in AtA are indeed zero, but
+                                  // not that the others have the
+                                  // correct value
+  Assert(AtA.el(0, 0) == 0, ExcMessage("Wrong element in AtA!"));
+  Assert(AtA.el(0, 1) == 0, ExcMessage("Wrong element in AtA!"));
+  Assert(AtA.el(1, 0) == 0, ExcMessage("Wrong element in AtA!"));
+
+                                  // now also check the one nonzero
+                                  // element
+  if ((n_procs == 1) || (my_id == 1))
+    Assert(AtA.el(1, 1) == 2, ExcMessage("Wrong element in AtA!"));
+
+  deallog << "OK" << endl;
+}
+
+
+
+int main (int argc, char **argv)
+{
+  Utilities::System::MPI_InitFinalize mpi_initialization(argc, argv);
+
+  unsigned int myid = Utilities::System::get_this_mpi_process (MPI_COMM_WORLD);
+  deallog.push(Utilities::int_to_string(myid));
+
+  if (myid == 0)
+    {
+      std::ofstream logfile(output_file_for_mpi("trilinos_sparse_matrix_mmult_01").c_str());
+      deallog.attach(logfile);
+      deallog << std::setprecision(4);
+      deallog.depth_console(0);
+      deallog.threshold_double(1.e-10);
+
+      test();
+    }
+  else
+    test();
+
+}
diff --git a/tests/mpi/trilinos_sparse_matrix_mmult_01/ncpu_1/cmp/generic b/tests/mpi/trilinos_sparse_matrix_mmult_01/ncpu_1/cmp/generic
new file mode 100644 (file)
index 0000000..be8d055
--- /dev/null
@@ -0,0 +1,2 @@
+
+DEAL:0::OK
diff --git a/tests/mpi/trilinos_sparse_matrix_mmult_01/ncpu_2/cmp/generic b/tests/mpi/trilinos_sparse_matrix_mmult_01/ncpu_2/cmp/generic
new file mode 100644 (file)
index 0000000..be8d055
--- /dev/null
@@ -0,0 +1,2 @@
+
+DEAL:0::OK

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.