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: "Sebastien Loriot (GeometryFactory)" <>
- 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:43:01 +0200
- Authentication-results: mail3-smtp-sop.national.inria.fr; spf=None ; spf=Pass ; spf=None
- Ironport-phdr: 9a23:f47UVRMHYw3q9dXtFfYl6mtUPXoX/o7sNwtQ0KIMzox0I/nzrarrMEGX3/hxlliBBdydt6oazbKO+4nbGkU4qa6bt34DdJEeHzQksu4x2zIaPcieFEfgJ+TrZSFpVO5LVVti4m3peRMNQJW2aFLduGC94iAPERvjKwV1Ov71GonPhMiryuy+4ZLebxlJiTanfb9+MAi9oBnMuMURnYZsMLs6xAHTontPdeRWxGdoKkyWkh3h+Mq+/4Nt/jpJtf45+MFOTav1f6IjTbxFFzsmKHw65NfqtRbYUwSC4GYXX3gMnRpJBwjF6wz6Xov0vyDnuOdxxDWWMMvrRr0yRD+s7bpkSAXwhSkaKTA5/mHZhM9+gq1Vrx2upQBwzYHPbYGJN/dzZL/Rcc8USGdDWMtaSixPApm7b4sKF+cNM/tWoJfmqFQUthu+GQisC/3pyjRVgHH22rc60+UmEQHawQcsG9MAv2rTrNXwKqgSTfq1zanIwzrZc/xWxTD96JLIcxwkr/2HQLV9cdHVyUkoDQPFj1KQqZD/MzOQ0uQAqGeb7+96WuKuj24rsR1+oj+qxso1jITCm4wbylfB9SpjwYY1I8W1SEhhbtG/C5RQsz+VO5FqTcMlRmFkoCk6yroBuZ60eCgF1o4ryALYa/yCa4SI4xTjVPyQIThinn5ldqi/iw6u/kiv0OL8UtO40FFQripKlNTHq3MD1wTL58SZVvdw+l2t1DWP2gzJ9O1IP104mbDbJpMv2rI9lJkevV7dEiL5mUj6lrKae0Eq9+Sy6enqYrPrrYKGOYBukAHxKKEul9S/AesmNggOWHCW+eGm273i+U31WbtLjucrnqXAvpDXJcoWqrS2Aw9S1YYj5BK/ACm83NsEmnkHKUpJeBOBj4f3J1HDOO70Aeu7jli2kzpmx+rKMqP/DpnQNHTOkLnsca5460FGyQozyd5f54hTCrEEOP/zQUvxtNvCDh84NAy03/zrCNB51o4FVmKPB7WWMKLWsVOS+u0vJOyMaJcPuDnhM/gl++LujXghlFABcqmmx58XZGmlEfRnOEWWfWfsgswaEWoRpQo/TOnqiEWYXjJJZnayWbg85jAhB468A4fDXNPlvbqawS3uHoFKfnsUTReXAHLwfsOFXe0NYWScOIh6gzkcXP+gTYEmkhqhvQu/x7t8JffP4X4kssfo29FxouHSjhov7idcDsKH0mjLQXsnsHkPQmp85611qFdhy1qFmYx/mf1fCZQTyP5OVwomL4/yxudmDMrjG0iJKsyNT0ynRcngBDUZQdc4wttIaEF4TYbxxivf1janVudG34eAA4Y5p/qFjirBYv1lwnOD75EPylwvQ89BL2qj3/ct+A3aBoqPmEKcxf/zKfYsmRXV/WLG9lKg+VlCWVcpA6rAVHEbIEDRqIahvx6Qf/qVEb0idzB554uCJ69NMIC7iFxHQLL8Io2bbT/q3Wi3AhmMy/WHa4+4I2g=
No on linux. The precompiled demo on windows is available here:
https://www.cgal.org/demo/4.12/polyhedron_3.zip
with common dlls (to be put in the same directy as Polyhedron_3.exe
in the aforementioned archive):
https://www.cgal.org/demo/4.12/CGAL-demoDLLs.zip
Sebastien.
On 07/09/2018 11:40 AM, sergio campo wrote:
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) < <mailto:>> 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
<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
<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
<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
<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/
<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
<https://sympa.inria.fr/sympa/info/cgal-discuss>
- [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, sergio, 07/04/2018
- 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, Sebastien Loriot (GeometryFactory), 07/09/2018
- 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, sergio campo, 07/09/2018
- 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, Sebastien Loriot (GeometryFactory), 07/09/2018
- 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, Zesen Qian, 07/10/2018
- 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, Sebastien Loriot (GeometryFactory), 07/10/2018
- 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, sergio campo, 07/10/2018
- 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, Zesen Qian, 07/15/2018
- 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, sergio campo, 07/16/2018
- 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, sergio campo, 07/10/2018
- 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, Sebastien Loriot (GeometryFactory), 07/10/2018
- 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, Zesen Qian, 07/10/2018
- 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, Sebastien Loriot (GeometryFactory), 07/09/2018
- 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, sergio campo, 07/09/2018
- 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, Sebastien Loriot (GeometryFactory), 07/09/2018
Archive powered by MHonArc 2.6.18.