500_4D_Raster.cpp File Reference

[← Previous Example] [Next Example → 600_4DField.cpp ]. More...

#include <fiber/bundle/Bundle.hpp>
#include <fiber/finit/FinitAPI.h>
#include <fiber/field/ArrayRef.hpp>
#include <eagle/PhysicalSpace.hpp>
#include <memcore/Console.hpp>
#include <fiber/field/UniformCartesianArray.hpp>
#include <fiber/grid/types/RegularlyFragmentedGrid.hpp>

Functions
int	main ()

Detailed Description

[← Previous Example] [Next Example → 600_4DField.cpp ].

In FiberLib Tutorial

#include <fiber/bundle/Bundle.hpp>
#include <fiber/finit/FinitAPI.h>
 
#include <fiber/field/ArrayRef.hpp>
#include <eagle/PhysicalSpace.hpp>
#include <memcore/Console.hpp>
#include <fiber/field/UniformCartesianArray.hpp>
#include <fiber/grid/types/RegularlyFragmentedGrid.hpp>
 
 
using namespace Fiber;
using namespace std;
using namespace Eagle::PhysicalSpace;
 
 
// ######################################################
// LOAD
//#######################################################
int     main()
{
        Finit();
 
BundlePtr BP = new Bundle();
        BP->bindToNew("TestRaster4D.f5");
 
Grid&G = BP[0.0]["4Draster"];
Representation&Rep = G.makeCartesianRepresentation(3);
 
MultiIndex<3> NumberOfFragments = {1,1,1};
MultiIndex<4> FragDims4    = {17,23,1,5};
MultiIndex<4> FragDims             = {17,23,1};
 
        Rep[FIBER_POSITIONS] = new RegularlyFragmentedField<3> (NumberOfFragments);
RefPtr<RegularlyFragmentedField<3> > PositionsF = Rep(FIBER_POSITIONS);
 
        Rep["Field4D"] = new RegularlyFragmentedField<3> (*PositionsF);
RefPtr<RegularlyFragmentedField<3> > Field4DF = Rep("Field4D");
 
        Rep["Field3D"] = new RegularlyFragmentedField<3> (*PositionsF);
RefPtr<RegularlyFragmentedField<3> > Field3DF = Rep("Field3D");
 
RefPtr<UniformCartesianArray> Coordinates = new UniformCartesianArray(FragDims, {0,0,0}, {1,1,0});
RefPtr<MemArray<4,double > > Data4Ptr = new MemArray<4,double >(FragDims4);
RefPtr<MemArray<3,double > > Data3Ptr = new MemArray<3,double >(FragDims);
 
const RefPtr<FragmentID>&FragID = PositionsF->getFragmentID( MIndex(0,0,0) );
 
        PositionsF->setPersistentData(Coordinates,FragID);
        Field4DF->setPersistentData(Data4Ptr,FragID);
        Field3DF->setPersistentData(Data3Ptr,FragID);
 
RegularlyFragmentedGrid<3>::create( G, PositionsF);
 
        Console() << "Done";
 
        return 0;
}