Source author record

Klas Modin

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

15 published item(s)

preprint2020arXiv

A Multiscale Theory for Image Registration and Nonlinear Inverse Problems

In an influential paper, Tadmor, Nezzar and Vese (Multiscale Model. Simul. (2004)) introduced a hierarchical decomposition of an image as a sum of constituents of different scales. Here we construct analogous hierarchical expansions for diffeomorphisms, in the context of image registration, with the sum replaced by composition of maps. We treat this as a special case of a general framework for multiscale decompositions, applicable to a wide range of imaging and nonlinear inverse problems. As a paradigmatic example of the latter, we consider the Calderón inverse conductivity problem. We prove that we can simultaneously perform a numerical reconstruction and a multiscale decomposition of the unknown conductivity, driven by the inverse problem itself. We provide novel convergence proofs which work in the general abstract settings, yet are sharp enough to settle an open problem on the hierarchical decompostion of Tadmor, Nezzar and Vese for arbitrary functions in $L^2$. We also give counterexamples that show the optimality of our general results.

preprint2016arXiv

Discrete Variational Derivative Methods for the EPDiff equation

The aim of this paper is the derivation of structure preserving schemes for the solution of the EPDiff equation, with particular emphasis on the two dimensional case. We develop three different schemes based on the Discrete Variational Derivative Method (DVDM) on a rectangular domain discretized with a regular, structured, orthogonal grid. We present numerical experiments to support our claims: we investigate the preservation of energy and linear momenta, the reversibility, and the empirical convergence of the schemes. The quality of our schemes is finally tested by simulating the interaction of singular wave fronts.

preprint2015arXiv

Diffeomorphic density matching by optimal information transport

We address the following problem: given two smooth densities on a manifold, find an optimal diffeomorphism that transforms one density into the other. Our framework builds on connections between the Fisher-Rao information metric on the space of probability densities and right-invariant metrics on the infinite-dimensional manifold of diffeomorphisms. This optimal information transport, and modifications thereof, allows us to construct numerical algorithms for density matching. The algorithms are inherently more efficient than those based on optimal mass transport or diffeomorphic registration. Our methods have applications in medical image registration, texture mapping, image morphing, non-uniform random sampling, and mesh adaptivity. Some of these applications are illustrated in examples.

preprint2015arXiv

Geometry of discrete-time spin systems

Classical Hamiltonian spin systems are continuous dynamical systems on the symplectic phase space $(S^2)^n$. In this paper we investigate the underlying geometry of a time discretization scheme for classical Hamiltonian spin systems called the spherical midpoint method. As it turns out, this method displays a range of interesting geometrical features, that yield insights and sets out general strategies for geometric time discretizations of Hamiltonian systems on non-canonical symplectic manifolds. In particular, our study provides two new, completely geometric proofs that the discrete-time spin systems obtained by the spherical midpoint method preserve symplecticity. The study follows two paths. First, we introduce an extended version of the Hopf fibration to show that the spherical midpoint method can be seen as originating from the classical midpoint method on $T^*\mathbf{R}^{2n}$ for a collective Hamiltonian. Symplecticity is then a direct, geometric consequence. Second, we propose a new discretization scheme on Riemannian manifolds called the Riemannian midpoint method. We determine its properties with respect to isometries and Riemannian submersions and, as a special case, we show that the spherical midpoint method is of this type for a non-Euclidean metric. In combination with Kähler geometry, this provides another geometric proof of symplecticity.

preprint2014arXiv

Generalised Hunter-Saxton equations, optimal information transport, and factorisation of diffeomorphisms

We study geodesic equations for a family of right-invariant Riemannian metrics on the group of diffeomorphisms of a compact manifold. The metrics descend to Fisher's information metric on the space of smooth probability densities. The right reduced geodesic equations are higher-dimensional generalisations of the $μ$--Hunter--Saxton equation, used to model liquid crystals under influence of magnetic fields. Local existence and uniqueness results are established by proving smoothness of the geodesic spray. The descending property of the metrics is used to obtain a novel factorisation of diffeomorphisms. Analogous to the polar factorisation in optimal mass transport, this factorisation solves an optimal information transport problem. It can be seen as an infinite-dimensional version of $QR$ factorisation of matrices.

preprint2014arXiv

Integrability of Nonholonomically Coupled Oscillators

We study a family of nonholonomic mechanical systems. These systems consist of harmonic oscillators coupled through nonholonomic constraints. In particular, the family includes the so called contact oscillator, which has been used as a test problem for numerical methods for nonholonomic mechanics. Furthermore, the systems under study constitute simple models for continuously variable transmission gearboxes. The main result is that each system in the family is integrable reversible with respect to the canonical reversibility map on the cotangent bundle. This result may explain previous numerical observations, that some discretisations of the contact oscillator have favourable structure preserving properties.

preprint2013arXiv

Collective symplectic integrators

We construct symplectic integrators for Lie-Poisson systems. The integrators are standard symplectic (partitioned) Runge--Kutta methods. Their phase space is a symplectic vector space with a Hamiltonian action with momentum map $J$ whose range is the target Lie--Poisson manifold, and their Hamiltonian is collective, that is, it is the target Hamiltonian pulled back by $J$. The method yields, for example, a symplectic midpoint rule expressed in 4 variables for arbitrary Hamiltonians on $\mathfrak{so}(3)^*$. The method specializes in the case that a sufficiently large symmetry group acts on the fibres of $J$, and generalizes to the case that the vector space carries a bifoliation. Examples involving many classical groups are presented.

preprint2013arXiv

Geometric Generalisations of SHAKE and RATTLE

A geometric analysis of the Shake and Rattle methods for constrained Hamiltonian problems is carried out. The study reveals the underlying differential geometric foundation of the two methods, and the exact relation between them. In addition, the geometric insight naturally generalises Shake and Rattle to allow for a strictly larger class of constrained Hamiltonian systems than in the classical setting. In order for Shake and Rattle to be well defined, two basic assumptions are needed. First, a nondegeneracy assumption, which is a condition on the Hamiltonian, i.e., on the dynamics of the system. Second, a coisotropy assumption, which is a condition on the geometry of the constrained phase space. Non-trivial examples of systems fulfilling, and failing to fulfill, these assumptions are given.

preprint2012arXiv

Geodesic Warps by Conformal Mappings

In recent years there has been considerable interest in methods for diffeomorphic warping of images, with applications e.g.\ in medical imaging and evolutionary biology. The original work generally cited is that of the evolutionary biologist D'Arcy Wentworth Thompson, who demonstrated warps to deform images of one species into another. However, unlike the deformations in modern methods, which are drawn from the full set of diffeomorphism, he deliberately chose lower-dimensional sets of transformations, such as planar conformal mappings. In this paper we study warps of such conformal mappings. The approach is to equip the infinite dimensional manifold of conformal embeddings with a Riemannian metric, and then use the corresponding geodesic equation in order to obtain diffeomorphic warps. After deriving the geodesic equation, a numerical discretisation method is developed. Several examples of geodesic warps are then given. We also show that the equation admits totally geodesic solutions corresponding to scaling and translation, but not to affine transformations.

preprint2012arXiv

On Hodge decomposition and conformal variational problems

The main result is the identification of the orthogonal complement of the subalgebra of conformal vector field inside the algebra of all vector fields of a compact flat 2-manifold. As a fundamental tool, the complete Hodge decomposition for manifold with boundary is used. The identification allows the derivation of governing differential equations for variational problems on the space of conformal vector fields. Several examples are given. In addition, the paper also gives a review, in full detail, of already known vector field decompositions involving subalgebras of volume preserving and symplectic vector fields.

preprint2011arXiv

Geometric Integration of Hamiltonian Systems Perturbed by Rayleigh Damping

Explicit and semi-explicit geometric integration schemes for dissipative perturbations of Hamiltonian systems are analyzed. The dissipation is characterized by a small parameter $ε$, and the schemes under study preserve the symplectic structure in the case $ε=0$. In the case $0<ε\ll 1$ the energy dissipation rate is shown to be asymptotically correct by backward error analysis. Theoretical results on monotone decrease of the modified Hamiltonian function for small enough step sizes are given. Further, an analysis proving near conservation of relative equilibria for small enough step sizes is conducted. Numerical examples, verifying the analyses, are given for a planar pendulum and an elastic 3--D pendulum. The results are superior in comparison with a conventional explicit Runge-Kutta method of the same order.