diff --git a/CMakeLists_files.cmake b/CMakeLists_files.cmake index d98cdb999..5be7b04b9 100644 --- a/CMakeLists_files.cmake +++ b/CMakeLists_files.cmake @@ -107,11 +107,11 @@ list (APPEND PUBLIC_HEADER_FILES opm/autodiff/ImpesTPFAAD.hpp opm/autodiff/FullyImplicitBlackoilSolver.hpp opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp - opm/autodiff/FullyImplicitSolver.hpp - opm/autodiff/FullyImplicitSolver_impl.hpp opm/autodiff/NewtonIterationBlackoilCPR.hpp opm/autodiff/NewtonIterationBlackoilInterface.hpp opm/autodiff/NewtonIterationBlackoilSimple.hpp + opm/autodiff/NewtonSolver.hpp + opm/autodiff/NewtonSolver_impl.hpp opm/autodiff/LinearisedBlackoilResidual.hpp opm/autodiff/RateConverter.hpp opm/autodiff/RedistributeDataHandles.hpp diff --git a/opm/autodiff/FullyImplicitSolver.hpp b/opm/autodiff/NewtonSolver.hpp similarity index 90% rename from opm/autodiff/FullyImplicitSolver.hpp rename to opm/autodiff/NewtonSolver.hpp index cbd15487a..c2c4f5144 100644 --- a/opm/autodiff/FullyImplicitSolver.hpp +++ b/opm/autodiff/NewtonSolver.hpp @@ -18,8 +18,8 @@ along with OPM. If not, see . */ -#ifndef OPM_FULLYIMPLICITSOLVER_HEADER_INCLUDED -#define OPM_FULLYIMPLICITSOLVER_HEADER_INCLUDED +#ifndef OPM_NEWTONSOLVER_HEADER_INCLUDED +#define OPM_NEWTONSOLVER_HEADER_INCLUDED #include #include @@ -27,9 +27,9 @@ namespace Opm { - /// A fully implicit solver class suitable for general models. + /// A Newton solver class suitable for general fully-implicit models. template - class FullyImplicitSolver + class NewtonSolver { public: // --------- Types and enums --------- @@ -65,8 +65,8 @@ namespace Opm { /// Construct solver for a given model. /// \param[in] param parameters controlling nonlinear Newton process /// \param[in, out] model physical simulation model - explicit FullyImplicitSolver(const SolverParameter& param, - PhysicalModel& model); + explicit NewtonSolver(const SolverParameter& param, + PhysicalModel& model); /// Take a single forward step, after which the states will be modified /// according to the physical model. @@ -106,6 +106,6 @@ namespace Opm { }; } // namespace Opm -#include "FullyImplicitSolver_impl.hpp" +#include "NewtonSolver_impl.hpp" -#endif // OPM_FULLYIMPLICITSOLVER_HEADER_INCLUDED +#endif // OPM_NEWTONSOLVER_HEADER_INCLUDED diff --git a/opm/autodiff/FullyImplicitSolver_impl.hpp b/opm/autodiff/NewtonSolver_impl.hpp similarity index 85% rename from opm/autodiff/FullyImplicitSolver_impl.hpp rename to opm/autodiff/NewtonSolver_impl.hpp index bc001e820..64bb9c246 100644 --- a/opm/autodiff/FullyImplicitSolver_impl.hpp +++ b/opm/autodiff/NewtonSolver_impl.hpp @@ -20,16 +20,16 @@ along with OPM. If not, see . */ -#ifndef OPM_FULLYIMPLICITSOLVER_IMPL_HEADER_INCLUDED -#define OPM_FULLYIMPLICITSOLVER_IMPL_HEADER_INCLUDED +#ifndef OPM_NEWTONSOLVER_IMPL_HEADER_INCLUDED +#define OPM_NEWTONSOLVER_IMPL_HEADER_INCLUDED -#include +#include namespace Opm { template - FullyImplicitSolver::FullyImplicitSolver(const SolverParameter& param, - PhysicalModel& model) + NewtonSolver::NewtonSolver(const SolverParameter& param, + PhysicalModel& model) : param_(param), model_(model), newtonIterations_(0), @@ -38,13 +38,13 @@ namespace Opm } template - unsigned int FullyImplicitSolver::newtonIterations () const + unsigned int NewtonSolver::newtonIterations () const { return newtonIterations_; } template - unsigned int FullyImplicitSolver::linearIterations () const + unsigned int NewtonSolver::linearIterations () const { return linearIterations_; } @@ -52,7 +52,7 @@ namespace Opm template int - FullyImplicitSolver:: + NewtonSolver:: step(const double dt, ReservoirState& reservoir_state, WellState& well_state) @@ -128,7 +128,7 @@ namespace Opm template - void FullyImplicitSolver::SolverParameter:: + void NewtonSolver::SolverParameter:: reset() { // default values for the solver parameters @@ -141,7 +141,7 @@ namespace Opm } template - FullyImplicitSolver::SolverParameter:: + NewtonSolver::SolverParameter:: SolverParameter() { // set default values @@ -149,7 +149,7 @@ namespace Opm } template - FullyImplicitSolver::SolverParameter:: + NewtonSolver::SolverParameter:: SolverParameter( const parameter::ParameterGroup& param ) { // set default values @@ -172,9 +172,9 @@ namespace Opm template void - FullyImplicitSolver::detectNewtonOscillations(const std::vector>& residual_history, - const int it, const double relaxRelTol, - bool& oscillate, bool& stagnate) const + NewtonSolver::detectNewtonOscillations(const std::vector>& residual_history, + const int it, const double relaxRelTol, + bool& oscillate, bool& stagnate) const { // The detection of oscillation in two primary variable results in the report of the detection // of oscillation for the solver. @@ -209,8 +209,8 @@ namespace Opm template void - FullyImplicitSolver::stabilizeNewton(V& dx, V& dxOld, const double omega, - const RelaxType relax_type) const + NewtonSolver::stabilizeNewton(V& dx, V& dxOld, const double omega, + const RelaxType relax_type) const { // The dxOld is updated with dx. // If omega is equal to 1., no relaxtion will be appiled. diff --git a/opm/autodiff/SimulatorFullyImplicitBlackoil_impl.hpp b/opm/autodiff/SimulatorFullyImplicitBlackoil_impl.hpp index a860e5612..d25d0ed5f 100644 --- a/opm/autodiff/SimulatorFullyImplicitBlackoil_impl.hpp +++ b/opm/autodiff/SimulatorFullyImplicitBlackoil_impl.hpp @@ -24,7 +24,7 @@ #include #include -#include +#include #include #include #include @@ -236,7 +236,7 @@ namespace Opm typedef BlackoilModel Model; typedef typename Model::ModelParameter ModelParam; ModelParam modelParam( param_ ); - typedef FullyImplicitSolver Solver; + typedef NewtonSolver Solver; typedef typename Solver::SolverParameter SolverParam; SolverParam solverParam( param_ );