20 #ifndef LIBMESH_TRANSIENT_RB_EVALUATION_H 21 #define LIBMESH_TRANSIENT_RB_EVALUATION_H 24 #include "libmesh/rb_evaluation.h" 25 #include "libmesh/rb_temporal_discretization.h" 34 class TransientRBThetaExpansion;
80 virtual void clear()
override;
90 bool resize_error_bound_data=
true)
override;
158 const bool write_binary_data=
true)
override;
167 bool read_error_bound_data=
true,
168 const bool read_binary_data=
true)
override;
263 #endif // LIBMESH_TRANSIENT_RB_EVALUATION_H std::vector< std::vector< std::vector< Number > > > Mq_Mq_representor_innerprods
virtual void resize_data_structures(const unsigned int Nmax, bool resize_error_bound_data=true) override
Resize and clear the data vectors corresponding to the value of Nmax.
TransientRBEvaluation & operator=(const TransientRBEvaluation &)=delete
virtual Real uncached_compute_residual_dual_norm(const unsigned int N)
Compute the dual norm of the residual for the solution saved in RB_solution.
Number cached_Fq_term
Cached residual terms.
DenseMatrix< Number > RB_LHS_matrix
Cached data for subsequent solves.
virtual void clear_riesz_representors() override
Clear all the Riesz representors that are used to compute the RB residual (and hence error bound)...
std::vector< DenseMatrix< Number > > RB_M_q_vector
Dense matrices for the RB mass matrices.
void cache_online_residual_terms(const unsigned int N)
Helper function for caching the terms in the online residual assembly that do not change in time...
virtual Real compute_residual_dual_norm(const unsigned int N)
Compute the dual norm of the residual for the solution saved in RB_solution_vector.
DenseVector< Number > old_RB_solution
The RB solution at the previous time-level.
DenseMatrix< Number > cached_Aq_Aq_matrix
virtual Real compute_residual_dual_norm(const unsigned int N) override
Compute the dual norm of the residual for the solution saved in RB_solution_vector.
This class is part of the rbOOmit framework.
The libMesh namespace provides an interface to certain functionality in the library.
virtual Real get_error_bound_normalization() override
DenseVector< Number > cached_Fq_Mq_vector
DenseVector< Number > RB_RHS_save
std::vector< std::vector< std::unique_ptr< NumericVector< Number > > > > M_q_representor
Vector storing the mass matrix representors.
virtual void legacy_write_offline_data_to_files(const std::string &directory_name="offline_data", const bool write_binary_data=true) override
Write out all the data to text files in order to segregate the Offline stage from the Online stage...
bool _rb_solve_data_cached
Check that the data has been cached in case of using rb_solve_again.
DenseMatrix< Number > cached_Aq_Mq_matrix
std::vector< Real > initial_L2_error_all_N
Vector storing initial L2 error for all 1 <= N <= RB_size.
virtual Real residual_scaling_numer(Real alpha_LB)
Specifies the residual scaling on the numerator to be used in the a posteriori error bound...
DenseVector< Number > cached_Fq_Aq_vector
std::vector< std::vector< std::vector< Number > > > Fq_Mq_representor_innerprods
Vectors storing the residual representor inner products to be used in computing the residuals online...
std::vector< std::vector< Real > > RB_output_error_bounds_all_k
The error bounds for each RB output for all time-levels from the most recent rb_solve.
std::vector< Real > error_bound_all_k
The error bound data for all time-levels from the most recent rb_solve.
virtual void clear() override
Clear this TransientRBEvaluation object.
RBEvaluation Parent
The type of the parent.
virtual Real rb_solve_again()
If a solve has already been performed, then we cached some data and we can perform a new solve much m...
std::vector< DenseVector< Number > > RB_initial_condition_all_N
The RB initial conditions (i.e.
DenseMatrix< Number > cached_Mq_Mq_matrix
virtual Real rb_solve(unsigned int N)
Perform online solve with the N RB basis functions, for the set of parameters in current_params, where 0 <= N <= RB_size.
This class is part of the rbOOmit framework.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
virtual void legacy_read_offline_data_from_files(const std::string &directory_name="offline_data", bool read_error_bound_data=true, const bool read_binary_data=true) override
Read in the saved Offline reduced basis data to initialize the system for Online solves.
std::vector< DenseVector< Number > > RB_temporal_solution_data
Array storing the solution data at each time level from the most recent solve.
DenseMatrix< Number > RB_RHS_matrix
std::vector< std::vector< Number > > RB_outputs_all_k
The RB outputs for all time-levels from the most recent rb_solve.
virtual ~TransientRBEvaluation()
Define a class that encapsulates the details of a "generalized Euler" temporal discretization to be u...
std::vector< std::vector< std::vector< std::vector< Number > > > > Aq_Mq_representor_innerprods
TransientRBEvaluation(const Parallel::Communicator &comm_in)
Constructor.
virtual Real rb_solve(unsigned int N) override
Perform online solve with the N RB basis functions, for the set of parameters in current_params, where 0 <= N <= RB_size.
DenseMatrix< Number > RB_L2_matrix
Dense RB L2 matrix.