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