testParallelIntegralInvariantEstimator.cpp File Reference

#include <iostream>
#include <type_traits>
#include "DGtal/base/Common.h"
#include "DGtal/shapes/implicit/ImplicitBall.h"
#include "DGtal/shapes/GaussDigitizer.h"
#include "DGtal/topology/LightImplicitDigitalSurface.h"
#include "DGtal/topology/DigitalSurface.h"
#include "DGtal/graph/DepthFirstVisitor.h"
#include "DGtal/graph/GraphVisitorRange.h"
#include "DGtal/geometry/surfaces/estimation/IIGeometricFunctors.h"
#include "DGtal/geometry/surfaces/estimation/IntegralInvariantVolumeEstimator.h"
#include "DGtal/geometry/surfaces/estimation/ParallelIIEstimator.h"
#include "DGtal/kernel/domains/DomainSplitter.h"
#include "DGtal/io/boards/Board2D.h"

Include dependency graph for testParallelIntegralInvariantEstimator.cpp:

Go to the source code of this file.

Functions
bool	testCurvature2dP (double h)
int	main (int, char **)

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

Bastien DOIGNIES (basti.nosp@m.en.d.nosp@m.oigni.nosp@m.es@l.nosp@m.iris..nosp@m.cnrs.nosp@m..fr ) Laboratoire d'InfoRmatique en Image et Systèmes d'information - LIRIS (CNRS, UMR 5205), INSA-Lyon, France LAboratoire de MAthématiques - LAMA (CNRS, UMR 5127), Université de Savoie, France

Date

2014/06/26

Functions for testing class IntegralInvariantVolumeEstimator and IIGeometricFunctor.

This file is part of the DGtal library.

Definition in file testParallelIntegralInvariantEstimator.cpp.

Function Documentation

main()

int main	(	int	,
		char **
	)

Definition at line 196 of file testParallelIntegralInvariantEstimator.cpp.

{
    trace.beginBlock ( "Testing class ParrallelIIEstimator with IntegralInvariantVolumeEstimator in 2d" );
 
    bool res = testCurvature2dP( 0.4 );
    trace.emphase() << ( res ? "Passed." : "Error." ) << std::endl;
    trace.endBlock();
    return res ? 0 : 1;
}

References DGtal::Trace::beginBlock(), DGtal::Trace::emphase(), DGtal::Trace::endBlock(), testCurvature2dP(), and DGtal::trace.

testCurvature2dP()

bool testCurvature2dP

(

double

h

)

[exampleParallelII-type]

[exampleParallelII-type]

[exampleParallelII-construction]

[exampleParallelII-construction]

Definition at line 64 of file testParallelIntegralInvariantEstimator.cpp.

{
  typedef ImplicitBall<Z2i::Space> ImplicitShape;
  typedef GaussDigitizer<Z2i::Space, ImplicitShape> DigitalShape;
  typedef LightImplicitDigitalSurface<Z2i::KSpace,DigitalShape> Boundary;
  typedef DigitalSurface< Boundary > MyDigitalSurface;
  typedef DepthFirstVisitor< MyDigitalSurface > Visitor;
  typedef GraphVisitorRange< Visitor > VisitorRange;
  typedef VisitorRange::ConstIterator VisitorConstIterator;
 
  typedef functors::IICurvatureFunctor<Z2i::Space> MyIICurvatureFunctor;
  typedef IntegralInvariantVolumeEstimator< Z2i::KSpace, DigitalShape, MyIICurvatureFunctor > MyIICurvatureEstimator;
 
  //typedef AxisDomainSplitter<Z2i::Domain> Splitter;
  typedef RegularDomainSplitter<HyperRectDomain<Z2i::Space>> Splitter;
  typedef ParallelIIEstimator<MyIICurvatureEstimator, Splitter> MyIICurvatureEstimatorP;
 
  typedef MyIICurvatureEstimator::Quantity Value;
 
  static_assert(std::is_same_v<typename MyIICurvatureEstimator::Quantity, typename MyIICurvatureEstimatorP::Quantity>);
 
  double re = 10;
  double radius = 15;
 
  trace.beginBlock( "[PARALLEL] Shape initialisation ..." );
 
  ImplicitShape ishape( Z2i::RealPoint( 0, 0 ), radius );
  DigitalShape dshape;
  dshape.attach( ishape );
  dshape.init( Z2i::RealPoint( -20.0, -20.0 ), Z2i::RealPoint( 20.0, 20.0 ), h );
 
  Z2i::KSpace K;
  if ( !K.init( dshape.getLowerBound(), dshape.getUpperBound(), true ) )
  {
    trace.error() << "Problem with Khalimsky space" << std::endl;
    return false;
  }
 
  Z2i::KSpace::Surfel bel = Surfaces<Z2i::KSpace>::findABel( K, dshape, 10000 );
  Boundary boundary( K, dshape, SurfelAdjacency<Z2i::KSpace::dimension>( true ), bel );
  MyDigitalSurface surf ( boundary );
 
 
 
  trace.endBlock();
 
  trace.beginBlock( "Curvature estimator initialisation ...");
 
  // Visitor ranges are typically unique and single pass. We need
  // to create one for each estimator in this case
  VisitorRange range( new Visitor( surf, *surf.begin() ));
  VisitorConstIterator ibegin = range.begin();
  VisitorConstIterator iend = range.end();
  // Parallel iterations
  VisitorRange rangeP( new Visitor( surf, *surf.begin() ));
  VisitorConstIterator ibeginP = rangeP.begin();
  VisitorConstIterator iendP = rangeP.end();
 
  MyIICurvatureFunctor curvatureFunctor;
  curvatureFunctor.init( h, re );
 
  MyIICurvatureEstimator  curvatureEstimator (    curvatureFunctor);
  curvatureEstimator.attach( K, dshape );
  curvatureEstimator.setParams( re/h );
  curvatureEstimator.init( h, ibegin, iend );
 
  // Parallel version expects a domain splitter instance and a number of
  // threads. Subsequent arguments are forwarded to the underlying estimator.
  Splitter splitter;
  MyIICurvatureEstimatorP curvatureEstimatorP( splitter, 4, curvatureFunctor );
  curvatureEstimatorP.attach( K, dshape );
  curvatureEstimatorP.setParams( re/h );
  curvatureEstimatorP.init( h, ibeginP, iendP );
 
  trace.endBlock();
 
  trace.beginBlock( "Curvature estimator evaluation ...");
 
  std::vector< Value > results, resultsP;
  std::back_insert_iterator< std::vector< Value > > resultsIt ( results  );
  std::back_insert_iterator< std::vector< Value > > resultsItP( resultsP );
 
  curvatureEstimator .eval( ibegin , iend , resultsIt  );
  curvatureEstimatorP.eval( ibeginP, iendP, resultsItP );
 
  trace.endBlock();
 
  trace.beginBlock ( "Comparing results of integral invariant 2D curvature ..." );
  // Debugging
  Board2D board;
  board << dshape.getDomain();
  for(auto p: dshape.getDomain())
    if (dshape(p))
      board << p;
  board.saveSVG("debugging-para.svg");
 
  trace.info() << dshape.getDomain() << std::endl;
  trace.info() << "Splits: "<<std::endl;
  auto splits = splitter(dshape.getDomain(), 4);
  for(const auto &sd: splits)
    trace.info() << sd.domain << std::endl;
 
 
  auto rsize  = results.size();
  auto rsizeP = resultsP.size();
 
  if (rsize != rsizeP)
  {
    trace.error() << "Size mismatch between parallel and non-parallel versions: " << rsize << " / " << rsizeP;
    trace.endBlock();
    return false;
  }
 
  for ( unsigned int i = 0; i < rsize; ++i )
  {
    if (std::abs(results[i] - resultsP[i]) >= 1e-6)
    {
      trace.error() << "Result mismatch between parallel and non-parallel versions at linel " << i << ": " << results[i] << " / " << resultsP[i] << "\n";
      trace.endBlock();
      return false;
    }
  }
  trace.endBlock();
  return true;
}

References DGtal::DigitalSurface< TDigitalSurfaceContainer >::begin(), DGtal::Trace::beginBlock(), DGtal::Trace::endBlock(), DGtal::Trace::error(), DGtal::Surfaces< TKSpace >::findABel(), DGtal::Trace::info(), DGtal::KhalimskySpaceND< dim, TInteger >::init(), K, LibBoard::Board::saveSVG(), splitter, and DGtal::trace.

Referenced by main().