eim_classes.h

Go to the documentation of this file.

#ifndef EIM_CLASSES_H
#define EIM_CLASSES_H
 
// libMesh includes
#include "libmesh/rb_eim_construction.h"
#include "libmesh/rb_eim_evaluation.h"
#include "libmesh/auto_ptr.h" // libmesh_make_unique
 
// Example includes
#include "assembly.h"
 
#ifndef LIBMESH_HAVE_CXX14_MAKE_UNIQUE
using libMesh::make_unique;
#endif
 
// A simple subclass of RBEIMEvaluation. Overload
// evaluate_parametrized_function to define the
// function that we "empirically" interpolate.
class SimpleEIMEvaluation : public RBEIMEvaluation
{
public:
 
  SimpleEIMEvaluation(const libMesh::Parallel::Communicator & comm)
    : RBEIMEvaluation(comm)
  {
    attach_parametrized_function(&g_x);
    attach_parametrized_function(&g_y);
    attach_parametrized_function(&g_z);
  }
 
  virtual std::unique_ptr<RBTheta> build_eim_theta(unsigned int index)
  {
    return libmesh_make_unique<ThetaEIM>(*this, index);
  }
 
  Gx g_x;
  Gy g_y;
  Gz g_z;
};
 
// A simple subclass of RBEIMConstruction.
class SimpleEIMConstruction : public RBEIMConstruction
{
public:
 
  SimpleEIMConstruction (EquationSystems & es,
                         const std::string & name_in,
                         const unsigned int number_in)
    : Parent(es, name_in, number_in)
  {
  }
 
  typedef RBEIMConstruction Parent;
 
  virtual std::unique_ptr<ElemAssembly> build_eim_assembly(unsigned int index)
  {
    return libmesh_make_unique<AssemblyEIM>(*this, index);
  }
 
  virtual void init_data()
  {
    Parent::init_data();
 
    set_inner_product_assembly(eim_ip);
  }
 
  virtual void init_implicit_system()
  {
    this->add_variable ("L2_proj_var", libMesh::FIRST);
  }
 
  virtual void init_explicit_system()
  {
    Gx_var = get_explicit_system().add_variable ("x_comp_of_G", libMesh::FIRST);
    Gy_var = get_explicit_system().add_variable ("y_comp_of_G", libMesh::FIRST);
    Gz_var = get_explicit_system().add_variable ("z_comp_of_G", libMesh::FIRST);
  }
 
  unsigned int Gx_var;
  unsigned int Gy_var;
  unsigned int Gz_var;
 
  Ex6EIMInnerProduct eim_ip;
};
 
#endif