Fem.h

// This file is a part of the OpenSurgSim project.
// Copyright 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_FEM_H
#define SURGSIM_PHYSICS_FEM_H

#include "SurgSim/DataStructures/Vertices.h"
#include "SurgSim/Framework/Asset.h"
#include "SurgSim/Physics/FemElementStructs.h"

namespace SurgSim
{
namespace Physics
{

template <class VertexData, class Element>
class Fem : public SurgSim::DataStructures::Vertices<VertexData>, public SurgSim::Framework::Asset,
        public std::enable_shared_from_this<Fem<VertexData, Element>>
{
public:
    Fem();

    size_t addElement(std::shared_ptr<Element> element);

    size_t getNumElements() const;

    const std::vector<std::shared_ptr<Element>>& getElements() const;

    std::vector<std::shared_ptr<Element>>& getElements();

    std::shared_ptr<Element> getElement(size_t id) const;

    size_t addBoundaryCondition(size_t boundaryCondition);

    const std::vector<size_t>& getBoundaryConditions() const;

    std::vector<size_t>& getBoundaryConditions();

    size_t getBoundaryCondition(size_t id) const;

protected:
    template <class PlyType, class FemType>
    bool loadFemFile(const std::string& filename);

    std::vector<std::shared_ptr<Element>> m_elements;

    std::vector<size_t> m_boundaryConditions;
};

} // namespace Physics
} // namespace SurgSim

#include "SurgSim/Physics/Fem-inl.h"

#endif // SURGSIM_PHYSICS_FEM_H