app-bholger/gfem/upcastFormula.h

Go to the documentation of this file.

#pragma once

#include <memory>
#include "toolbox/sharedPointer.hh"
#include "formula/elementFormulaRCP.hh"
#include "formula/formula.hh"

namespace concepts {
namespace gfem {

template<class ScalarT, class ScalarT2>
class UpcastFormula : public concepts::ElementFormula<ScalarT> {
public:
  //typedef typename Realtype<ScalarT>::type G;
  typedef Real G;

  UpcastFormula( const concepts::ElementFormula<ScalarT2>& other ) 
   : formula(other.clone()) 
  {
  }

  virtual ScalarT operator() (const ElementWithCell<G>& elm, const Real p,
      const Real t = 0.0) const
  {
    return formula->operator()(elm, p, t);
  }
  virtual ScalarT operator() (const ElementWithCell<G>& elm, const Real2d& p,
      const Real t = 0.0) const 
  {
    return formula->operator()(elm, p, t);
  }

  virtual ScalarT operator() (const ElementWithCell<G>& elm, const Real3d& p,
      const Real t = 0.0) const
  {
    return formula->operator()(elm, p, t);
  }

  virtual UpcastFormula<ScalarT, ScalarT2>* clone() const {
    return new UpcastFormula<ScalarT, ScalarT2>(*formula);
  }

private:
  std::auto_ptr< concepts::ElementFormula<ScalarT2> > formula;
};

template<class ScalarT>
class FormulaToElemFormula : public concepts::ElementFormula<ScalarT> {
public:
  typedef Real G;

  FormulaToElemFormula(RCP< concepts::Formula<ScalarT> > form) 
  : form(form)
  {

  }

  ~FormulaToElemFormula() {

  }

  FormulaToElemFormula<ScalarT>* clone() const {
    return new FormulaToElemFormula<ScalarT>(*this);
  }

  ScalarT operator() (const ElementWithCell<G>& elm, const Real p,
      const Real t = 0.0) const 
  {
    return form->operator()(p, t);
  }
  ScalarT operator() (const ElementWithCell<G>& elm, const Real2d& p,
      const Real t = 0.0) const 
  {
    return form->operator()(p, t);  
  }
  ScalarT operator() (const ElementWithCell<G>& elm, const Real3d& p,
      const Real t = 0.0) const 
  {
    return form->operator()(p, t);
  }

private:
  RCP<Formula<ScalarT> > form;
};


} //namespace
} //namespace