gsPeriodicParametrization.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsModeling/gsParametrization.h>
 
namespace gismo
{
 
template <class T>
class gsPeriodicParametrization : public gsParametrization<T>
{
 
public:
 
    typedef memory::shared_ptr<gsPeriodicParametrization<T> > uPtr;
 
    class FlatMesh
    {
        typedef typename gsMesh<T>::VertexHandle VertexHandle;
    public:
 
        FlatMesh(const gsMesh<T>& unfolded)
            : m_unfolded(unfolded)
        {}
 
        gsMesh<T> createRestrictedFlatMesh() const;
 
    protected:
 
        real_t correspondingV(const VertexHandle& v0,
                              const VertexHandle& v1,
                              real_t u) const;
 
        void addThreeFlatTrianglesOneOut(gsMesh<T>& mesh,
                                         const VertexHandle& v0,
                                         const VertexHandle& v1,
                                         const VertexHandle& v2) const;
 
        void addThreeFlatTrianglesTwoOut(gsMesh<T>& mesh,
                                         const VertexHandle& v0,
                                         const VertexHandle& v1,
                                         const VertexHandle& v2) const;
 
        void addOneFlatTriangleNotIntersectingBoundary(gsMesh<T>& mesh,
                                                       const VertexHandle& v0,
                                                       const VertexHandle& v1,
                                                       const VertexHandle& v2) const;
 
    protected: // members
        gsHalfEdgeMesh<T> m_unfolded;
    };
 
public:
 
    gsPeriodicParametrization(const gsMesh<T>& mesh,
                              const gsMatrix<T>& verticesV0,
                              const gsMatrix<T>& paramsV0,
                              const gsMatrix<T>& verticesV1,
                              const gsMatrix<T>& paramsV1,
                              const gsOptionList &list = gsParametrization<T>::defaultOptions())
        : gsParametrization<T>(mesh, list),
          m_paramsV0(paramsV0), m_paramsV1(paramsV1),
          m_indicesV0(this->indices(verticesV0)),
          m_indicesV1(this->indices(verticesV1))
    {
        GISMO_ASSERT(this->m_paramsV0.rows() == 1, "one row expected in paramsV0");
        GISMO_ASSERT(this->m_paramsV1.rows() == 1, "one row expected in paramsV1");
 
        GISMO_ASSERT(this->m_paramsV0.cols() == static_cast<index_t>(m_indicesV0.size()),
                     "paramsV0 and valuesV0 are required to have the same number of cols");
        GISMO_ASSERT(this->m_paramsV1.cols() == static_cast<index_t>(m_indicesV1.size()),
                     "paramsV1 and valuesV1 are required to have the same number of cols");
    }
 
    void restrictMatrices(gsMatrix<T>& uv, const gsMatrix<T>& xyz,
                          real_t uMin = 0, real_t uMax = 1) const;
 
    virtual void compute() = 0;
 
protected:
 
    void initParameterPoints();
 
protected: // members
 
    const gsMatrix<T> m_paramsV0;         
    const gsMatrix<T> m_paramsV1;         
    const std::vector<size_t> m_indicesV0; 
    const std::vector<size_t> m_indicesV1; 
 
};
 
} // namespace gismo
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsPeriodicParametrization.hpp)
#endif