-
Notifications
You must be signed in to change notification settings - Fork 28
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Add identity parameter-parameter dependence
Adds a parameter dependence that simply outputs the given parameter value.
- Loading branch information
1 parent
c3586c1
commit 82a0f84
Showing
3 changed files
with
146 additions
and
4 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
141 changes: 141 additions & 0 deletions
141
src/libcadet/model/paramdep/IdentityParameterDependence.cpp
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,141 @@ | ||
| // ============================================================================= | ||
| // CADET | ||
| // | ||
| // Copyright ┬® 2008-2022: The CADET Authors | ||
| // Please see the AUTHORS and CONTRIBUTORS file. | ||
| // | ||
| // All rights reserved. This program and the accompanying materials | ||
| // are made available under the terms of the GNU Public License v3.0 (or, at | ||
| // your option, any later version) which accompanies this distribution, and | ||
| // is available at http://www.gnu.org/licenses/gpl.html | ||
| // ============================================================================= | ||
|
|
||
| #include "model/paramdep/ParameterDependenceBase.hpp" | ||
| #include "ParamReaderHelper.hpp" | ||
| #include "cadet/Exceptions.hpp" | ||
| #include "SimulationTypes.hpp" | ||
| #include "cadet/ParameterId.hpp" | ||
|
|
||
| #include <sstream> | ||
| #include <vector> | ||
| #include <iomanip> | ||
|
|
||
| namespace cadet | ||
| { | ||
|
|
||
| namespace model | ||
| { | ||
|
|
||
| /** | ||
| * @brief Defines an identity parameter dependence | ||
| */ | ||
| class IdentityParameterStateDependence : public ParameterStateDependenceBase | ||
| { | ||
| public: | ||
|
|
||
| IdentityParameterStateDependence() { } | ||
| virtual ~IdentityParameterStateDependence() CADET_NOEXCEPT { } | ||
|
|
||
| static const char* identifier() { return "IDENTITY"; } | ||
| virtual const char* name() const CADET_NOEXCEPT { return IdentityParameterStateDependence::identifier(); } | ||
|
|
||
| virtual int jacobianElementsPerRowLiquid() const CADET_NOEXCEPT { return 0; } | ||
| virtual int jacobianElementsPerRowCombined() const CADET_NOEXCEPT { return 0; } | ||
|
|
||
| virtual void analyticJacobianLiquidAdd(const ColumnPosition& colPos, double param, double const* y, int comp, double factor, int offset, int row, linalg::DoubleSparseMatrix& jac) const { } | ||
| virtual void analyticJacobianCombinedAddLiquid(const ColumnPosition& colPos, double param, double const* yLiquid, double const* ySolid, int comp, double factor, int offset, int row, linalg::DoubleSparseMatrix& jac) const { } | ||
| virtual void analyticJacobianCombinedAddSolid(const ColumnPosition& colPos, double param, double const* yLiquid, double const* ySolid, int bnd, double factor, int offset, int row, linalg::DoubleSparseMatrix& jac) const { } | ||
|
|
||
| CADET_PARAMETERSTATEDEPENDENCE_BOILERPLATE | ||
|
|
||
| protected: | ||
|
|
||
| virtual bool configureImpl(IParameterProvider& paramProvider, UnitOpIdx unitOpIdx, ParticleTypeIdx parTypeIdx, const std::string& name) | ||
| { | ||
| return true; | ||
| } | ||
|
|
||
| template <typename StateType, typename ParamType> | ||
| typename DoubleActivePromoter<StateType, ParamType>::type liquidParameterImpl(const ColumnPosition& colPos, const ParamType& param, StateType const* y, int comp) const | ||
| { | ||
| return param; | ||
| } | ||
|
|
||
| template <typename RowIterator> | ||
| void analyticJacobianLiquidAddImpl(const ColumnPosition& colPos, double param, double const* y, int comp, double factor, int offset, RowIterator jac) const { } | ||
|
|
||
| template <typename StateType, typename ParamType> | ||
| typename DoubleActivePromoter<StateType, ParamType>::type combinedParameterLiquidImpl(const ColumnPosition& colPos, const ParamType& param, StateType const* yLiquid, StateType const* ySolid, int comp) const | ||
| { | ||
| return param; | ||
| } | ||
|
|
||
| template <typename RowIterator> | ||
| void analyticJacobianCombinedAddLiquidImpl(const ColumnPosition& colPos, double param, double const* yLiquid, double const* ySolid, int comp, double factor, int offset, RowIterator jac) const { } | ||
|
|
||
| template <typename StateType, typename ParamType> | ||
| typename DoubleActivePromoter<StateType, ParamType>::type combinedParameterSolidImpl(const ColumnPosition& colPos, const ParamType& param, StateType const* yLiquid, StateType const* ySolid, int bnd) const | ||
| { | ||
| return param; | ||
| } | ||
|
|
||
| template <typename RowIterator> | ||
| void analyticJacobianCombinedAddSolidImpl(const ColumnPosition& colPos, double param, double const* yLiquid, double const* ySolid, int bnd, double factor, int offset, RowIterator jac) const { } | ||
| }; | ||
|
|
||
|
|
||
| /** | ||
| * @brief Defines an identity parameter dependence | ||
| */ | ||
| class IdentityParameterParameterDependence : public ParameterParameterDependenceBase | ||
| { | ||
| public: | ||
|
|
||
| IdentityParameterParameterDependence() { } | ||
| virtual ~IdentityParameterParameterDependence() CADET_NOEXCEPT { } | ||
|
|
||
| static const char* identifier() { return "IDENTITY"; } | ||
| virtual const char* name() const CADET_NOEXCEPT { return IdentityParameterParameterDependence::identifier(); } | ||
|
|
||
| CADET_PARAMETERPARAMETERDEPENDENCE_BOILERPLATE | ||
|
|
||
| protected: | ||
|
|
||
| virtual bool configureImpl(IParameterProvider& paramProvider, UnitOpIdx unitOpIdx, ParticleTypeIdx parTypeIdx, BoundStateIdx bndIdx, const std::string& name) | ||
| { | ||
| return true; | ||
| } | ||
|
|
||
| template <typename ParamType> | ||
| ParamType getValueImpl(const IModel& model, const ColumnPosition& colPos, int comp, int parType, int bnd) const | ||
| { | ||
| return 0.0; | ||
| } | ||
|
|
||
| template <typename ParamType> | ||
| ParamType getValueImpl(const IModel& model, const ColumnPosition& colPos, int comp, int parType, int bnd, const ParamType& val) const | ||
| { | ||
| return val; | ||
| } | ||
|
|
||
| }; | ||
|
|
||
|
|
||
| namespace paramdep | ||
| { | ||
| void registerIdentityParamDependence(std::unordered_map<std::string, std::function<model::IParameterStateDependence*()>>& paramDeps) | ||
| { | ||
| paramDeps[IdentityParameterStateDependence::identifier()] = []() { return new IdentityParameterStateDependence(); }; | ||
| paramDeps["NONE"] = []() { return new IdentityParameterStateDependence(); }; | ||
| } | ||
|
|
||
| void registerIdentityParamDependence(std::unordered_map<std::string, std::function<model::IParameterParameterDependence*()>>& paramDeps) | ||
| { | ||
| paramDeps[IdentityParameterParameterDependence::identifier()] = []() { return new IdentityParameterParameterDependence(); }; | ||
| paramDeps["NONE"] = []() { return new IdentityParameterParameterDependence(); }; | ||
| } | ||
| } // namespace paramdep | ||
|
|
||
| } // namespace model | ||
|
|
||
| } // namespace cadet |