When doing reinit() on a distributed vector, it is not enough to ensure that the...

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

diff --git a/deal.II/lac/source/trilinos_vector.cc b/deal.II/lac/source/trilinos_vector.cc
index 5197e98cb44c0dd1acc569b43f3ba4b9c930af0d..b53d8bed257890b8709954b09bc17c7e585fc08b 100755 (executable)
--- a/deal.II/lac/source/trilinos_vector.cc
+++ b/deal.II/lac/source/trilinos_vector.cc
@@ -133,7 +133,7 @@ namespace TrilinosWrappers
                                        // generate the vector.
       if (allow_different_maps == false)
         {
-         if (local_range() != v.local_range())
+         if (vector->Map().SameAs(v.vector->Map()) == false)
            {
              vector.reset();
              communicator.reset (Utilities::Trilinos::
@@ -147,11 +147,12 @@ namespace TrilinosWrappers
            }
          else if (fast == false)
            {
+                                              // old and new vectors
+                                              // have exactly the
+                                              // same map, i.e. size
+                                              // and parallel
+                                              // distribution
              int ierr;
-             Assert (vector->Map().SameAs(v.vector->Map()) == true,
-                     ExcMessage ("The Epetra maps in the assignment operator ="
-                                 " do not match, even though the local_range "
-                                 " seems to be the same. Check vector setup!"));
              ierr = vector->GlobalAssemble (last_action);      
              AssertThrow (ierr == 0, ExcTrilinosError(ierr));