raster_cp.h

// -*-C++-*-

#ifndef __RASTER_CP_H
#define __RASTER_CP_H

#include "raster.h"

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

    ~Raster_CubicPixel ();

    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;


protected:
    double Value (int basisidx, const Point &p);
    
    double Value (const IVector &basisgrd, const Point &p);

private:
    RVector ibgrid;            
    RVector iggrid;            
    RGenericSparseMatrix *G;   
    RGenericSparseMatrix *GI;  
};

#endif // !__RASTER_CP_H