bform.hh

Go to the documentation of this file.

 
#ifndef aglowav2bform_hh
#define aglowav2bform_hh
 
#ifdef __GUNG__
#pragma interface
#endif
 
#include "operator/bilinearForm.hh"
#include "aglowav2/xfy.hh"
 
namespace aglowav2 {
 
  // *********************************************************** AglowavBF00 **
 
  template<class F = concepts::Real>
  class AglowavBF00 : public concepts::BilinearForm<F> {
  private:
    typedef struct Nfld {
      struct Nfld*        lnk;
      const Haar3d000<F>* elm;
      F*                  valXY;
      F*                  valYX;
    } Nfld;
 
    struct Cell {
      F* mtsttrl;
      F* mtrltst;
 
      const Haar3d000<F>* elm;
 
      Cell() : mtsttrl(0), mtrltst(0), elm(0) {}
    };
 
    virtual AglowavBF00* clone() const { return new AglowavBF00(bf_, XFY_); }
  public:
    AglowavBF00(concepts::BilinearForm<F>& bf, F00<F>& XFY);
    ~AglowavBF00();
 
    void operator()(const concepts::Element<F>& elmX,
        const concepts::Element<F>& elmY,
        concepts::ElementMatrix<F>& em);
    void operator()(const Haar3d000<F>& elmX, const Haar3d000<F>& elmY,
        concepts::ElementMatrix<F>& em);
 
    uint memory() const {return mem_;}
 
  private:
    void nearfield_(const Haar3d000<F>& elmX, const Haar3d000<F>& elmY,
        F* valXY, uint ldXY, F* valYX, uint ldYX);
 
    void updateChld_(const Haar3d000<F>& tst, const Haar3d000<F>& trl,
         F* valtsttrl, uint ldtsttrl, F* valtrltst, uint ldtrltst);
 
    void updateElm_(const Haar3d000<F>& tst, const Haar3d000<F>& trl,
        F* valtsttrl, uint ldtsttrl, F* valtrltst, uint ldtrltst);
 
    concepts::BilinearForm<F>& bf_;
    F00<F>& XFY_;
    Nfld** nfld_;
    uint mem_;
 
    Nfld*** nfldLast_;
 
    F* valXY_;
    F* vXY_;
    F* valYX_;
    F* vYX_;
    F* foo_;
 
    Cell* cache_;
    Cell* buf_;
    bool reuse_;
  };
 
  // *********************************************************** WavIdentity **
 
  template<class F = concepts::Real>
  class WavIdentity : public concepts::BilinearForm<F> {
  public:
    void operator()(const concepts::Element<F>& elmX,
        const concepts::Element<F>& elmY,
        concepts::ElementMatrix<F>& em);
    void operator()(const Haar3d000<F>& elmX, const Haar3d000<F>& elmY,
        concepts::ElementMatrix<F>& em);
    virtual WavIdentity* clone() const { return new WavIdentity(); }
  };
 
} // namespace aglowav2
 
#endif // aglowav2bform_hh