Geogram  Version 1.9.1
A programming library of geometric algorithms
optimal_transport_2d.h
Go to the documentation of this file.
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 #ifndef H_EXPLORAGRAM_OPTIMAL_TRANSPORT_OPTIMAL_TRANSPORT_2D_H
41 #define H_EXPLORAGRAM_OPTIMAL_TRANSPORT_OPTIMAL_TRANSPORT_2D_H
42 
46 
53 namespace GEO {
54 
79  Mesh* omega,
80  index_t nb_points,
81  const double* points,
82  double* centroids,
83  Mesh* RVD=nullptr,
84  bool verbose=false,
85  index_t nb_air_particles = 0,
86  const double* air_particles = nullptr,
87  index_t air_particles_stride = 0,
88  double air_fraction = 0.0,
89  const double* weights_in = nullptr,
90  double* weights_out = nullptr,
91  index_t nb_iter = 1000
92  );
93 
94 
111  public:
120  Mesh* mesh,
121  const std::string& delaunay = "BPOW2d",
122  bool BRIO=false
123  );
124 
129 
133  void get_RVD(Mesh& M) override;
134 
138  void compute_Laguerre_centroids(double* centroids) override;
139 
145  double total_mesh_mass() const;
146 
147  protected:
151  void call_callback_on_RVD() override;
152 
153  public:
162  };
163 
164  /*********************************************************************/
165 
166 }
167 
168 #endif
Internal representation of polygons for GenericVoronoiDiagram.
Represents a mesh.
Definition: mesh.h:2701
Computes semi-discrete optimal transport maps.
void get_RVD(Mesh &M) override
Computes a mesh with the restricted Voronoi diagram.
GEOGen::Polygon work_
Used by clipping operations.
void compute_Laguerre_centroids(double *centroids) override
Computes the centroids of the Laguerre cells.
OptimalTransportMap2d(Mesh *mesh, const std::string &delaunay="BPOW2d", bool BRIO=false)
OptimalTransportMap2d constructor.
GEOGen::Polygon clipped_
Used by clipping operations.
void call_callback_on_RVD() override
Calls the callback for each intersection between a Laguerre cell and a simplex of the background mesh...
~OptimalTransportMap2d() override
OptimalTransportMap destructor.
double total_mesh_mass() const
Gets the total mass of the mesh.
Computes semi-discrete optimal transport maps.
#define EXPLORAGRAM_API
Linkage declaration for exploragram symbols.
Definition: defs.h:18
Included by all headers in exploragram.
Internal representation of polygons for GenericVoronoiDiagram.
Global Vorpaline namespace.
Definition: basic.h:55
void EXPLORAGRAM_API compute_Laguerre_centroids_2d(Mesh *omega, index_t nb_points, const double *points, double *centroids, Mesh *RVD=nullptr, bool verbose=false, index_t nb_air_particles=0, const double *air_particles=nullptr, index_t air_particles_stride=0, double air_fraction=0.0, const double *weights_in=nullptr, double *weights_out=nullptr, index_t nb_iter=1000)
Computes the centroids of the Laguerre cells that correspond to optimal transport in 2D.
geo_index_t index_t
The type for storing and manipulating indices.
Definition: numeric.h:329
Base class for semi-discrete optimal transport.