GCC Code Coverage Report

Directory: ./
File: examples/geogram/intersect/main.cpp
Date: 2026-09-07 02:28:19
Exec Total Coverage
Lines: 82 87 94.3%
Functions: 1 1 100.0%
Branches: 137 346 39.6%
Line Branch Exec Source
1 /*
2 * Copyright (c) 2000-2022 Inria
3 * All rights reserved.
4 *
5 * Redistribution and use in source and binary forms, with or without
6 * modification, are permitted provided that the following conditions are met:
7 *
8 * * Redistributions of source code must retain the above copyright notice,
9 * this list of conditions and the following disclaimer.
10 * * Redistributions in binary form must reproduce the above copyright notice,
11 * this list of conditions and the following disclaimer in the documentation
12 * and/or other materials provided with the distribution.
13 * * Neither the name of the ALICE Project-Team nor the names of its
14 * contributors may be used to endorse or promote products derived from this
15 * software without specific prior written permission.
16 *
17 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
18 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
19 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
20 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
21 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
22 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
23 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
24 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
25 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
26 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
27 * POSSIBILITY OF SUCH DAMAGE.
28 *
29 * Contact: Bruno Levy
30 *
31 * https://www.inria.fr/fr/bruno-levy
32 *
33 * Inria,
34 * Domaine de Voluceau,
35 * 78150 Le Chesnay - Rocquencourt
36 * FRANCE
37 *
38 */
39
40 #include <geogram/basic/common.h>
41 #include <geogram/basic/logger.h>
42 #include <geogram/basic/command_line.h>
43 #include <geogram/basic/command_line_args.h>
44 #include <geogram/basic/stopwatch.h>
45 #include <geogram/basic/file_system.h>
46 #include <geogram/mesh/mesh.h>
47 #include <geogram/mesh/mesh_io.h>
48 #include <geogram/mesh/mesh_repair.h>
49 #include <geogram/mesh/mesh_surface_intersection.h>
50
51 /*
52 * Computes auto-intersections in a mesh.
53 */
54
55
56 13 int main(int argc, char** argv) {
57 using namespace GEO;
58
59 // Needs to be called once.
60 13 GEO::initialize(GEO::GEOGRAM_INSTALL_ALL);
61
2/4
✓ Branch 2 taken 13 times.
✗ Branch 3 not taken.
✓ Branch 5 taken 13 times.
✗ Branch 6 not taken.
 26 Stopwatch W_total("Total time");
62
63
64 try {
65
66 std::vector<std::string> filenames;
67
68
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 CmdLine::import_arg_group("standard");
69
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 CmdLine::import_arg_group("algo");
70
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::declare_arg(
71
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
 26 "post",true,"glue triangle and unglue non-manifold edges"
72 );
73
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::declare_arg(
74
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
 26 "verbose",true,"makes intersection algorithm more chatty"
75 );
76
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::declare_arg(
77
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
 26 "fine_verbose",true,"makes intersection algorithm even more chatty"
78 );
79
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::declare_arg(
80
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
 26 "Delaunay",true,"use Delaunay triangulation to remesh intersections"
81 );
82
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::declare_arg(
83
2/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
 26 "detect_intersecting_neighbors",true,
84
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 "test also neighboring triangles for intersection"
85 );
86
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::declare_arg(
87
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
 26 "remove_internal_shells",true,"remove internal shells"
88 );
89
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::declare_arg(
90
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
 26 "simplify_coplanar_facets",true,"simplify coplanar facets"
91 );
92
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::declare_arg(
93
2/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
 26 "coplanar_angle_tolerance",0.0,
94
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 "maximum angle (in degrees) between coplanar facets"
95 );
96
4/10
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✓ Branch 7 taken 13 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✓ Branch 10 taken 13 times.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
 26 CmdLine::declare_arg("expr","","Region classification expression");
97
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::declare_arg(
98
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
 26 "monster_threshold",100000,"monster threshold"
99 );
100
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::declare_arg(
101
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
 26 "dry_run",false,"Do not insert triangulations in global mesh"
102 );
103
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::declare_arg(
104
2/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
 26 "save_skeleton",false,
105
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 "Save skeleton of intersection in skeleton.geogram"
106 );
107
108 13 if(
109
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 !CmdLine::parse(
110
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 13 times.
✗ Branch 4 not taken.
 26 argc, argv, filenames, "inputfile <outputfile|none>"
111 )
112 ) {
113 return 1;
114 }
115
116
117 std::string output_filename =
118
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 13 times.
✓ Branch 3 taken 13 times.
✗ Branch 4 not taken.
 13 filenames.size() >= 2 ? filenames[1] : std::string("out.meshb");
119
120
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 Logger::div("Data I/O");
121
122
2/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
 13 Logger::out("I/O") << "Output = " << output_filename << std::endl;
123
124
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 Mesh A;
125
126
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 MeshIOFlags load_flags;
127 13 if(
128
3/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 13 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 13 times.
✗ Branch 6 not taken.
 26 FileSystem::extension(output_filename) == "mesh" ||
129
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 13 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 13 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
 26 FileSystem::extension(output_filename) == "meshb"
130 ) {
131 load_flags.set_attributes(MESH_ALL_ATTRIBUTES);
132 }
133
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
✓ Branch 4 taken 13 times.
 13 if(!mesh_load(filenames[0],A, load_flags)) {
134 return 1;
135 }
136
137 {
138
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 Logger::div("Intersect");
139
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 MeshSurfaceIntersection I(A);
140
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 I.set_verbose(CmdLine::get_arg_bool("verbose"));
141
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 I.set_fine_verbose(CmdLine::get_arg_bool("fine_verbose"));
142
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 I.set_delaunay(CmdLine::get_arg_bool("Delaunay"));
143 I.set_detect_intersecting_neighbors(
144
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 CmdLine::get_arg_bool("detect_intersecting_neighbors")
145 );
146 I.set_radial_sort(
147
2/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
 13 CmdLine::get_arg_bool("remove_internal_shells") ||
148
2/12
✗ Branch 0 not taken.
✓ Branch 1 taken 13 times.
✗ Branch 3 not taken.
✗ Branch 4 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
✗ Branch 10 not taken.
✓ Branch 11 taken 13 times.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
 26 CmdLine::get_arg_bool("simplify_coplanar_facets") ||
149
2/12
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✗ Branch 4 not taken.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 13 times.
✗ Branch 10 not taken.
✗ Branch 11 not taken.
✗ Branch 12 not taken.
✗ Branch 13 not taken.
 13 CmdLine::get_arg("expr") != ""
150 );
151
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 I.set_monster_threshold(
152
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 CmdLine::get_arg_uint("monster_threshold")
153 );
154 I.set_dry_run(
155
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 CmdLine::get_arg_bool("dry_run")
156 );
157
158
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 Mesh skel;
159
3/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
 26 if(CmdLine::get_arg_bool("save_skeleton")) {
160 I.set_build_skeleton(&skel);
161 }
162
163
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 I.intersect();
164
165
3/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
 26 if(CmdLine::get_arg("expr") != "") {
166 ✗ I.classify(CmdLine::get_arg("expr"));
167
3/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 13 times.
✗ Branch 7 not taken.
 26 } else if(CmdLine::get_arg_bool("remove_internal_shells")) {
168
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 I.remove_internal_shells();
169 }
170
171
3/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 13 times.
✗ Branch 7 not taken.
 26 if(CmdLine::get_arg_bool("simplify_coplanar_facets")) {
172
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 I.simplify_coplanar_facets(
173
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 26 CmdLine::get_arg_double("coplanar_angle_tolerance")
174 );
175 }
176
177
3/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✓ Branch 7 taken 13 times.
 26 if(CmdLine::get_arg_bool("save_skeleton")) {
178 ✗ mesh_save(skel,"skeleton.geogram");
179 }
180
181 13 }
182
183
3/6
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 13 times.
✗ Branch 7 not taken.
 26 if(CmdLine::get_arg_bool("post")) {
184
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 Logger::div("Post");
185
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 Stopwatch W_post("Post");
186
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 mesh_repair(
187 A,
188 MeshRepairMode(
189 MESH_REPAIR_COLOCATE | MESH_REPAIR_DUP_F
190 ),
191 0.0
192 );
193 13 }
194
195
196
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
✓ Branch 6 taken 13 times.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
✗ Branch 9 not taken.
 26 Logger::div("Data I/O");
197
198
1/2
✓ Branch 0 taken 13 times.
✗ Branch 1 not taken.
 13 if(output_filename != "none") {
199
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 MeshIOFlags save_flags;
200 13 if(
201
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 13 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 13 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
 26 FileSystem::extension(output_filename) == "mesh" ||
202
3/8
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 13 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 13 times.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
✗ Branch 8 not taken.
 26 FileSystem::extension(output_filename) == "meshb"
203 ) {
204 save_flags.set_attributes(MESH_ALL_ATTRIBUTES);
205 }
206
1/2
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
 13 mesh_save(A, output_filename, save_flags);
207 }
208
209
210 26 }
211 ✗ catch(const std::exception& e) {
212 ✗ std::cerr << "Received an exception: " << e.what() << std::endl;
213 return 1;
214 ✗ }
215
2/4
✓ Branch 1 taken 13 times.
✗ Branch 2 not taken.
✓ Branch 4 taken 13 times.
✗ Branch 5 not taken.
 13 Logger::div("Total time");
216
217 13 return 0;
218 13 }
219