raster_bl.h

// -*-C++-*-

#ifndef __RASTER_BL_H
#define __RASTER_BL_H

#include "raster.h"

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

    virtual ~Raster_Blob();

    virtual double Value (const Point &p, int i) const;

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

    virtual void NodeValues (int node, RVector &nv) const;
    RVector NodeValues (int node) const
    { RVector tmp(slen); NodeValues (node, tmp); return tmp; }

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

    void Map_MeshToGrid (const RVector &mvec, RVector &gvec) const;
    void Map_MeshToGrid (const CVector &mvec, CVector &gvec) const;

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

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

    void Map_GridToMesh (const RVector &gvec, RVector &mvec) const;
    void Map_GridToMesh (const CVector &gvec, CVector &mvec) const;

    void Map_SolToMesh (const RVector &svec, RVector &mvec) const;
    void Map_SolToMesh (const CVector &svec, CVector &mvec) 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_SolToBasis (const RVector &svec, RVector &bvec) const;
    void Map_SolToBasis (const CVector &svec, CVector &bvec) const;

    void Map_BasisToSol (const RVector &bvec, RVector &svec) const;
    void Map_BasisToSol (const CVector &bvec, CVector &svec) const;
    
    void Map_SolToGrid (const RVector &svec, RVector &gvec) const;
    void Map_SolToGrid (const CVector &svec, CVector &gvec) const;

    void Map_GridToSol (const RVector &gvec, RVector &svec) const;
    void Map_GridToSol (const CVector &gvec, CVector &svec) const;

protected:
    double ComputeBasisScale ();

    void ComputeNodeValues ();

    RCompRowMatrix *nodevals;
    // Sparse matrix of size nlen x slen
    // containing nodal basis values for each basis function
    // Initialised by ComputeNodeValues()

    double sup;               // blob support radius in units of grid spacing
    RVector igrid;            // inverse grid spacing in x, y (, z)
    int npad;                 // number of padding layers around bounding box
    Point bbmin_pad, bbmax_pad; // bounding box including padding
};

#endif // !__RASTER_BL_H