]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix compressing TrilinosWrappers::SparsityPattern with empty column map
authorDaniel Arndt <arndtd@ornl.gov>
Thu, 23 May 2019 03:49:28 +0000 (23:49 -0400)
committerDaniel Arndt <arndtd@ornl.gov>
Thu, 23 May 2019 03:57:39 +0000 (23:57 -0400)
doc/news/changes/minor/20190522ArndtAnselmann [new file with mode: 0644]
source/lac/trilinos_sparsity_pattern.cc

diff --git a/doc/news/changes/minor/20190522ArndtAnselmann b/doc/news/changes/minor/20190522ArndtAnselmann
new file mode 100644 (file)
index 0000000..7ee6432
--- /dev/null
@@ -0,0 +1,3 @@
+Fixed: TrilinosWrappers::SparsityPattern can deal with empty column maps.
+<br>
+(Daniel Arndt, Mathias Anselmann, 2019/05/22)
index 84403175c62462ab295ac1c9429016c7bd444ea4..c4482795050d74cdd0d6c0b1df9ef15eb96b41e5 100644 (file)
@@ -791,7 +791,8 @@ namespace TrilinosWrappers
     if (nonlocal_graph.get() != nullptr)
       {
         if (nonlocal_graph->IndicesAreGlobal() == false &&
-            nonlocal_graph->RowMap().NumMyElements() > 0)
+            nonlocal_graph->RowMap().NumMyElements() > 0 &&
+            column_space_map->NumMyElements() > 0)
           {
             // Insert dummy element at (row, column) that corresponds to row 0
             // in local index counting.
@@ -813,6 +814,7 @@ namespace TrilinosWrappers
             AssertThrow(ierr == 0, ExcTrilinosError(ierr));
           }
         Assert(nonlocal_graph->RowMap().NumMyElements() == 0 ||
+                 column_space_map->NumMyElements() == 0 ||
                  nonlocal_graph->IndicesAreGlobal() == true,
                ExcInternalError());
 

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.