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

Added more (inactive) debug output.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-02-08 13:43:22 +01:00
parent e998753236
commit 46f29380b3
1 changed files with 8 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
opm/polymer/polymermodel.cpp
View File

@ -176,7 +176,7 @@ residual_c(double c, void *data)
+ rhor*((1.0 - porosity)/porosity)*(ads - ads0) + rhor*((1.0 - porosity)/porosity)*(ads - ads0)
+ p->dtpv*(p->outflux*ff*mc + p->influx_polymer); + p->dtpv*(p->outflux*ff*mc + p->influx_polymer);
#ifdef EXTRA_DEBUG_OUTPUT #ifdef EXTRA_DEBUG_OUTPUT
std::cout << "res(c) = " << res << std::endl; std::cout << "c = " << c << " s = " << s << " c-residual = " << res << std::endl;
#endif #endif
return res; return res;
} }
@ -241,7 +241,7 @@ get_parameters_c(struct PolymerSolverData *d, int cell, NonlinearSolverCtrl* ctr
double src = d->source[cell]; double src = d->source[cell];
p.influx = src > 0 ? -src : 0.0; p.influx = src > 0 ? -src : 0.0;
p.influx_polymer = src > 0 ? -src*compute_mc(d->inflow_c, d->props, d->polydata) : 0.0; p.influx_polymer = src > 0 ? -src*compute_mc(d->inflow_c, d->props, d->polydata) : 0.0;
p.outflux = d->source[cell] <= 0 ? -d->source[cell] : 0.0; p.outflux = src <= 0 ? -src : 0.0;
p.porosity = d->porosity[cell]; p.porosity = d->porosity[cell];
p.psdata = d; p.psdata = d;
p.cell = cell; p.cell = cell;
@ -272,8 +272,12 @@ get_parameters_c(struct PolymerSolverData *d, int cell, NonlinearSolverCtrl* ctr
} }
} }
#ifdef EXTRA_DEBUG_OUTPUT #ifdef EXTRA_DEBUG_OUTPUT
std::cout << "in: " << p.influx << " in(polymer): " << p.influx_polymer std::cout << "Cell: " << cell
<< " out: " << p.outflux << std::endl; << " in: " << p.influx
<< " in(polymer): " << p.influx_polymer
<< " out: " << p.outflux << '\n'
<< " c0: " << p.c0
<< " s0: " << p.s0 << std::endl;
#endif #endif
return p; return p;
} }