CGAL users discussion list

["Sebastien Loriot (GeometryFactory)" <>]

No idea what you did, if I open the polyhedron demo, load your model
call orient polygon soup and remove self intersections, everything is 
done in less than a minute on my computer and the output model does not
contain any self-intersection.

Do you have intermediate writing to "string" with a rounding of the 
coordinates of the point not using the max precision 
(std::setprecision(17) for the ostream)?

Sebastien.

On 07/04/2018 11:33 AM, sergio wrote:
> This is may application that i compile in windows in Release with the showed
> CMakeList.txt and cgal_vtk_test.cpp files. When I try to repair this file:
>
> Number intersection: 290
> https://drive.google.com/file/d/1uqcbZqgtbiiXYjKnzzb5d-OXm4JAHqdx/view?usp=sharing
>
> 1. my application build in Release generate the following file:
>
> Number intersections: 179
> Add some visible holes and some deformations.
> Repairing time:..................201.128 seconds.
>
> https://drive.google.com/file/d/13Vfw0pP5Ymw0oo4RnfRjKTI_aTndb6P4/view?usp=sharing
>
>
> 2. my application compiled in Debug, it generate the following file:
>
> Number intersections: 4
> Repairing time:..................1072.96 seconds.
>
> https://drive.google.com/file/d/1x1Z7qaNcnfGC8a24qzJI-DZEF3qqeLqq/view?usp=sharing
>
>
> with warnings:
>
> Expression : !"Returning no output. Filling hole with extra triangles is not
> successful!"
> File       :
> C:\dev\CGAL-4.12\include\CGAL/Polygon_mesh_processing/internal/Hole_filling/Triangulate_hole_polyline.h
> Line       : 1103
> Explanation:
> Refer to the bug-reporting instructions at
> https://www.cgal.org/bug_report.html
>
> The problem is that in debug build it takes around 17 min for reparacion and
> in Release it  only takes 3 minutes, but in relese build it modify the
> original mesh and add a big hole(pleas check links to files)
>
> What's wrong?
>
>
>
>
>
> -------------------------------------------------------------------------------
> CMakeList.txt
> -------------------------------------------------------------------------------
>
> cmake_minimum_required(VERSION 2.8.4)
>
> project( cgal_vtk_test )
>
> # Find CGAL
> find_package(CGAL REQUIRED COMPONENTS Core) # If the dependency is required,
> use REQUIRED option - if it's not found CMake will issue an error
> include( ${CGAL_USE_FILE} )
>
> # Find VTK
> find_package(VTK REQUIRED)
> include(${VTK_USE_FILE})
>
>
> # Setup your executable
> #include_directories (BEFORE "include")
>
> include( CGAL_CreateSingleSourceCGALProgram )
> create_single_source_cgal_program( "cgal_vtk_test.cpp" ) # This will create
> an executable target with name 'cgal_vtk_test'
>
> # Add VTK link libraries to your executable target
> target_link_libraries(cgal_vtk_test ${VTK_LIBRARIES})
> #target_link_libraries(cgal_vtk_test)
>
> -------------------------------------------------------------------------------
> cgal_vtk_test.cpp
> -------------------------------------------------------------------------------
>
> #include <fstream>
> #include <iomanip>
> #include <ctime>
>
> //
> ----------------------------------------------------------------------------
> // CGAL
> // ----------------------------------------------------------------------
>
>
>
> #include <limits>
>
> #include <CGAL/Polygon_mesh_processing/repair.h>
> #include <CGAL/Polygon_mesh_processing/internal/repair_extra.h>
> #include <CGAL/Polygon_mesh_processing/corefinement.h>
>
>
> #include <CGAL/Polygon_mesh_processing/polygon_soup_to_polygon_mesh.h>
> #include <CGAL/Polyhedron_3.h>
> #include <CGAL/Surface_mesh/Surface_mesh.h>
> #include <CGAL/Polygon_mesh_processing/corefinement.h>
> #include <CGAL/IO/print_wavefront.h>
> #include <CGAL/IO/OBJ_reader.h>
> #include <CGAL/IO/print_wavefront.h>
> #include <CGAL/assertions.h>
>
> #include <CGAL/IO/STL_reader.h>
> #include <CGAL/IO/STL_writer.h>
>
>
> typedef CGAL::Exact_predicates_inexact_constructions_kernel Kernel;
> typedef CGAL::Surface_mesh<Kernel::Point_3>             Mesh;
> typedef CGAL::Polyhedron_3<Kernel> Polyhedron;
>
> void addStringToFileName(std::string  & fileName, const std::string & str)
> {
>         size_t lastindex = fileName.find_last_of(".");
>         std::string name = fileName.substr(0, lastindex);
>
>         fileName = name + str;
> }
>
> bool repairOBJFile(std::string & fileName)
> {
>         clock_t begin = clock();
>
>         std::cout << std::endl;
>         std::cout << std::endl;
>
>         std::cout << 
> "..........................................................."
> << std::endl;
>         std::cout << ".                        Repairing                       
> ."
> << std::endl;
>         std::cout << 
> "..........................................................."
> << std::endl;
>
>         std::ifstream input(fileName);
>         if (!input)
>         {
>                 std::cerr << "Cannot open file " << std::endl;
>                 return false;
>         }
>
>         std::vector<Kernel::Point_3> points;
>         std::vector< std::vector<std::size_t> > polygons;
>         if (!CGAL::read_OBJ(input, points, polygons))
>         {
>                 std::cerr << "Error parsing the OFF file " << std::endl;
>                 return false;
>         }
>         //Mesh mesh;
>         Polyhedron mesh;
>         CGAL::Polygon_mesh_processing::orient_polygon_soup(points, polygons);
>         CGAL::Polygon_mesh_processing::polygon_soup_to_polygon_mesh(points,
> polygons, mesh);
>
>         bool solved =
> CGAL::Polygon_mesh_processing::remove_self_intersections(mesh);
>
>         if (!solved)
>         {
>                 std::cout << "Some self-intersection could not be fixed" << 
> std::endl;
>                 //return false;
>         }
>
>         addStringToFileName(fileName, "_fixed.obj");
>         std::ofstream ofs(fileName);
>         CGAL::print_polyhedron_wavefront(ofs, mesh);
>
>         clock_t end = clock();
>         double elapsed_secs = double(end - begin) / CLOCKS_PER_SEC;
>
>         std::cout << std::endl;
>         std::cout << std::endl;
>
>         std::cout << "Repairing time:.................." << elapsed_secs << "
> seconds."<< std::endl;
>
>         return true;
> }
>   
> int main(int argc, char* argv[])
> {
>         std::string fileName = argv[1];
>         repairOBJFile(fileName);
>         return 0;
> }
>
>
>
>
> --
> Sent from: http://cgal-discuss.949826.n4.nabble.com/