rb_eim_evaluation.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_EVALUATION_H
#define LIBMESH_RB_EIM_EVALUATION_H
 
// libMesh includes
#include "libmesh/point.h"
#include "libmesh/rb_evaluation.h"
#include "libmesh/replicated_mesh.h"
#include "libmesh/rb_theta_expansion.h"
 
// C++ includes
#include <memory>
 
namespace libMesh
{
 
class RBParameters;
class RBParametrizedFunction;
class Elem;
class RBTheta;
 
class RBEIMEvaluation : public RBEvaluation
{
public:
 
  RBEIMEvaluation (const libMesh::Parallel::Communicator & comm_in);
 
  virtual ~RBEIMEvaluation ();
 
  typedef RBEvaluation Parent;
 
  virtual void clear() override;
 
  virtual void resize_data_structures(const unsigned int Nmax,
                                      bool resize_error_bound_data=true) override;
 
  void attach_parametrized_function(RBParametrizedFunction * pf);
 
 
  unsigned int get_n_parametrized_functions() const;
 
  ReplicatedMesh & get_interpolation_points_mesh();
 
  Number evaluate_parametrized_function(unsigned int var_index,
                                        const Point & p,
                                        const Elem & elem);
 
  virtual Real rb_solve(unsigned int N) override;
 
  void rb_solve(DenseVector<Number> & EIM_rhs);
 
  virtual Real get_error_bound_normalization() override;
 
  void initialize_eim_theta_objects();
 
  std::vector<std::unique_ptr<RBTheta>> & get_eim_theta_objects();
 
  virtual std::unique_ptr<RBTheta> build_eim_theta(unsigned int index);
 
  virtual void legacy_write_offline_data_to_files(const std::string & directory_name = "offline_data",
                                                  const bool write_binary_data=true) override;
 
  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;
 
  //----------- PUBLIC DATA MEMBERS -----------//
 
  DenseMatrix<Number> interpolation_matrix;
 
  std::vector<Point> interpolation_points;
 
  std::vector<unsigned int> interpolation_points_var;
 
  std::vector<Elem *> interpolation_points_elem;
 
private:
 
  void legacy_write_out_interpolation_points_elem(const std::string & directory_name);
 
  void legacy_read_in_interpolation_points_elem(const std::string & directory_name);
 
  std::vector<RBParametrizedFunction *> _parametrized_functions;
 
  std::vector<std::unique_ptr<RBTheta>> _rb_eim_theta_objects;
 
  RBThetaExpansion _empty_rb_theta_expansion;
 
  RBParameters _previous_parameters;
 
  unsigned int _previous_N;
 
  Real _previous_error_bound;
 
  ReplicatedMesh _interpolation_points_mesh;
 
};
 
}
 
#endif // LIBMESH_RB_EIM_EVALUATION_H