070-RegularResample.cpp File Reference

[← Previous Example] [Next Example → 130-PerformanceCreator.cpp ]. More...

#include <fiber/bundle/Bundle.hpp>
#include <fiber/finit/FinitAPI.h>
#include <fiber/field/UniformCartesianArray.hpp>
#include <fiber/baseop/Resample.hpp>

Functions
int	main ()

Detailed Description

[← Previous Example] [Next Example → 130-PerformanceCreator.cpp ].

In FiberLib Tutorial

#include <fiber/bundle/Bundle.hpp>
#include <fiber/finit/FinitAPI.h>
#include <fiber/field/UniformCartesianArray.hpp>
#include <fiber/baseop/Resample.hpp>
 
using namespace Fiber;
using namespace Eagle;
using namespace Eagle::PhysicalSpace;
 
 
int     main()
{
        // Initialize I/O layers
        Finit();
 
RefPtr<Cache>&theGlobalCache = Cache::MemCache();
const   int kB = 1024, MB = kB*kB;
        theGlobalCache->setMaxCost( 3*MB );
 
BundlePtr       BP;
 
Grid&InputGrid = BP[1.0]["DemoGrid"];
 
Representation&myCartesianRepresentation = InputGrid.makeCartesianRepresentation(3);
        BP->bindTo("070_RegularSample_LowRes.f5");
 
//MultiIndex<3> DemoGridSize = {8,16,12};
 
MultiIndex<3> DemoGridSize = {3,2,2};
 
        Verbose(0) << "Creating low-resolution Grid of dimension " << DemoGridSize
                   << " bound to " << BP->boundUrl();
 
RefPtr<UniformCartesianArray>
        CoordsP = new UniformCartesianArray( {-1.0,-1.0,-1.0} , DemoGridSize, point{+1.0, +1.0, +1.0} );
 
RefPtr<MemBase> M = CoordsP;
        myCartesianRepresentation.setPositions( new Field( M ) );
 
//      myCartesianRepresentation.setPositions( new Field( CoordsP ) );
 
 
#if 1
Ref<MemArray<3, double>> ScalarData( DemoGridSize );
        myCartesianRepresentation[ "scalardata"] = new Field( ScalarData );
 
const UniformCartesianArray&Coords= *CoordsP;
 
        for(auto&value : ScalarData->myChunk()->get_vector())
                value = -1.0;
 
        for(auto I : ScalarData->Size() )
        {
        const point&P = Coords[I];
 
        double  x = P.x(), y = P.y(), z = P.z();
 
                (*ScalarData)[ I ] += x*x + y*y - z*z;
        }
#endif
 
 
#if     1
        Verbose(0) << "#################### CREATE OUTPUT GRID AND ASSOCIATE WITH FILE ###################";
        {
        BundlePtr       OutputBP;
        Grid&OutputGrid = OutputBP[0.1]["ResampledDemoGrid"];
 
                OutputBP->bindTo("070_RegularSample_HighRes.f5");
 
                Verbose(0) << "#################### RESAMPLE AND SAVE VIA DESTRUCTION ###################";
 
        RefPtr<Field>
                ResampledField = ResampleRegularGrid3D(OutputGrid, InputGrid,
                                                       "scalardata",
                                                       {2.0, 2.0, 2.0} );
        }
 
#endif
 
        Verbose(0) << "Demo done.";
        return true;
}