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