#include <iostream>
#include <vector>
#include <algorithm>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/SpaceND.h"
#include "DGtal/geometry/tools/GenericLatticeConvexHull.h"
#include "DGtal/geometry/tools/AffineGeometry.h"
#include "DGtalCatch.h"

Include dependency graph for testGenericLatticeConvexHull.cpp:

Functions
std::mt19937	g (rd())

template<typename Point >
std::vector< Point >	randomPointsInBall (int nb, int R)

template<typename Point >
std::vector< Point >	makeRandomLatticePointsFromDirVectors (int nb, const vector< Point > &V)

	SCENARIO ("GenericLatticeConvexHull< ConvexHullIntegralKernel< 2 > > unit tests", "[genquickhull][integral_kernel][2d]")

	SCENARIO ("GenericLatticeConvexHull< ConvexHullIntegralKernel< 3 > > unit tests", "[genquickhull][integral_kernel][3d]")

Variables
std::random_device	rd

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: Jacques-Olivier Lachaud (jacqu.nosp@m.es-o.nosp@m.livie.nosp@m.r.la.nosp@m.chaud.nosp@m.@uni.nosp@m.v-sav.nosp@m.oie..nosp@m.fr ) Laboratory of Mathematics (CNRS, UMR 5127), University of Savoie, France

Date: 2025/09/04

Functions for testing class GenericLatticeConvexHull.

This file is part of the DGtal library.

Definition in file testGenericLatticeConvexHull.cpp.

Function Documentation

◆ g()

std::mt19937 g ( rd() )

Referenced by makeRandomLatticePointsFromDirVectors().

◆ makeRandomLatticePointsFromDirVectors()

template<typename Point >

std::vector< Point > makeRandomLatticePointsFromDirVectors	(	int	nb,
		const vector< Point > &	V
	)

Definition at line 65 of file testGenericLatticeConvexHull.cpp.

{
  std::uniform_int_distribution<int> U(-10, 10);
  vector< Point > P;
  int n = V[0].size();
  int m = V.size();
  Point A;
  for ( auto i = 0; i < n; i++ )
    A[ i ] = U( g );
  P.push_back( A );
  for ( auto k = 0; k < nb; k++ )
    {
      Point B = A;
      for ( auto i = 0; i < m; i++ )
        {
          int l = U( g );
          B += l * V[ i ];
        }
      P.push_back( B );
    }
  std::shuffle( P.begin(), P.end(), g );
  return P;
}

References g().

Referenced by SCENARIO().

◆ randomPointsInBall()

template<typename Point >

std::vector< Point > randomPointsInBall	(	int	nb,
		int	R
	)

Examples: geometry/tools/checkLatticeBallQuickHull.cpp.

Definition at line 49 of file testGenericLatticeConvexHull.cpp.

{
  std::vector< Point > V;
  Point c = Point::diagonal( R );
  double R2 = (double) R * (double) R;
  for ( int i = 0; i < nb; ) {
    Point p;
    for ( DGtal::Dimension k = 0; k < Point::dimension; ++k )
      p[ k ] = rand() % (2*R);
    if ( ( p - c ).squaredNorm() < R2 ) { V.push_back( p ); i++; }
  }
  return V;
}

References DGtal::R.

◆ SCENARIO() [1/2]

SCENARIO	(	"GenericLatticeConvexHull< ConvexHullIntegralKernel< 2 > > unit tests"	,
		""	[genquickhull][integral_kernel][2d]
	)

Definition at line 93 of file testGenericLatticeConvexHull.cpp.

{
  typedef GenericLatticeConvexHull< 2, int > CvxHull;
  typedef SpaceND< 2, int >                Space;      
  typedef Space::Point                     Point;
 
  
  GIVEN( "Given a set { } " ) {
    std::vector<Point> V = { };
    CvxHull hull;
    bool ok = hull.compute( V );
    std::cout << hull << std::endl;
    THEN( "Everything went fine and dim=-1." ) {
      REQUIRE( ok );
      REQUIRE( hull.affine_dimension == -1 );
      REQUIRE( hull.count() == 0 );
    }
  }
  GIVEN( "Given a set { (2,1), (2,1) } " ) {
    std::vector<Point> V
      = { Point(2,1), Point(2,1) };
    CvxHull hull;
    bool ok = hull.compute( V );
    std::cout << hull << std::endl;
    THEN( "Everything went fine and dim=0." ) {
      REQUIRE( ok );
      REQUIRE( hull.affine_dimension == 0 );
      REQUIRE( hull.count() == 1 );
    }
  }
  GIVEN( "Given a set { (0,0), (-4,-1), (8, 2), (16,4) } " ) {
    std::vector<Point> V
      = { Point(0,0), Point(-4,-1), Point(8,2), Point(16,4) };
    CvxHull hull;
    bool ok = hull.compute( V );
    std::cout << hull << std::endl;
    THEN( "Everything went fine and dim=1." ) {
      REQUIRE( ok );
      REQUIRE( hull.affine_dimension == 1 );
      REQUIRE( hull.count() == 6 );
    }
  }
  GIVEN( "Given a set { (0,0), (-4,-1), (-3,5), (7,3), (5, -2), (2, 2) } " ) {
    std::vector<Point> V
      = { Point(0,0), Point(-4,-1), Point(-3,5), Point(7,3), Point(5, -2), Point(2,2) };
    CvxHull hull;
    bool ok = hull.compute( V );
    std::cout << hull << std::endl;
    THEN( "Everything went fine and dim=2." ) {
      REQUIRE( ok );
      REQUIRE( hull.affine_dimension == 2 );
      REQUIRE( hull.count() > V.size() );
    }
  }
 
}

References GIVEN(), and REQUIRE().

◆ SCENARIO() [2/2]

SCENARIO	(	"GenericLatticeConvexHull< ConvexHullIntegralKernel< 3 > > unit tests"	,
		""	[genquickhull][integral_kernel][3d]
	)

Definition at line 150 of file testGenericLatticeConvexHull.cpp.

{
  typedef GenericLatticeConvexHull< 3, int > CvxHull;
  typedef SpaceND< 3, int >                Space;      
  typedef Space::Point                     Point;
 
  std::vector< Point > V1 = { Point{ 5,-2, 1 } };
  std::vector< Point > X1 = makeRandomLatticePointsFromDirVectors( 20, V1 );
  std::vector< Point > V2 = { Point{ 5,-2, 1 }, Point{ -3, 4, 2 } };
  std::vector< Point > X2 = makeRandomLatticePointsFromDirVectors( 20, V2 );
  std::vector< Point > V3 = { Point{ 5,-2, 1 }, Point{ -3, 4, 2 }, Point{ 1, -7, 11 } };
  std::vector< Point > X3 = makeRandomLatticePointsFromDirVectors( 20, V3 );
  
  GIVEN( "Given a 1-d lattice set in Z3 " ) {
    CvxHull hull;
    bool ok = hull.compute( X1 );
    std::cout << hull << std::endl;
    THEN( "Everything went fine and dim=1." ) {
      CAPTURE( hull.positions );
      REQUIRE( ok );
      REQUIRE( hull.affine_dimension == 1 );
      REQUIRE( hull.count() >= 2 );
    }
  }
  GIVEN( "Given a 2-d lattice set in Z3 " ) {
    CvxHull hull;
    bool ok = hull.compute( X2 );
    std::cout << hull << std::endl;
    THEN( "Everything went fine and dim=2." ) {
      REQUIRE( ok );
      REQUIRE( hull.affine_dimension == 2 );
      REQUIRE( hull.count() >= 3 );
    }
  }
  GIVEN( "Given a 3-d lattice set in Z3 " ) {
    CvxHull hull;
    bool ok = hull.compute( X3 );
    std::cout << hull << std::endl;
    THEN( "Everything went fine and dim=3." ) {
      REQUIRE( ok );
      REQUIRE( hull.affine_dimension == 3 );
      REQUIRE( hull.count() >= 4 );
    }
  }
}

References CAPTURE(), GIVEN(), makeRandomLatticePointsFromDirVectors(), and REQUIRE().

Variable Documentation

◆ rd

std::random_device rd

Definition at line 44 of file testGenericLatticeConvexHull.cpp.

Functions

Variables