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