BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Martin Rumpf

Martin Rumpf appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical AnalysisComputer Visionmath.APGraphicsmath.OC

Catalog footprint

What is connected

22works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

Finite Element Approximation of Large-Scale Isometric Deformations of Parametrized Surfaces

In this paper, the numerical approximation of isometric deformations of thin elastic shells is discussed. To this end, for a thin shell represented by a parametrized surface, it is shown how to transform the stored elastic energy for an isometric deformation such that the highest order term is quadratic. For this reformulated model, existence of optimal isometric deformations is shown. A finite element approximation is obtained using the Discrete Kirchhoff Triangle (DKT) approach and the convergence of discrete minimizers to a continuous minimizer is demonstrated. In that respect, this paper generalizes the results by Bartels for the approximation of bending isometries of plates. A Newton scheme is derived to numerically simulate large bending isometries of shells. The proven convergence properties are experimentally verified and characteristics of isometric deformations are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Two-scale elastic shape optimization for additive manufacturing

In this paper, a two-scale approach for elastic shape optimization of fine-scale structures in additive manufacturing is investigated. To this end, a free material optimization is performed on the macro-scale using elasticity tensors in a set of microscopically realizable tensors. A database of these realizable tensors and their cost values is obtained with a shape and topology optimization on microscopic cells, working within a fixed set of elasticity tensors samples. This microscopic optimization takes into account manufacturability constraints via predefined material bridges to neighbouring cells at the faces of the microscopic fundamental cell. For the actual additive manufacturing on a chosen fine-scale, a piece-wise constant elasticity tensor ansatz on grid cells of a macroscopic mesh is applied. The macroscopic optimization is performed in an efficient online phase, whereas the associated cell-wise optimal material patterns are retrieved from the database that was computed offline. For that, the set of admissible realizable elasticity tensors is parametrized using tensor product cubic B-splines over the unit square matching the precomputed samples. This representation is then efficiently used in an interior point method for the free material optimization on the macro-scale.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence of the Time Discrete Metamorphosis Model on Hadamard Manifolds

Continuous image morphing is a classical task in image processing. The metamorphosis model proposed by Trouvé, Younes and coworkers casts this problem in the frame of Riemannian geometry and geodesic paths between images. The associated metric in the space of images incorporates dissipation caused by a viscous flow transporting image intensities and its variations along motion paths. In many applications, images are maps from the image domain into a manifold (e.g. in diffusion tensor imaging (DTI) the manifold of symmetric positive definite matrices with a suitable Riemannian metric). In this paper, we propose a generalized metamorphosis model for manifold-valued images, where the range space is a finite-dimensional Hadamard manifold. A corresponding time discrete version was presented by Neumayer et al. based on the general variational time discretization proposed by Berkels et al. Here, we prove the Mosco--convergence of the time discrete metamorphosis functional to the proposed manifold-valued metamorphosis model, which implies the convergence of time discrete geodesic paths to a geodesic path in the (time continuous) metamorphosis model. In particular, the existence of geodesic paths is established. In fact, images as maps into Hadamard manifold are not only relevant in applications, but it is also shown that the joint convexity of the distance function which characterizes Hadamard manifolds is a crucial ingredient to establish existence of the metamorphosis model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Geometric optimization using nonlinear rotation-invariant coordinates

Geometric optimization problems are at the core of many applications in geometry processing. The choice of a representation fitting an optimization problem can considerably simplify solving the problem. We consider the Nonlinear Rotation-Invariant Coordinates (NRIC) that represent the nodal positions of a discrete triangular surface with fixed combinatorics as a vector that stacks all edge lengths and dihedral angles of the mesh. It is known that this representation associates a unique vector to an equivalence class of nodal positions that differ by a rigid body motion. Moreover, integrability conditions that ensure the existence of nodal positions that match a given vector of edge lengths and dihedral angles have been established. The goal of this paper is to develop the machinery needed to use the NRIC for solving geometric optimization problems. First, we use the integrability conditions to derive an implicit description of the space of discrete surfaces as a submanifold of an Euclidean space and a corresponding description of its tangent spaces. Secondly, we reformulate the integrability conditions using quaternions and provide explicit formulas for their first and second derivatives facilitating the use of Hessians in NRIC-based optimization problems. Lastly, we introduce a fast and robust algorithm that reconstructs nodal positions from almost integrable NRIC. We demonstrate the benefits of this approach on a collection of geometric optimization problems. Comparisons to alternative approaches indicate that NRIC-based optimization is particularly effective for problems involving near-isometric deformations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Image Morphing in Deep Feature Spaces: Theory and Applications

This paper combines image metamorphosis with deep features. To this end, images are considered as maps into a high-dimensional feature space and a structure-sensitive, anisotropic flow regularization is incorporated in the metamorphosis model proposed by Miller, Trouvé, Younes and coworkers. For this model a variational time discretization of the Riemannian path energy is presented and the existence of discrete geodesic paths minimizing this energy is demonstrated. Furthermore, convergence of discrete geodesic paths to geodesic paths in the time continuous model is investigated. The spatial discretization is based on a finite difference approximation in image space and a stable spline approximation in deformation space, the fully discrete model is optimized using the iPALM algorithm. Numerical experiments indicate that the incorporation of semantic deep features is superior to intensity-based approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On Material Optimisation for Nonlinearly Elastic Plates and Shells

This paper investigates the optimal distribution of hard and soft material on elastic plates. In the class of isometric deformations stationary points of a Kirchhoff plate functional with incorporated material hardness function are investigated and a compliance cost functional is taken into account. Under symmetry assumptions on the material distribution and the load it is shown that cylindrical solutions are stationary points. Furthermore, it is demonstrated that the optimal design of cylindrically deforming, clamped rectangular plates is non trivial, i.e. with a material distribution which is not just depending on one axial direction on the plate. Analytical results are complemented with numerical optimization results using a suitable finite element discretization and a phase field description of the material phases. Finally, using numerical methods an outlook on the optimal design of non isometrically deforming plates and shells is given.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Generalized optimal transport with singular sources

We present a generalized optimal transport model in which the mass-preserving constraint for the $L^2$-Wasserstein distance is relaxed by introducing a source term in the continuity equation. The source term is also incorporated in the path energy by means of its squared $L^2$-norm in time of a functional with linear growth in space. This extension of the original transport model enables local density modulation, which is a desirable feature in applications such as image warping and blending. A key advantage of the use of a functional with linear growth in space is that it allows for singular sources and sinks, which can be supported on points or lines. On a technical level, the $L^2$-norm in time ensures a disintegration of the source in time, which we use to obtain the well-posedness of the model and the existence of geodesic paths. Furthermore, a numerical scheme based on the proximal splitting approach (Papadakis et al., 2014) is presented. We compare our model with the corresponding model involving the $L^2(L^2)$-norm of the source, which merges the metamorphosis approach and the optimal transport approaches in imaging. Selected numerical test cases show strikingly different behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimization of the branching pattern in coherent phase transitions

Branching can be observed at the austenite-martensite interface of martensitic phase transformations. For a model problem, Kohn and Müller studied a branching pattern with optimal scaling of the energy with respect to its parameters. Here, we present finite element simulations that suggest a topologically different class of branching patterns and derive a novel, low dimensional family of patterns. After a geometric optimization within this family, the resulting pattern bears a striking resemblance to our simulation. The novel microstructure admits the same scaling exponents but results in a significantly lower upper energy bound.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic Dominance Constraints in Elastic Shape Optimization

This paper deals with shape optimization for elastic materials under stochastic loads. It transfers the paradigm of stochastic dominance, which allows for flexible risk aversion via comparison with benchmark random variables, from finite-dimensional stochastic programming to shape optimization. Rather than handling risk aversion in the objective, this enables risk aversion by including dominance constraints that single out subsets of nonanticipative shapes which compare favorably to a chosen stochastic benchmark. This new class of stochastic shape optimization problems arises by optimizing over such feasible sets. The analytical description is built on risk-averse cost measures. The underlying cost functional is of compliance type plus a perimeter term, in the implementation shapes are represented by a phase field which permits an easy estimate of a regularized perimeter. The analytical description and the numerical implementation of dominance constraints are built on risk-averse measures for the cost functional. A suitable numerical discretization is obtained using finite elements both for the displacement and the phase field function. Different numerical experiments demonstrate the potential of the proposed stochastic shape optimization model and in particular the impact of high variability of forces or probabilities in the different realizations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

3D Adaptive Central Schemes: part I Algorithms for Assembling the Dual Mesh

Central schemes are frequently used for incompressible and compressible flow calculations. The present paper is the first in a forthcoming series where a new approach to a 2nd order accurate Finite Volume scheme operating on cartesian grids is discussed. Here we start with an adaptively refined cartesian primal grid in 3D and present a construction technique for the staggered dual grid based on $L^{\infty}$-Voronoi cells. The local refinement constellation on the primal grid leads to a finite number of uniquely defined local patterns on a primal cell. Assembling adjacent local patterns forms the dual grid. All local patterns can be analysed in advance. Later, running the numerical scheme on staggered grids, all necessary geometric information can instantly be retrieved from lookup-tables. The new scheme is compared to established ones in terms of algorithmical complexity and computational effort.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A $BV$ Functional and its Relaxation for Joint Motion Estimation and Image Sequence Recovery

The estimation of motion in an image sequence is a fundamental task in image processing. Frequently, the image sequence is corrupted by noise and one simultaneously asks for the underlying motion field and a restored sequence. In smoothly shaded regions of the restored image sequence the brightness constancy assumption along motion paths leads to a pointwise differential condition on the motion field. At object boundaries which are edge discontinuities both for the image intensity and for the motion field this condition is no longer well defined. In this paper a total-variation type functional is discussed for joint image restoration and motion estimation. This functional turns out not to be lower semicontinuous, and in particular fine-scale oscillations may appear around edges. By the general theory of vector valued $BV$ functionals its relaxation leads to the appearance of a singular part of the energy density, which can be determined by the solution of a local minimization problem at edges. Based on bounds for the singular part of the energy and under appropriate assumptions on the local intensity variation one can exclude the existence of microstructures and obtain a model well-suited for simultaneous image restoration and motion estimation. Indeed, the relaxed model incorporates a generalized variational formulation of the brightness constancy assumption. The analytical findings are related to ambiguity problems in motion estimation such as the proper distinction between foreground and background motion at object edges.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A generalized model for optimal transport of images including dissipation and density modulation

In this paper the optimal transport and the metamorphosis perspectives are combined. For a pair of given input images geodesic paths in the space of images are defined as minimizers of a resulting path energy. To this end, the underlying Riemannian metric measures the rate of transport cost and the rate of viscous dissipation. Furthermore, the model is capable to deal with strongly varying image contrast and explicitly allows for sources and sinks in the transport equations which are incorporated in the metric related to the metamorphosis approach by Trouvé and Younes. In the non-viscous case with source term existence of geodesic paths is proven in the space of measures. The proposed model is explored on the range from merely optimal transport to strongly dissipative dynamics. For this model a robust and effective variational time discretization of geodesic paths is proposed. This requires to minimize a discrete path energy consisting of a sum of consecutive image matching functionals. These functionals are defined on corresponding pairs of intensity functions and on associated pairwise matching deformations. Existence of time discrete geodesics is demonstrated. Furthermore, a finite element implementation is proposed and applied to instructive test cases and to real images. In the non-viscous case this is compared to the algorithm proposed by Benamou and Brenier including a discretization of the source term. Finally, the model is generalized to define discrete weighted barycentres with applications to textures and objects.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Nested Variational Time Discretization for Parametric Anisotropic Willmore Flow

A variational time discretization of anisotropic Willmore flow combined with a spatial discretization via piecewise affine finite elements is presented. Here, both the energy and the metric underlying the gradient flow are anisotropic, which in particular ensures that Wulff shapes are invariant up to scaling under the gradient flow. In each time step of the gradient flow a nested optimization problem has to be solved. Thereby, an outer variational problem reflects the time discretization of the actual Willmore flow and involves an approximate anisotropic $L^2$-distance between two consecutive time steps and a fully implicit approximation of the anisotropic Willmore energy. The anisotropic mean curvature needed to evaluate the energy integrand is replaced by the time discrete, approximate speed from an inner, fully implicit variational scheme for anisotropic mean curvature motion. To solve the nested optimization problem a Newton method for the associated Lagrangian is applied. Computational results for the evolution of curves underline the robustness of the new scheme, in particular with respect to large time steps.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Posteriori Error Control for the Binary Mumford-Shah Model

The binary Mumford-Shah model is a widespread tool for image segmentation and can be considered as a basic model in shape optimization with a broad range of applications in computer vision, ranging from basic segmentation and labeling to object reconstruction. This paper presents robust a posteriori error estimates for a natural error quantity, namely the area of the non properly segmented region. To this end, a suitable strictly convex and non-constrained relaxation of the originally non-convex functional is investigated and Repin's functional approach for a posteriori error estimation is used to control the numerical error for the relaxed problem in the $L^2$-norm. In combination with a suitable cut out argument, a fully practical estimate for the area mismatch is derived. This estimate is incorporated in an adaptive meshing strategy. Two different adaptive primal-dual finite element schemes, and the most frequently used finite difference discretization are investigated and compared. Numerical experiments show qualitative and quantitative properties of the estimates and demonstrate their usefulness in practical applications.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A posteriori error estimates for sequential laminates in shape optimization

A posteriori error estimates are derived in the context of two-dimensional structural elastic shape optimization under the compliance objective. It is known that the optimal shape features are microstructures that can be constructed using sequential lamination. The descriptive parameters explicitly depend on the stress. To derive error estimates the dual weighted residual approach for control problems in PDE constrained optimization is employed, involving the elastic solution and the microstructure parameters. Rigorous estimation of interpolation errors ensures robustness of the estimates while local approximations are used to obtain fully practical error indicators. Numerical results show sharply resolved interfaces between regions of full and intermediate material density.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bézier curves in the space of images

Bézier curves are a widespread tool for the design of curves in Euclidian space. This paper generalizes the notion of Bézier curves to the infinite-dimensional space of images. To this end the space of images is equipped with a Riemannian metric which measures the cost of image transport and intensity variation in the sense of the metamorphosis model by Miller and Younes. Bézier curves are then computed via the Riemannian version of de Casteljau's algorithm, which is based on a hierarchical scheme of convex combination along geodesic curves. Geodesics are approximated using a variational discretization of the Riemannian path energy. This leads to a generalized de Casteljau method to compute suitable discrete Bézier curves in image space. Selected test cases demonstrate qualitative properties of the approach. Furthermore, a Bézier approach for the modulation of face interpolation and shape animation via image sketches is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On Isogeometric Subdivision Methods for PDEs on Surfaces

Subdivision surfaces are proven to be a powerful tool in geometric modeling and computer graphics, due to the great flexibility they offer in capturing irregular topologies. This paper discusses the robust and efficient implementation of an isogeometric discretization approach to partial differential equations on surfaces using subdivision methodology. Elliptic equations with the Laplace-Beltrami and the surface bi-Laplacian operator as well as the associated eigenvalue problems are considered. Thereby, efficiency relies on the proper choice of a numerical quadrature scheme which preserves the expected higher order consistency. A particular emphasis is on the robustness of the approach in the vicinity of extraordinary vertices. In this paper, the focus is on Loop's subdivision scheme on triangular meshes. Based on a series of numerical experiments, different quadrature schemes are compared and a mid-edge quadrature, which is easy-to-implement via lookup tables, turns out to be a preferable choice due to its robustness and efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Time Discrete Geodesic Paths in the Space of Images

In this paper the space of images is considered as a Riemannian manifold using the metamorphosis approach, where the underlying Riemannian metric simultaneously measures the cost of image transport and intensity variation. A robust and effective variational time discretization of geodesics paths is proposed. This requires to minimize a discrete path energy consisting of a sum of consecutive image matching functionals over a set of image intensity maps and pairwise matching deformations. For square-integrable input images the existence of discrete, connecting geodesic paths defined as minimizers of this variational problem is shown. Furthermore, $Γ$-convergence of the underlying discrete path energy to the continuous path energy is proved. This includes a diffeomorphism property for the induced transport and the existence of a square-integrable weak material derivative in space and time. A spatial discretization via finite elements combined with an alternating descent scheme in the set of image intensity maps and the set of matching deformations is presented to approximate discrete geodesic paths numerically. Computational results underline the efficiency of the proposed approach and demonstrate important qualitative properties.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cortical Surface Co-Registration based on MRI Images and Photos

Brain shift, i.e. the change in configuration of the brain after opening the dura mater, is a key problem in neuronavigation. We present an approach to co-register intra-operative microscope images with pre-operative MRI to adapt and optimize intra-operative neuronavigation. The tools are a robust classification of sulci on MRI extracted cortical surfaces, guided user marking of most prominent sulci on a microscope image, and the actual variational registration method with a fidelity energy for 3D deformations of the cortical surface combined with a higher order, linear elastica type prior energy. Furthermore, the actual registration is validated on an artificial testbed with known ground truth deformation and on real data of a neuro clinical patient.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variational time discretization of geodesic calculus

We analyze a variational time discretization of geodesic calculus on finite- and certain classes of infinite-dimensional Riemannian manifolds. We investigate the fundamental properties of discrete geodesics, the associated discrete logarithm, discrete exponential maps, and discrete parallel transport, and we prove convergence to their continuous counterparts. The presented analysis is based on the direct methods in the calculus of variation, on $Γ$-convergence, and on weighted finite element error estimation. The convergence results of the discrete geodesic calculus are experimentally confirmed for a basic model on a two-dimensional Riemannian manifold. This provides a theoretical basis for applications in computer vision: In particular, discrete geodesics offer an effective tool for shape morphing and the computation of a distance between shapes, the discrete logarithm allows a linear representation of strongly nonlinear shape variability, the discrete exponential map provides a robust tool for shape extrapolation, and the discrete parallel transport can be used to transfer geometric details from one shape to another. The basic operations are exemplarily illustrated for two different shape spaces, the space of viscous volumetric objects and the space of discrete viscous shells.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Discrete geodesic calculus in the space of viscous fluidic objects

Based on a local approximation of the Riemannian distance on a manifold by a computationally cheap dissimilarity measure, a time discrete geodesic calculus is developed, and applications to shape space are explored. The dissimilarity measure is derived from a deformation energy whose Hessian reproduces the underlying Riemannian metric, and it is used to define length and energy of discrete paths in shape space. The notion of discrete geodesics defined as energy minimizing paths gives rise to a discrete logarithmic map, a variational definition of a discrete exponential map, and a time discrete parallel transport. This new concept is applied to a shape space in which shapes are considered as boundary contours of physical objects consisting of viscous material. The flexibility and computational efficiency of the approach is demonstrated for topology preserving shape morphing, the representation of paths in shape space via local shape variations as path generators, shape extrapolation via discrete geodesic flow, and the transfer of geometric features.

0 citations0 reviewsNo code linkedOpen access