wmuhofvgeometry.h

Go to the documentation of this file.

#ifndef WMUHOFVGEOMETRY_H
#define WMUHOFVGEOMETRY_H
 
// Wm includes
#include "lib/wmunstructuredpatch.h"
#include "lib/wmunstructuredgeometry.h"
#include "compute/wmcomputemachine.h"
 
class WmUHOFVGeometry
{
public:
    WmUHOFVGeometry(const WmUnstructuredPatch& patch,
                    const int facePolynomialOrder,
                    const int volumePolynomialOrder);
 
    ~WmUHOFVGeometry();
 
    int getNumElements() const
    {
        return _numElements;
    }
 
    int getElementGeometrySizeInBytes() const
    {
        return _numElements * _elementGeometryBlockSize * sizeof(real);
    }
 
    int getElementGeometryBlockSize() const
    {
        return _elementGeometryBlockSize;
    }
    int getElementGeometryFaceBlockSize() const
    {
        return _elementFaceGeometryBlockSize;
    }
    int getElementGeometryInternalBlockSize() const
    {
        return _elementInternalGeometryBlockSize;
    }
 
    int getNumFaceQuadraturePoints() const
    {
        return _numFaceQuadraturePoints;
    }
    int getNumVolumeQuadraturePoints() const
    {
        return _numVolumeQuadraturePoints;
    }
 
    const cl_mem* getElementGeometryMemPtr() const
    {
        return &_mem_elementGeometry;
    }
    const cl_mem* getVertexDxGeometryMemPtr() const
    {
        return &_mem_vertexDxs;
    }
 
    const real* getElementGeometry() const;
 
protected:
    void generateGeometryRange(const int startIndex, const int endIndex);
    void generateReconstructionGeometryRange(const int startIndex, const int endIndex);
 
    void getDxOnElement(const real centroid[3],
                        const real v_tC[][3],
                        const real normalizedPosition[3],
                        real dx[3]) const;
    void getDxOnFace(const real centroid[3],
                     const real f_tC[][3],
                     const real normalizedFacePosition[3],
                     real dx[3]) const;
 
    real generateFaceJacobian(const real f_tC[][3],
                              const real* normalizedFacePosition) const;
    real generateVolumeJacobian(const real v_tC[][3],
                                const real* normalizedPosition) const;
 
    void convertFaceTransformCoordinates(const real v_tC[][3],
                                         const int fIndex,
                                         real f_tC[][3]) const;
 
    void printElement(const int i) const;
 
    const WmUnstructuredPatch& _patch;
    const WmUnstructuredGeometry& _uG;
 
    int _numElements;
    int _numReconstructionElements;
 
    int _elementFaceGeometryBlockSize;
    int _elementGeometryBlockSize;
    int _elementInternalGeometryBlockSize;
 
    real* _elementGeometry;
    real* _vertexDxs;
 
    int _facePolynomialOrder;
    int _volumePolynomialOrder;
 
    int _numFaceQuadraturePoints;
    int _numVolumeQuadraturePoints;
 
    cl_mem _mem_elementGeometry;
    cl_mem _mem_vertexDxs;
 
private:
    //    WmUHOFVGeometry(){}
};
#endif // WMUHOFVGEOMETRY_H