71 const ReservoirProperties& rp,
72 const std::vector<double>& orig_sat,
73 const std::vector<double>& cap_press)
76 cap_press_(cap_press),
79 std::pair<double, double> vols = poreSatVolumes(grid, rp, orig_sat);
80 orig_satvol_ = vols.second;
81 int num_cells = orig_sat.size();
82 cp_new_.resize(num_cells);
83 sat_.resize(num_cells);
87 double operator()(
double dp)
const
89 std::transform(cap_press_.begin(), cap_press_.end(), cp_new_.begin(),
90 [dp](
const double& input) { return input + dp; });
92 std::pair<double, double> vols = poreSatVolumes(grid_, rp_, sat_);
93 return (vols.second - orig_satvol_)/vols.first;
96 const std::vector<double>& lastSaturations()
const
102 void computeSaturations()
const
104 int num_cells = grid_.numberOfCells();
105 for (
int c = 0; c < num_cells; ++c) {
106 sat_[c] = rp_.saturationFromCapillaryPressure(c, cp_new_[c]);
110 const GridInterface& grid_;
111 const ReservoirProperties& rp_;
112 const std::vector<double>& cap_press_;
114 mutable std::vector<double> cp_new_;
115 mutable std::vector<double> sat_;