rb_eim_construction.h

Go to the documentation of this file.

// rbOOmit: An implementation of the Certified Reduced Basis method.
// Copyright (C) 2009, 2010 David J. Knezevic
 
// This file is part of rbOOmit.
 
// rbOOmit is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
 
// rbOOmit is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.
 
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
#ifndef LIBMESH_RB_EIM_CONSTRUCTION_H
#define LIBMESH_RB_EIM_CONSTRUCTION_H
 
// rbOOmit includes
#include "libmesh/rb_construction.h"
#include "libmesh/rb_assembly_expansion.h"
#include "libmesh/rb_eim_assembly.h"
 
// libMesh includes
#include "libmesh/mesh_function.h"
#include "libmesh/coupling_matrix.h"
 
// C++ includes
 
namespace libMesh
{
 
class RBEIMConstruction : public RBConstruction
{
public:
 
  enum BEST_FIT_TYPE { PROJECTION_BEST_FIT, EIM_BEST_FIT };
 
  RBEIMConstruction (EquationSystems & es,
                     const std::string & name,
                     const unsigned int number);
 
  virtual ~RBEIMConstruction ();
 
  typedef RBEIMConstruction sys_type;
 
  typedef RBConstruction Parent;
 
  virtual void clear() override;
 
  virtual void process_parameters_file (const std::string & parameters_filename) override;
 
  void set_best_fit_type_flag (const std::string & best_fit_type_string);
 
  virtual void print_info() override;
 
  virtual void initialize_rb_construction(bool skip_matrix_assembly=false,
                                          bool skip_vector_assembly=false) override;
 
  virtual Real train_reduced_basis(const bool resize_rb_eval_data=true) override;
 
  virtual Real truth_solve(int plot_solution) override;
 
  virtual Real compute_best_fit_error();
 
  virtual void init_context_with_sys(FEMContext & c, System & sys);
 
  virtual void init_explicit_system() = 0;
 
  virtual void init_implicit_system() = 0;
 
  Number evaluate_mesh_function(unsigned int var_number,
                                Point p);
 
  void set_point_locator_tol(Real point_locator_tol);
 
  Real get_point_locator_tol() const;
 
  virtual void initialize_eim_assembly_objects();
 
  std::vector<std::unique_ptr<ElemAssembly>> & get_eim_assembly_objects();
 
  virtual std::unique_ptr<ElemAssembly> build_eim_assembly(unsigned int bf_index) = 0;
 
  ExplicitSystem & get_explicit_system();
 
  virtual void load_basis_function(unsigned int i) override;
 
  virtual void load_rb_solution() override;
 
  void set_explicit_sys_subvector(NumericVector<Number>& dest,
                                  unsigned int var,
                                  NumericVector<Number>& source);
 
  void get_explicit_sys_subvector(NumericVector<Number>& dest,
                                  unsigned int var,
                                  NumericVector<Number>& localized_source);
 
  void init_dof_map_between_systems();
 
  void plot_parametrized_functions_in_training_set(const std::string & pathname);
 
  //----------- PUBLIC DATA MEMBERS -----------//
 
  BEST_FIT_TYPE best_fit_type_flag;
 
protected:
 
  virtual void init_data() override;
 
  virtual void enrich_RB_space() override;
 
  virtual void update_system() override;
 
  virtual void update_RB_system_matrices() override;
 
  virtual Real get_RB_error_bound() override;
 
  void initialize_parametrized_functions_in_training_set();
 
  bool _parametrized_functions_in_training_set_initialized;
 
  std::vector<std::unique_ptr<NumericVector<Number>>> _parametrized_functions_in_training_set;
 
private:
 
  std::unique_ptr<MeshFunction> _mesh_function;
 
  std::unique_ptr<NumericVector<Number>> _ghosted_meshfunction_vector;
 
  RBAssemblyExpansion _empty_rb_assembly_expansion;
 
  std::vector<std::unique_ptr<ElemAssembly>> _rb_eim_assembly_objects;
 
  std::string _explicit_system_name;
 
  std::vector<std::vector<dof_id_type>> _dof_map_between_systems;
 
  std::vector<std::unique_ptr<NumericVector<Number>>> _matrix_times_bfs;
 
  Real _point_locator_tol;
};
 
} // namespace libMesh
 
#endif // LIBMESH_RB_EIM_CONSTRUCTION_H