Skip to content
Snippets Groups Projects
Commit 635777c6 authored by Elias Pipping's avatar Elias Pipping
Browse files

[Cleanup] Continue renaming: alpha -> logState

parent 5d190b0f
No related branches found
No related tags found
No related merge requests found
...@@ -8,7 +8,7 @@ template <class SingletonVectorType, class VectorType> ...@@ -8,7 +8,7 @@ template <class SingletonVectorType, class VectorType>
class DieterichStateUpdater class DieterichStateUpdater
: public StateUpdater<SingletonVectorType, VectorType> { : public StateUpdater<SingletonVectorType, VectorType> {
public: public:
DieterichStateUpdater(SingletonVectorType _alpha_initial, DieterichStateUpdater(SingletonVectorType _logState_initial,
Dune::BitSetVector<1> const &_nodes, double _L); Dune::BitSetVector<1> const &_nodes, double _L);
virtual void nextTimeStep(); virtual void nextTimeStep();
...@@ -18,7 +18,7 @@ class DieterichStateUpdater ...@@ -18,7 +18,7 @@ class DieterichStateUpdater
private: private:
SingletonVectorType logState_o; SingletonVectorType logState_o;
SingletonVectorType alpha; SingletonVectorType logState;
Dune::BitSetVector<1> const &nodes; Dune::BitSetVector<1> const &nodes;
double const L; double const L;
double tau; double tau;
...@@ -26,13 +26,13 @@ class DieterichStateUpdater ...@@ -26,13 +26,13 @@ class DieterichStateUpdater
template <class SingletonVectorType, class VectorType> template <class SingletonVectorType, class VectorType>
DieterichStateUpdater<SingletonVectorType, VectorType>::DieterichStateUpdater( DieterichStateUpdater<SingletonVectorType, VectorType>::DieterichStateUpdater(
SingletonVectorType _alpha_initial, Dune::BitSetVector<1> const &_nodes, SingletonVectorType _logState_initial, Dune::BitSetVector<1> const &_nodes,
double _L) double _L)
: alpha(_alpha_initial), nodes(_nodes), L(_L) {} : logState(_logState_initial), nodes(_nodes), L(_L) {}
template <class SingletonVectorType, class VectorType> template <class SingletonVectorType, class VectorType>
void DieterichStateUpdater<SingletonVectorType, VectorType>::nextTimeStep() { void DieterichStateUpdater<SingletonVectorType, VectorType>::nextTimeStep() {
logState_o = alpha; logState_o = logState;
} }
template <class SingletonVectorType, class VectorType> template <class SingletonVectorType, class VectorType>
...@@ -47,14 +47,14 @@ void DieterichStateUpdater<SingletonVectorType, VectorType>::solve( ...@@ -47,14 +47,14 @@ void DieterichStateUpdater<SingletonVectorType, VectorType>::solve(
for (size_t i = 0; i < nodes.size(); ++i) for (size_t i = 0; i < nodes.size(); ++i)
if (nodes[i][0]) { if (nodes[i][0]) {
double const V = velocity_field[i].two_norm(); double const V = velocity_field[i].two_norm();
alpha[i] = state_update_dieterich_euler(tau, V / L, logState_o[i]); logState[i] = state_update_dieterich_euler(tau, V / L, logState_o[i]);
} }
} }
template <class SingletonVectorType, class VectorType> template <class SingletonVectorType, class VectorType>
void DieterichStateUpdater<SingletonVectorType, VectorType>::extractLogState( void DieterichStateUpdater<SingletonVectorType, VectorType>::extractLogState(
SingletonVectorType &logState) { SingletonVectorType &_logState) {
logState = alpha; _logState = logState;
} }
#endif #endif
...@@ -7,7 +7,7 @@ ...@@ -7,7 +7,7 @@
template <class SingletonVectorType, class VectorType> template <class SingletonVectorType, class VectorType>
class RuinaStateUpdater : public StateUpdater<SingletonVectorType, VectorType> { class RuinaStateUpdater : public StateUpdater<SingletonVectorType, VectorType> {
public: public:
RuinaStateUpdater(SingletonVectorType _alpha_initial, RuinaStateUpdater(SingletonVectorType _logState_initial,
Dune::BitSetVector<1> const &_nodes, double _L); Dune::BitSetVector<1> const &_nodes, double _L);
virtual void nextTimeStep(); virtual void nextTimeStep();
...@@ -17,7 +17,7 @@ class RuinaStateUpdater : public StateUpdater<SingletonVectorType, VectorType> { ...@@ -17,7 +17,7 @@ class RuinaStateUpdater : public StateUpdater<SingletonVectorType, VectorType> {
private: private:
SingletonVectorType logState_o; SingletonVectorType logState_o;
SingletonVectorType alpha; SingletonVectorType logState;
Dune::BitSetVector<1> const &nodes; Dune::BitSetVector<1> const &nodes;
double const L; double const L;
double tau; double tau;
...@@ -25,13 +25,13 @@ class RuinaStateUpdater : public StateUpdater<SingletonVectorType, VectorType> { ...@@ -25,13 +25,13 @@ class RuinaStateUpdater : public StateUpdater<SingletonVectorType, VectorType> {
template <class SingletonVectorType, class VectorType> template <class SingletonVectorType, class VectorType>
RuinaStateUpdater<SingletonVectorType, VectorType>::RuinaStateUpdater( RuinaStateUpdater<SingletonVectorType, VectorType>::RuinaStateUpdater(
SingletonVectorType _alpha_initial, Dune::BitSetVector<1> const &_nodes, SingletonVectorType _logState_initial, Dune::BitSetVector<1> const &_nodes,
double _L) double _L)
: alpha(_alpha_initial), nodes(_nodes), L(_L) {} : logState(_logState_initial), nodes(_nodes), L(_L) {}
template <class SingletonVectorType, class VectorType> template <class SingletonVectorType, class VectorType>
void RuinaStateUpdater<SingletonVectorType, VectorType>::nextTimeStep() { void RuinaStateUpdater<SingletonVectorType, VectorType>::nextTimeStep() {
logState_o = alpha; logState_o = logState;
} }
template <class SingletonVectorType, class VectorType> template <class SingletonVectorType, class VectorType>
...@@ -45,14 +45,14 @@ void RuinaStateUpdater<SingletonVectorType, VectorType>::solve( ...@@ -45,14 +45,14 @@ void RuinaStateUpdater<SingletonVectorType, VectorType>::solve(
for (size_t i = 0; i < nodes.size(); ++i) for (size_t i = 0; i < nodes.size(); ++i)
if (nodes[i][0]) { if (nodes[i][0]) {
double const V = velocity_field[i].two_norm(); double const V = velocity_field[i].two_norm();
alpha[i] = state_update_ruina(tau, V * tau / L, logState_o[i]); logState[i] = state_update_ruina(tau, V * tau / L, logState_o[i]);
} }
} }
template <class SingletonVectorType, class VectorType> template <class SingletonVectorType, class VectorType>
void RuinaStateUpdater<SingletonVectorType, VectorType>::extractLogState( void RuinaStateUpdater<SingletonVectorType, VectorType>::extractLogState(
SingletonVectorType &logState) { SingletonVectorType &_logState) {
logState = alpha; _logState = logState;
} }
#endif #endif
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment