From 1d9b0f2530fcf65ed202cc9f12043a199d06fa13 Mon Sep 17 00:00:00 2001
From: DavidSCN <dav.schneider@tum.de>
Date: Wed, 3 Feb 2021 13:49:27 +0100
Subject: [PATCH] Add an experimental cell loop

---
 include/mf_nh_operator.h | 23 ++++++++++++++---------
 1 file changed, 14 insertions(+), 9 deletions(-)

diff --git a/include/mf_nh_operator.h b/include/mf_nh_operator.h
index 11b19c7..3196791 100644
--- 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)