020-CreateBundleAndSave.cpp File Reference

[← Previous Example] [Next Example → 021_CreateMultiBundlesAndSave.cpp ]. More...

#include <fiber/bundle/Bundle.hpp>
#include <fiber/finit/FinitAPI.h>
#include <fiber/field/ArrayRef.hpp>
#include <eagle/PhysicalSpace.hpp>

Functions
int	main ()

Detailed Description

[← Previous Example] [Next Example → 021_CreateMultiBundlesAndSave.cpp ].

In FiberLib Tutorial

#include <fiber/bundle/Bundle.hpp>
#include <fiber/finit/FinitAPI.h>
 
#include <fiber/field/ArrayRef.hpp>
#include <eagle/PhysicalSpace.hpp>
 
 
using namespace Fiber;
using namespace Eagle::PhysicalSpace;
 
 
int     main()
{
        // Initialize I/O layers
        Finit(); 
 
 
        // 
        // Create a bundle object with a grid at T=1.0, named "myGrid", in 
        // a three-dimensional cartesian representation, using the default 
        // coordinate system.
        //
BundlePtr      BP;
Representation&myCartesianRepresentation = BP[1.0]["myGrid"].makeCartesianRepresentation(3); 
 
        //
        // Define the size of the grid to be used hereforth 
        //
const MultiIndex<3> GridDims = MIndex(11,13,17);
 
        // 
        // Create a coordinate field and fill it with values 
        //
RefPtr<Field>&Positions = myCartesianRepresentation[ FIBER_POSITIONS ];
        Positions = new Field();
 
        {
        ArrayRef<point, 3>
                Coordinates(GridDims); 
 
                for(MultiIndex<3> I : GridDims)
                {
                point&P = Coordinates[ I ]; 
                        P.x() = I[0]*0.3; 
                        P.y() = I[1]*0.2; 
                        P.z() = I[2]*0.1;
                } 
 
                Positions->setPersistentData(Coordinates);
        } 
        Assert( Positions->isUnfragmented() );
 
 
        // 
        // Define a scalar field called "values" and set numbers 
        // to it using some formula based on the indices, just to 
        // get something non-trivial here.
        //
RefPtr<Field>&Values = myCartesianRepresentation["values"];
        Values = new Field();
        {
        ArrayRef<double, 3>
                GridValues(GridDims); 
 
                for(MultiIndex<3> I : GridDims)
                {
                double& value = GridValues[ I ]; 
                        value = (I[0]*0.3)*(I[0]*0.3) +
                                (I[1]*0.2)*(I[0]*0.2) - 
                                (I[2]*0.1)*(I[2]*0.1);
                }
 
                Values->setPersistentData( GridValues );
        } 
        Assert( Values->isUnfragmented() );
 
        // 
        // Save the bundle to a file. 
        //
        BP.save("020_CreateBundleAndSave.f5");
        return 0;
}