diff --git a/opm/polymer/SimulatorCompressiblePolymer.cpp b/opm/polymer/SimulatorCompressiblePolymer.cpp index 59fd08c9b..544543d78 100644 --- a/opm/polymer/SimulatorCompressiblePolymer.cpp +++ b/opm/polymer/SimulatorCompressiblePolymer.cpp @@ -126,8 +126,6 @@ namespace Opm WellsManager& wells_manager_; const Wells* wells_; const PolymerInflowInterface& polymer_inflow_; - const std::vector& src_; - const FlowBoundaryConditions* bcs_; const double* gravity_; // Solvers CompressibleTpfaPolymer psolver_; @@ -190,8 +188,6 @@ namespace Opm wells_manager_(wells_manager), wells_(wells_manager.c_wells()), polymer_inflow_(polymer_inflow), - src_(src), - bcs_(bcs), gravity_(gravity), psolver_(grid, props, rock_comp_props, poly_props, linsolver, param.getDefault("nl_pressure_residual_tolerance", 0.0), diff --git a/opm/polymer/SimulatorPolymer.cpp b/opm/polymer/SimulatorPolymer.cpp index 547c97189..827b80ddd 100644 --- a/opm/polymer/SimulatorPolymer.cpp +++ b/opm/polymer/SimulatorPolymer.cpp @@ -203,7 +203,6 @@ namespace Opm polymer_inflow_(polymer_inflow), src_(src), bcs_(bcs), - gravity_(gravity), psolver_(grid, props, rock_comp_props, poly_props, linsolver, param.getDefault("nl_pressure_residual_tolerance", 0.0), param.getDefault("nl_pressure_change_tolerance", 1.0), diff --git a/opm/polymer/TransportSolverTwophasePolymer.cpp b/opm/polymer/TransportSolverTwophasePolymer.cpp index e6179bf07..767a7fafd 100644 --- a/opm/polymer/TransportSolverTwophasePolymer.cpp +++ b/opm/polymer/TransportSolverTwophasePolymer.cpp @@ -929,7 +929,7 @@ namespace Opm tmp_x[1]=x[1]; } res_eq.computeJacobiRes(tmp_x, dres_s_dsdc, dres_c_dsdc); - double dFx_dx,dFx_dy,dFy_dx,dFy_dy; + double dFx_dx(0),dFx_dy(0),dFy_dx(0),dFy_dy(0); double det = dFx_dx*dFy_dy - dFy_dx*dFx_dy; if(use_sc){ dFx_dx=(dres_s_dsdc[0]-tmp_x[1]*dres_s_dsdc[1]/tmp_x[0]);