]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Make tests work again
authorPeter Munch <peterrmuench@gmail.com>
Mon, 5 Apr 2021 05:58:49 +0000 (07:58 +0200)
committerPeter Munch <peterrmuench@gmail.com>
Mon, 5 Apr 2021 05:58:49 +0000 (07:58 +0200)
include/deal.II/multigrid/mg_transfer_global_coarsening.templates.h

index af7b5e691c988cd8413604c2b11f09978b4706bb..5e1cc25cca93de090e3af553fdd7a959bc17710a 100644 (file)
@@ -2184,7 +2184,7 @@ MGTwoLevelTransfer<dim, LinearAlgebra::distributed::Vector<Number>>::prolongate(
                     weights[i];
               else
                 for (unsigned int i = 0; i < scheme.dofs_per_cell_fine; ++i)
-                  dst.local_element(indices_fine[i]) =
+                  this->vec_fine.local_element(indices_fine[i]) =
                     this->vec_coarse.local_element(indices_coarse[i]);
 
               indices_fine += scheme.dofs_per_cell_fine;
@@ -2264,7 +2264,8 @@ MGTwoLevelTransfer<dim, LinearAlgebra::distributed::Vector<Number>>::prolongate(
                       evaluation_data_fine[i][v] * weights[i];
                 else
                   for (unsigned int i = 0; i < scheme.dofs_per_cell_fine; ++i)
-                    dst.local_element(indices[i]) = evaluation_data_fine[i][v];
+                    this->vec_fine.local_element(indices[i]) =
+                      evaluation_data_fine[i][v];
 
                 indices += scheme.dofs_per_cell_fine;
 
@@ -2279,10 +2280,10 @@ MGTwoLevelTransfer<dim, LinearAlgebra::distributed::Vector<Number>>::prolongate(
                                             // in do_restrict_add does not work
 
   if (schemes.size() > 0 && schemes.front().fine_element_is_continuous)
-    {
-      this->vec_fine.compress(VectorOperation::add);
-      dst.copy_locally_owned_data_from(this->vec_fine);
-    }
+    this->vec_fine.compress(VectorOperation::add);
+
+  if (true || schemes.size() > 0 && schemes.front().fine_element_is_continuous)
+    dst.copy_locally_owned_data_from(this->vec_fine);
 }
 
 
@@ -2352,7 +2353,8 @@ MGTwoLevelTransfer<dim, LinearAlgebra::distributed::Vector<Number>>::
               else
                 for (unsigned int i = 0; i < scheme.dofs_per_cell_fine; ++i)
                   distribute_local_to_global(indices_coarse[i],
-                                             src.local_element(indices_fine[i]),
+                                             this->vec_fine.local_element(
+                                               indices_fine[i]),
                                              this->vec_coarse);
 
               indices_fine += scheme.dofs_per_cell_fine;
@@ -2401,7 +2403,8 @@ MGTwoLevelTransfer<dim, LinearAlgebra::distributed::Vector<Number>>::
                       this->vec_fine.local_element(indices[i]) * weights[i];
                 else
                   for (unsigned int i = 0; i < scheme.dofs_per_cell_fine; ++i)
-                    evaluation_data_fine[i][v] = src.local_element(indices[i]);
+                    evaluation_data_fine[i][v] =
+                      this->vec_fine.local_element(indices[i]);
 
                 indices += scheme.dofs_per_cell_fine;
 

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.