From 607b939208eeeaf0397a855d8ef730b7fa8721d2 Mon Sep 17 00:00:00 2001
From: Max Luebke <mluebke@gfz-potsdam.de>
Date: Thu, 13 Apr 2023 17:26:18 +0200
Subject: [PATCH] fix: remove unnecessary RoundToZero function from
 DiffusionModule

---
 include/poet/DiffusionModule.hpp |  2 --
 src/DiffusionModule.cpp          | 11 -----------
 2 files changed, 13 deletions(-)

diff --git a/include/poet/DiffusionModule.hpp b/include/poet/DiffusionModule.hpp
index 52d0199d0..63c514189 100644
--- a/include/poet/DiffusionModule.hpp
+++ b/include/poet/DiffusionModule.hpp
@@ -86,8 +86,6 @@ private:
 
   void initialize(poet::DiffusionParams args);
 
-  void RoundToZero(double *field, uint32_t cell_count) const;
-
   Grid &grid;
   uint8_t dim;
 
diff --git a/src/DiffusionModule.cpp b/src/DiffusionModule.cpp
index b5a60ac25..bca4728d1 100644
--- a/src/DiffusionModule.cpp
+++ b/src/DiffusionModule.cpp
@@ -182,11 +182,6 @@ void DiffusionModule::simulate(double dt) {
     } else {
       tug::diffusion::ADI_2D(this->diff_input, in_field, in_alpha.data());
     }
-
-    // TODO: do not use hardcoded index for O, H and charge
-    if (i > 2) {
-      this->RoundToZero(in_field, this->n_cells_per_prop);
-    }
   }
 
   std::cout << " done!\n";
@@ -195,12 +190,6 @@ void DiffusionModule::simulate(double dt) {
 
   transport_t += sim_a_transport - sim_b_transport;
 }
-inline void DiffusionModule::RoundToZero(double *field,
-                                         uint32_t cell_count) const {
-  for (uint32_t i = 0; i < cell_count; i++) {
-    field[i] = ((int32_t)(field[i] / ZERO_MULTIPLICATOR)) * ZERO_MULTIPLICATOR;
-  }
-}
 
 void DiffusionModule::end() {
   // R["simtime_transport"] = transport_t;