From e53a3fd3f4ce6004998f63d85181072e5df1b204 Mon Sep 17 00:00:00 2001
From: jakobtorben <Jakob.Torben@sintef.no>
Date: Tue, 8 Oct 2024 18:40:40 +0200
Subject: [PATCH] Remove unused well apply

---
 opm/simulators/wells/BlackoilWellModel.hpp    |  3 --
 .../wells/BlackoilWellModel_impl.hpp          | 42 -------------------
 2 files changed, 45 deletions(-)

diff --git a/opm/simulators/wells/BlackoilWellModel.hpp b/opm/simulators/wells/BlackoilWellModel.hpp
index 51615cda5..b8e39b6ea 100644
--- a/opm/simulators/wells/BlackoilWellModel.hpp
+++ b/opm/simulators/wells/BlackoilWellModel.hpp
@@ -297,9 +297,6 @@ template<class Scalar> class WellContributions;
                 return this->computeWellBlockAveragePressures();
             }
 
-            // subtract Binv(D)rw from r;
-            void apply( BVector& r) const;
-
             // subtract B*inv(D)*C * x from A*x
             void apply(const BVector& x, BVector& Ax) const;
 
diff --git a/opm/simulators/wells/BlackoilWellModel_impl.hpp b/opm/simulators/wells/BlackoilWellModel_impl.hpp
index 456b625e7..363b1a7bf 100644
--- a/opm/simulators/wells/BlackoilWellModel_impl.hpp
+++ b/opm/simulators/wells/BlackoilWellModel_impl.hpp
@@ -1734,48 +1734,6 @@ namespace Opm {
 
     }
 
-
-    template<typename TypeTag>
-    void
-    BlackoilWellModel<TypeTag>::
-    apply(BVector& r) const
-    {
-        BVector r_local;
-        for (auto& well : well_container_) {
-            auto cells = well->cells();
-            r_local.resize(cells.size());
-
-            if (this->param_.nonlinear_solver_ == "nldd") {
-                // transfer global cells index to local subdomain cells index
-
-                auto domain_cells = domains_cells_[well_domain_.at(well->name())];
-
-                // Assuming domain_cells is sorted
-                for (size_t i = 0; i < cells.size(); ++i) {
-                    auto it = std::lower_bound(domain_cells.begin(), domain_cells.end(), cells[i]);
-                    if (it != domain_cells.end() && *it == cells[i]) {
-                        // Found the cell, get the index
-                        auto local_index = std::distance(domain_cells.begin(), it);
-                        cells[i] = local_index;
-                    } else {
-                        std::cerr << "Cell value " << cells[i] << " not found in domain_cells." << std::endl;
-                    }
-                }
-            }
-
-            for (size_t i = 0; i < cells.size(); ++i) {
-                r_local[i] = r[cells[i]];
-            }
-
-            well->apply(r_local);
-
-            for (size_t i = 0; i < cells.size(); ++i) {
-                r[cells[i]] = r_local[i];
-            }
-        }
-    }
-
-
     // Ax = A x - C D^-1 B x
     template<typename TypeTag>
     void