wmunstructuredreconstruction.h

Go to the documentation of this file.

 
#ifndef WMUNSTRUCTUREDRECONSTRUCTION_H_
#define WMUNSTRUCTUREDRECONSTRUCTION_H_
 
// Wm includes
#include "lib/wmunstructuredpatch.h"
#include "lib/geometry/unstructured_geometry/wmunstructuredgeometry.h"
 
class WmUnstructuredStencilConstructor;
 
class WmUnstructuredReconstruction
{
public:
    WmUnstructuredReconstruction(const WmUnstructuredPatch& patch,
                                 int numReconstructionElements,
                                 int numInternalElements,
                                 int spatialOrder,
                                 int stencilSize = -1,
                                 int weightOrder = 1);
 
    ~WmUnstructuredReconstruction();
 
    int getSpatialOrder() const
    {
        return _spatialOrder;
    }
 
    int getLSize() const
    {
        return _lSize;
    }
 
    int getLInverseSize() const
    {
        return _lSize;
    }
 
    int getBSize() const
    {
        return _bSize;
    }
 
    int getStencilSize() const
    {
        return _stencilSize;
    }
 
    int getWeightsSize() const
    {
        return _wSize;
    }
 
    int getNumDifferentials() const
    {
        return _numDifferentials;
    }
 
    const real* getB() const
    {
        return _b;
    }
 
    const int* getStencils() const
    {
        return _stencils;
    }
 
    const real* getStencilsDxs() const
    {
        return _stencildxs;
    }
 
    const real* getWeights() const
    {
        return _w;
    }
 
    const real* getLInverse() const
    {
        return _lInverse;
    }
 
protected:
    void generateStencils(const int startIndex, const int endIndex);
 
    void generateB(const int startIndex, const int endIndex);
 
    void generateWeights(const int startIndex, const int endIndex);
 
    void generateL(const int startIndex, const int endIndex);
 
    void invertL(const int startIndex, const int endIndex);
 
    void setBComponent(const int elementIndex, const int* index, const real& val);
 
    real getBComponent(const int elementIndex, const int* index) const;
 
    void setLComponent(const int elementIndex, const int* index, const real& val);
 
    real getLComponent(const int elementIndex, const int* index) const;
 
    void setLInverseComponent(const int elementIndex, const int* index, const real& val);
 
    real getLInverseComponent(const int elementIndex, const int* index) const;
 
    void constructStencil(const int elementIndex,
                          const int stencilLength,
                          int* stencil,
                          real* stencildxs);
 
    int convertIndices(const int* indices) const;
 
    const WmUnstructuredPatch& _patch;
 
    int _weightOrder;
 
    int _numDimsPerPrimitive;
 
    int _numDimsPerNode;
 
    int _numReconstructionElements;
 
    int _numInternalElements;
 
    int _wSize;
 
    int _bSize;
 
    int _lSize;
 
    int _spatialOrder;
 
    int _stencilSize;
 
    int _numDifferentials;
 
    WmUnstructuredStencilConstructor* _usc;
 
    real* _lInverse;
 
    real* _l;
 
    real* _b;
 
    real* _w;
 
    int* _stencils;
 
    real* _stencildxs;
 
private:
};
 
#endif // WMUNSTRUCTUREDRECONSTRUCTION_H_