Skip to main content
Log in

Complex line bundle Laplacians

  • Original Article
  • Published:
The Visual Computer Aims and scope Submit manuscript

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.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10
Fig. 11

Similar content being viewed by others

References

  1. Alexa, M., Wardetzky, M.: Discrete Laplacians on general polygonal meshes. ACM Trans. Graph. 30(4), 102 (2011)

    Article  Google Scholar 

  2. 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)

    Google Scholar 

  3. 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)

    Article  Google Scholar 

  4. Diaz-Gutierrez, P., Eppstein, D., Gopi, M.: Curvature aware fundamental cycles. Comput. Graph. Forum 28(7), 2015–2024 (2009)

    Article  Google Scholar 

  5. Dodziuk, J.: Finite-difference approach to the Hodge theory of harmonic forms. Am. J. Math. 98(1), 79–104 (1976)

    Article  MathSciNet  MATH  Google Scholar 

  6. Dong, S., Bremer, P., Garland, M., Pascucci, V., Hart, J.: Spectral surface quadrangulation. ACM Trans. Graph. 25(3), 1057–1066 (2006)

    Article  Google Scholar 

  7. 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)

    Chapter  Google Scholar 

  8. Erickson, J., Whittlesey, K.: Greedy optimal homotopy and homology generators. In: Symposium on Discrete Algorithms, pp. 1038–1046. SIAM, Philadelphia (2005)

    Google Scholar 

  9. Frankel, T.: The Geometry of Physics: An Introduction. Cambridge University Press, Cambridge (2011)

    Book  MATH  Google Scholar 

  10. Gross, P.W., Kotiuga, P.R.: Electromagnetic Theory and Computation: A Topological Approach. Cambridge University Press, Cambridge (2004)

    Book  MATH  Google Scholar 

  11. Hatcher, A.: Algebraic Topology. Cambridge University Press, Cambridge (2002)

    MATH  Google Scholar 

  12. 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)

    Article  MathSciNet  MATH  Google Scholar 

  13. 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)

    Chapter  Google Scholar 

  14. 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

  15. 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)

    Article  Google Scholar 

  16. Lévy, B.: Laplace–Beltrami eigenfunctions: towards an algorithm that understands geometry. In: International Conference on Shape Modeling and Applications (2006)

    Google Scholar 

  17. do Carmo, M.P.: Riemannian Geometry. Birkhäuser, Boston (1992)

    Book  MATH  Google Scholar 

  18. Mémoli, F.: Spectral Gromov–Wasserstein distances for shape matching. In: International Conference on Computer Vision, pp. 256–263 (2009)

    Google Scholar 

  19. 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)

    Google Scholar 

  20. Ovsjanikov, M., Sun, J., Guibas, L.: Global intrinsic symmetries of shapes. Comput. Graph. Forum 27(5), 1341–1348 (2008)

    Article  Google Scholar 

  21. Peinecke, N., Wolter, F.-E.: Mass density Laplace-spectra for image recognition. In: Cyberworlds, pp. 409–416. IEEE Press, New York (2007)

    Google Scholar 

  22. Peinecke, N., Wolter, F.-E., Reuter, M.: Laplace-spectra as fingerprints for image recognition. Comput. Aided Des. 39(6), 460–476 (2007)

    Article  Google Scholar 

  23. Reuter, M.: Hierarchical shape segmentation and registration via topological features of Laplace–Beltrami eigenfunctions. J. Comput. Vis. 89, 287–308 (2010)

    Article  Google Scholar 

  24. 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)

    Article  Google Scholar 

  25. 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)

    Google Scholar 

  26. 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)

    Google Scholar 

  27. 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)

    Article  Google Scholar 

  28. 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)

    Article  Google Scholar 

  29. Ruggeri, M., Patanè, G., Spagnuolo, M., Saupe, D.: Spectral-driven isometry-invariant matching of 3D shapes. J. Comput. Vis. 89, 248–265 (2010)

    Article  Google Scholar 

  30. 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)

    Google Scholar 

  31. Solin, P.: Partial Differential Equations and the Finite Element Method. Wiley, New York (2005)

    Book  Google Scholar 

  32. 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)

    Article  Google Scholar 

  33. Taubin, G.: A signal processing approach to fair surface design. In: Computer Graphics and Interactive Techniques, pp. 351–358. ACM, New York (1995)

    Google Scholar 

  34. Vais, A., Berger, B., Wolter, F.-E.: Spectral computations on nontrivial line bundles. Comput. Graph. 36(5), 398–409 (2012)

    Article  Google Scholar 

  35. Vallet, B., Lévy, B.: Spectral geometry processing with manifold harmonics. Comput. Graph. Forum 27(2), 251–260 (2008)

    Article  Google Scholar 

  36. de Verdière, E.C., Lazarus, F.: Optimal system of loops on an orientable surface. Discrete Comput. Geom. 33(3), 507–534 (2005)

    Article  MathSciNet  MATH  Google Scholar 

  37. Wardetzky, M., Mathur, S., Kälberer, F., Grinspun, E.: Discrete Laplace operators: no free lunch. In: Symposium on Geometry Processing, pp. 33–37 (2007)

    Google Scholar 

  38. 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

    Google Scholar 

  39. 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

  40. 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

  41. 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)

    Chapter  Google Scholar 

  42. Xu, G.: Discrete Laplace–Beltrami operators and their convergence. Comput. Aided Geom. Des. 21(8), 767–784 (2004)

    Article  MATH  Google Scholar 

  43. Zhang, H., Van Kaick, O., Dyer, R.: Spectral mesh processing. Comput. Graph. Forum 29(6), 1865–1894 (2010)

    Article  Google Scholar 

Download references

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

Authors

Corresponding author

Correspondence to Alexander Vais.

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.

Fig. 12
figure 12

Fundamental group of the torus acting on its universal covering space by deck transformations

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 mk (M,∂M) where m denotes the dimension of the manifold.

figure a

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

Reprints 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

Download citation

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00371-012-0737-5

Keywords

Navigation