aglowav2/bform.hh

Go to the documentation of this file.

/* Bilinear form for the agglomerated wavelets.
 */

#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