mixed finite elements for variational surface modeling · 2019-03-06 · previous work simple...
TRANSCRIPT
![Page 1: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/1.jpg)
Mixed Finite Elements for Variational Surface Modeling Alec Jacobson Elif Tosun Olga Sorkine Denis Zorin New York University
![Page 2: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/2.jpg)
Produce high-quality surfaces Via energy minimization Or solving Partial Differential Equations
Motivation
Laplacian energy Laplacian gradient energy
Biharmonic equation Triharmonic equation
![Page 3: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/3.jpg)
Motivation Obtain different boundary conditions
Region
Curve
Point
![Page 4: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/4.jpg)
Motivation Higher-order equations on mesh (i.e. piecewise
linear elements) Dealing with higher-order derivatives not
straightforward
![Page 5: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/5.jpg)
Previous work Simple domains, analytic boundaries
[Bloor and Wilson 1990] Model shaped minimization of curvature variation energy
[Moreton and Séquin 1992] Interpolate curve networks, local quadratic fits and finite differences
[Welch and Witkin 1994] Uniform-weight discrete Laplacian
[Taubin 1995] Cotangent-weight discrete Laplacian
[Pinkall and Polthier 1993], [Wardetzky et al. 2007],
[Reuter et al. 2009] Wilmore flow, using FEM with aux variables
Position and co-normal specification on boundary [Clarenz et al. 2004]
Linear systems for k-harmonic equations Uses discretized Laplacian operator
[Botsch and Kobbelt 04]
![Page 6: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/6.jpg)
Previous work Simple domains, analytic boundaries
[Bloor and Wilson 1990] Model shaped minimization of curvature variation energy
[Moreton and Séquin 1992] Interpolate curve networks, local quadratic fits and finite differences
[Welch and Witkin 1994] Uniform-weight discrete Laplacian
[Taubin 1995] Cotangent-weight discrete Laplacian
[Pinkall and Polthier 1993], [Wardetzky et al. 2007],
[Reuter et al. 2009] Wilmore flow, using FEM with aux variables
Position and co-normal specification on boundary [Clarenz et al. 2004]
Linear systems for k-harmonic equations Uses discretized Laplacian operator
[Botsch and Kobbelt 04]
![Page 7: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/7.jpg)
Standard Finite Element Method Requires at least C1 elements for fourth order
Can’t use standard triangle meshes High order surfaces exist, (e.g. Argyris triangle)
Require many extra degrees of freedom Not popular due to complexity
Low order, C0, workarounds E.g. mixed elements
![Page 8: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/8.jpg)
Discrete Geometric Discretization Derive mesh analog of geometric quantity
E.g. Laplace-Beltrami operator integrated over vertex area Re-expressed using only first-order
derivatives Use average value as energy of vertex area
Used often in geometric modeling No obvious way to connect to continuous
case
![Page 9: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/9.jpg)
Mixed Elements Introduce additional
variable to convert high order problem to low order
![Page 10: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/10.jpg)
Use Langrange multipliers to enforce constraint
Constraint structure also makes certain boundary types easier
Mixed Elements Introduce additional
variable to convert high order problem to low order
![Page 11: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/11.jpg)
Our original higher order problem
Introduce an additional variable
Two second order problems Can use just linear elements
Curve Fixed boundary curve Specified tangents:
Mixed Elements
![Page 12: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/12.jpg)
Discretize with piecewise linear approximations for variables
Mixed Elements
Discrete Laplacian Mass matrix
![Page 13: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/13.jpg)
Discretize with piecewise linear approximations for variables
Mixed Elements
Discrete Laplacian Mass matrix
![Page 14: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/14.jpg)
Matrix form, curve boundary conditions
Diagonalized, lumped mass matrices eliminate auxiliary variable
Mixed Elements
Discrete Laplacian Mass matrix Neumann matrix
Where and
![Page 15: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/15.jpg)
Curve Fixed boundary curve Specified tangents:
Point Single fixed points on
surface
Boundary Conditions
![Page 16: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/16.jpg)
Boundary Conditions Region
Fixed part of mesh outside solved region
![Page 17: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/17.jpg)
Use Lagrangian to enforce region condition
Discretize with piecewise linear approximations for variables
May also eliminate aux. variable
Mixed Elements
Discrete Laplacian
Mass matrix
![Page 18: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/18.jpg)
Boundary Conditions
Difference in right-hand side Curve conditions don’t require lumped mass
matrix But we use it in practice, for speed and numerical
accuracy Equivalent to [Botsch and Kobbelt, 2004]
Specified tangents ≈ parameter for continuity control
Region: Curve:
![Page 19: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/19.jpg)
Boundary Conditions Region
Fixed part of mesh outside solved region
![Page 20: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/20.jpg)
Boundary Conditions Convert high order problem to low order problem
Use Langrange multipliers to enforce constraint
![Page 21: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/21.jpg)
Convert high order problem to low order problem
Use Langrange multipliers to enforce constraint
Notice similarity to Lagrangian for biharmonic
Boundary Conditions
![Page 22: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/22.jpg)
Discretization, formulation works the same way
Eliminate auxiliary variables Leaving system with only
Mixed Elements
Discrete Laplacian
Mass matrix
where
![Page 23: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/23.jpg)
Curve Fixed boundary curve Specified tangents and curvatures: ,
Leads to singular systems
Boundary Conditions
![Page 24: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/24.jpg)
Boundary Conditions Curve Region
Fixed boundary curve and one ring into interior
Specified curvatures:
![Page 25: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/25.jpg)
Experimental Results Tested convergence of our systems Randomly generated domains of varying
irregularity One vertex placed randomly in each square of grid Parameter controlled variation from regular
Connected using Triangle Library Control minimal interior angles
![Page 26: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/26.jpg)
Specify boundary conditions using analytic target functions: Try to reproduce original function by solving system:
Measure error between analytic target and our mixed FEM approximation
Experimental Results
![Page 27: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/27.jpg)
Experimental Results Nearly optimal convergence for biharmonic
![Page 28: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/28.jpg)
Experimental Results Boundary conditions perform differently for
triharmonic
![Page 29: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/29.jpg)
Applications Filling in holes: Laplacian energy
input
![Page 30: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/30.jpg)
Applications Filling in holes: Laplacian energy
input region constraint
![Page 31: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/31.jpg)
Applications Filling in holes: Laplacian energy
input region constraint manipulating tangent controls
![Page 32: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/32.jpg)
Applications Filling in holes: Laplacian energy
manipulating tangent controls
![Page 33: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/33.jpg)
Applications Filling in holes: Laplacian gradient energy
input region constraint
manipulating curvature controls
![Page 34: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/34.jpg)
Applications
Specifying tangents in Laplacian energy around regions
![Page 35: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/35.jpg)
Applications
![Page 36: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/36.jpg)
Applications Biharmonic Triharmonic
![Page 37: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/37.jpg)
Summary Technique for discretizing energies or PDEs
Reduce to low order by introducing variables Use constraints to enforce region boundary conditions Lump mass matrix
Convergence for fourth- and sixth-order PDEs
![Page 38: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/38.jpg)
Summary Technique for discretizing energies or PDEs
Reduce to low order by introducing variables Use constraints to enforce region boundary conditions Lump mass matrix
Convergence for fourth- and sixth-order PDEs Future work
Improve convergence of triharmonic solution Explore using non-flat metric
![Page 39: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/39.jpg)
Acknowledgement of funding This work was supported in part by an NSF award
IIS- 0905502.
![Page 40: Mixed Finite Elements for Variational Surface Modeling · 2019-03-06 · Previous work Simple domains, analytic boundaries [Bloor and Wilson 1990] Model shaped minimization of curvature](https://reader034.vdocument.in/reader034/viewer/2022050323/5f7d3f00ab8e2e3a483fd3b3/html5/thumbnails/40.jpg)
Mixed Finite Elements for Variational Surface Modeling
Alec Jacobson ([email protected]) Elif Tosun Olga Sorkine Denis Zorin