BundleBaseSpaceFunctions.hpp

#ifndef __FIBER_BASEOP_BUNDLEFUNCTIONS_HPP
#define __FIBER_BASEOP_BUNDLEFUNCTIONS_HPP
 
#include "gridopDllApi.h"
 
#include <grid/Grid.hpp>
 
#include <grid/CartesianChart.hpp>
#include <bundle/Slice.hpp>
 
namespace Fiber
{
 
inline MemCore::RefPtr<Field> GetCartesianPositions( MemCore::RefPtr<Grid> grid )
{
        if (!grid)
                return NullPtr();
 
        return grid->CartesianPositions(); 
/*
        MemCore::RefPtr<Skeleton> skel = grid->findVertices(); 
        if(!skel){ puts("GetCartesianPositions() no skeleton found"); return NullPtr();} 
 
        MemCore::RefPtr< Chart > chart = grid->makeChart( typeid(Fiber::CartesianChart3D) ); 
        MemCore::RefPtr<Representation> rep = (*skel)( chart ); 
        if(!rep){ puts("GetCartesianPositions() no representation found"); return NullPtr();} 
 
        return (*rep)("Positions"); 
*/
}
 
 
inline MemCore::RefPtr<Field> GetFragments( MemCore::RefPtr<Grid> grid )
{
        if (!grid)
                return NullPtr();
 
MemCore::RefPtr<Skeleton> source_skel = grid->findVertices(); 
        if(!source_skel){ puts("GetFragments() no skeleton found"); return NullPtr();} 
 
MemCore::RefPtr<Skeleton> skel = (*grid)( SkeletonID(3, 2) ); 
        if(!skel){ puts("GetFragments() no skeleton found"); return NullPtr();} 
 
MemCore::RefPtr<Representation> rep = (*skel)( source_skel ); 
        if(!rep){ puts("GetFragments() no representation found"); return NullPtr();} 
 
        assert(rep); 
        return (*rep)("Positions"); 
}
 
inline MemCore::RefPtr<Field> GetFragmentsBBs( MemCore::RefPtr<Grid> grid )
{
        if (!grid)
                return NullPtr();
 
        MemCore::RefPtr<Skeleton> skel = (*grid)( SkeletonID(3, 2) ); 
        if(!skel){ puts("GetFragmentsBBs() no skeleton found"); return NullPtr();} 
 
        MemCore::RefPtr< Chart > chart = grid->makeChart( typeid(CartesianChart3D) ); 
        MemCore::RefPtr<Representation> rep = (*skel)( chart ); 
        if(!rep){ puts("GetFragmentsBBs() no representation found"); return NullPtr();} 
 
        assert(rep); 
        return (*rep)("Positions"); 
}
 
 
inline MemCore::RefPtr<Field> GetUniGridMappers( MemCore::RefPtr<Slice> slice, MemCore::RefPtr<Grid> source_grid )
{
        if (!slice)
                return NullPtr();
 
        MemCore::RefPtr< Grid > unigrid = (*slice)( "UniMappers" ); 
        if(!unigrid){ puts("GetUniGridMappers() no grid named 'UniMappers' found!"); return NullPtr(); } 
        MemCore::RefPtr<Skeleton> uniskel = (*unigrid)( SkeletonID(3, 1) ); 
        if(!uniskel){ puts("GetUniGridMappers() no skel(dim 3,index depth 1) in unigrid found!"); return NullPtr(); } 
        MemCore::RefPtr<Skeleton> fragskel = (*source_grid)( SkeletonID(3, 2) ); 
        if(!fragskel){ puts("GetUniGridMappers() no skel(dim 3,index depth 2) in source_grid found!"); return NullPtr(); } 
        MemCore::RefPtr<Representation> rep = (*fragskel)( uniskel );
        if(!rep){ puts("GetUniGridMappers() no representation for unimappers found!"); return NullPtr(); } 
 
        return (*rep)("Positions");
}
 
} // namespace Fiber
 
 
#endif // __BUNDLEFUNCTIONS_HPP