You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
78 lines
2.7 KiB
78 lines
2.7 KiB
1 year ago
|
// This file is part of libigl, a simple c++ geometry processing library.
|
||
|
|
//
|
||
|
|
// Copyright (C) 2018 Alec Jacobson <alecjacobson@gmail.com>
|
||
|
|
//
|
||
|
|
// This Source Code Form is subject to the terms of the Mozilla Public License
|
||
|
|
// v. 2.0. If a copy of the MPL was not distributed with this file, You can
|
||
|
|
// obtain one at http://mozilla.org/MPL/2.0/.
|
||
|
|
#ifndef IGL_HEAT_GEODESICS_H
|
||
|
|
#define IGL_HEAT_GEODESICS_H
|
||
|
|
#include "igl_inline.h"
|
||
|
|
#include "min_quad_with_fixed.h"
|
||
|
|
#include <Eigen/Sparse>
|
||
|
|
#include <Eigen/Sparse>
|
||
|
|
namespace igl
|
||
|
|
{
|
||
|
|
/// Precomputation data for heat_geodesics_solve
|
||
|
|
///
|
||
|
|
/// \fileinfo
|
||
|
|
template <typename Scalar>
|
||
|
|
struct HeatGeodesicsData
|
||
|
|
{
|
||
|
|
/// Gradient operator
|
||
|
|
Eigen::SparseMatrix<Scalar> Grad;
|
||
|
|
/// Divergence operator
|
||
|
|
Eigen::SparseMatrix<Scalar> Div;
|
||
|
|
/// Number of gradient components
|
||
|
|
int ng;
|
||
|
|
/// List of boundary vertex indices
|
||
|
|
Eigen::VectorXi b;
|
||
|
|
/// Cached solvers for Dirichet, Neumann problems
|
||
|
|
min_quad_with_fixed_data<Scalar> Dirichlet,Neumann,Poisson;
|
||
|
|
/// Whether to use intrinsic Delaunay Laplacian
|
||
|
|
bool use_intrinsic_delaunay = false;
|
||
|
|
};
|
||
|
|
/// Precompute factorized solvers for computing a fast approximation of
|
||
|
|
/// geodesic distances on a mesh (V,F). [Crane et al. 2013]
|
||
|
|
///
|
||
|
|
/// @param[in] V #V by dim list of mesh vertex positions
|
||
|
|
/// @param[in] F #F by 3 list of mesh face indices into V
|
||
|
|
/// @param[out] data precomputation data (see heat_geodesics_solve)
|
||
|
|
///
|
||
|
|
/// \fileinfo
|
||
|
|
template < typename DerivedV, typename DerivedF, typename Scalar >
|
||
|
|
IGL_INLINE bool heat_geodesics_precompute(
|
||
|
|
const Eigen::MatrixBase<DerivedV> & V,
|
||
|
|
const Eigen::MatrixBase<DerivedF> & F,
|
||
|
|
HeatGeodesicsData<Scalar> & data);
|
||
|
|
/// \overload
|
||
|
|
/// @param[in] t "heat" parameter (smaller --> more accurate, less stable)
|
||
|
|
///
|
||
|
|
/// \fileinfo
|
||
|
|
template < typename DerivedV, typename DerivedF, typename Scalar >
|
||
|
|
IGL_INLINE bool heat_geodesics_precompute(
|
||
|
|
const Eigen::MatrixBase<DerivedV> & V,
|
||
|
|
const Eigen::MatrixBase<DerivedF> & F,
|
||
|
|
const Scalar t,
|
||
|
|
HeatGeodesicsData<Scalar> & data);
|
||
|
|
/// Compute fast approximate geodesic distances using precomputed data from a
|
||
|
|
/// set of selected source vertices (gamma).
|
||
|
|
///
|
||
|
|
/// @param[in] data precomputation data (see heat_geodesics_precompute)
|
||
|
|
/// @param[in] gamma #gamma list of indices into V of source vertices
|
||
|
|
/// @param[out] D #V list of distances to gamma
|
||
|
|
///
|
||
|
|
/// \fileinfo
|
||
|
|
template < typename Scalar, typename Derivedgamma, typename DerivedD>
|
||
|
|
IGL_INLINE void heat_geodesics_solve(
|
||
|
|
const HeatGeodesicsData<Scalar> & data,
|
||
|
|
const Eigen::MatrixBase<Derivedgamma> & gamma,
|
||
|
|
Eigen::PlainObjectBase<DerivedD> & D);
|
||
|
|
}
|
||
|
|
|
||
|
|
#ifndef IGL_STATIC_LIBRARY
|
||
|
|
#include "heat_geodesics.cpp"
|
||
|
|
#endif
|
||
|
|
|
||
|
|
#endif
|