Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] remove_self_intersections changes the mesh and the reparation is less when it is build in Release than when i build my application in Debug

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] remove_self_intersections changes the mesh and the reparation is less when it is build in Release than when i build my application in Debug


Chronological Thread 
  • From: sergio campo <>
  • To:
  • Subject: Re: [cgal-discuss] remove_self_intersections changes the mesh and the reparation is less when it is build in Release than when i build my application in Debug
  • Date: Mon, 9 Jul 2018 11:40:19 +0200
  • Authentication-results: mail2-smtp-roc.national.inria.fr; spf=None ; spf=Pass ; spf=None
  • Ironport-phdr: 9a23:bfJpPBVynTuvlF25ZI6gngIWRTXV8LGtZVwlr6E/grcLSJyIuqrYbBGFt8tkgFKBZ4jH8fUM07OQ7/i+HzRYqb+681k6OKRWUBEEjchE1ycBO+WiTXPBEfjxciYhF95DXlI2t1uyMExSBdqsLwaK+i764jEdAAjwOhRoLerpBIHSk9631+ev8JHPfglEnjWwba9zIRmssQndqtQdjJd/JKo21hbHuGZDdf5MxWNvK1KTnhL86dm18ZV+7SleuO8v+tBZX6nicKs2UbJXDDI9M2Ao/8LrrgXMTRGO5nQHTGoblAdDDhXf4xH7WpfxtTb6tvZ41SKHM8D6Uaw4VDK/5KptVRTmijoINyQh/W7VhMx+jKxVrhG8qRJh34HZe5uaOOZkc67HYd8WWWhMU8BMXCJBGIO8aI4PAvIaPelGtYn9qFoOrRyjDgawHuzvzydIhnjr1qA9zeQuCx3J3A0mH9IBvnTUq8/1NKAJXOC6yanH1zTDb/dM1Tjh74jIdwksrPeRVrx+dsrRzFMgFwLDjliIr4zqJSia2f8Js2eF9eZvSeWvi2s/pwF+pTivwcgsio3XiY4P1l/E8iB5zJ40Jd2+Uk57YMSrHIFetyGAL4d5XswiTHtsuCogzb0Go5G7cS4Xw5ok3x7Sc+KLf5SM7x75V+ucIS10iGx5dL+8nRq/8Eitx+vhXceuyllKtDBKktzUu3ANyRPT7s+HR+N4/ki72DaP0xnf5f9ZLkwpjKbbJZ8sz74qmpoctkTDGSD2mEHog6OMakok/e2o5/zmYrXguJCcK5d5hh/iPqkqgMCyAuQ1PhIQU2SH+umwzrzu8E3hTLVPlPI2k63ZsJ7AJcQco660GxFV3Zw/5Ba+CTem39oYnX4dLFJCYx+HgI3pNEvPIPD8F/uwn1OskDJzy/DcIrLhGonNLmTEkLr5Ybl95FRTyA4qwd9C5pJUEa0OIO/oWk/qr9HYFR84Mwmsw+n9Etl914UeWXiOAqCDKq/Sv0WItaoSJfKRbtoVpCrlMKpioOX/iGcw314bZ6igm5UNL2upG+xvZESfb32ridgIFSIGvxE1Uffx22GFBDVcbnL3U6Mn7SwgE6qnC53CT8ajmu+vxiC+S4ZbfG9HDRiHGHLwcM3QRf4HaSSdLMh7lRQLULGgT8kq0hT451yy8KZuMueBon5QjpnkztUgv7SCxyF3ziR9CoGm60/ISmh1mm0SQDpvhfJwpEV8zhGI1q0q2qUER+wW3OtAV0IBDbCZ1/ZzUomgVQfIf9PPQ1GjEI3/XGMBC+kpytpLWH5TXtWviheZgnivCr4R0rGHXdk6rP+a0H/2KMJwjX3B0ft5gg==

Thanks for your fast answer! :). I set std::setprecision(17)   with the same results. Have you execute the demo with Windows? In my case I executed the program in Windows 10.

On Mon, Jul 9, 2018 at 10:27 AM, Sebastien Loriot (GeometryFactory) <> wrote:
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/


--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://sympa.inria.fr/sympa/info/cgal-discuss






Archive powered by MHonArc 2.6.18.

Top of Page