testGeodesicsInHeat.cpp File Reference

#include <iostream>
#include "DGtal/base/Common.h"
#include "ConfigTest.h"
#include "DGtalCatch.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/dec/PolygonalCalculus.h"
#include "DGtal/dec/GeodesicsInHeat.h"
#include "DGtal/shapes/SurfaceMesh.h"
#include "DGtal/shapes/MeshHelpers.h"

Go to the source code of this file.

Functions
	TEST_CASE ("Testing GeodesicsInHeat")

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@l.nosp@m.iris..nosp@m.cnrs.nosp@m..fr ) Laboratoire d'InfoRmatique en Image et Systemes d'information - LIRIS (CNRS, UMR 5205), CNRS, France

Date

2021/09/21

Functions for testing class GeodesicsInHeat.

This file is part of the DGtal library.

Definition in file testGeodesicsInHeat.cpp.

Function Documentation

TEST_CASE()

TEST_CASE

(

"Testing GeodesicsInHeat"

)

Definition at line 47 of file testGeodesicsInHeat.cpp.

{
  typedef SurfaceMesh< RealPoint,RealPoint > Mesh;
  std::vector<RealPoint> positions = { RealPoint( 0, 0, 0 ) ,
    RealPoint( 1, 0, 0 ) ,
    RealPoint( 0, 1, 0 ) ,
    RealPoint( 1, 1, 0 ) ,
    RealPoint( 0, 0, 1 ) ,
    RealPoint( 1, 0, 1 ) ,
    RealPoint( 0, 1, 1 ) ,
    RealPoint( 1, 1, 1 ) ,
    RealPoint( 1, 0, 2 ) ,
    RealPoint( 0, 0, 2 ) };
  std::vector<Mesh::Vertices> faces = { { 1, 0, 2, 3 },
    { 0, 1, 5, 4 } ,
    { 1, 3, 7, 5 } ,
    { 3, 2, 6, 7 } ,
    { 2, 0, 4, 6 } ,
    { 4, 5, 8, 9 } };
  
  Mesh box(positions.cbegin(), positions.cend(),
           faces.cbegin(), faces.cend());
  
  PolygonalCalculus<RealPoint,RealVector> boxCalculus(box);
  
  
  SECTION("Construction and basic operators")
  {
    GeodesicsInHeat<PolygonalCalculus<RealPoint,RealVector>> heat(boxCalculus);
    REQUIRE( heat.isValid() == false );
    
    trace.beginBlock("init solvers");
    heat.init(0.1);
    trace.endBlock();
    REQUIRE( heat.isValid()  );
    
    heat.addSource(0);
    GeodesicsInHeat<PolygonalCalculus<RealPoint,RealVector>>::Vector d = heat.compute();
    REQUIRE( d.size() == positions.size() );
    REQUIRE( d[5] == Approx(1.444608) );
    heat.clearSource();
    auto sources = heat.source();
    REQUIRE(sources.sum() == 0);
  }
}

References DGtal::Trace::beginBlock(), DGtal::Trace::endBlock(), REQUIRE(), SECTION(), and DGtal::trace.