parallelIIShortcuts.cpp

Go to the documentation of this file.

#include <iostream>
#include "ConfigExamples.h"
 
// Helpers
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
 
#include "DGtal/helpers/Shortcuts.h"
#include "DGtal/helpers/ShortcutsGeometry.h"

 
using namespace DGtal;
 
// Using standard 3D digital space.
typedef Shortcuts<Z3i::KSpace> SH3;
typedef ShortcutsGeometry<Z3i::KSpace> SHG3;
 
int main()
{

    auto params = SH3::defaultParameters() | SHG3::defaultParameters();
    params( "polynomial", "goursat" )( "gridstep", 0.125 );
    auto implicit_shape  = SH3::makeImplicitShape3D  ( params );
    auto digitized_shape = SH3::makeDigitizedImplicitShape3D( implicit_shape, params );
    auto K         = SH3::getKSpace( params );
    auto surface   = SH3::makeDigitalSurface( digitized_shape, K, params );
    auto surfels   = SH3::getSurfelRange( surface, params );

    trace.info()<< "Input vol domain: "<< digitized_shape->getDomain() << std::endl;
    //Sequential
    trace.beginBlock("Single thread");
    auto curv      = SHG3::getIIMeanCurvatures( digitized_shape, surfels, params );
    trace.endBlock();
 
    //Parallel
    trace.beginBlock("4 threads on default axis");
    auto curv_par4  = SHG3::getIIMeanCurvatures( digitized_shape, surfels,
                                                params( "ii-thread-number", 4 ));
    trace.endBlock();
 
    //Parallel
    trace.beginBlock("8 threads on axis 0");
    auto curv_par8_0 = SHG3::getIIMeanCurvatures( digitized_shape, surfels,
                                                params( "ii-thread-number", 8 )
                                                ( "ii-split-axis", 0 ) );
    trace.endBlock();
 
    //Parallel
    trace.beginBlock("8 threads on axis 1");
    auto curv_par8_1  = SHG3::getIIMeanCurvatures( digitized_shape, surfels,
                                                params( "ii-thread-number", 8 )
                                                ( "ii-split-axis", 1 ) );
    trace.endBlock();
 
    //Parallel
    trace.beginBlock("8 threads on axis 2");
    auto curv_par8_2  = SHG3::getIIMeanCurvatures( digitized_shape, surfels,
                                                params( "ii-thread-number", 8 )
                                                ( "ii-split-axis", 2 ) );
    trace.endBlock();
 
    //Parallel
    trace.beginBlock("16 threads on axis 2");
    auto curv_par16  = SHG3::getIIMeanCurvatures( digitized_shape, surfels,
                                                 params( "ii-thread-number", 16 )
                                                 ( "ii-split-axis", 2 ) );
    trace.endBlock();
 
    return 0;
}
//                                                                           //