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/auto_ptr.h"
#include "libmesh/point.h"
#include "libmesh/rb_evaluation.h"
#include "libmesh/elem.h"
#include "libmesh/serial_mesh.h"

// C++ includes

namespace libMesh
{

class RBParameters;
class RBParametrizedFunction;

// ------------------------------------------------------------
// RBEIMEvaluation class definition

class RBEIMEvaluation : public RBEvaluation
{
public:

  RBEIMEvaluation (const libMesh::Parallel::Communicator &comm_in
                   LIBMESH_CAN_DEFAULT_TO_COMMWORLD);

  virtual ~RBEIMEvaluation ();

  typedef RBEvaluation Parent;

  virtual void clear();

  virtual void resize_data_structures(const unsigned int Nmax,
                                      bool resize_error_bound_data=true);

  void attach_parametrized_function(RBParametrizedFunction* pf);


  unsigned int get_n_parametrized_functions() const;

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

  void rb_solve(DenseVector<Number>& EIM_rhs);

  void initialize_eim_theta_objects();

  std::vector<RBTheta*> get_eim_theta_objects();

  virtual UniquePtr<RBTheta> build_eim_theta(unsigned int index);

  virtual void write_offline_data_to_files(const std::string& directory_name = "offline_data",
                                           const bool write_binary_data=true);

  virtual void read_offline_data_from_files(const std::string& directory_name = "offline_data",
                                            bool read_error_bound_data=true,
                                            const bool read_binary_data=true);

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

  Point extra_interpolation_point;
  unsigned int extra_interpolation_point_var;
  Elem* extra_interpolation_point_elem;

  DenseVector<Number> extra_interpolation_matrix_row;

private:

  void write_out_interpolation_points_elem(const std::string& directory_name);

  void read_in_interpolation_points_elem(const std::string& directory_name);

  std::vector<RBParametrizedFunction*> _parametrized_functions;

  std::vector<RBTheta*> _rb_eim_theta_objects;

  RBThetaExpansion _empty_rb_theta_expansion;

  RBParameters _previous_parameters;

  unsigned int _previous_N;

  Real _previous_error_bound;

  SerialMesh _interpolation_points_mesh;

};

}

#endif // LIBMESH_RB_EIM_EVALUATION_H