Reaktoro
A unified framework for modeling chemically reactive systems
OptimumSolverSimplex.hpp
1 // Reaktoro is a unified framework for modeling chemically reactive systems.
2 //
3 // Copyright (C) 2014-2018 Allan Leal
4 //
5 // This library is free software; you can redistribute it and/or
6 // modify it under the terms of the GNU Lesser General Public
7 // License as published by the Free Software Foundation; either
8 // version 2.1 of the License, or (at your option) any later version.
9 //
10 // This library is distributed in the hope that it will be useful,
11 // but WITHOUT ANY WARRANTY; without even the implied warranty of
12 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 // Lesser General Public License for more details.
14 //
15 // You should have received a copy of the GNU Lesser General Public License
16 // along with this library. If not, see <http://www.gnu.org/licenses/>.
17 
18 #pragma once
19 
20 // C++ includes
21 #include <memory>
22 
23 // Reaktoro includes
24 #include <Reaktoro/Optimization/OptimumSolverBase.hpp>
25 
26 namespace Reaktoro {
27 
28 // Forward declarations
29 struct OptimumOptions;
30 struct OptimumProblem;
31 struct OptimumResult;
32 struct OptimumState;
33 
36 {
37 public:
40 
43 
45  virtual ~OptimumSolverSimplex();
46 
49 
54  auto feasible(const OptimumProblem& problem, OptimumState& state, const OptimumOptions& options) -> OptimumResult;
55 
60  auto simplex(const OptimumProblem& problem, OptimumState& state, const OptimumOptions& options) -> OptimumResult;
61 
65  virtual auto solve(const OptimumProblem& problem, OptimumState& state) -> OptimumResult;
66 
71  virtual auto solve(const OptimumProblem& problem, OptimumState& state, const OptimumOptions& options) -> OptimumResult;
72 
76  virtual auto dxdp(VectorConstRef dgdp, VectorConstRef dbdp) -> Vector;
77 
79  virtual auto clone() const -> OptimumSolverBase*;
80 
81 private:
82  struct Impl;
83 
84  std::unique_ptr<Impl> pimpl;
85 };
86 
87 } // namespace Reaktoro
The class that implements the simplex algorithm for linear programming problems.
Definition: OptimumSolverSimplex.hpp:36
auto feasible(const OptimumProblem &problem, OptimumState &state, const OptimumOptions &options) -> OptimumResult
Find a feasible point for the linear optimisation problem.
Definition: OptimumSolverSimplex.cpp:478
The base class for all optimization algorithms.
Definition: OptimumSolverBase.hpp:36
virtual auto clone() const -> OptimumSolverBase *
Return a clone of this instance.
Definition: OptimumSolverSimplex.cpp:506
OptimumSolverSimplex()
Construct a default OptimumSolverSimplex instance.
Definition: OptimumSolverSimplex.cpp:461
A type that describes the non-linear constrained optimisation problem.
Definition: OptimumProblem.hpp:49
The namespace containing all components of the Reaktoro library.
Definition: ChemicalScalar.hpp:24
auto operator=(OptimumSolverSimplex other) -> OptimumSolverSimplex &
Assign a copy of an OptimumSolverSimplex instance.
Definition: OptimumSolverSimplex.cpp:472
virtual auto dxdp(VectorConstRef dgdp, VectorConstRef dbdp) -> Vector
Return the sensitivity dx/dp of the solution x with respect to a vector of parameters p.
Definition: OptimumSolverSimplex.cpp:499
A type that describes the options of a optimisation calculation.
Definition: OptimumOptions.hpp:187
virtual ~OptimumSolverSimplex()
Destroy this OptimumSolverSimplex instance.
Definition: OptimumSolverSimplex.cpp:469
A type that describes the state of an optimum solution.
Definition: OptimumState.hpp:29
A type that describes the result of an optimisation calculation.
Definition: OptimumResult.hpp:24
auto simplex(const OptimumProblem &problem, OptimumState &state, const OptimumOptions &options) -> OptimumResult
Solve the linear optimisation problem with starting from a feasible point.
Definition: OptimumSolverSimplex.cpp:483
virtual auto solve(const OptimumProblem &problem, OptimumState &state) -> OptimumResult
Solve the linear optimisation problem by finding a feasible point and then applying a simplex algorit...
Definition: OptimumSolverSimplex.cpp:488
Eigen::Ref< const Eigen::VectorXd > VectorConstRef
< Alias to Eigen type Ref<VectorXd>.
Definition: Matrix.hpp:31