Home
Publications
Presentations
Support
Download
Examples
Developers
Class Docs
Namespaces
Classes
Files
Mailing Lists
Gallery
Wiki

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

  virtual void print_info();

  virtual void initialize_rb_construction();

  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 AutoPtr<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;

  AutoPtr< 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

---

Site Created By: libMesh Developers
Last modified: February 05 2013 19:54:48 UTC

Hosted By: