rb_scm_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_SCM_CONSTRUCTION_H
#define LIBMESH_RB_SCM_CONSTRUCTION_H

// Configuration data
#include "libmesh/libmesh_config.h"

// Currently, the RBSCMConstruction is only usable
// if SLEPc is enabled.
#if defined(LIBMESH_HAVE_SLEPC) && (LIBMESH_HAVE_GLPK)

// rbOOmit includes
#include "libmesh/rb_construction_base.h"
#include "libmesh/rb_scm_evaluation.h"

// libMesh includes
#include "libmesh/condensed_eigen_system.h"

// C++ includes

namespace libMesh
{

// ------------------------------------------------------------
// RBSCMConstruction class definition

class RBSCMConstruction : public RBConstructionBase<CondensedEigenSystem>
{
public:

  RBSCMConstruction (EquationSystems& es,
                     const std::string& name_in,
                     const unsigned int number_in);

  virtual ~RBSCMConstruction ();

  typedef RBSCMConstruction sys_type;

  typedef RBConstructionBase<CondensedEigenSystem> Parent;

  virtual void clear ();

  void set_rb_scm_evaluation(RBSCMEvaluation& rb_scm_eval_in);

  RBSCMEvaluation& get_rb_scm_evaluation();

  RBThetaExpansion& get_rb_theta_expansion();

  virtual void resize_SCM_vectors ();

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

  virtual void print_info();

  virtual void set_eigensolver_properties(int ) { }

  void set_RB_system_name(const std::string& new_name)
  { RB_system_name = new_name; }

  Real get_SCM_training_tolerance() const                         { return SCM_training_tolerance; }
  void set_SCM_training_tolerance(Real SCM_training_tolerance_in) { this->SCM_training_tolerance = SCM_training_tolerance_in; }

  virtual void perform_SCM_greedy();

  virtual void attach_deflation_space() {}

protected:

  virtual void add_scaled_symm_Aq(unsigned int q_a, Number scalar);

  virtual void load_matrix_B() ;

  virtual void compute_SCM_bounding_box();

  virtual void evaluate_stability_constant();

  virtual void enrich_C_J(unsigned int new_C_J_index);

  virtual std::pair<unsigned int,Real> compute_SCM_bounds_on_training_set();

  Number B_inner_product(const NumericVector<Number>& v, const NumericVector<Number>& w) const;

  Number Aq_inner_product(unsigned int q,
                          const NumericVector<Number>& v,
                          const NumericVector<Number>& w);

  virtual Real SCM_greedy_error_indicator(Real LB, Real UB) { return fabs(UB-LB)/fabs(UB); }

  //----------- PROTECTED DATA MEMBERS -----------//

  Real SCM_training_tolerance;

  std::string RB_system_name;

private:

  RBSCMEvaluation* rb_scm_eval;

};

} // namespace libMesh

#endif // LIBMESH_HAVE_SLEPC && LIBMESH_HAVE_GLPK

#endif // LIBMESH_RB_SCM_CONSTRUCTION_H