app-bholger/gfem/gfemAdvectionIgnore.h

Go to the documentation of this file.

#pragma once

#include "formula/elementFormula.hh"
#include "operator/bilinearForm.hh"
#include "gfemQuad.h"
#include "gfemQuadFast.h"
#include "riesz.h"

namespace concepts {
namespace gfem {

class GfemAdvectionIgnore : public hp2D::Advection<Cmplx> 
{
public:
    typedef concepts::RCP<const concepts::ElementFormula<Cmplx> >
      ElementFormulaRCP;

    typedef concepts::RCP<const concepts::ElementFormula<Point< Cmplx, 2 > > >
      ElementFormulaPoint2dRCP;

    GfemAdvectionIgnore(ElementFormulaRCP frm1, 
        ElementFormulaRCP frm2)
      : hp2D::Advection<Cmplx>(frm1, frm2, false)
    { }

    GfemAdvectionIgnore(ElementFormulaPoint2dRCP frm)
      : hp2D::Advection<Cmplx>(frm, false)
    { }

    virtual GfemAdvectionIgnore* clone() const { 
      return new GfemAdvectionIgnore(this->frm_);
    }

    virtual void operator()(
          const concepts::Element<Real>& elmX,
          const concepts::Element<Real>& elmY,
          concepts::ElementMatrix<Cmplx>& em) 
    {
      const hp2D::Quad<Real>* elmYH = dynamic_cast<const hp2D::Quad<Real>* >(&elmY);
      const hp2D::Quad<Real>* elmXH = dynamic_cast<const hp2D::Quad<Real>* >(&elmX);
      if (elmYH && elmXH) { 
        hp2D::Advection<Cmplx>::operator()(*elmXH, *elmYH, em);              
        return;
      }

      // FIXME: this should not be necessary ... change concepts
      uint n = elmX.T().m();
      uint m = elmX.T().m();
      em.resize(n, m);
      em.zeros();

      // simply ignore all non-quad elements
      return;
    } 

    virtual ~GfemAdvectionIgnore() {}
};

} // namespace
} // namespace