raster_px.h

// -*-C++-*-

#ifndef __RASTER_PX_H
#define __RASTER_PX_H

#include "raster.h"

// =========================================================================
class STOASTLIB Raster_Pixel: public Raster {
public:
    Raster_Pixel (const IVector &_bdim, const IVector &_gdim, Mesh *mesh,
        RDenseMatrix *bb = 0);

    ~Raster_Pixel ();

    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;

    const RGenericSparseMatrix &Mesh2BasisMatrix() const { return *C; }

    const RGenericSparseMatrix &Basis2MeshMatrix() const { return *CI; }

private:
    bool grid_is_basis;       
    RGenericSparseMatrix *G;  
    RGenericSparseMatrix *GI; 
    RGenericSparseMatrix *C;  
    RGenericSparseMatrix *CI; 
    RCompRowMatrix *D;        
};
    
#endif // !__RASTER_PX_H