parallelIIShortcuts.cpp File Reference

#include <iostream>
#include "ConfigExamples.h"
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/helpers/Shortcuts.h"
#include "DGtal/helpers/ShortcutsGeometry.h"
#include "DGtal/io/viewers/PolyscopeViewer.h"
#include "DGtal/io/colormaps/GradientColorMap.h"

Include dependency graph for parallelIIShortcuts.cpp:

Go to the source code of this file.

Typedefs
typedef Shortcuts< Z3i::KSpace >	SH3
typedef ShortcutsGeometry< Z3i::KSpace >	SHG3

Functions
int	main ()

Detailed Description

This program is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.

Author

David Coeurjolly (david.nosp@m..coe.nosp@m.urjol.nosp@m.ly@c.nosp@m.nrs.f.nosp@m.r) Laboratoire d'InfoRmatique en Image et Systèmes d'information - LIRIS (CNRS, UMR 5205), INSA-Lyon, France

Date

2026/06/08

This file is part of the DGtal library.

Definition in file parallelIIShortcuts.cpp.

Typedef Documentation

SH3

typedef Shortcuts<Z3i::KSpace> SH3

Definition at line 52 of file parallelIIShortcuts.cpp.

SHG3

typedef ShortcutsGeometry<Z3i::KSpace> SHG3

Definition at line 53 of file parallelIIShortcuts.cpp.

Function Documentation

main()

int main

(

void

)

[Parallel-instantiation]

[Parallel-instantiation]

[Parallel-run]

[Parallel-run]

Definition at line 56 of file parallelIIShortcuts.cpp.

{
 
    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;
}

References DGtal::Display3D< Space, KSpace >::addQuantity(), DGtal::Trace::beginBlock(), DGtal::Shortcuts< TKSpace >::defaultParameters(), DGtal::ShortcutsGeometry< TKSpace >::defaultParameters(), DGtal::Display3D< Space, KSpace >::draw(), DGtal::Trace::endBlock(), DGtal::ShortcutsGeometry< TKSpace >::getIIMeanCurvatures(), DGtal::Shortcuts< TKSpace >::getKSpace(), DGtal::Shortcuts< TKSpace >::getSurfelRange(), DGtal::Trace::info(), K, DGtal::Shortcuts< TKSpace >::makeDigitalSurface(), DGtal::Shortcuts< TKSpace >::makeDigitizedImplicitShape3D(), DGtal::Shortcuts< TKSpace >::makeImplicitShape3D(), DGtal::PolyscopeViewer< Space, KSpace >::show(), splitter, surface, and DGtal::trace.