From: heister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Sun, 16 Jun 2013 21:20:41 +0000 (+0000)
Subject: add test
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=2a4acfe5e20c7685fec72b3e4877c0754716985f;p=dealii-svn.git

add test

git-svn-id: https://svn.dealii.org/branches/branch_unify_linear_algebra@29830 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/tests/gla/block_mat_02.cc b/tests/gla/block_mat_02.cc
new file mode 100644
index 0000000000..73794854e6
--- /dev/null
+++ b/tests/gla/block_mat_02.cc
@@ -0,0 +1,245 @@
+//---------------------------------------------------------------------------
+//    $Id: ghost_01.cc 28440 2013-02-17 14:35:08Z heister $
+//    Version: $Name$ 
+//
+//    Copyright (C) 2004, 2005, 2010 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.
+//
+//---------------------------------------------------------------------------
+
+
+// test problem in distribute_sparsity_pattern (block version)
+// now Trilinos and PETSc produce the same block matrix
+
+
+#define USE_PETSC_LA
+
+#include "../tests.h"
+#include <deal.II/lac/generic_linear_algebra.h>
+#include <deal.II/base/index_set.h>
+#include <fstream>
+#include <iostream>
+#include <iomanip>
+#include <vector>
+#include <deal.II/lac/full_matrix.h>
+#include <deal.II/lac/constraint_matrix.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_dgp.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/dofs/dof_tools.h>
+#include <deal.II/dofs/dof_renumbering.h>
+#include <deal.II/lac/sparsity_tools.h>
+#include <deal.II/lac/block_sparsity_pattern.h>
+
+#include "gla.h"
+
+template <class LA, int dim> 
+void test ()
+{
+  unsigned int myid = Utilities::MPI::this_mpi_process (MPI_COMM_WORLD);
+  unsigned int numproc = Utilities::MPI::n_mpi_processes (MPI_COMM_WORLD);
+  
+  if (myid==0)
+    deallog << "numproc=" << numproc << std::endl;
+
+
+  parallel::distributed::Triangulation<dim>
+    triangulation (MPI_COMM_WORLD,
+                     typename Triangulation<dim>::MeshSmoothing
+                     (Triangulation<dim>::smoothing_on_refinement |
+                      Triangulation<dim>::smoothing_on_coarsening));
+
+  GridGenerator::hyper_cube(triangulation);
+  triangulation.refine_global(2);
+
+  const FESystem<dim> stokes_fe(FE_Q<dim>(1),
+				1,FE_DGP<dim>(0),1);
+  
+  
+  DoFHandler<dim> stokes_dof_handler(triangulation);
+  std::vector<unsigned int> stokes_sub_blocks (1+1,0);
+  stokes_sub_blocks[1] = 1;
+  stokes_dof_handler.distribute_dofs (stokes_fe);
+  DoFRenumbering::component_wise (stokes_dof_handler, stokes_sub_blocks);
+
+  std::vector<types::global_dof_index> stokes_dofs_per_block (2);
+  DoFTools::count_dofs_per_block (stokes_dof_handler, stokes_dofs_per_block,
+                                    stokes_sub_blocks);
+
+  const unsigned int n_u = stokes_dofs_per_block[0],
+		     n_p = stokes_dofs_per_block[1];
+
+  std::vector<IndexSet> stokes_partitioning, stokes_relevant_partitioning;
+
+  IndexSet stokes_index_set = stokes_dof_handler.locally_owned_dofs();
+  stokes_partitioning.push_back(stokes_index_set.get_view(0,n_u));
+  stokes_partitioning.push_back(stokes_index_set.get_view(n_u,n_u+n_p));
+
+  IndexSet stokes_relevant_set;
+  DoFTools::extract_locally_relevant_dofs (stokes_dof_handler,
+					   stokes_relevant_set);
+  stokes_relevant_partitioning.push_back(stokes_relevant_set.get_view(0,n_u));
+  stokes_relevant_partitioning.push_back(stokes_relevant_set.get_view(n_u,n_u+n_p));
+
+  stokes_relevant_set.print(deallog.get_file_stream());
+  
+  ConstraintMatrix cm(stokes_relevant_set);
+      
+  DoFTools::make_hanging_node_constraints (stokes_dof_handler,
+                                               cm);
+
+  //boundary conditions ?
+  cm.close();
+
+      
+  BlockCompressedSimpleSparsityPattern sp(stokes_relevant_partitioning);
+
+  Table<2,DoFTools::Coupling> coupling (1+1, 1+1);
+    coupling[0][0] = DoFTools::always;
+    coupling[1][1] = DoFTools::always;
+    coupling[0][1] = DoFTools::always;
+
+    DoFTools::make_sparsity_pattern (stokes_dof_handler,
+        coupling, sp,
+        cm, false,
+        Utilities::MPI::this_mpi_process(MPI_COMM_WORLD));
+
+    SparsityTools::distribute_sparsity_pattern(sp,
+					       stokes_dof_handler.n_locally_owned_dofs_per_processor(),
+					       MPI_COMM_WORLD, stokes_relevant_set);
+
+    
+  sp.compress();
+
+  typename LA::MPI::BlockSparseMatrix stokes_matrix;
+  stokes_matrix.reinit (stokes_partitioning, sp, MPI_COMM_WORLD);
+
+  stokes_matrix.print(deallog.get_file_stream());
+
+
+				   // done
+  if (myid==0)
+    deallog << "OK" << std::endl;
+}
+
+template <int dim> 
+void test_LA_Trilinos ()
+{
+  typedef LA_Trilinos LA;
+  
+  unsigned int myid = Utilities::MPI::this_mpi_process (MPI_COMM_WORLD);
+  unsigned int numproc = Utilities::MPI::n_mpi_processes (MPI_COMM_WORLD);
+  
+  if (myid==0)
+    deallog << "numproc=" << numproc << std::endl;
+
+
+  parallel::distributed::Triangulation<dim>
+    triangulation (MPI_COMM_WORLD,
+                     typename Triangulation<dim>::MeshSmoothing
+                     (Triangulation<dim>::smoothing_on_refinement |
+                      Triangulation<dim>::smoothing_on_coarsening));
+
+  GridGenerator::hyper_cube(triangulation);
+  triangulation.refine_global(2);
+
+  const FESystem<dim> stokes_fe(FE_Q<dim>(1),
+				1,FE_DGP<dim>(0),1);
+  
+  
+  DoFHandler<dim> stokes_dof_handler(triangulation);
+  std::vector<unsigned int> stokes_sub_blocks (1+1,0);
+  stokes_sub_blocks[1] = 1;
+  stokes_dof_handler.distribute_dofs (stokes_fe);
+  DoFRenumbering::component_wise (stokes_dof_handler, stokes_sub_blocks);
+
+  std::vector<types::global_dof_index> stokes_dofs_per_block (2);
+  DoFTools::count_dofs_per_block (stokes_dof_handler, stokes_dofs_per_block,
+                                    stokes_sub_blocks);
+
+  const unsigned int n_u = stokes_dofs_per_block[0],
+		     n_p = stokes_dofs_per_block[1];
+
+  std::vector<IndexSet> stokes_partitioning, stokes_relevant_partitioning;
+
+  IndexSet stokes_index_set = stokes_dof_handler.locally_owned_dofs();
+  stokes_partitioning.push_back(stokes_index_set.get_view(0,n_u));
+  stokes_partitioning.push_back(stokes_index_set.get_view(n_u,n_u+n_p));
+
+  IndexSet stokes_relevant_set;
+  DoFTools::extract_locally_relevant_dofs (stokes_dof_handler,
+					   stokes_relevant_set);
+  stokes_relevant_partitioning.push_back(stokes_relevant_set.get_view(0,n_u));
+  stokes_relevant_partitioning.push_back(stokes_relevant_set.get_view(n_u,n_u+n_p));
+  stokes_relevant_set.print(deallog.get_file_stream());
+  
+  ConstraintMatrix cm(stokes_relevant_set);
+      
+  DoFTools::make_hanging_node_constraints (stokes_dof_handler,
+                                               cm);
+
+  //boundary conditions ?
+  cm.close();
+  TrilinosWrappers::BlockSparsityPattern
+    //typename LA::MPI::CompressedBlockSparsityPattern
+    sp(stokes_partitioning, MPI_COMM_WORLD);
+  
+
+    Table<2,DoFTools::Coupling> coupling (1+1, 1+1);
+
+    coupling[0][0] = DoFTools::always;
+    coupling[1][1] = DoFTools::always;
+    coupling[0][1] = DoFTools::always;
+
+    DoFTools::make_sparsity_pattern (stokes_dof_handler,
+                                     coupling, sp,
+                                     cm, false,
+                                     Utilities::MPI::
+                                     this_mpi_process(MPI_COMM_WORLD));
+    
+  sp.compress();
+
+  //sp.print(deallog.get_file_stream());
+
+  LA::MPI::BlockSparseMatrix stokes_matrix;
+  stokes_matrix.reinit (sp);
+
+  stokes_matrix.print(deallog.get_file_stream());
+
+  
+				   // done
+  if (myid==0)
+    deallog << "OK" << std::endl;
+}
+
+
+int main (int argc, char **argv)
+{
+  Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
+  MPILogInitAll log(__FILE__);
+
+  {
+    deallog.push("PETSc");
+    test<LA_PETSc,2>();
+    deallog.pop();	
+    deallog.push("Trilinos");
+    test<LA_Trilinos,2>();
+    deallog.pop();
+    deallog.push("TrilinosAlt");
+    test_LA_Trilinos<2>();
+    deallog.pop();	
+  }
+  
+  // compile, don't run
+  //if (myid==9999)
+    //  test<LA_Dummy>();
+  
+
+}
diff --git a/tests/gla/block_mat_02/ncpu_1/cmp/generic b/tests/gla/block_mat_02/ncpu_1/cmp/generic
new file mode 100644
index 0000000000..d64fb73026
--- /dev/null
+++ b/tests/gla/block_mat_02/ncpu_1/cmp/generic
@@ -0,0 +1,769 @@
+
+DEAL:0:PETSc::numproc=1
+{[0,40]}
+Block (0, 0)
+(0,0) 0.00000
+(0,1) 0.00000
+(0,2) 0.00000
+(0,3) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(1,2) 0.00000
+(1,3) 0.00000
+(1,4) 0.00000
+(1,5) 0.00000
+(2,0) 0.00000
+(2,1) 0.00000
+(2,2) 0.00000
+(2,3) 0.00000
+(2,6) 0.00000
+(2,7) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(3,4) 0.00000
+(3,5) 0.00000
+(3,6) 0.00000
+(3,7) 0.00000
+(3,8) 0.00000
+(4,1) 0.00000
+(4,3) 0.00000
+(4,4) 0.00000
+(4,5) 0.00000
+(4,9) 0.00000
+(4,10) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,5) 0.00000
+(5,7) 0.00000
+(5,8) 0.00000
+(5,9) 0.00000
+(5,10) 0.00000
+(5,13) 0.00000
+(6,2) 0.00000
+(6,3) 0.00000
+(6,6) 0.00000
+(6,7) 0.00000
+(6,15) 0.00000
+(6,16) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,5) 0.00000
+(7,6) 0.00000
+(7,7) 0.00000
+(7,8) 0.00000
+(7,15) 0.00000
+(7,16) 0.00000
+(7,17) 0.00000
+(8,3) 0.00000
+(8,5) 0.00000
+(8,7) 0.00000
+(8,8) 0.00000
+(8,10) 0.00000
+(8,13) 0.00000
+(8,16) 0.00000
+(8,17) 0.00000
+(8,21) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(9,9) 0.00000
+(9,10) 0.00000
+(9,11) 0.00000
+(9,12) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,8) 0.00000
+(10,9) 0.00000
+(10,10) 0.00000
+(10,11) 0.00000
+(10,12) 0.00000
+(10,13) 0.00000
+(10,14) 0.00000
+(11,9) 0.00000
+(11,10) 0.00000
+(11,11) 0.00000
+(11,12) 0.00000
+(12,9) 0.00000
+(12,10) 0.00000
+(12,11) 0.00000
+(12,12) 0.00000
+(12,13) 0.00000
+(12,14) 0.00000
+(13,5) 0.00000
+(13,8) 0.00000
+(13,10) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(13,14) 0.00000
+(13,17) 0.00000
+(13,21) 0.00000
+(13,22) 0.00000
+(14,10) 0.00000
+(14,12) 0.00000
+(14,13) 0.00000
+(14,14) 0.00000
+(14,21) 0.00000
+(14,22) 0.00000
+(15,6) 0.00000
+(15,7) 0.00000
+(15,15) 0.00000
+(15,16) 0.00000
+(15,18) 0.00000
+(15,19) 0.00000
+(16,6) 0.00000
+(16,7) 0.00000
+(16,8) 0.00000
+(16,15) 0.00000
+(16,16) 0.00000
+(16,17) 0.00000
+(16,18) 0.00000
+(16,19) 0.00000
+(16,20) 0.00000
+(17,7) 0.00000
+(17,8) 0.00000
+(17,13) 0.00000
+(17,16) 0.00000
+(17,17) 0.00000
+(17,19) 0.00000
+(17,20) 0.00000
+(17,21) 0.00000
+(17,23) 0.00000
+(18,15) 0.00000
+(18,16) 0.00000
+(18,18) 0.00000
+(18,19) 0.00000
+(19,15) 0.00000
+(19,16) 0.00000
+(19,17) 0.00000
+(19,18) 0.00000
+(19,19) 0.00000
+(19,20) 0.00000
+(20,16) 0.00000
+(20,17) 0.00000
+(20,19) 0.00000
+(20,20) 0.00000
+(20,21) 0.00000
+(20,23) 0.00000
+(21,8) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(21,17) 0.00000
+(21,20) 0.00000
+(21,21) 0.00000
+(21,22) 0.00000
+(21,23) 0.00000
+(21,24) 0.00000
+(22,13) 0.00000
+(22,14) 0.00000
+(22,21) 0.00000
+(22,22) 0.00000
+(22,23) 0.00000
+(22,24) 0.00000
+(23,17) 0.00000
+(23,20) 0.00000
+(23,21) 0.00000
+(23,22) 0.00000
+(23,23) 0.00000
+(23,24) 0.00000
+(24,21) 0.00000
+(24,22) 0.00000
+(24,23) 0.00000
+(24,24) 0.00000
+Block (0, 1)
+(0,0) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(2,0) 0.00000
+(2,2) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(4,1) 0.00000
+(4,4) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,6) 0.00000
+(6,2) 0.00000
+(6,8) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,8) 0.00000
+(7,9) 0.00000
+(8,3) 0.00000
+(8,6) 0.00000
+(8,9) 0.00000
+(8,12) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,6) 0.00000
+(10,7) 0.00000
+(11,5) 0.00000
+(12,5) 0.00000
+(12,7) 0.00000
+(13,6) 0.00000
+(13,7) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(14,7) 0.00000
+(14,13) 0.00000
+(15,8) 0.00000
+(15,10) 0.00000
+(16,8) 0.00000
+(16,9) 0.00000
+(16,10) 0.00000
+(16,11) 0.00000
+(17,9) 0.00000
+(17,11) 0.00000
+(17,12) 0.00000
+(17,14) 0.00000
+(18,10) 0.00000
+(19,10) 0.00000
+(19,11) 0.00000
+(20,11) 0.00000
+(20,14) 0.00000
+(21,12) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(21,15) 0.00000
+(22,13) 0.00000
+(22,15) 0.00000
+(23,14) 0.00000
+(23,15) 0.00000
+(24,15) 0.00000
+Block (1, 0)
+Block (1, 1)
+(0,0) 0.00000
+(1,1) 0.00000
+(2,2) 0.00000
+(3,3) 0.00000
+(4,4) 0.00000
+(5,5) 0.00000
+(6,6) 0.00000
+(7,7) 0.00000
+(8,8) 0.00000
+(9,9) 0.00000
+(10,10) 0.00000
+(11,11) 0.00000
+(12,12) 0.00000
+(13,13) 0.00000
+(14,14) 0.00000
+(15,15) 0.00000
+DEAL:0:PETSc::OK
+DEAL:0:Trilinos::numproc=1
+{[0,40]}
+Block (0, 0)
+(0,0) 0.00000
+(0,1) 0.00000
+(0,2) 0.00000
+(0,3) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(1,2) 0.00000
+(1,3) 0.00000
+(1,4) 0.00000
+(1,5) 0.00000
+(2,0) 0.00000
+(2,1) 0.00000
+(2,2) 0.00000
+(2,3) 0.00000
+(2,6) 0.00000
+(2,7) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(3,4) 0.00000
+(3,5) 0.00000
+(3,6) 0.00000
+(3,7) 0.00000
+(3,8) 0.00000
+(4,1) 0.00000
+(4,3) 0.00000
+(4,4) 0.00000
+(4,5) 0.00000
+(4,9) 0.00000
+(4,10) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,5) 0.00000
+(5,7) 0.00000
+(5,8) 0.00000
+(5,9) 0.00000
+(5,10) 0.00000
+(5,13) 0.00000
+(6,2) 0.00000
+(6,3) 0.00000
+(6,6) 0.00000
+(6,7) 0.00000
+(6,15) 0.00000
+(6,16) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,5) 0.00000
+(7,6) 0.00000
+(7,7) 0.00000
+(7,8) 0.00000
+(7,15) 0.00000
+(7,16) 0.00000
+(7,17) 0.00000
+(8,3) 0.00000
+(8,5) 0.00000
+(8,7) 0.00000
+(8,8) 0.00000
+(8,10) 0.00000
+(8,13) 0.00000
+(8,16) 0.00000
+(8,17) 0.00000
+(8,21) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(9,9) 0.00000
+(9,10) 0.00000
+(9,11) 0.00000
+(9,12) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,8) 0.00000
+(10,9) 0.00000
+(10,10) 0.00000
+(10,11) 0.00000
+(10,12) 0.00000
+(10,13) 0.00000
+(10,14) 0.00000
+(11,9) 0.00000
+(11,10) 0.00000
+(11,11) 0.00000
+(11,12) 0.00000
+(12,9) 0.00000
+(12,10) 0.00000
+(12,11) 0.00000
+(12,12) 0.00000
+(12,13) 0.00000
+(12,14) 0.00000
+(13,5) 0.00000
+(13,8) 0.00000
+(13,10) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(13,14) 0.00000
+(13,17) 0.00000
+(13,21) 0.00000
+(13,22) 0.00000
+(14,10) 0.00000
+(14,12) 0.00000
+(14,13) 0.00000
+(14,14) 0.00000
+(14,21) 0.00000
+(14,22) 0.00000
+(15,6) 0.00000
+(15,7) 0.00000
+(15,15) 0.00000
+(15,16) 0.00000
+(15,18) 0.00000
+(15,19) 0.00000
+(16,6) 0.00000
+(16,7) 0.00000
+(16,8) 0.00000
+(16,15) 0.00000
+(16,16) 0.00000
+(16,17) 0.00000
+(16,18) 0.00000
+(16,19) 0.00000
+(16,20) 0.00000
+(17,7) 0.00000
+(17,8) 0.00000
+(17,13) 0.00000
+(17,16) 0.00000
+(17,17) 0.00000
+(17,19) 0.00000
+(17,20) 0.00000
+(17,21) 0.00000
+(17,23) 0.00000
+(18,15) 0.00000
+(18,16) 0.00000
+(18,18) 0.00000
+(18,19) 0.00000
+(19,15) 0.00000
+(19,16) 0.00000
+(19,17) 0.00000
+(19,18) 0.00000
+(19,19) 0.00000
+(19,20) 0.00000
+(20,16) 0.00000
+(20,17) 0.00000
+(20,19) 0.00000
+(20,20) 0.00000
+(20,21) 0.00000
+(20,23) 0.00000
+(21,8) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(21,17) 0.00000
+(21,20) 0.00000
+(21,21) 0.00000
+(21,22) 0.00000
+(21,23) 0.00000
+(21,24) 0.00000
+(22,13) 0.00000
+(22,14) 0.00000
+(22,21) 0.00000
+(22,22) 0.00000
+(22,23) 0.00000
+(22,24) 0.00000
+(23,17) 0.00000
+(23,20) 0.00000
+(23,21) 0.00000
+(23,22) 0.00000
+(23,23) 0.00000
+(23,24) 0.00000
+(24,21) 0.00000
+(24,22) 0.00000
+(24,23) 0.00000
+(24,24) 0.00000
+Block (0, 1)
+(0,0) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(2,0) 0.00000
+(2,2) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(4,1) 0.00000
+(4,4) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,6) 0.00000
+(6,2) 0.00000
+(6,8) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,8) 0.00000
+(7,9) 0.00000
+(8,3) 0.00000
+(8,6) 0.00000
+(8,9) 0.00000
+(8,12) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,6) 0.00000
+(10,7) 0.00000
+(11,5) 0.00000
+(12,5) 0.00000
+(12,7) 0.00000
+(13,6) 0.00000
+(13,7) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(14,7) 0.00000
+(14,13) 0.00000
+(15,8) 0.00000
+(15,10) 0.00000
+(16,8) 0.00000
+(16,9) 0.00000
+(16,10) 0.00000
+(16,11) 0.00000
+(17,9) 0.00000
+(17,11) 0.00000
+(17,12) 0.00000
+(17,14) 0.00000
+(18,10) 0.00000
+(19,10) 0.00000
+(19,11) 0.00000
+(20,11) 0.00000
+(20,14) 0.00000
+(21,12) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(21,15) 0.00000
+(22,13) 0.00000
+(22,15) 0.00000
+(23,14) 0.00000
+(23,15) 0.00000
+(24,15) 0.00000
+Block (1, 0)
+Block (1, 1)
+(0,0) 0.00000
+(1,1) 0.00000
+(2,2) 0.00000
+(3,3) 0.00000
+(4,4) 0.00000
+(5,5) 0.00000
+(6,6) 0.00000
+(7,7) 0.00000
+(8,8) 0.00000
+(9,9) 0.00000
+(10,10) 0.00000
+(11,11) 0.00000
+(12,12) 0.00000
+(13,13) 0.00000
+(14,14) 0.00000
+(15,15) 0.00000
+DEAL:0:Trilinos::OK
+DEAL:0:TrilinosAlt::numproc=1
+{[0,40]}
+Block (0, 0)
+(0,0) 0.00000
+(0,1) 0.00000
+(0,2) 0.00000
+(0,3) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(1,2) 0.00000
+(1,3) 0.00000
+(1,4) 0.00000
+(1,5) 0.00000
+(2,0) 0.00000
+(2,1) 0.00000
+(2,2) 0.00000
+(2,3) 0.00000
+(2,6) 0.00000
+(2,7) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(3,4) 0.00000
+(3,5) 0.00000
+(3,6) 0.00000
+(3,7) 0.00000
+(3,8) 0.00000
+(4,1) 0.00000
+(4,3) 0.00000
+(4,4) 0.00000
+(4,5) 0.00000
+(4,9) 0.00000
+(4,10) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,5) 0.00000
+(5,7) 0.00000
+(5,8) 0.00000
+(5,9) 0.00000
+(5,10) 0.00000
+(5,13) 0.00000
+(6,2) 0.00000
+(6,3) 0.00000
+(6,6) 0.00000
+(6,7) 0.00000
+(6,15) 0.00000
+(6,16) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,5) 0.00000
+(7,6) 0.00000
+(7,7) 0.00000
+(7,8) 0.00000
+(7,15) 0.00000
+(7,16) 0.00000
+(7,17) 0.00000
+(8,3) 0.00000
+(8,5) 0.00000
+(8,7) 0.00000
+(8,8) 0.00000
+(8,10) 0.00000
+(8,13) 0.00000
+(8,16) 0.00000
+(8,17) 0.00000
+(8,21) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(9,9) 0.00000
+(9,10) 0.00000
+(9,11) 0.00000
+(9,12) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,8) 0.00000
+(10,9) 0.00000
+(10,10) 0.00000
+(10,11) 0.00000
+(10,12) 0.00000
+(10,13) 0.00000
+(10,14) 0.00000
+(11,9) 0.00000
+(11,10) 0.00000
+(11,11) 0.00000
+(11,12) 0.00000
+(12,9) 0.00000
+(12,10) 0.00000
+(12,11) 0.00000
+(12,12) 0.00000
+(12,13) 0.00000
+(12,14) 0.00000
+(13,5) 0.00000
+(13,8) 0.00000
+(13,10) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(13,14) 0.00000
+(13,17) 0.00000
+(13,21) 0.00000
+(13,22) 0.00000
+(14,10) 0.00000
+(14,12) 0.00000
+(14,13) 0.00000
+(14,14) 0.00000
+(14,21) 0.00000
+(14,22) 0.00000
+(15,6) 0.00000
+(15,7) 0.00000
+(15,15) 0.00000
+(15,16) 0.00000
+(15,18) 0.00000
+(15,19) 0.00000
+(16,6) 0.00000
+(16,7) 0.00000
+(16,8) 0.00000
+(16,15) 0.00000
+(16,16) 0.00000
+(16,17) 0.00000
+(16,18) 0.00000
+(16,19) 0.00000
+(16,20) 0.00000
+(17,7) 0.00000
+(17,8) 0.00000
+(17,13) 0.00000
+(17,16) 0.00000
+(17,17) 0.00000
+(17,19) 0.00000
+(17,20) 0.00000
+(17,21) 0.00000
+(17,23) 0.00000
+(18,15) 0.00000
+(18,16) 0.00000
+(18,18) 0.00000
+(18,19) 0.00000
+(19,15) 0.00000
+(19,16) 0.00000
+(19,17) 0.00000
+(19,18) 0.00000
+(19,19) 0.00000
+(19,20) 0.00000
+(20,16) 0.00000
+(20,17) 0.00000
+(20,19) 0.00000
+(20,20) 0.00000
+(20,21) 0.00000
+(20,23) 0.00000
+(21,8) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(21,17) 0.00000
+(21,20) 0.00000
+(21,21) 0.00000
+(21,22) 0.00000
+(21,23) 0.00000
+(21,24) 0.00000
+(22,13) 0.00000
+(22,14) 0.00000
+(22,21) 0.00000
+(22,22) 0.00000
+(22,23) 0.00000
+(22,24) 0.00000
+(23,17) 0.00000
+(23,20) 0.00000
+(23,21) 0.00000
+(23,22) 0.00000
+(23,23) 0.00000
+(23,24) 0.00000
+(24,21) 0.00000
+(24,22) 0.00000
+(24,23) 0.00000
+(24,24) 0.00000
+Block (0, 1)
+(0,0) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(2,0) 0.00000
+(2,2) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(4,1) 0.00000
+(4,4) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,6) 0.00000
+(6,2) 0.00000
+(6,8) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,8) 0.00000
+(7,9) 0.00000
+(8,3) 0.00000
+(8,6) 0.00000
+(8,9) 0.00000
+(8,12) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,6) 0.00000
+(10,7) 0.00000
+(11,5) 0.00000
+(12,5) 0.00000
+(12,7) 0.00000
+(13,6) 0.00000
+(13,7) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(14,7) 0.00000
+(14,13) 0.00000
+(15,8) 0.00000
+(15,10) 0.00000
+(16,8) 0.00000
+(16,9) 0.00000
+(16,10) 0.00000
+(16,11) 0.00000
+(17,9) 0.00000
+(17,11) 0.00000
+(17,12) 0.00000
+(17,14) 0.00000
+(18,10) 0.00000
+(19,10) 0.00000
+(19,11) 0.00000
+(20,11) 0.00000
+(20,14) 0.00000
+(21,12) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(21,15) 0.00000
+(22,13) 0.00000
+(22,15) 0.00000
+(23,14) 0.00000
+(23,15) 0.00000
+(24,15) 0.00000
+Block (1, 0)
+Block (1, 1)
+(0,0) 0.00000
+(1,1) 0.00000
+(2,2) 0.00000
+(3,3) 0.00000
+(4,4) 0.00000
+(5,5) 0.00000
+(6,6) 0.00000
+(7,7) 0.00000
+(8,8) 0.00000
+(9,9) 0.00000
+(10,10) 0.00000
+(11,11) 0.00000
+(12,12) 0.00000
+(13,13) 0.00000
+(14,14) 0.00000
+(15,15) 0.00000
+DEAL:0:TrilinosAlt::OK
diff --git a/tests/gla/block_mat_02/ncpu_2/cmp/generic b/tests/gla/block_mat_02/ncpu_2/cmp/generic
new file mode 100644
index 0000000000..ba1c48878a
--- /dev/null
+++ b/tests/gla/block_mat_02/ncpu_2/cmp/generic
@@ -0,0 +1,786 @@
+
+DEAL:0:PETSc::numproc=2
+{[0,17], [21,22], [25,34], [37,38]}
+Block (0, 0)
+(0,0) 0.00000
+(0,1) 0.00000
+(0,2) 0.00000
+(0,3) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(1,2) 0.00000
+(1,3) 0.00000
+(1,4) 0.00000
+(1,5) 0.00000
+(2,0) 0.00000
+(2,1) 0.00000
+(2,2) 0.00000
+(2,3) 0.00000
+(2,6) 0.00000
+(2,7) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(3,4) 0.00000
+(3,5) 0.00000
+(3,6) 0.00000
+(3,7) 0.00000
+(3,8) 0.00000
+(4,1) 0.00000
+(4,3) 0.00000
+(4,4) 0.00000
+(4,5) 0.00000
+(4,9) 0.00000
+(4,10) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,5) 0.00000
+(5,7) 0.00000
+(5,8) 0.00000
+(5,9) 0.00000
+(5,10) 0.00000
+(5,13) 0.00000
+(6,2) 0.00000
+(6,3) 0.00000
+(6,6) 0.00000
+(6,7) 0.00000
+(6,15) 0.00000
+(6,16) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,5) 0.00000
+(7,6) 0.00000
+(7,7) 0.00000
+(7,8) 0.00000
+(7,15) 0.00000
+(7,16) 0.00000
+(7,17) 0.00000
+(8,3) 0.00000
+(8,5) 0.00000
+(8,7) 0.00000
+(8,8) 0.00000
+(8,10) 0.00000
+(8,13) 0.00000
+(8,16) 0.00000
+(8,17) 0.00000
+(8,21) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(9,9) 0.00000
+(9,10) 0.00000
+(9,11) 0.00000
+(9,12) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,8) 0.00000
+(10,9) 0.00000
+(10,10) 0.00000
+(10,11) 0.00000
+(10,12) 0.00000
+(10,13) 0.00000
+(10,14) 0.00000
+(11,9) 0.00000
+(11,10) 0.00000
+(11,11) 0.00000
+(11,12) 0.00000
+(12,9) 0.00000
+(12,10) 0.00000
+(12,11) 0.00000
+(12,12) 0.00000
+(12,13) 0.00000
+(12,14) 0.00000
+(13,5) 0.00000
+(13,8) 0.00000
+(13,10) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(13,14) 0.00000
+(13,17) 0.00000
+(13,21) 0.00000
+(13,22) 0.00000
+(14,10) 0.00000
+(14,12) 0.00000
+(14,13) 0.00000
+(14,14) 0.00000
+(14,21) 0.00000
+(14,22) 0.00000
+Block (0, 1)
+(0,0) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(2,0) 0.00000
+(2,2) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(4,1) 0.00000
+(4,4) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,6) 0.00000
+(6,2) 0.00000
+(6,8) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,8) 0.00000
+(7,9) 0.00000
+(8,3) 0.00000
+(8,6) 0.00000
+(8,9) 0.00000
+(8,12) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,6) 0.00000
+(10,7) 0.00000
+(11,5) 0.00000
+(12,5) 0.00000
+(12,7) 0.00000
+(13,6) 0.00000
+(13,7) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(14,7) 0.00000
+(14,13) 0.00000
+Block (1, 0)
+Block (1, 1)
+(0,0) 0.00000
+(1,1) 0.00000
+(2,2) 0.00000
+(3,3) 0.00000
+(4,4) 0.00000
+(5,5) 0.00000
+(6,6) 0.00000
+(7,7) 0.00000
+DEAL:0:PETSc::OK
+DEAL:0:Trilinos::numproc=2
+{[0,17], [21,22], [25,34], [37,38]}
+Block (0, 0)
+(0,0) 0.00000
+(0,1) 0.00000
+(0,2) 0.00000
+(0,3) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(1,2) 0.00000
+(1,3) 0.00000
+(1,4) 0.00000
+(1,5) 0.00000
+(2,0) 0.00000
+(2,1) 0.00000
+(2,2) 0.00000
+(2,3) 0.00000
+(2,6) 0.00000
+(2,7) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(3,4) 0.00000
+(3,5) 0.00000
+(3,6) 0.00000
+(3,7) 0.00000
+(3,8) 0.00000
+(4,1) 0.00000
+(4,3) 0.00000
+(4,4) 0.00000
+(4,5) 0.00000
+(4,9) 0.00000
+(4,10) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,5) 0.00000
+(5,7) 0.00000
+(5,8) 0.00000
+(5,9) 0.00000
+(5,10) 0.00000
+(5,13) 0.00000
+(6,2) 0.00000
+(6,3) 0.00000
+(6,6) 0.00000
+(6,7) 0.00000
+(6,15) 0.00000
+(6,16) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,5) 0.00000
+(7,6) 0.00000
+(7,7) 0.00000
+(7,8) 0.00000
+(7,15) 0.00000
+(7,16) 0.00000
+(7,17) 0.00000
+(8,3) 0.00000
+(8,5) 0.00000
+(8,7) 0.00000
+(8,8) 0.00000
+(8,10) 0.00000
+(8,13) 0.00000
+(8,16) 0.00000
+(8,17) 0.00000
+(8,21) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(9,9) 0.00000
+(9,10) 0.00000
+(9,11) 0.00000
+(9,12) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,8) 0.00000
+(10,9) 0.00000
+(10,10) 0.00000
+(10,11) 0.00000
+(10,12) 0.00000
+(10,13) 0.00000
+(10,14) 0.00000
+(11,9) 0.00000
+(11,10) 0.00000
+(11,11) 0.00000
+(11,12) 0.00000
+(12,9) 0.00000
+(12,10) 0.00000
+(12,11) 0.00000
+(12,12) 0.00000
+(12,13) 0.00000
+(12,14) 0.00000
+(13,5) 0.00000
+(13,8) 0.00000
+(13,10) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(13,14) 0.00000
+(13,17) 0.00000
+(13,21) 0.00000
+(13,22) 0.00000
+(14,10) 0.00000
+(14,12) 0.00000
+(14,13) 0.00000
+(14,14) 0.00000
+(14,21) 0.00000
+(14,22) 0.00000
+Block (0, 1)
+(0,0) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(2,0) 0.00000
+(2,2) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(4,1) 0.00000
+(4,4) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,6) 0.00000
+(6,2) 0.00000
+(6,8) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,8) 0.00000
+(7,9) 0.00000
+(8,3) 0.00000
+(8,6) 0.00000
+(8,9) 0.00000
+(8,12) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,6) 0.00000
+(10,7) 0.00000
+(11,5) 0.00000
+(12,5) 0.00000
+(12,7) 0.00000
+(13,6) 0.00000
+(13,7) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(14,7) 0.00000
+(14,13) 0.00000
+Block (1, 0)
+Block (1, 1)
+(0,0) 0.00000
+(1,1) 0.00000
+(2,2) 0.00000
+(3,3) 0.00000
+(4,4) 0.00000
+(5,5) 0.00000
+(6,6) 0.00000
+(7,7) 0.00000
+DEAL:0:Trilinos::OK
+DEAL:0:TrilinosAlt::numproc=2
+{[0,17], [21,22], [25,34], [37,38]}
+Block (0, 0)
+(0,0) 0.00000
+(0,1) 0.00000
+(0,2) 0.00000
+(0,3) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(1,2) 0.00000
+(1,3) 0.00000
+(1,4) 0.00000
+(1,5) 0.00000
+(2,0) 0.00000
+(2,1) 0.00000
+(2,2) 0.00000
+(2,3) 0.00000
+(2,6) 0.00000
+(2,7) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(3,4) 0.00000
+(3,5) 0.00000
+(3,6) 0.00000
+(3,7) 0.00000
+(3,8) 0.00000
+(4,1) 0.00000
+(4,3) 0.00000
+(4,4) 0.00000
+(4,5) 0.00000
+(4,9) 0.00000
+(4,10) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,5) 0.00000
+(5,7) 0.00000
+(5,8) 0.00000
+(5,9) 0.00000
+(5,10) 0.00000
+(5,13) 0.00000
+(6,2) 0.00000
+(6,3) 0.00000
+(6,6) 0.00000
+(6,7) 0.00000
+(6,15) 0.00000
+(6,16) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,5) 0.00000
+(7,6) 0.00000
+(7,7) 0.00000
+(7,8) 0.00000
+(7,15) 0.00000
+(7,16) 0.00000
+(7,17) 0.00000
+(8,3) 0.00000
+(8,5) 0.00000
+(8,7) 0.00000
+(8,8) 0.00000
+(8,10) 0.00000
+(8,13) 0.00000
+(8,16) 0.00000
+(8,17) 0.00000
+(8,21) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(9,9) 0.00000
+(9,10) 0.00000
+(9,11) 0.00000
+(9,12) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,8) 0.00000
+(10,9) 0.00000
+(10,10) 0.00000
+(10,11) 0.00000
+(10,12) 0.00000
+(10,13) 0.00000
+(10,14) 0.00000
+(11,9) 0.00000
+(11,10) 0.00000
+(11,11) 0.00000
+(11,12) 0.00000
+(12,9) 0.00000
+(12,10) 0.00000
+(12,11) 0.00000
+(12,12) 0.00000
+(12,13) 0.00000
+(12,14) 0.00000
+(13,5) 0.00000
+(13,8) 0.00000
+(13,10) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(13,14) 0.00000
+(13,17) 0.00000
+(13,21) 0.00000
+(13,22) 0.00000
+(14,10) 0.00000
+(14,12) 0.00000
+(14,13) 0.00000
+(14,14) 0.00000
+(14,21) 0.00000
+(14,22) 0.00000
+Block (0, 1)
+(0,0) 0.00000
+(1,0) 0.00000
+(1,1) 0.00000
+(2,0) 0.00000
+(2,2) 0.00000
+(3,0) 0.00000
+(3,1) 0.00000
+(3,2) 0.00000
+(3,3) 0.00000
+(4,1) 0.00000
+(4,4) 0.00000
+(5,1) 0.00000
+(5,3) 0.00000
+(5,4) 0.00000
+(5,6) 0.00000
+(6,2) 0.00000
+(6,8) 0.00000
+(7,2) 0.00000
+(7,3) 0.00000
+(7,8) 0.00000
+(7,9) 0.00000
+(8,3) 0.00000
+(8,6) 0.00000
+(8,9) 0.00000
+(8,12) 0.00000
+(9,4) 0.00000
+(9,5) 0.00000
+(10,4) 0.00000
+(10,5) 0.00000
+(10,6) 0.00000
+(10,7) 0.00000
+(11,5) 0.00000
+(12,5) 0.00000
+(12,7) 0.00000
+(13,6) 0.00000
+(13,7) 0.00000
+(13,12) 0.00000
+(13,13) 0.00000
+(14,7) 0.00000
+(14,13) 0.00000
+Block (1, 0)
+Block (1, 1)
+(0,0) 0.00000
+(1,1) 0.00000
+(2,2) 0.00000
+(3,3) 0.00000
+(4,4) 0.00000
+(5,5) 0.00000
+(6,6) 0.00000
+(7,7) 0.00000
+DEAL:0:TrilinosAlt::OK
+
+{[2,3], [5,8], 10, [12,24], [27,28], [31,40]}
+Block (0, 0)
+(15,6) 0.00000
+(15,7) 0.00000
+(15,15) 0.00000
+(15,16) 0.00000
+(15,18) 0.00000
+(15,19) 0.00000
+(16,6) 0.00000
+(16,7) 0.00000
+(16,8) 0.00000
+(16,15) 0.00000
+(16,16) 0.00000
+(16,17) 0.00000
+(16,18) 0.00000
+(16,19) 0.00000
+(16,20) 0.00000
+(17,7) 0.00000
+(17,8) 0.00000
+(17,13) 0.00000
+(17,16) 0.00000
+(17,17) 0.00000
+(17,19) 0.00000
+(17,20) 0.00000
+(17,21) 0.00000
+(17,23) 0.00000
+(18,15) 0.00000
+(18,16) 0.00000
+(18,18) 0.00000
+(18,19) 0.00000
+(19,15) 0.00000
+(19,16) 0.00000
+(19,17) 0.00000
+(19,18) 0.00000
+(19,19) 0.00000
+(19,20) 0.00000
+(20,16) 0.00000
+(20,17) 0.00000
+(20,19) 0.00000
+(20,20) 0.00000
+(20,21) 0.00000
+(20,23) 0.00000
+(21,8) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(21,17) 0.00000
+(21,20) 0.00000
+(21,21) 0.00000
+(21,22) 0.00000
+(21,23) 0.00000
+(21,24) 0.00000
+(22,13) 0.00000
+(22,14) 0.00000
+(22,21) 0.00000
+(22,22) 0.00000
+(22,23) 0.00000
+(22,24) 0.00000
+(23,17) 0.00000
+(23,20) 0.00000
+(23,21) 0.00000
+(23,22) 0.00000
+(23,23) 0.00000
+(23,24) 0.00000
+(24,21) 0.00000
+(24,22) 0.00000
+(24,23) 0.00000
+(24,24) 0.00000
+Block (0, 1)
+(15,8) 0.00000
+(15,10) 0.00000
+(16,8) 0.00000
+(16,9) 0.00000
+(16,10) 0.00000
+(16,11) 0.00000
+(17,9) 0.00000
+(17,11) 0.00000
+(17,12) 0.00000
+(17,14) 0.00000
+(18,10) 0.00000
+(19,10) 0.00000
+(19,11) 0.00000
+(20,11) 0.00000
+(20,14) 0.00000
+(21,12) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(21,15) 0.00000
+(22,13) 0.00000
+(22,15) 0.00000
+(23,14) 0.00000
+(23,15) 0.00000
+(24,15) 0.00000
+Block (1, 0)
+Block (1, 1)
+(8,8) 0.00000
+(9,9) 0.00000
+(10,10) 0.00000
+(11,11) 0.00000
+(12,12) 0.00000
+(13,13) 0.00000
+(14,14) 0.00000
+(15,15) 0.00000
+{[2,3], [5,8], 10, [12,24], [27,28], [31,40]}
+Block (0, 0)
+(15,15) 0.00000
+(15,16) 0.00000
+(15,18) 0.00000
+(15,19) 0.00000
+(15,6) 0.00000
+(15,7) 0.00000
+(16,15) 0.00000
+(16,16) 0.00000
+(16,17) 0.00000
+(16,18) 0.00000
+(16,19) 0.00000
+(16,20) 0.00000
+(16,6) 0.00000
+(16,7) 0.00000
+(16,8) 0.00000
+(17,16) 0.00000
+(17,17) 0.00000
+(17,19) 0.00000
+(17,20) 0.00000
+(17,21) 0.00000
+(17,23) 0.00000
+(17,7) 0.00000
+(17,8) 0.00000
+(17,13) 0.00000
+(18,15) 0.00000
+(18,16) 0.00000
+(18,18) 0.00000
+(18,19) 0.00000
+(19,15) 0.00000
+(19,16) 0.00000
+(19,17) 0.00000
+(19,18) 0.00000
+(19,19) 0.00000
+(19,20) 0.00000
+(20,16) 0.00000
+(20,17) 0.00000
+(20,19) 0.00000
+(20,20) 0.00000
+(20,21) 0.00000
+(20,23) 0.00000
+(21,17) 0.00000
+(21,20) 0.00000
+(21,21) 0.00000
+(21,22) 0.00000
+(21,23) 0.00000
+(21,24) 0.00000
+(21,8) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(22,21) 0.00000
+(22,22) 0.00000
+(22,23) 0.00000
+(22,24) 0.00000
+(22,13) 0.00000
+(22,14) 0.00000
+(23,17) 0.00000
+(23,20) 0.00000
+(23,21) 0.00000
+(23,22) 0.00000
+(23,23) 0.00000
+(23,24) 0.00000
+(24,21) 0.00000
+(24,22) 0.00000
+(24,23) 0.00000
+(24,24) 0.00000
+Block (0, 1)
+(15,8) 0.00000
+(15,10) 0.00000
+(16,8) 0.00000
+(16,9) 0.00000
+(16,10) 0.00000
+(16,11) 0.00000
+(17,9) 0.00000
+(17,11) 0.00000
+(17,12) 0.00000
+(17,14) 0.00000
+(18,10) 0.00000
+(19,10) 0.00000
+(19,11) 0.00000
+(20,11) 0.00000
+(20,14) 0.00000
+(21,12) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(21,15) 0.00000
+(22,13) 0.00000
+(22,15) 0.00000
+(23,14) 0.00000
+(23,15) 0.00000
+(24,15) 0.00000
+Block (1, 0)
+Block (1, 1)
+(8,8) 0.00000
+(9,9) 0.00000
+(10,10) 0.00000
+(11,11) 0.00000
+(12,12) 0.00000
+(13,13) 0.00000
+(14,14) 0.00000
+(15,15) 0.00000
+{[2,3], [5,8], 10, [12,24], [27,28], [31,40]}
+Block (0, 0)
+(15,15) 0.00000
+(15,16) 0.00000
+(15,18) 0.00000
+(15,19) 0.00000
+(15,6) 0.00000
+(15,7) 0.00000
+(16,15) 0.00000
+(16,16) 0.00000
+(16,17) 0.00000
+(16,18) 0.00000
+(16,19) 0.00000
+(16,20) 0.00000
+(16,6) 0.00000
+(16,7) 0.00000
+(16,8) 0.00000
+(17,16) 0.00000
+(17,17) 0.00000
+(17,19) 0.00000
+(17,20) 0.00000
+(17,21) 0.00000
+(17,23) 0.00000
+(17,7) 0.00000
+(17,8) 0.00000
+(17,13) 0.00000
+(18,15) 0.00000
+(18,16) 0.00000
+(18,18) 0.00000
+(18,19) 0.00000
+(19,15) 0.00000
+(19,16) 0.00000
+(19,17) 0.00000
+(19,18) 0.00000
+(19,19) 0.00000
+(19,20) 0.00000
+(20,16) 0.00000
+(20,17) 0.00000
+(20,19) 0.00000
+(20,20) 0.00000
+(20,21) 0.00000
+(20,23) 0.00000
+(21,17) 0.00000
+(21,20) 0.00000
+(21,21) 0.00000
+(21,22) 0.00000
+(21,23) 0.00000
+(21,24) 0.00000
+(21,8) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(22,21) 0.00000
+(22,22) 0.00000
+(22,23) 0.00000
+(22,24) 0.00000
+(22,13) 0.00000
+(22,14) 0.00000
+(23,17) 0.00000
+(23,20) 0.00000
+(23,21) 0.00000
+(23,22) 0.00000
+(23,23) 0.00000
+(23,24) 0.00000
+(24,21) 0.00000
+(24,22) 0.00000
+(24,23) 0.00000
+(24,24) 0.00000
+Block (0, 1)
+(15,8) 0.00000
+(15,10) 0.00000
+(16,8) 0.00000
+(16,9) 0.00000
+(16,10) 0.00000
+(16,11) 0.00000
+(17,9) 0.00000
+(17,11) 0.00000
+(17,12) 0.00000
+(17,14) 0.00000
+(18,10) 0.00000
+(19,10) 0.00000
+(19,11) 0.00000
+(20,11) 0.00000
+(20,14) 0.00000
+(21,12) 0.00000
+(21,13) 0.00000
+(21,14) 0.00000
+(21,15) 0.00000
+(22,13) 0.00000
+(22,15) 0.00000
+(23,14) 0.00000
+(23,15) 0.00000
+(24,15) 0.00000
+Block (1, 0)
+Block (1, 1)
+(8,8) 0.00000
+(9,9) 0.00000
+(10,10) 0.00000
+(11,11) 0.00000
+(12,12) 0.00000
+(13,13) 0.00000
+(14,14) 0.00000
+(15,15) 0.00000
+