A Simple Discretization of the Vector Dirichlet Energy

Total Page:16

File Type:pdf, Size:1020Kb

A Simple Discretization of the Vector Dirichlet Energy DOI:10.1111/cgf.14070 Eurographics Symposium on Geometry Processing 2020 Volume 39 (2020), Number 5 Q. Huang and A. Jacobson (Guest Editors) A Simple Discretization of the Vector Dirichlet Energy Oded Stein1, Max Wardetzky2, Alec Jacobson3 and Eitan Grinspun3;1 1Columbia University, USA 2University of Göttingen, Germany 3University of Toronto, Canada Abstract We present a simple and concise discretization of the covariant derivative vector Dirichlet energy for triangle meshes in 3D using Crouzeix-Raviart finite elements. The discretization is based on linear discontinuous Galerkin elements, and is simple to implement, without compromising on quality: there are two degrees of freedom for each mesh edge, and the sparse Dirichlet energy matrix can be constructed in a single pass over all triangles using a short formula that only depends on the edge lengths, reminiscent of the scalar cotangent Laplacian. Our vector Dirichlet energy discretization can be used in a variety of applications, such as the calculation of Killing fields, parallel transport of vectors, and smooth vector field design. Experiments suggest convergence and suitability for applications similar to other discretizations of the vector Dirichlet energy. 1. Introduction The covariant derivative r generalizes the gradient of scalar func- tions to vector fields defined on surfaces. As the gradient does for scalar functions, the covariant derivative measures the infinitesimal change of a vector field in every direction. As with the gradient’s 1 R 2 scalar Dirichlet energy, Escalar(u) := 2 W kruk dx for a smooth scalar function u and a smooth surface W, the covariant derivative has a corresponding vector Dirichlet energy, Z 1 2 E(u) := krukF dx , (1) transporting the 2 W red vector (enlarged) denoising a vector field by smoothing across the surface where u is a smooth vector-valued function on W, and k·kF is the Frobenius norm. Much like the scalar Dirichlet energy Escalar does Figure 1: Parallel transporting a vector across a surface with our for scalar functions, the vector Dirichlet energy E measures the Crouzeix-Raviart discretization of vector heat flow (left), and de- smoothness of a vector field. noising a vector field with our discretization of the vector diffusion equation (right). Just as Escalar is useful for scalar data processing, E has many uses in vector field processing. While E has been employed scalar Our discretization is based on Crouzeix-Raviart finite elements. in many geometry processing methods and applications, E has, in Its most notable features are its simplicity and ease of implementa- comparison, seen less usage in practice. A key reason for the wide tion, without sacrificing any quality: adoption of Escalar is the existence of a simple, reliable, and read- ily available finite element discretization: the cotangent Laplacian. • the real, sparse Dirichlet energy matrix can be built in a single There are a variety of existing discretizations of E, but they can pass over all faces using a simple formula depending only on not be implemented using quite as few lines of code as the cotan- each (individual) triangle’s edge lengths; gent Laplacian. The matrix for the cotangent Laplacian can be con- • the degrees of freedom directly correspond to the edges of the 1 structed using only the expression Lcotani j = 2 cotai j + cotbi j mesh, making the results easy to understand and visualize; for each vertex pair i; j (where ai j;bi j are the two angles opposite • no preprocessing or intermediate data structures beyond simple the edge i j), as well as summation [MDSB03, (5)]. We advocate for matrices are required; a discretization of E that aspires to the simplicity of the cotangent • applications, performance and convergence are similar to other Laplacian. discretizations of E. c 2020 The Author(s) Computer Graphics Forum c 2020 The Eurographics Association and John Wiley & Sons Ltd. Published by John Wiley & Sons Ltd. https://www.eg.org https://diglib.eg.org 82 O. Stein, M. Wardetzky, A. Jacobson & E. Grinspun / A Simple Discretization of the Vector Dirichlet Energy The sparse 2m×2m-matrix L of the discrete energy E for a mesh The covariant derivative on a smooth surface is this vector gradi- 3 with m edges ent of R restricted to the tangent space of the surface, where it 1: for all faces f , consecutive edges ei;e j 2 f do becomes dependent on the surface’s curvature (while this defini- 3 2: li = len(ei); l j = len(e j); lk = len(ek) // ek 3rd edge tion only works for surfaces embedded inR , it generalizes to other 3: s = orientation(ei) · orientation(e j) // can be 1 or −1 kinds of surfaces [Lee97, Section 5]). 2 2 2 2 4: a = s(li + l j − lk ) =(4lil j area( f )) Operators like ru, which consume a vector field and return a 2 2 2 5: b = s(li + l j − lk )=(lil j) vector field, have an associated inner product : that generalizes the 2 2 scalar product · . The operation : for X;Y represented as matrices 6: L(ei;ei) += li =area( f ); L(ei +m;ei +m) += li =area( f ) can be written as X : Y := tr(X|Y), where X| is the transpose of // L(a;b) x expands to L(a;b) = x; L(b;a) = x an operator X that takes the surface’s metric into account [Lee97]. 7: L(ei;e j) a; L(ei + m;e j + m) a X|Y is sometimes also referred to as composing the adjoint of X 8: L(ei;e j + m) b; L(e j;ei + m) −b with Y. ∗ We demonstrate our discretization for three established applica- The covariant derivative r has an adjoint operator r , similar tions of the vector Dirichlet energy—smoothing vector fields, cre- to how the gradient of a scalar function is adjoint to the divergence ating Killing vector fields, and efficient parallel transport of vectors of a vector field. The adjoint covariant derivative is defined by in- (see Figure 1). We study the convergence of our method and two tegration by parts over a surface W, previous methods, and we find that our method converges at the Z Z Z ru : X dx = u · r∗X dx + u · X(n) dG same rate as previous methods. W W ¶W (2) The way we arrive at our discretization parallels the discretiza- n unit boundary normal at ¶W , tion of Escalar by the cotangent Laplacian for applications in scalar where X is an operator that takes a vector field and returns a vector geometry processing. Desbrun et al. [DMSB99] adopted the well- field, and u is a vector field. This adjoint operator is used to define established cotangent Laplacian [Mac49], demonstrated its utility, a Laplace operator for vector fields, the Bochner Laplacian (also and spurred its adoption in geometry processing. Similarly, we referred to as the connection Laplacian) [Pet06, pp. 209], employ a discretization introduced by Stein et al. [SJWG20] as part D u = r∗ru . of a routine to process scalar functions, and demonstrate its utility B n for vector field processing. Our discretization generalizes popular In R , the Bochner Laplacian is (up to sign) the usual vector field and extensively studied methods from numerical analysis of finite Laplacian defined by applying the scalar Laplacian in each coor- 2 elements for flat domains in R , which simplifies mathematical dinate. It appears, for example, in the Navier-Stokes equations for study of the method. fluid mechanics [Fan19]. 2.2. The vector Dirichlet energy 2. Theoretical background The Bochner Laplacian DB defines a vector Laplace equation, In this section we provide a brief overview of the covariant deriva- DBu = 0 , (3) tive and the vector Dirichlet energy. Readers who are already fa- which shares many properties with the scalar Laplace equation miliar with these topics, or who want to get straight to the appli- Du = 0, such as smoothness of the associated flow [BGV96, Chap- cations, can safely skip this section. The discretization and matrix ter 2]. In general, the right-hand side can be any function, but we implementation follow in Section 3. will assume that it is zero for simplicity. The Bochner Laplacian has an associated energy, constructed us- 2.1. The covariant derivative ing integration by parts. Let u;v be smooth vector fields on the sur- face W. Then we can write (3), by multiplying with a test function For a vector field u, the covariant derivative ru generalizes the v and integrating, as gradient of a scalar function, r f [Lee97, Section 4]. Z Z Z We will be working on a smooth, orientable surface W. The DBu · v dx = ru : rv dx − rnu · v dG W W ¶W (4) covariant derivative is an operator r : T ( ) × T ( ) !T ( ), W W W n unit boundary normal at ¶W , where T (W) is the space of tangent vector fields on the smooth surface W. It is usually written as rvu for two vector fields v;u, where we used (2). By calculus of variation [Eva92, Section 8], the and when the subscript v is omitted, ru is interpreted as the op- solution u of (4) is the minimizer of the vector Dirichlet energy, erator that takes v as an argument and returns rvu. In the flat Eu- 3 Z clidean space R , with a vector field u written as a column vector 1 2 | E(u) := krukF dx , (5) u = ux uy uz , the covariant derivative is simply componen- 2 W twise differentiation, 0 1 2 ¶xux ¶yux ¶zux where kXkF := X : X is the Frobenius norm. The boundary term in ru = @¶xuy ¶yuy ¶zuyA . (4) implies that minimizers will fulfill the Neumann boundary con- ¶xuz ¶yuz ¶zuz dition, rnu = 0 at ¶W, if no explicit boundary condition is applied c 2020 The Author(s) Computer Graphics Forum c 2020 The Eurographics Association and John Wiley & Sons Ltd.
Recommended publications
  • HARMONIC FUNCTIONS for BOUNDED and UNBOUNDED P(X)
    UNIVERSITY OF JYVASKYL¨ A¨ UNIVERSITAT¨ JYVASKYL¨ A¨ DEPARTMENT OF MATHEMATICS INSTITUT FUR¨ MATHEMATIK AND STATISTICS UND STATISTIK REPORT 130 BERICHT 130 EXISTENCE AND UNIQUENESS OF p(x)-HARMONIC FUNCTIONS FOR BOUNDED AND UNBOUNDED p(x) JUKKA KEISALA JYVASKYL¨ A¨ 2011 UNIVERSITY OF JYVASKYL¨ A¨ UNIVERSITAT¨ JYVASKYL¨ A¨ DEPARTMENT OF MATHEMATICS INSTITUT FUR¨ MATHEMATIK AND STATISTICS UND STATISTIK REPORT 130 BERICHT 130 EXISTENCE AND UNIQUENESS OF p(x)-HARMONIC FUNCTIONS FOR BOUNDED AND UNBOUNDED p(x) JUKKA KEISALA JYVASKYL¨ A¨ 2011 Editor: Pekka Koskela Department of Mathematics and Statistics P.O. Box 35 (MaD) FI{40014 University of Jyv¨askyl¨a Finland ISBN 978-951-39-4299-1 ISSN 1457-8905 Copyright c 2011, Jukka Keisala and University of Jyv¨askyl¨a University Printing House Jyv¨askyl¨a2011 Contents 0 Introduction 3 0.1 Notation and prerequisities . 6 1 Constant p, 1 < p < 1 8 1.1 The direct method of calculus of variations . 10 1.2 Dirichlet energy integral . 10 2 Infinity harmonic functions, p ≡ 1 13 2.1 Existence of solutions . 14 2.2 Uniqueness of solutions . 18 2.3 Minimizing property and related topics . 22 3 Variable p(x), with 1 < inf p(x) < sup p(x) < +1 24 4 Variable p(x) with p(·) ≡ +1 in a subdomain 30 4.1 Approximate solutions uk ...................... 32 4.2 Passing to the limit . 39 5 One-dimensional case, where p is continuous and sup p(x) = +1 45 5.1 Discussion . 45 5.2 Preliminary results . 47 5.3 Measure of fp(x) = +1g is positive . 48 6 Appendix 51 References 55 0 Introduction In this licentiate thesis we study the Dirichlet boundary value problem ( −∆ u(x) = 0; if x 2 Ω; p(x) (0.1) u(x) = f(x); if x 2 @Ω: Here Ω ⊂ Rn is a bounded domain, p :Ω ! (1; 1] a measurable function, f : @Ω ! R the boundary data, and −∆p(x)u(x) is the p(x)-Laplace operator, which is written as p(x)−2 −∆p(x)u(x) = − div jru(x)j ru(x) for finite p(x).
    [Show full text]
  • Berestycki, Introduction to the Gaussian Free Field and Liouville
    Introduction to the Gaussian Free Field and Liouville Quantum Gravity Nathanaël Berestycki [Draft lecture notes: July 19, 2016] 1 Contents 1 Definition and properties of the GFF6 1.1 Discrete case∗ ...................................6 1.2 Green function..................................8 1.3 GFF as a stochastic process........................... 10 1.4 GFF as a random distribution: Dirichlet energy∗ ................ 13 1.5 Markov property................................. 15 1.6 Conformal invariance............................... 16 1.7 Circle averages.................................. 16 1.8 Thick points.................................... 18 1.9 Exercises...................................... 19 2 Liouville measure 22 2.1 Preliminaries................................... 23 2.2 Convergence and uniform integrability in the L2 phase............ 24 2.3 The GFF viewed from a Liouville typical point................. 25 2.4 General case∗ ................................... 27 2.5 The phase transition in Gaussian multiplicative chaos............. 29 2.6 Conformal covariance............................... 29 2.7 Random surfaces................................. 32 2.8 Exercises...................................... 32 3 The KPZ relation 34 3.1 Scaling exponents; KPZ theorem........................ 34 3.2 Applications of KPZ to exponents∗ ....................... 36 3.3 Proof in the case of expected Minkowski dimension.............. 37 3.4 Sketch of proof of Hausdorff KPZ using circle averages............ 39 3.5 Proof of multifractal spectrum of Liouville
    [Show full text]
  • Properly-Weighted Graph Laplacian for Semi-Supervised Learning
    PROPERLY-WEIGHTED GRAPH LAPLACIAN FOR SEMI-SUPERVISED LEARNING JEFF CALDER Department of Mathematics, University of Minnesota DEJAN SLEPČEV Department of Mathematical Sciences, Carnegie Mellon University Abstract. The performance of traditional graph Laplacian methods for semi-supervised learning degrades substantially as the ratio of labeled to unlabeled data decreases, due to a degeneracy in the graph Laplacian. Several approaches have been proposed recently to address this, however we show that some of them remain ill-posed in the large-data limit. In this paper, we show a way to correctly set the weights in Laplacian regularization so that the estimator remains well posed and stable in the large-sample limit. We prove that our semi-supervised learning algorithm converges, in the infinite sample size limit, to the smooth solution of a continuum variational problem that attains the labeled values continuously. Our method is fast and easy to implement. 1. Introduction For many applications of machine learning, such as medical image classification and speech recognition, labeling data requires human input and is expensive [13], while unlabeled data is relatively cheap. Semi-supervised learning aims to exploit this dichotomy by utilizing the geometric or topological properties of the unlabeled data, in conjunction with the labeled data, to obtain better learning algorithms. A significant portion of the semi-supervised literature is on transductive learning, whereby a function is learned only at the unlabeled points, and not as a parameterized function on an ambient space. In the transductive setting, graph based algorithms, such the graph Laplacian-based learning pioneered by [53], are widely used and have achieved great success [3,26,27,44–48,50, 52].
    [Show full text]
  • A Simple Discretization of the Vector Dirichlet Energy
    Eurographics Symposium on Geometry Processing 2020 Volume 39 (2020), Number 5 Q. Huang and A. Jacobson (Guest Editors) A Simple Discretization of the Vector Dirichlet Energy Oded Stein1, Max Wardetzky2, Alec Jacobson3 and Eitan Grinspun3;1 1Columbia University, USA 2University of Göttingen, Germany 3University of Toronto, Canada Abstract We present a simple and concise discretization of the covariant derivative vector Dirichlet energy for triangle meshes in 3D using Crouzeix-Raviart finite elements. The discretization is based on linear discontinuous Galerkin elements, and is simple to implement, without compromising on quality: there are two degrees of freedom for each mesh edge, and the sparse Dirichlet energy matrix can be constructed in a single pass over all triangles using a short formula that only depends on the edge lengths, reminiscent of the scalar cotangent Laplacian. Our vector Dirichlet energy discretization can be used in a variety of applications, such as the calculation of Killing fields, parallel transport of vectors, and smooth vector field design. Experiments suggest convergence and suitability for applications similar to other discretizations of the vector Dirichlet energy. 1. Introduction The covariant derivative r generalizes the gradient of scalar func- tions to vector fields defined on surfaces. As the gradient does for scalar functions, the covariant derivative measures the infinitesimal change of a vector field in every direction. As with the gradient’s 1 R 2 scalar Dirichlet energy, Escalar(u) := 2 W kruk dx for a smooth scalar function u and a smooth surface W, the covariant derivative has a corresponding vector Dirichlet energy, Z 1 2 E(u) := krukF dx , (1) transporting the 2 W red vector (enlarged) denoising a vector field by smoothing across the surface where u is a smooth vector-valued function on W, and k·kF is the Frobenius norm.
    [Show full text]
  • The Laplacian in Its Different Guises
    U.U.D.M. Project Report 2019:8 The Laplacian in its different guises Oskar Tegby Examensarbete i matematik, 15 hp Handledare: Martin Strömqvist Examinator: Jörgen Östensson Februari 2019 Department of Mathematics Uppsala University Contents 1 p = 2: The Laplacian 3 1.1 In the Dirichlet problem . 3 1.2 In vector calculus . 6 1.3 As a minimizer . 8 1.4 In the complex plane . 13 1.5 As the mean value property . 16 1.6 Viscosity solutions . 17 1.7 In stochastic processes . 18 2 2 < p < and p = : The p-Laplacian and -Laplacian 20 2.1 Ded∞uction of th∞e p-Laplacian and -Laplaci∞an . 20 2.2 As a minimizer . ∞. 22 2.3 Uniqueness of its solutions . 24 2.4 Existence of its solutions . 25 2.4.1 Method 1: Weak solutions and Sobolev spaces . 25 2.4.2 Method 2: Viscosity solutions and Perron’s method . 27 2.5 In the complex plane . 28 2.5.1 The nonlinear Cauchy-Riemann equations . 28 2.5.2 Quasiconformal maps . 29 2.6 In the asymptotic expansion . 36 2.7 As a Tug of War game with or without drift . 41 2.8 Uniqueness of the solutions . 44 1 Introduction The Laplacian n ∂2u ∆u := u = ∇ · ∇ ∂x2 i=1 i � is a scalar operator which gives the divergence of the gradient of a scalar field. It is prevalent in the famous Dirichlet problem, whose importance cannot be overstated. It entails finding the solution u to the problem ∆u = 0 in Ω �u = g on ∂Ω. The Dirichlet problem is of fundamental importance in mathematics and its applications, and the efforts to solve the problem has led to many revolutionary ideas and important advances in mathematics.
    [Show full text]
  • Arxiv:1703.06851V1
    SYMMETRIES AND CONSERVATION LAWS OF A NONLINEAR SIGMA MODEL WITH GRAVITINO JÜRGEN JOST, ENNO KEßLER, JÜRGEN TOLKSDORF, RUIJUN WU, MIAOMIAO ZHU Abstract. We show that the action functional of the nonlinear sigma model with gravitino considered in [18] is invariant under rescaled conformal transformations, super Weyl transfor- mations and diffeomorphisms. We give a careful geometric explanation how a variation of the metric leads to the corresponding variation of the spinors. In particular cases and despite using only commutative variables, the functional possesses a degenerate super symmetry. The corre- sponding conservation laws lead to a geometric interpretation of the energy-momentum tensor and supercurrent as holomorphic sections of appropriate bundles. 1. Introduction The main motivation for the introduction of the two-dimensional nonlinear supersymmetric sigma model in quantum field theory, or more specifically super gravity and super string theory, are its symmetries, see for instance [3, 9]. Furthermore, as argued in [20], the functional is determined by its symmetries together with suitable bounds on the order of its Euler–Lagrange equations. While super symmetric models are usually formulated using anti-commutative vari- ables, in [18] an analogue of the two-dimensional nonlinear supersymmetric sigma model using only commutative variables was introduced. Here we would like to give a detailed geometric account of the symmetries of the purely commutative model considered in [18]. In physics as well as in geometry, symmetries of a system always have important implications. In particular, differentiable symmetries of a system give rise to conservation laws—this is ba- sically what the famous Noether’s theorem states, see e.g.
    [Show full text]
  • Regularity in the Calculus of Variations
    Regularity in the Calculus of Variations Connor Mooney 1 Contents 1 Introduction 3 2 Scalar Equations 4 2.1 Gradient Bounds Inherited From the Boundary . 7 2.2 De Giorgi-Nash-Moser . 8 2.3 Density Estimates and H¨olderRegularity . 10 2.4 Harnack Inequalities . 11 2.5 Boundary Estimates . 13 2.6 Dirichlet Problem for Minimal Surfaces . 16 2.7 Interior Gradient Estimate for Minimal Surfaces . 18 2.8 Optimal Gradient Bound . 22 2.9 Concluding Remarks on Minimal Surfaces . 24 2.10 Stability in the Calculus of Variations . 25 3 Elliptic Systems 27 3.1 Basics . 27 3.1.1 Existence . 27 3.1.2 Quasiconvexity . 28 3.1.3 Lower Semicontinuity . 29 3.1.4 Rank One Convexity . 31 3.2 Polyconvexity, Quasiconvexity, and Rank One Convexity . 32 3.2.1 Quadratic Functions . 32 3.3 Sverak's Example . 34 3.4 Linear Theory . 35 3.5 Partial Regularity . 37 3.6 Caccioppoli inequality for quasiconvex functionals . 39 3.7 Counterexamples . 40 3.8 Reverse H¨olderInequality . 42 3.9 Special Structure: Radial . 43 3.10 Quasiconformal Mappings . 44 2 1 Introduction In these notes we outline the regularity theory for minimizers in the calculus of variations. The first part of the notes deals with the scalar case, with emphasis on the minimal surface equation. The second deals with vector mappings, which have different regularity properties due to the loss of the maximum principle. These notes are based on parts of courses given by Prof. Ovidiu Savin in Fall 2011 and Spring 2014. 3 2 Scalar Equations Scalar divergence-form equations generally arise by taking minimizers of interesting energy functionals.
    [Show full text]
  • Calculus of Variations Lecture Notes Riccardo Cristoferi May 9 2016
    Calculus of Variations Lecture Notes Riccardo Cristoferi May 9 2016 Contents Introduction 7 Chapter 1. Some classical problems1 1.1. The brachistochrone problem1 1.2. The hanging cable problem3 1.3. Minimal surfaces of revolution4 1.4. The isoperimetric problem5 N Chapter 2. Minimization on R 9 2.1. The one dimensional case9 2.2. The multidimensional case 11 2.3. Local representation close to non degenerate points 12 2.4. Convex functions 14 2.5. On the notion of differentiability 15 Chapter 3. First order necessary conditions for one dimensional scalar functions 17 3.1. The first variation - C2 theory 17 3.2. The first variation - C1 theory 24 3.3. Lagrangians of some special form 29 3.4. Solution of the brachistochrone problem 32 3.5. Problems with free ending points 37 3.6. Isoperimetric problems 38 3.7. Solution of the hanging cable problem 42 3.8. Broken extremals 44 Chapter 4. First order necessary conditions for general functions 47 4.1. The Euler-Lagrange equation 47 4.2. Natural boundary conditions 49 4.3. Inner variations 51 4.4. Isoperimetric problems 53 4.5. Holomic constraints 54 Chapter 5. Second order necessary conditions 57 5.1. Non-negativity of the second variation 57 5.2. The Legendre-Hadamard necessary condition 58 5.3. The Weierstrass necessary condition 59 Chapter 6. Null lagrangians 63 Chapter 7. Lagrange multipliers and eigenvalues of the Laplace operator 67 Chapter 8. Sufficient conditions 75 8.1. Coercivity of the second variation 75 6 CONTENTS 8.2. Jacobi conjugate points 78 8.3.
    [Show full text]
  • Gaussian Free Field, Liouville Quantum Gravity and Gaussian Multiplicative
    Gaussian free field, Liouville quantum gravity and Gaussian multiplicative chaos Nathanaël Berestycki Ellen Powell [Draft lecture notes: March 16, 2021] 1 Contents 1 Definition and properties of the GFF7 1.1 Discrete case...................................7 1.2 Continuous Green function............................ 11 1.3 GFF as a stochastic process........................... 14 1.4 Integration by parts and Dirichlet energy.................... 17 1.5 Reminders about functions spaces........................ 18 1.6 GFF as a random distribution.......................... 22 1.7 Itô’s isometry for the GFF............................ 25 1.8 Cameron–Martin space of the Dirichlet GFF.................. 26 1.9 Markov property................................. 28 1.10 Conformal invariance............................... 29 1.11 Circle averages.................................. 30 1.12 Thick points.................................... 31 1.13 Exercises...................................... 34 2 Liouville measure 36 2.1 Preliminaries................................... 37 2.2 Convergence and uniform integrability in the L2 phase............ 38 2.3 Weak convergence to Liouville measure..................... 40 2.4 The GFF viewed from a Liouville typical point................. 41 2.5 General case.................................... 44 2.6 The phase transition for the Liouville measure................. 46 2.7 Conformal covariance............................... 46 2.8 Random surfaces................................. 48 2.9 Exercises.....................................
    [Show full text]
  • A Discrete Laplace-Beltrami Operator for Simplicial Surfaces
    A discrete Laplace-Beltrami operator for simplicial surfaces Alexander I. Bobenko Boris A. Springborn February 24, 2006 Institut f¨ur Mathematik, Technische Universit¨at Berlin, Strasse des 17. Juni 136, 10623 Berlin, Germany [email protected] [email protected] Abstract We define a discrete Laplace-Beltrami operator for simplicial sur- faces (Definition 16). It depends only on the intrinsic geometry of the surface and its edge weights are positive. Our Laplace operator is sim- ilar to the well known finite-elements Laplacian (the so called “cotan formula”) except that it is based on the intrinsic Delaunay triangula- tion of the simplicial surface. This leads to new definitions of discrete harmonic functions, discrete mean curvature, and discrete minimal sur- faces. The definition of the discrete Laplace-Beltrami operator depends on the existence and uniqueness of Delaunay tessellations in piecewise flat surfaces. While the existence is known, we prove the uniqueness. Using Rippa’s Theorem we show that, as claimed, Musin’s harmonic index provides an optimality criterion for Delaunay triangulations, and this can be used to prove that the edge flipping algorithm terminates also in the setting of piecewise flat surfaces. Keywords: Laplace operator, Delaunay triangulation, Dirichlet energy, arXiv:math/0503219v3 [math.DG] 24 Feb 2006 simplicial surfaces, discrete differential geometry 1 Dirichlet energy of piecewise linear functions Let S be a simplicial surface in 3-dimensional Euclidean space, i.e. a geomet- 3 ric simplicial complex in Ê whose carrier S is a 2-dimensional submanifold, Research for this article was supported by the DFG Research Unit 565 “Polyhedral Surfaces” and the DFG Research Center Matheon “Mathematics for key technologies” in Berlin.
    [Show full text]
  • Discrete Differential Geometry: an Applied Introduction
    DISCRETE DIFFERENTIAL GEOMETRY: AN APPLIED INTRODUCTION Keenan Crane Last updated: February 25, 2021 Contents Chapter 1. Introduction 3 1.1. Disclaimer 6 1.2. Copyright 6 1.3. Acknowledgements 6 Chapter 2. Combinatorial Surfaces 7 2.1. Abstract Simplicial Complex 8 2.2. Anatomy of a Simplicial Complex: Star, Closure, and Link 10 2.3. Simplicial Surfaces 14 2.4. Adjacency Matrices 16 2.5. Halfedge Mesh 17 2.6. Written Exercises 21 2.7. Coding Exercises 26 Chapter 3. A Quick and Dirty Introduction to Differential Geometry 28 3.1. The Geometry of Surfaces 28 3.2. Derivatives and Tangent Vectors 31 3.3. The Geometry of Curves 34 3.4. Curvature of Surfaces 37 3.5. Geometry in Coordinates 41 Chapter 4. A Quick and Dirty Introduction to Exterior Calculus 45 4.1. Exterior Algebra 46 4.2. Examples of Wedge and Star in Rn 52 4.3. Vectors and 1-Forms 54 4.4. Differential Forms and the Wedge Product 58 4.5. Hodge Duality 62 4.6. Differential Operators 67 4.7. Integration and Stokes’ Theorem 73 4.8. Discrete Exterior Calculus 77 Chapter 5. Curvature of Discrete Surfaces 84 5.1. Vector Area 84 5.2. Area Gradient 87 5.3. Volume Gradient 89 5.4. Other Definitions 91 5.5. Gauss-Bonnet 94 5.6. Numerical Tests and Convergence 95 Chapter 6. The Laplacian 100 6.1. Basic Properties 100 6.2. Discretization via FEM 103 6.3. Discretization via DEC 107 1 CONTENTS 2 6.4. Meshes and Matrices 110 6.5. The Poisson Equation 112 6.6.
    [Show full text]
  • On the Existance of Minimizers of the Variable Exponent Dirichlet Energy Integral
    COMMUNICATIONS ON Website: http://AIMsciences.org PURE AND APPLIED ANALYSIS Volume 5, Number 3, September 2006 pp. 415{422 ON THE EXISTANCE OF MINIMIZERS OF THE VARIABLE EXPONENT DIRICHLET ENERGY INTEGRAL Peter A. HastÄ oÄ Department of Mathematics and Statistics P.O. Box 68, FI-00014 University of Helsinki, Finland (Communicated by Stefano Bianchini) Abstract. In this note we consider the Dirichlet energy integral in the variable exponent case under minimal assumptions on the exponent. First we show that the Dirichlet energy integral always has a minimizer if the boundary values are in L1. Second, we give an example which shows that if the so-called \jump- condition", known to be su±cient, is violated, then a minimizer need not exist for unbounded boundary values. 1. Introduction. After decades of sporadic research e®orts, the ¯eld of variable exponent function spaces entered into a phase of great activity starting in the late 1990's. The researchers were motivated both by inherent interest in developing techniques that work in non-translation invariant, but otherwise very natural, func- tion spaces, and by interesting applications (see, e.g., [3, 24]). For a survey of the present state of the ¯eld, the reader is referred to the article [8] which also includes a comprehensive bibliography of over a hundred titles on the subject. Two questions can be singled out as having received most attention: integral operators (Hardy{ Littlewood, Riesz, other singular integrals, see, e.g., [6, 7, 22]) and p(x)-Laplacian equations with Dirichlet boundary values. The latter is the subject matter of this paper.
    [Show full text]