32template <NonlinearSolverTag NLTag>
48template <ODESystemTag ODETag, NonlinearSolverTag NLTag>
85 void assemble(std::vector<GlobalVector*>
const& x_new_timestep,
86 std::vector<GlobalVector*>
const& x_prev,
87 int const process_id)
override;
91 std::vector<GlobalIndexType>
94 return _ode.getIndicesOfResiduumWithoutInitialCompensation();
104 int const process_id)
override;
106 void applyKnownSolutions(
GlobalVector& x)
const override;
112 int const process_id)
override
114 _ode.updateConstraints(lower, upper, process_id);
117 bool isLinear()
const override {
return _ode.isLinear(); }
121 _ode.preIteration(iter, x);
126 return _ode.postIteration(x);
131 const int process_id)
const override
133 return _ode.getMatrixSpecifications(process_id);
144 std::vector<NumLib::IndexValueVector<Index>>
const* _known_solutions =
152 std::size_t _Jac_id = 0u;
153 std::size_t _M_id = 0u;
154 std::size_t _K_id = 0u;
155 std::size_t _b_id = 0u;
188 void assemble(std::vector<GlobalVector*>
const& x_new_timestep,
189 std::vector<GlobalVector*>
const& x_prev,
190 int const process_id)
override;
194 std::vector<GlobalIndexType>
197 return _ode.getIndicesOfResiduumWithoutInitialCompensation();
202 _mat_trans->computeA(*_M, *_K, A);
207 _mat_trans->computeRhs(*_M, *_K, *_b, x_prev, rhs);
211 int const process_id)
override;
213 void applyKnownSolutions(
GlobalVector& x)
const override;
218 bool isLinear()
const override {
return _ode.isLinear(); }
222 _ode.preIteration(iter, x);
227 return _ode.postIteration(x);
232 const int process_id)
const override
234 return _ode.getMatrixSpecifications(process_id);
245 std::vector<NumLib::IndexValueVector<Index>>
const* _known_solutions =
252 std::size_t _M_id = 0u;
253 std::size_t _K_id = 0u;
254 std::size_t _b_id = 0u;
Global vector based on Eigen vector.
virtual TimeDiscretization & getTimeDiscretization()=0
Exposes the used time discretization scheme.
GlobalMatrix * _Jac
the Jacobian of the residual
std::unique_ptr< MatTrans > _mat_trans
the object used to compute the matrix/vector for the nonlinear solver
void preIteration(const unsigned iter, GlobalVector const &x) override
ODE & _ode
ode the ODE being wrapped
MathLib::MatrixSpecifications getMatrixSpecifications(const int process_id) const override
std::vector< GlobalIndexType > getIndicesOfResiduumWithoutInitialCompensation() const override
void updateConstraints(GlobalVector &lower, GlobalVector &upper, int const process_id) override
TimeDisc & _time_disc
the time discretization to being used
GlobalMatrix * _M
Matrix .
GlobalMatrix * _K
Matrix .
GlobalVector * _b
Matrix .
bool isLinear() const override
TimeDisc & getTimeDiscretization() override
Exposes the used time discretization scheme.
IterationResult postIteration(GlobalVector const &x) override
MathLib::MatrixVectorTraits< GlobalMatrix >::Index Index
void getA(GlobalMatrix &A) const override
GlobalMatrix * _K
Matrix .
void preIteration(const unsigned iter, GlobalVector const &x) override
GlobalMatrix * _M
Matrix .
bool isLinear() const override
ODE & _ode
ode the ODE being wrapped
MathLib::MatrixVectorTraits< GlobalMatrix >::Index Index
std::vector< GlobalIndexType > getIndicesOfResiduumWithoutInitialCompensation() const override
MathLib::MatrixSpecifications getMatrixSpecifications(const int process_id) const override
IterationResult postIteration(GlobalVector const &x) override
GlobalVector * _b
Matrix .
TimeDisc & getTimeDiscretization() override
Exposes the used time discretization scheme.
TimeDisc & _time_disc
the time discretization to being used
void getRhs(GlobalVector const &x_prev, GlobalVector &rhs) const override
std::unique_ptr< MatTrans > _mat_trans
the object used to compute the matrix/vector for the nonlinear solver
ODESystemTag
Tag used to specify the type of ODE.
IterationResult
Status flags telling the NonlinearSolver if an iteration succeeded.
@ FirstOrderImplicitQuasilinear