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
{

// ------------------------------------------------------------
// RBEIMConstruction class definition

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();

  virtual void process_parameters_file (const std::string& parameters_filename);

  void set_best_fit_type_flag (const std::string& best_fit_type_string);

  virtual void print_info();

  virtual void initialize_rb_construction(bool skip_matrix_assembly=false,
                                          bool skip_vector_assembly=false);

  virtual Real train_reduced_basis(const std::string& directory_name = "offline_data",
                                   const bool resize_rb_eval_data=true);

  virtual Real truth_solve(int plot_solution);

  virtual Real compute_best_fit_error();

  virtual void init_context(FEMContext &c);

  Number evaluate_mesh_function(unsigned int var_number,
                                Point p);

  virtual void initialize_eim_assembly_objects();

  std::vector<ElemAssembly*> get_eim_assembly_objects();

  virtual UniquePtr<ElemAssembly> build_eim_assembly(unsigned int bf_index) = 0;

  //----------- PUBLIC DATA MEMBERS -----------//

  BEST_FIT_TYPE best_fit_type_flag;


  CouplingMatrix _coupling_matrix;

protected:

  virtual void init_data();

  virtual void enrich_RB_space();

  virtual void update_system();

  virtual void update_RB_system_matrices();

  virtual Real get_RB_error_bound();

  virtual bool greedy_termination_test(Real training_greedy_error, int count);

  void initialize_parametrized_functions_in_training_set();

  bool _parametrized_functions_in_training_set_initialized;

  std::vector< NumericVector<Number>* > _parametrized_functions_in_training_set;

private:

  MeshFunction* _mesh_function;

  bool _performing_extra_greedy_step;

  UniquePtr< NumericVector<Number> > _ghosted_meshfunction_vector;

  RBAssemblyExpansion _empty_rb_assembly_expansion;

  std::vector<ElemAssembly*> _rb_eim_assembly_objects;

};

} // namespace libMesh

#endif // LIBMESH_RB_EIM_CONSTRUCTION_H