FemPlyReaderDelegate.h

// This file is a part of the OpenSurgSim project.
// Copyright 2014-2015, SimQuest Solutions Inc.
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//     http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.

#ifndef SURGSIM_PHYSICS_FEMPLYREADERDELEGATE_H
#define SURGSIM_PHYSICS_FEMPLYREADERDELEGATE_H

#include <array>
#include <memory>

#include "SurgSim/DataStructures/PlyReaderDelegate.h"
#include "SurgSim/Physics/Fem.h"

namespace SurgSim
{
namespace Physics
{

class FemPlyReaderDelegate : public SurgSim::DataStructures::PlyReaderDelegate
{
public:
    FemPlyReaderDelegate();

protected:
    // \return Name of the element (1/2/3D), which this delegate processes.
    virtual std::string getElementName() const = 0;

    bool registerDelegate(SurgSim::DataStructures::PlyReader* reader) override;
    bool fileIsAcceptable(const SurgSim::DataStructures::PlyReader& reader) override;

    virtual void endParseFile() = 0;

    virtual void* beginVertices(const std::string& elementName, size_t vertexCount);

    virtual void processVertex(const std::string& elementName) = 0;

    virtual void endVertices(const std::string& elementName);

    void* beginFemElements(const std::string& elementName, size_t elementCount);

    virtual void processFemElement(const std::string& elementName) = 0;

    void endFemElements(const std::string& elementName);

    void* beginMaterials(const std::string& elementName, size_t materialCount);

    void endMaterials(const std::string& elementName);

    void* beginBoundaryConditions(const std::string& elementName, size_t boundaryConditionCount);

    virtual void processBoundaryCondition(const std::string& elementName) = 0;

protected:
    struct Vertex6DData
    {
        double x = std::numeric_limits<double>::signaling_NaN();
        double y = std::numeric_limits<double>::signaling_NaN();
        double z = std::numeric_limits<double>::signaling_NaN();
        int64_t overrun1 = -1; 
        double thetaX = std::numeric_limits<double>::signaling_NaN();
        double thetaY = std::numeric_limits<double>::signaling_NaN();
        double thetaZ = std::numeric_limits<double>::signaling_NaN();
        int64_t overrun2 = -1; 
    } m_vertexData;

    bool m_hasBoundaryConditions;

    unsigned int m_boundaryConditionData;

    struct MaterialData
    {
        double massDensity = std::numeric_limits<double>::signaling_NaN();
        double poissonRatio = std::numeric_limits<double>::signaling_NaN();
        double youngModulus = std::numeric_limits<double>::signaling_NaN();
        int64_t overrun = -1; 
    } m_materialData;

    bool m_hasMaterial;
    bool m_hasPerElementMaterial;

    struct ElementData
    {
        unsigned int type = 0;   // “LinearBeam”, “CorotationalTetrahedron”…
        int64_t overrun1 = -1; 

        unsigned int* indices = nullptr;
        unsigned int vertexCount = 0;
        int64_t overrun2 = -1; 
        double massDensity = std::numeric_limits<double>::signaling_NaN();
        double poissonRatio = std::numeric_limits<double>::signaling_NaN();
        double youngModulus = std::numeric_limits<double>::signaling_NaN();
    } m_elementData;
    bool m_hasRotationDOF;
};

} // namespace Physics
} // namespace SurgSim

#endif // SURGSIM_PHYSICS_FEMPLYREADERDELEGATE_H