Add an experimental cell loop

davidscn · Feb 3, 2021 · 1d9b0f2 · 1d9b0f2
1 parent 15fee79
commit 1d9b0f2
Showing 1 changed file with 14 additions and 9 deletions.
diff --git a/include/mf_nh_operator.h b/include/mf_nh_operator.h
@@ -672,19 +672,24 @@ NeoHookOperator<dim, fe_degree, n_q_points_1d, Number>::vmult_add(
   // https://www.dealii.org/developer/doxygen/deal.II/matrix__free_8h_source.html#l00109
 
   // 1. make sure ghosts are updated
-  if (mask & MFMask::MPI)
-    src.update_ghost_values();
+  //  if (mask & MFMask::MPI)
+  //    src.update_ghost_values();
+
+  data_current->cell_loop(&NeoHookOperator::local_apply_cell<mask>,
+                          this,
+                          dst,
+                          src);
 
   // 2. loop over all locally owned cell blocks
-  local_apply_cell<mask>(*data_current,
-                         dst,
-                         src,
-                         std::make_pair<unsigned int, unsigned int>(
-                           0, data_current->n_cell_batches()));
+  //  local_apply_cell<mask>(*data_current,
+  //                         dst,
+  //                         src,
+  //                         std::make_pair<unsigned int, unsigned int>(
+  //                           0, data_current->n_cell_batches()));
 
   // 3. communicate results with MPI
-  if (mask & MFMask::MPI)
-    dst.compress(VectorOperation::add);
+  //  if (mask & MFMask::MPI)
+  //    dst.compress(VectorOperation::add);
 
   // 4. constraints
   if (mask & MFMask::RW)