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

WIP debugging.

Browse Source
This commit is contained in:
Kai Bao 2015-10-05 17:22:46 +02:00
parent 0fee0059a6
commit ca5dd0dca4
3 changed files with 55 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

41
opm/autodiff/BlackoilMultiSegmentModel_impl.hpp
View File

@ -913,8 +913,11 @@ namespace Opm {
break;
}
}
std::cout << " xw.bhp() is " << xw.bhp()[w] << std::endl;
std::cout << " ctrl_index " << ctrl_index << " nwc " << nwc << std::endl;
if (int debug = 0) {
std::cout << " xw.bhp() is " << xw.bhp()[w] << std::endl;
std::cout << " ctrl_index " << ctrl_index << " nwc " << nwc << std::endl;
}
if (ctrl_index != nwc) {
// Constraint number ctrl_index was broken, switch to it.
@ -1116,6 +1119,38 @@ namespace Opm {
start_segment += nseg;
}
if (int debug = 0) {
std::cout << "output bhp_targets " << std::endl;
for (int i = 0; i < nw; ++i) {
std::cout << i << " " << bhp_targets(i) << std::endl;
}
std::cout << " bhp_well_elems " << std::endl;
std::cout << " the size of bhp_well_elems is " << bhp_well_elems.size() << std::endl;
for (int i = 0; i < bhp_well_elems.size(); ++i) {
std::cout << " bhp_well_elems " << i << " is " << bhp_well_elems[i] << std::endl;
}
std::cout << " rate_well_elems " << std::endl;
std::cout << " the size of rate_well_elems is " << rate_well_elems.size() << std::endl;
for (int i = 0; i < rate_well_elems.size(); ++i) {
std::cout << " rate_well_elems " << i << " is " << rate_well_elems[i] << std::endl;
}
std::cout << " bhp_top_elems " << std::endl;
std::cout << " the size of bhp_top_elems is " << bhp_top_elems.size() << std::endl;
for (int i = 0; i < bhp_top_elems.size(); ++i) {
std::cout << " bhp_top_elems " << i << " is " << bhp_top_elems[i] << std::endl;
}
std::cout << " rate_top_elems " << std::endl;
std::cout << " the size of the rate_top_elems " << rate_top_elems.size() << std::endl;
for (int i = 0; i < rate_top_elems.size(); ++i) {
std::cout << " rate_top_elems " << i << " is " << rate_top_elems[i] << std::endl;
}
std::cout << " the others_elems " << std::endl;
std::cout << " the size of others_elems is " << others_elems.size() << std::endl;
for(int i = 0; i < others_elems.size(); ++i) {
std::cout << "others_elems " << i << " is " << others_elems[i] << std::endl;
}
}
// for each segment: 1, if the segment is the top segment, then control equation
// 2, if the segment is not the top segment, then the pressure equation
@ -1476,7 +1511,7 @@ namespace Opm {
// TODO: handling the THP control related.
}
if (int debug = 1) {
if (int debug = 0) {
std::cout << " output all the well state informations after updateWellState()" << std::endl;
const int np = well_state.numberOfPhases();
const int nw = well_state.numberOfWells();

2
opm/autodiff/SimulatorFullyImplicitBlackoilMultiSegment_impl.hpp
View File

@ -119,7 +119,7 @@ namespace Opm
}
// for DEBUGGING OUTPUT
if (int debug = 1) {
if (int debug = 0) {
std::cout << " the number of the wells from EclipseState " << wells_ecl.size() << std::endl;
for (size_t i = 0; i < wells_ecl.size(); ++i) {
std::cout << " well name " << wells_ecl[i]->name() << std::endl;

32
opm/autodiff/WellStateMultiSegment.hpp
View File

@ -361,7 +361,7 @@ namespace Opm
}
}
if (int debug = 1) {
if (int debug = 0) {
std::cout << " output all the well state informations after initialization " << std::endl;
const int np = numberOfPhases();
const int nw = numberOfWells();
@ -408,25 +408,25 @@ namespace Opm
}
std::cout << std::endl;
}
}
std::cout << " output all the information from the wellMap " << std::endl;
std::cout << " output all the information from the wellMap " << std::endl;
for (WellMapType::const_iterator iter = wellMap().begin(); iter != wellMap().end(); ++iter) {
std::cout << " well name : " << iter->first << std::endl;
const MapentryType &wellmapInfo = iter->second;
std::cout << " well number : " << wellmapInfo.well_number << " start segment " << wellmapInfo.start_segment
<< " number of segment : " << wellmapInfo.number_of_segments << std::endl;
std::cout << " start perforation : " << wellmapInfo.start_perforation << " number of perforations : " << wellmapInfo.number_of_perforations << std::endl;
const int nseg = wellmapInfo.number_of_segments;
std::cout << " start performation ofr each segment and number of perforation that each segment has" << std::endl;
for (int i = 0; i < nseg; ++i) {
std::cout << " segment " << i << " start perforation " << wellmapInfo.start_perforation_segment[i]
<< " number of perforations " << wellmapInfo.number_of_perforations_segment[i] << std::endl;
for (WellMapType::const_iterator iter = wellMap().begin(); iter != wellMap().end(); ++iter) {
std::cout << " well name : " << iter->first << std::endl;
const MapentryType &wellmapInfo = iter->second;
std::cout << " well number : " << wellmapInfo.well_number << " start segment " << wellmapInfo.start_segment
<< " number of segment : " << wellmapInfo.number_of_segments << std::endl;
std::cout << " start perforation : " << wellmapInfo.start_perforation << " number of perforations : " << wellmapInfo.number_of_perforations << std::endl;
const int nseg = wellmapInfo.number_of_segments;
std::cout << " start performation ofr each segment and number of perforation that each segment has" << std::endl;
for (int i = 0; i < nseg; ++i) {
std::cout << " segment " << i << " start perforation " << wellmapInfo.start_perforation_segment[i]
<< " number of perforations " << wellmapInfo.number_of_perforations_segment[i] << std::endl;
}
}
std::cout << " output the well state right after intialization is DONE! " << std::endl;
std::cin.ignore();
}
std::cout << " output the well state right after intialization is DONE! " << std::endl;
std::cin.ignore();
}
std::vector<double>& segPhaseRates() { return segphaserates_; }