TangentialVectors.hpp

#ifndef __BASEOP_TANGENTIALVECTORS_HPP
#define __BASEOP_TANGENTIALVECTORS_HPP

#include "gridopDllApi.h"
#include <grid/Representation.hpp>
#include <grid/Grid.hpp>
#include <eagle/PhysicalSpace.hpp>


namespace Fiber
{

extern gridop_API const char TangentialVectorFieldName[];

extern gridop_API const char CurveParameterFieldName[];

extern gridop_API RefPtr<MemArray<1, Eagle::tvector3> >
        ComputeTangentialVectors(const MemArray<1, Eagle::point3>&Vertices,
                                 const MemArray<1, std::vector<index_t> >&EdgesArray,
                                 const MemBase::Creator_t&Crec);

extern gridop_API RefPtr<MemArray<1, Eagle::tvector3> >
        getTangentialVectors(Grid&G, const string&FieldName = TangentialVectorFieldName);


struct  gridop_API TangentialVectorsAndCurveParameter
{
        RefPtr<MemArray<1, Eagle::tvector3> >   TangentialVectors;
        double                                  max_len;
        RefPtr<MemArray<1, double > >           CurveParameter;

        void compute(const MemArray<1, Eagle::point3>&Vertices,
                     const MemArray<1, std::vector<index_t> >&EdgesArray,
                     const MemBase::Creator_t&Crec = NullPtr() );


        TangentialVectorsAndCurveParameter(const MemArray<1, Eagle::point3>&Vertices,
                                           const MemArray<1, std::vector<index_t> >&EdgesArray)
        {
                compute(Vertices, EdgesArray);
        } 

        struct  gridop_API Exception : std::exception
        {}; 

        struct  gridop_API NoEdges : Exception {};
        struct  gridop_API NoVertices : Exception {};
        struct  gridop_API NoVertexField : Exception {}; 
        struct  gridop_API NoEdgesAsVertices : NoEdges {};
        struct  gridop_API NoEdgePositions : NoEdges {}; 
        struct  gridop_API NoVertexFieldCreator : NoVertices {};
        struct  gridop_API NoVertexFieldData : NoVertices {};

        TangentialVectorsAndCurveParameter(Grid&G, 
                                           const string&TangentFieldName = TangentialVectorFieldName,
                                           const string&LengthFieldName = CurveParameterFieldName);
};


} // namespace Fiber

#endif // __BASEOP_TANGENTIALVECTORS_HPP