From 9c1c87a04051c4e1702bbcf146700f74c8cb7848 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= Date: Fri, 9 Oct 2015 11:10:46 +0200 Subject: [PATCH] Remember to initialize the data members. --- opm/autodiff/WellStateMultiSegment.hpp | 20 ++++++++++++++++++++ 1 file changed, 20 insertions(+) diff --git a/opm/autodiff/WellStateMultiSegment.hpp b/opm/autodiff/WellStateMultiSegment.hpp index 220235082..5ce218bb5 100644 --- a/opm/autodiff/WellStateMultiSegment.hpp +++ b/opm/autodiff/WellStateMultiSegment.hpp @@ -84,16 +84,36 @@ namespace Opm nseg += wells[iw]->numberOfSegments(); } + bhp().resize(nw); + thp().resize(nw); top_segment_loc_.resize(nw); + temperature().resize(nw, 273.15 + 20); // standard temperature for now // deciding to add the following variables temporarily // TODO: making it better later nseg_ = nseg; nperf_ = nperf; + wellRates().resize(nw * np, 0.0); + + currentControls().resize(nw); + for(int iw = 0; iw < nw; ++iw) { + currentControls()[iw] = well_controls_get_current(wells[iw]->wellControls()); + } + + for (int iw = 0; iw < nw; ++iw) { + assert((wells[iw]->wellType() == INJECTOR) || (wells[iw]->wellType() == PRODUCER)); + } + int start_segment = 0; int start_perforation = 0; + perfPhaseRates().clear(); + perfPhaseRates().resize(nperf * np, 0.0); + + perfPress().clear(); + perfPress().resize(nperf, -1.0e100); + segphaserates_.clear(); segphaserates_.resize(nseg * np, 0.0);