raster_px2.h

// -*-C++-*-

#ifndef __RASTER_PX2_H
#define __RASTER_PX2_H

#include "raster.h"

void Tetsplit (Mesh *mesh, int el, int cut_orient, double cut_pos);
// exported temporarily

// =========================================================================
class STOASTLIB Raster_Pixel2: public Raster {
public:
    Raster_Pixel2 (const IVector &_bdim, const IVector &_gdim, Mesh *mesh,
        RDenseMatrix *bb = 0, double _map_tol=1e-10);

    ~Raster_Pixel2 ();

    double Value_nomask (const Point &p, int i, bool is_solidx=true) const;

    void Map_GridToBasis (const RVector &gvec, RVector &bvec) const;

    void Map_GridToBasis (const CVector &gvec, CVector &bvec) const;

    void Map_BasisToGrid (const RVector &bvec, RVector &gvec) const;

    void Map_BasisToGrid (const CVector &bvec, CVector &gvec) const;

    void Map_MeshToBasis (const RVector &mvec, RVector &bvec) const;

    void Map_BasisToMesh (const RVector &bvec, RVector &mvec) const;

    void Map_BasisToSol (const RVector &bvec, RVector &svec) const;

    void Map_SolToBasis (const RVector &svec, RVector &bvec) const;

    void Map_MeshToSol (const RVector &mvec, RVector &svec) const;

    void Map_SolToMesh (const RVector &svec, RVector &mvec) const;


private:
    bool grid_is_basis;       
    RCompRowMatrix *D;        
    RCompRowMatrix *Buu;      
    RCompRowMatrix *Bvv;      
    RCompRowMatrix *Buv;      

    RPrecon_IC *Buu_precon;   
    RPrecon_IC *Bvv_precon;   
    double map_tol;           

    RCompRowMatrix *CreateMixedMassmat () const;
    RCompRowMatrix *CreateMixedMassmat_tet4 () const;

    RCompRowMatrix *CreatePixelMassmat () const;

    int SutherlandHodgman (int el, int xgrid, int ygrid, Point *clip_poly,
        int npoly) const;
};
    
#endif // !__RASTER_PX2_H