065-F5moveSkeleton.cpp File Reference

[← Previous Example] [Next Example → 066-StorageTransformationWithPositions.cpp ]. More...

#include <algorithm>
#include <fiber/bundle/Bundle.hpp>
#include <fiber/bundle/StorageTransformations.hpp>
#include <fiber/finit/FinitAPI.h>
#include <memcore/Console.hpp>

Functions
int	main (int argc, const char **argv)

Detailed Description

[← Previous Example] [Next Example → 066-StorageTransformationWithPositions.cpp ].

In FiberLib Tutorial

 
#include <algorithm>
#include <fiber/bundle/Bundle.hpp>
#include <fiber/bundle/StorageTransformations.hpp>
#include <fiber/finit/FinitAPI.h>
#include <memcore/Console.hpp>
 
using namespace Fiber;
using namespace Eagle::PhysicalSpace;
 
int     main(int argc, const char**argv)
{
        Finit();
 
string InFile = "065_F5moveSkeleton_InputFile.f5";
 
        if (argc>1)
        {
        BundlePtr BP = new Bundle() ;
                BP->bindTo(InFile);
                return 0;
        }
 
        {
        BundlePtr BPIn =  new Bundle ;
        Grid&myGrid = BPIn[0.0]["TestGrid"];
        RefPtr<Field> PositionsF = myGrid[FIBER_POSITIONS];
 
 
        Ref<MemArray<1,point> > Mem(100);
 
                std::fill(Mem->myChunk()->get_vector().begin(),
                          Mem->myChunk()->get_vector().end(),
                          point{1.,0.,2.} );
 
                Console() << "Mem size: " << Mem->nElements();
 
                PositionsF->setPersistentData(Mem);
 
                {
                Ref<MemArray<1,double>> Mem(100);
                        std::fill(Mem->myChunk()->get_vector().begin(),
                                  Mem->myChunk()->get_vector().end(), 42.0);
 
                        myGrid["Values"]->setPersistentData(Mem);
                }
 
                BPIn->save(InFile);
                Console() << "Save file " << InFile;
        }
 
string OutFile = "065_F5moveSkeleton_OutFile.f5";
        {
        BundlePtr BP = new Bundle() ;
                BP->bindToNew(OutFile);
 
 
 
 
        bool success;
        RefPtr<LoaderParameters> ld = new LoaderParameters(InFile);
                BP = Bundle::load(success, ld , BP  );
        BP->sync();
 
 
        RefPtr<Grid> myGrid = BP[0.0]("TestGrid");
 
        if (!myGrid)
        {
                Console() <<"No Grid" ;
                return 1;
        }
 
        const int PointSkeletonDims = 3;
 
        /* change root skeleton-id (e.g point) to refinement
           skeleton (that is: index depth, dims AND REFINEMENTLEVEL!) */
        SkeletonID mySkeletonIDCoarse(PointSkeletonDims,0);
        RefPtr<Skeleton> mySkeletonFine = (*myGrid)(mySkeletonIDCoarse);
 
                if (!mySkeletonFine)
                {
                        Console() << "Cannot find root (finest)  skeleton" ;
                        return false;
                }
                mySkeletonFine->setIsotropicRefinement(10);
 
                /* note: use index depth 2 for fragments instead of
                   index depth 0 which correspond incorrectly to point*/
 
        }
 
 
 
 
        return 0;
}