21template <
typename Solutions,
typename Vector>
30 for (
auto const& bc : *known_solutions)
32 for (std::size_t i = 0; i < bc.ids.size(); ++i)
47 TimeDiscretizedODESystem(
const int process_id,
ODE& ode,
50 _time_disc(time_discretization),
54 _ode.getMatrixSpecifications(process_id), _Jac_id);
56 _ode.getMatrixSpecifications(process_id), _M_id);
58 _ode.getMatrixSpecifications(process_id), _K_id);
60 _ode.getMatrixSpecifications(process_id), _b_id);
75 assemble(std::vector<GlobalVector*>
const& x_new_timestep,
76 std::vector<GlobalVector*>
const& x_prev,
81 auto const t = _time_disc.getCurrentTime();
82 auto const dt = _time_disc.getCurrentTimeIncrement();
83 auto const& x_curr = *x_new_timestep[process_id];
90 _ode.preAssemble(t, dt, x_curr);
91 _ode.assembleWithJacobian(t, dt, x_new_timestep, x_prev, process_id, *_M,
106 double const dt = _time_disc.getCurrentTimeIncrement();
107 _mat_trans->computeResidual(*_M, *_K, *_b, dt, x_new_timestep, x_prev, res);
114 _mat_trans->computeJacobian(*_Jac, Jac);
120 int const process_id)
123 _ode.getKnownSolutions(_time_disc.getCurrentTime(), x, process_id);
138 if (!_known_solutions)
144 std::vector<IndexType> ids;
145 for (
auto const& bc : *_known_solutions)
147 ids.insert(end(ids), begin(bc.ids), end(bc.ids));
151 std::vector<double> values(ids.size(), 0);
157 TimeDiscretizedODESystem(
const int process_id,
ODE& ode,
160 _time_disc(time_discretization),
164 ode.getMatrixSpecifications(process_id), _M_id);
166 ode.getMatrixSpecifications(process_id), _K_id);
168 ode.getMatrixSpecifications(process_id), _b_id);
182 assemble(std::vector<GlobalVector*>
const& x_new_timestep,
183 std::vector<GlobalVector*>
const& x_prev,
184 int const process_id)
188 auto const t = _time_disc.getCurrentTime();
189 auto const dt = _time_disc.getCurrentTimeIncrement();
190 auto const& x_curr = *x_new_timestep[process_id];
196 _ode.preAssemble(t, dt, x_curr);
197 _ode.assemble(t, dt, x_new_timestep, x_prev, process_id, *_M, *_K, *_b);
207 int const process_id)
210 _ode.getKnownSolutions(_time_disc.getCurrentTime(), x, process_id);
226 if (!_known_solutions)
232 std::vector<IndexType> ids;
233 std::vector<double> values;
234 for (
auto const& bc : *_known_solutions)
236 ids.insert(end(ids), begin(bc.ids), end(bc.ids));
237 values.insert(end(values), begin(bc.values), end(bc.values));
Global vector based on Eigen vector.
virtual void releaseMatrix(GlobalMatrix const &A)=0
virtual GlobalMatrix & getMatrix(std::size_t &id)=0
Get an uninitialized matrix with the given id.
virtual GlobalVector & getVector(std::size_t &id)=0
Get an uninitialized vector with the given id.
virtual void releaseVector(GlobalVector const &x)=0
std::unique_ptr< MatrixTranslator< ODETag > > createMatrixTranslator(TimeDiscretization const &timeDisc)
@ FirstOrderImplicitQuasilinear
void finalizeAssembly(PETScMatrix &A)
void applyKnownSolution(EigenMatrix &A, EigenVector &b, EigenVector &, const std::vector< EigenMatrix::IndexType > &vec_knownX_id, const std::vector< double > &vec_knownX_x)
void setVector(PETScVector &v, std::initializer_list< double > values)
void applyKnownSolutions(std::vector< Solutions > const *const known_solutions, Vector &x)
Applies known solutions to the solution vector x.
static NUMLIB_EXPORT MatrixProvider & provider
static NUMLIB_EXPORT VectorProvider & provider