ResInsight/ApplicationCode/GeoMech/GeoMechVisualization/RivFemPartGeometryGenerator.h

#pragma once

#include "cvfBase.h"
#include "cvfObject.h"
#include "cvfArray.h"

#include "RigFemPart.h"
#include "cvfStructGrid.h"

namespace cvf
{

class DrawableGeo;
class ScalarMapper;

}

class RigFemPartScalarDataAccess;

class RivFemPartTriangleToElmMapper : public cvf::Object
{
public:
    size_t              triangleCount() const               { return m_trianglesToElementIndex.size();}

    int                 elementIndex(size_t triangleIdx) const { return m_trianglesToElementIndex[triangleIdx]; }
    char                elementFace(size_t triangleIdx) const  { return m_trianglesToElmFace[triangleIdx]; }

    // Interface for building the mappings
    std::vector<int>&   triangleToElmIndexMap() { return m_trianglesToElementIndex; }
    std::vector<char>&  triangleToElmFaceMap()  { return m_trianglesToElmFace; }
 
private:
    std::vector<int>    m_trianglesToElementIndex;
    std::vector<char>   m_trianglesToElmFace;
};

//==================================================================================================
//
// 
//
//==================================================================================================
class RivFemPartGeometryGenerator : public cvf::Object
{
public:
    explicit RivFemPartGeometryGenerator(const RigFemPart* part);
    ~RivFemPartGeometryGenerator();

    // Setup methods

    void                        setElementVisibility(const cvf::UByteArray* cellVisibility);

    // Access, valid after generation is done

    const RigFemPart*           activePart() { return m_part.p(); }

    // Generated geometry

    cvf::ref<cvf::DrawableGeo>  generateSurface();
    cvf::ref<cvf::DrawableGeo>  createMeshDrawable();
    cvf::ref<cvf::DrawableGeo>  createOutlineMeshDrawable(double creaseAngle);

    const std::vector<size_t>&  quadVerticesToNodeIdxMapping() const   { return m_quadVerticesToNodeIdx;}
    const std::vector<size_t>&  quadVerticesToGlobalElmNodeIdx() const { return m_quadVerticesToGlobalElmNodeIdx;}
    const std::vector<size_t>&  quadVerticesToGlobalElmFaceNodeIdx() const { return m_quadVerticesToGlobalElmFaceNodeIdx; }

    RivFemPartTriangleToElmMapper* triangleToElementMapper() { return m_triangleMapper.p();}

    static cvf::ref<cvf::DrawableGeo> createMeshDrawableFromSingleElement(const RigFemPart* grid, size_t elementIndex);

private:
    static cvf::ref<cvf::UIntArray> 
                                lineIndicesFromQuadVertexArray(const cvf::Vec3fArray* vertexArray);
    void                        computeArrays();

private:
    // Input
    cvf::cref<RigFemPart>       m_part;                     // The part being processed
    cvf::cref<cvf::UByteArray>  m_elmVisibility;

    // Created arrays
    cvf::ref<cvf::Vec3fArray>   m_quadVertices;
    //cvf::ref<cvf::Vec3fArray> m_triangleVertices; // If needed, we will do it like this, I think
    std::vector<size_t>         m_quadVerticesToNodeIdx;
    std::vector<size_t>         m_quadVerticesToGlobalElmNodeIdx;
    std::vector<size_t>         m_quadVerticesToGlobalElmFaceNodeIdx;

        // Mappings
    cvf::ref<RivFemPartTriangleToElmMapper> m_triangleMapper;

};