Abstract
In the present work, we extend the theoretical and numerical discussion of the well-known Laplace–Beltrami operator by equipping the underlying manifolds with additional structure provided by vector bundles. Focusing on the particular class of flat complex line bundles, we examine a whole family of Laplacians including the Laplace–Beltrami operator as a special case.
To demonstrate that our proposed approach is numerically feasible, we describe a robust and efficient finite-element discretization, supplementing the theoretical discussion with first numerical spectral decompositions of those Laplacians.
Our method is based on the concept of introducing complex phase discontinuities into the finite element basis functions across a set of homology generators of the given manifold. More precisely, given an m-dimensional manifold M and a set of n generators that span the relative homology group H m−1(M,∂M), we have the freedom to choose n phase shifts, one for each generator, resulting in a n-dimensional family of Laplacians with associated spectra and eigenfunctions. The spectra and absolute magnitudes of the eigenfunctions are not influenced by the exact location of the paths, depending only on their corresponding homology classes.
Employing our discretization technique, we provide and discuss several interesting computational examples highlighting special properties of the resulting spectral decompositions. We examine the spectrum, the eigenfunctions and their zero sets which depend continuously on the choice of phase shifts.
Similar content being viewed by others
References
Alexa, M., Wardetzky, M.: Discrete Laplacians on general polygonal meshes. ACM Trans. Graph. 30(4), 102 (2011)
Bronstein, M.M., Kokkinos, I.: Scale-invariant heat kernel signatures for non-rigid shape recognition. In: Computer Vision and Pattern Recognition, pp. 1704–1711. IEEE Press, New York (2010)
Dey, T.K., Li, K., Sun, J., Cohen-Steiner, D.: Computing geometry-aware handle and tunnel loops in 3D models. ACM Trans. Graph. 27(3), 1–9 (2008)
Diaz-Gutierrez, P., Eppstein, D., Gopi, M.: Curvature aware fundamental cycles. Comput. Graph. Forum 28(7), 2015–2024 (2009)
Dodziuk, J.: Finite-difference approach to the Hodge theory of harmonic forms. Am. J. Math. 98(1), 79–104 (1976)
Dong, S., Bremer, P., Garland, M., Pascucci, V., Hart, J.: Spectral surface quadrangulation. ACM Trans. Graph. 25(3), 1057–1066 (2006)
Dziuk, G.: Finite elements for the Beltrami operator on arbitrary surfaces. In: Hildebrandt, S., Leis, R. (eds.) Partial Differential Equations and Calculus of Variations. Lecture Notes in Mathematics, vol. 1357, pp. 142–155. Springer, Berlin (1988)
Erickson, J., Whittlesey, K.: Greedy optimal homotopy and homology generators. In: Symposium on Discrete Algorithms, pp. 1038–1046. SIAM, Philadelphia (2005)
Frankel, T.: The Geometry of Physics: An Introduction. Cambridge University Press, Cambridge (2011)
Gross, P.W., Kotiuga, P.R.: Electromagnetic Theory and Computation: A Topological Approach. Cambridge University Press, Cambridge (2004)
Hatcher, A.: Algebraic Topology. Cambridge University Press, Cambridge (2002)
Hernandez, V., Roman, J.E., Vidal, V.: SLEPc: a scalable and flexible toolkit for the solution of eigenvalue problems. ACM Trans. Math. Softw. 31(3), 351–362 (2005)
Hildebrandt, K., Schulz, C., von Tycowicz, C., Polthier, K.: Eigenmodes of surface energies for shape analysis. In: Mourrain, B., Schaefer, S., Xu, G. (eds.) Advances in Geometric Modeling and Processing. Lecture Notes in Computer Science, vol. 6130, pp. 296–314. Springer, Berlin (2010)
Hou, T., Qin, H.: Robust dense registration of partial nonrigid shapes. IEEE Trans. Vis. Comput. Graph., preprint (2011). http://doi.ieeecomputersociety.org/10.1109/TVCG.2011.267
Kotiuga, P.: An algorithm to make cuts for magnetic scalar potentials in tetrahedral meshes based on the finite element method. IEEE Trans. Magn. 25(5), 4129–4131 (1989)
Lévy, B.: Laplace–Beltrami eigenfunctions: towards an algorithm that understands geometry. In: International Conference on Shape Modeling and Applications (2006)
do Carmo, M.P.: Riemannian Geometry. Birkhäuser, Boston (1992)
Mémoli, F.: Spectral Gromov–Wasserstein distances for shape matching. In: International Conference on Computer Vision, pp. 256–263 (2009)
Niethammer, M., Reuter, M., Wolter, F.-E., Bouix, S., Peinecke, N., Ko, M.S., Shenton, M.E.: Global medical shape analysis using the Laplace–Beltrami-spectrum. In: Medical Image Computing and Computer Assisted Intervention (2007)
Ovsjanikov, M., Sun, J., Guibas, L.: Global intrinsic symmetries of shapes. Comput. Graph. Forum 27(5), 1341–1348 (2008)
Peinecke, N., Wolter, F.-E.: Mass density Laplace-spectra for image recognition. In: Cyberworlds, pp. 409–416. IEEE Press, New York (2007)
Peinecke, N., Wolter, F.-E., Reuter, M.: Laplace-spectra as fingerprints for image recognition. Comput. Aided Des. 39(6), 460–476 (2007)
Reuter, M.: Hierarchical shape segmentation and registration via topological features of Laplace–Beltrami eigenfunctions. J. Comput. Vis. 89, 287–308 (2010)
Reuter, M., Biasotti, S., Giorgi, D., Patanè, G., Spagnuolo, M.: Discrete Laplace–Beltrami operators for shape analysis and segmentation. Comput. Graph. 33(3), 381–390 (2009)
Reuter, M., Niethammer, M., Wolter, F.-E., Bouix, S., Shenton, M.E.: Global medical shape analysis using the volumetric Laplace spectrum. In: Cyberworlds, pp. 417–426. IEEE Press, New York (2007)
Reuter, M., Wolter, F.-E., Peinecke, N.: Laplace-spectra as fingerprints for shape matching. In: Symposium on Solid and Physical Modeling, pp. 101–106 (2005)
Reuter, M., Wolter, F.-E., Peinecke, N.: Laplace–Beltrami spectra as shape DNA of surfaces and solids. Comput. Aided Des. 38(4), 342–366 (2006)
Reuter, M., Wolter, F.-E., Shenton, M.E., Niethammer, M.: Laplace–Beltrami eigenvalues and topological features of eigenfunctions for statistical shape analysis. Comput. Aided Des. 41(10), 739–755 (2009)
Ruggeri, M., Patanè, G., Spagnuolo, M., Saupe, D.: Spectral-driven isometry-invariant matching of 3D shapes. J. Comput. Vis. 89, 248–265 (2010)
Rustamov, R.M.: Laplace–Beltrami eigenfunctions for deformation invariant shape representation. In: Belyaev, A., Garland, M. (eds.) Symposium on Geometry Processing, pp. 225–233 (2007)
Solin, P.: Partial Differential Equations and the Finite Element Method. Wiley, New York (2005)
Sun, J., Ovsjanikov, M., Guibas, L.: A concise and provably informative multi-scale signature based on heat diffusion. Comput. Graph. Forum 28(5), 1383–1392 (2009)
Taubin, G.: A signal processing approach to fair surface design. In: Computer Graphics and Interactive Techniques, pp. 351–358. ACM, New York (1995)
Vais, A., Berger, B., Wolter, F.-E.: Spectral computations on nontrivial line bundles. Comput. Graph. 36(5), 398–409 (2012)
Vallet, B., Lévy, B.: Spectral geometry processing with manifold harmonics. Comput. Graph. Forum 27(2), 251–260 (2008)
de Verdière, E.C., Lazarus, F.: Optimal system of loops on an orientable surface. Discrete Comput. Geom. 33(3), 507–534 (2005)
Wardetzky, M., Mathur, S., Kälberer, F., Grinspun, E.: Discrete Laplace operators: no free lunch. In: Symposium on Geometry Processing, pp. 33–37 (2007)
Wolter, F.-E., Friese, K.-I.: Local and global geometric methods for analysis interrogation, reconstruction, modification and design of shape. In: Computer Graphics International, pp. 137–151 (2000). Also available as Welfenlab report No. 3. ISSN 1866-7996
Wolter, F.-E., Howind, T., Altschaffel, T., Reuter, M., Peinecke, N.: Laplace-Spektra—Anwendungen in Gestalt- und Bildkognition. Available as Welfenlab Report No. 7. ISSN 1866-7996
Wolter, F.-E., Peinecke, N., Reuter, M.: Verfahren zur Charakterisierung Von Objekten/A Method for the Characterization of Objects (Surfaces, Solids and Images). German Patent Application, June 2005 (pending), US Patent US2009/0169050 A1, July 2, 2009, 2006
Xin, S., He, Y., Fu, C., Wang, D., Lin, S., Chu, W., Cheng, J., Gu, X., Lui, L.: Euclidean geodesic loops on high-genus surfaces applied to the morphometry of vestibular systems. In: Medical Image Computing and Computer-Assisted Intervention—MICCAI 2011, pp. 384–392 (2011)
Xu, G.: Discrete Laplace–Beltrami operators and their convergence. Comput. Aided Geom. Des. 21(8), 767–784 (2004)
Zhang, H., Van Kaick, O., Dyer, R.: Spectral mesh processing. Comput. Graph. Forum 29(6), 1865–1894 (2010)
Acknowledgements
The authors would like to thank H. Thielhelm for valuable comments and suggestions. We also thank the AIM@SHAPE Repository for making accessible the three-dimensional models used in this paper.
Author information
Authors and Affiliations
Corresponding author
Appendix: Algebraic topology
Appendix: Algebraic topology
Let M be a path-connected topological space. A continuous map γ:[0,1]→M with γ(0)=γ(1)=p is called a loop with base point p. Two loops are called homotopic if one can be gradually deformed into the other. The set of equivalence classes of loops with a fixed base point is called the first fundamental group of M, denoted by π 1(M,p).
A covering of M is a space \(\widetilde{M}\) with a surjective local homeomorphism \(\rho:\widetilde{M} \to M\). If \(\widetilde{M}\) is simply connected, then ρ is called a universal covering of M. A deck transformation is a homeomorphism \(h: \widetilde{M} \to \widetilde{M}\) such that ρ∘h=ρ. Using these notions, another geometric interpretation of the first fundamental group is given by the fact that it is isomorphic to the deck transformation group of the universal cover of M.
To visualize this construction, consider a topological torus M=S 1×S 1 as shown in Fig. 12. Its fundamental group is generated by two loops, denoted by a and b. Cutting the torus along these loops creates a rectangle, a so-called fundamental domain. The universal covering is the Euclidean plane \(\widetilde{M}=\mathbb{R}^{2}\), which is tiled with copies of the fundamental domain. The fundamental group acts on the points \(\widetilde{q} \in\widetilde{M}\) by translation, sending \(\widetilde{q}\) to the corresponding point in another copy of the fundamental domain.
Now, assume M is a manifold represented by a singular simplicial complex and let R be an arbitrary ring. A k-chain is a formal linear combination of oriented k-simplices with coefficients in the ring R. The set of k-chains form a group C k under addition. The boundary operator ∂ k :C k →C k−1 is a linear operator that maps any oriented simplex to the chain consisting of its appropriately signed oriented boundary simplices. A chain α∈C k is called closed, or cycle, if ∂α=0 and it is called exact if it can be written as α=∂γ for some γ∈C k+1. Any exact chain is closed as a consequence of the fact that the boundary of a boundary is empty, i.e. ∂ 2=0. Therefore, the sequence of chain groups C k with the boundary operators in between form a chain complex.
Now let Z k be the group of closed k-chains and let B k be the group of exact k-cycles. Two cycles α,β∈Z k are called homologous, if α−β∈B k . The kth homology groups are the quotient groups H k (M,R):=Z k /B k induced by this equivalence relation.
The set of homomorphisms from H k to R form the k-th cohomology group H k(M,R).
For manifolds M with non-empty boundary ∂M, we will also need the so called relative homology groups H k (M,∂M) which are obtained by modifying the homology equivalence relation to treat any chain on the boundary as zero. In this relative homology, two relative k-cycles a,β are called homologous if α−β=∂γ+δ for a (k+1)-chain γ and some k-chain δ contained in the boundary ∂M.
The well-known Lefschetz duality theorem guarantees the existence of an isomorphism between the cohomology group H k(M) and the relative homology group H m−k (M,∂M) where m denotes the dimension of the manifold.
For the annulus above, the cycle α is a generator of H 1(M). It is also a generator of the first fundamental group π 1(M). The blue path β is not a cycle, since its boundary consists of two points. However, it is a relative cycle since these points lie on the boundary ∂M. In fact, β is a generator of H 1(M,∂M).
Rights and permissions
About this article
Cite this article
Vais, A., Berger, B. & Wolter, FE. Complex line bundle Laplacians. Vis Comput 29, 345–357 (2013). https://doi.org/10.1007/s00371-012-0737-5
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00371-012-0737-5