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

Removed unused code.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-04-12 15:11:30 +02:00
parent ad298a39d4
commit 346647ca37
1 changed files with 1 additions and 25 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
examples/polymer_reorder.cpp
View File

@ -472,17 +472,6 @@ main(int argc, char** argv)
src[num_cells - 1] = -flow_per_sec; src[num_cells - 1] = -flow_per_sec;
} }
TransportSource* tsrc = create_transport_source(2, 2);
double ssrc[] = { 1.0, 0.0 };
double ssink[] = { 0.0, 1.0 };
double zdummy[] = { 0.0, 0.0 };
for (int cell = 0; cell < num_cells; ++cell) {
if (src[cell] > 0.0) {
append_transport_source(cell, 2, 0, src[cell], ssrc, zdummy, tsrc);
} else if (src[cell] < 0.0) {
append_transport_source(cell, 2, 0, src[cell], ssink, zdummy, tsrc);
}
}
std::vector<double> reorder_src = src; std::vector<double> reorder_src = src;
// Boundary conditions. // Boundary conditions.
@ -643,18 +632,7 @@ main(int argc, char** argv)
ptime += pt; ptime += pt;
// Process transport sources (to include bdy terms). // Process transport sources (to include bdy terms).
if (use_reorder) { Opm::computeTransportSource(*grid->c_grid(), src, state.faceflux(), 1.0, reorder_src);
Opm::computeTransportSource(*grid->c_grid(), src, state.faceflux(), 1.0, reorder_src);
} else {
clear_transport_source(tsrc);
for (int cell = 0; cell < num_cells; ++cell) {
if (src[cell] > 0.0) {
append_transport_source(cell, 2, 0, src[cell], ssrc, zdummy, tsrc);
} else if (src[cell] < 0.0) {
append_transport_source(cell, 2, 0, src[cell], ssink, zdummy, tsrc);
}
}
}
// Find inflow rate. // Find inflow rate.
const double current_time = simtimer.currentTime(); const double current_time = simtimer.currentTime();
@ -765,6 +743,4 @@ main(int argc, char** argv)
outputState(*grid->c_grid(), state, simtimer.currentStepNum(), output_dir); outputState(*grid->c_grid(), state, simtimer.currentStepNum(), output_dir);
outputWaterCut(watercut, output_dir); outputWaterCut(watercut, output_dir);
} }
destroy_transport_source(tsrc);
} }