061-RegularFragmentTopologyWithHigherTopology.cpp File Reference

[← Previous Example] [Next Example → 065-F5moveSkeleton.cpp ]. More...

#include <fiber/bundle/Bundle.hpp>
#include <fiber/finit/FinitAPI.h>
#include <fiber/field/ArrayRef.hpp>
#include <fiber/field/RegularlyFragmentedField.hpp>
#include <fiber/field/UniformCartesianArray.hpp>
#include <fiber/grid/types/RegularlyFragmentedGrid.hpp>
#include <fiber/field/Cell.hpp>

Functions
int	main ()

Detailed Description

[← Previous Example] [Next Example → 065-F5moveSkeleton.cpp ].

In FiberLib Tutorial

#include <fiber/bundle/Bundle.hpp>
#include <fiber/finit/FinitAPI.h>
#include <fiber/field/ArrayRef.hpp>
#include <fiber/field/RegularlyFragmentedField.hpp>
#include <fiber/field/UniformCartesianArray.hpp>
#include <fiber/grid/types/RegularlyFragmentedGrid.hpp>
#include <fiber/field/Cell.hpp>
 
 
using namespace Fiber;
using namespace Eagle::PhysicalSpace;
 
int     main()
{
        // Initialize I/O layers (required for standalone programs)
        Finit();
 
BundlePtr      BP;
Grid&myGrid = BP[0.0]["DemoGrid"]; 
 
        //BP->bindTo("061_RegularlyFragmentedFieldWithHigherTopo.f5");
 
 
SkeletonID VertexSkeletonID(3,0); // Dimension 1 and index depth 0
 
Skeleton&VertexSkeleton = myGrid[VertexSkeletonID];
 
RefPtr<Chart> Chart = myGrid.makeCartesianChart();
 
Representation&VertexRep = VertexSkeleton[Chart];
 
 
 
SkeletonID CellSkeletonID(3,1); // Dimension 2 and index depth 1
 
Skeleton&CellSkeleton = myGrid[CellSkeletonID];
 
 
 
Representation&RepCellToVertex = CellSkeleton[VertexSkeleton];
 
MultiIndex<3> Fragments = {2,1,1};
 
RefPtr<RegularlyFragmentedField<3> > VertexF  = RegularlyFragmentedField<3>::newField(Fragments, RefPtr<Field>() );
 
        RegularlyFragmentedGrid<3>::create( myGrid, VertexF);
 
 
 
//RefPtr<RegularlyFragmentedField<3> > VertexF = new RegularlyFragmentedField<3> ( Fragments );
        VertexRep [FIBER_POSITIONS] = VertexF;
 
 
 
 
//RefPtr<Field> VertexF = VertexRep[FIBER_POSITIONS];
 
 
 
point start = {0,0,0};
point end   = {2,1,0};
MultiIndex<3> dims = {3,2,1};
 
RefPtr<UniformCartesianArray> VertexMem1 = new UniformCartesianArray(start, dims, end);
 
RefPtr<FragmentID> ID1 = VertexF->getFragmentID ({0,0,0});
 
        VertexF->setPersistentData(VertexMem1,ID1);
 
 
 
        start = {2,1,0};
        end   = {4,2,0};
        dims = {3,2,1};
 
RefPtr<UniformCartesianArray> VertexMem2 = new UniformCartesianArray(start, dims, end);
 
RefPtr<FragmentID> ID2 = VertexF->getFragmentID ({1,0,0});
 
        VertexF->setPersistentData(VertexMem2,ID2);
 
 
#if 1
 
RefPtr<RegularlyFragmentedField<3> > CellToVertexF = new RegularlyFragmentedField<3> ( Fragments );
        RepCellToVertex [FIBER_POSITIONS] = CellToVertexF;
 
        RegularlyFragmentedGrid<3>::create( myGrid, CellToVertexF);
 
//RefPtr<Field> CellToVertexF = RepCellToVertex["Temperature"];
 
RefPtr<MemArray<3,float > > CellMem1 = new MemArray<3,float > ({2,1,1});
 
        (*CellMem1)[{0,0,0}] = 10.2;
        (*CellMem1)[{1,0,0}] = 11.2;
 
RefPtr<FragmentID> IDa = CellToVertexF->getFragmentID ({0,0,0});
 
        CellToVertexF->setPersistentData(CellMem1,IDa);
 
 
RefPtr<MemArray<3,float > > CellMem2 = new MemArray<3,float > ({2,1,1});
 
        (*CellMem2)[{0,0,0}] = 20.2;
        (*CellMem2)[{1,0,0}] = 21.2;
 
RefPtr<FragmentID> IDb = CellToVertexF->getFragmentID ({1,0,0});
 
        CellToVertexF->setPersistentData(CellMem2,IDb);
 
#endif
        BP->save("061_RegularlyFragmentedFieldWithHigherTopo.f5");
 
        return 0;
}