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"
 
// libMesh includes
#include "libmesh/condensed_eigen_system.h"
 
// C++ includes
 
namespace libMesh
{
 
// Forward declarations
class RBSCMEvaluation;
 
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 () override;
 
  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