OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Remove unused well apply

Browse Source
This commit is contained in:
jakobtorben 2024-10-08 18:40:40 +02:00
parent 08bbbbf9fb
commit e53a3fd3f4
2 changed files with 0 additions and 45 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
opm/simulators/wells/BlackoilWellModel.hpp
View File

@ -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;

42
opm/simulators/wells/BlackoilWellModel_impl.hpp
View File

@ -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