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"
 
 
#ifdef DGTAL_WITH_POLYSCOPE_VIEWER
// Visualization
#include "DGtal/io/viewers/PolyscopeViewer.h"
#include "DGtal/io/colormaps/GradientColorMap.h"
#endif
 
 
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();
 
#ifdef DGTAL_WITH_POLYSCOPE_VIEWER
    PolyscopeViewer viewer;
 
    std::string objectName = "Surfels";
    viewer.draw(surfels, objectName); // Draws the object independently
    viewer.addQuantity(objectName, "Mean curvature", curv_par8_0);
 
    AxisDomainSplitter<Z3i::Domain> splitter(0);
    AxisDomainSplitter<Z3i::Domain>::SplitDomainsInfo splits = splitter(digitized_shape->getDomain(), 8);
    HueShadeColorMap<unsigned int> cmap(0,(unsigned int)splits.size());
    for(auto i=0; i< splits.size(); ++i)
    {
        viewer << cmap(i);
        viewer << splits[i].domain;
    }
    viewer.show();
#endif
 
    return 0;
}
//                                                                           //