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

Olivier Verdier

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

ResearcherUnclaimed source record
math.NAmath-phmath.MPComputer Visionmath.APmath.DGmath.SGeess.IVMachine Learningmath.CTmath.DSmath.OCmath.RAmath.RT

Catalog footprint

What is connected

20works
14topics
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

20 published item(s)

preprint2019arXiv

Invariant connections, Lie algebra actions, and foundations of numerical integration on manifolds

Motivated by numerical integration on manifolds, we relate the algebraic properties of invariant connections to their geometric properties. Using this perspective, we generalize some classical results of Cartan and Nomizu to invariant connections on algebroids. This has fundamental consequences for the theory of numerical integrators, giving a characterization of the spaces on which Butcher and Lie-Butcher series methods, which generalize Runge-Kutta methods, may be applied.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Spatiotemporal PET reconstruction using ML-EM with learned diffeomorphic deformation

Patient movement in emission tomography deteriorates reconstruction quality because of motion blur. Gating the data improves the situation somewhat: each gate contains a movement phase which is approximately stationary. A standard method is to use only the data from a few gates, with little movement between them. However, the corresponding loss of data entails an increase of noise. Motion correction algorithms have been implemented to take into account all the gated data, but they do not scale well, especially not in 3D. We propose a novel motion correction algorithm which addresses the scalability issue. Our approach is to combine an enhanced ML-EM algorithm with deep learning based movement registration. The training is unsupervised, and with artificial data. We expect this approach to scale very well to higher resolutions and to 3D, as the overall cost of our algorithm is only marginally greater than that of a standard ML-EM algorithm. We show that we can significantly decrease the noise corresponding to a limited number of gates.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The ML-EM algorithm in continuum: sparse measure solutions

Linear inverse problems $A μ= δ$ with Poisson noise and non-negative unknown $μ\geq 0$ are ubiquitous in applications, for instance in Positron Emission Tomography (PET) in medical imaging. The associated maximum likelihood problem is routinely solved using an expectation-maximisation algorithm (ML-EM). This typically results in images which look spiky, even with early stopping. We give an explanation for this phenomenon. We first regard the image $μ$ as a measure. We prove that if the measurements $δ$ are not in the cone $\{A μ, μ\geq 0\}$, which is typical of short exposure times, likelihood maximisers as well as ML-EM cluster points must be sparse, i.e., typically a sum of point masses. On the other hand, in the long exposure regime, we prove that cluster points of ML-EM will be measures without singular part. Finally, we provide concentration bounds for the probability to be in the sparse case.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Full affine equivariance and weak natural transformations in numerical analysis - the case of B-Series

Many algorithms in numerical analysis are affine equivariant: they are immune to changes of affine coordinates. This is because those algorithms are defined using affine invariant constructions. There is, however, a crucial ingredient missing: most algorithms are in fact defined regardless of the underlying dimension. As a result, they are also invariant with respect to non-invertible affine transformation from spaces of different dimensions. We formulate this property precisely: these algorithms fall short of being natural transformations between affine functors. We give a precise definition of what we call a weak natural transformation between functors, and illustrate the point using examples coming from numerical analysis, in particular B-Series.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Integrators on homogeneous spaces: Isotropy choice and connections

We consider numerical integrators of ODEs on homogeneous spaces (spheres, affine spaces, hyperbolic spaces). Homogeneous spaces are equipped with a built-in symmetry. A numerical integrator respects this symmetry if it is equivariant. One obtains homogeneous space integrators by combining a Lie group integrator with an isotropy choice. We show that equivariant isotropy choices combined with equivariant Lie group integrators produce equivariant homogeneous space integrators. Moreover, we show that the RKMK, Crouch--Grossman or commutator-free methods are equivariant. To show this, we give a novel description of Lie group integrators in terms of stage trees and motion maps, which unifies the known Lie group integrators. We then proceed to study the equivariant isotropy maps of order zero, which we call connections, and show that they can be identified with reductive structures and invariant principal connections. We give concrete formulas for connections in standard homogeneous spaces of interest, such as Stiefel, Grassmannian, isospectral, and polar decomposition manifolds. Finally, we show that the space of matrices of fixed rank possesses no connection.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Aromatic Butcher Series

We show that without other further assumption than affine equivariance and locality, a numerical integrator has an expansion in a generalized form of Butcher series (B-series) which we call aromatic B-series. We obtain an explicit description of aromatic B-series in terms of elementary differentials associated to aromatic trees, which are directed graphs generalizing trees. We also define a new class of integrators, the class of aromatic Runge-Kutta methods, that extends the class of Runge-Kutta methods, and have aromatic B-series expansion but are not B-series methods. Finally, those results are partially extended to the case of more general affine group equivariance.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High order semi-Lagrangian methods for the incompressible Navier-Stokes equations

We propose a class of semi-Lagrangian methods of high approximation order in space and time, based on spectral element space discretizations and exponential integrators of Runge-Kutta type. We discuss the extension of these methods to the Navier-Stokes equations, and their implementation using projections. Semi-Lagrangian methods up to order three are implemented and tested on various examples. The good performance of the methods for convection-dominated problems is demonstrated with numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multi-symplectic discretisation of wave map equations

We present a new multi-symplectic formulation of constrained Hamiltonian partial differential equations, and we study the associated local conservation laws. A multi-symplectic discretisation based on this new formulation is exemplified by means of the Euler box scheme. When applied to the wave map equation, this numerical scheme is explicit, preserves the constraint and can be seen as a generalisation of the Shake algorithm for constrained mechanical systems. Furthermore, numerical experiments show excellent conservation properties of the numerical solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A classification of volume preserving generating forms in R^3

In earlier work, Lomeli and Meiss used a generalization of the symplectic approach to study volume preserving generating differential forms. In particular, for the $\mathbb{R}^3$ case, the first to differ from the symplectic case, they derived thirty-six one-forms that generate exact volume preserving maps. Xue and Zanna had studied these differential forms in connection with the numerical solution of divergence-free differential equations: can such forms be used to devise new volume preserving integrators or to further understand existing ones? As a partial answer to this question, Xue and Zanna showed how six of the generating volume form were naturally associated to consistent, first order, volume preserving numerical integrators. In this paper, we investigate and classify the remaining cases. The main result is the reduction of the thirty-six cases to five essentially different cases, up to variable relabeling and adjunction. We classify these five cases, identifying two novel classes and associating the other three to volume preserving vector fields under a Hamiltonian or Lagrangian representation. We demonstrate how these generating form lead to consistent volume preserving schemes for volume preserving vector fields in $\mathbb{R}^3$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Reductions of Operator Pencils

We study problems associated with an operator pencil, i.e., a pair of operators on Banach spaces. Two natural problems to consider are linear constrained differential equations and the description of the generalized spectrum. The main tool to tackle either of those problems is the reduction of the pencil. There are two kinds of natural reduction operations associated to a pencil, which are conjugate to each other. Our main result is that those two kinds of reductions commute, under some mild assumptions that we investigate thoroughly. Each reduction exhibits moreover a pivot operator. The invertibility of all the pivot operators of all possible successive reductions corresponds to the notion of regular pencil in the finite dimensional case, and to the inf-sup condition for saddle point problems on Hilbert spaces. Finally, we show how to use the reduction and the pivot operators to describe the generalized spectrum of the pencil.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Reduction and Normal Forms of Matrix Pencils

Matrix pencils, or pairs of matrices, may be used in a variety of applications. In particular, a pair of matrices (E,A) may be interpreted as the differential equation E x' + A x = 0. Such an equation is invariant by changes of variables, or linear combination of the equations. This change of variables or equations is associated to a group action. The invariants corresponding to this group action are well known, namely the Kronecker indices and divisors. Similarly, for another group action corresponding to the weak equivalence, a complete set of invariants is also known, among others the strangeness. We show how to define those invariants in a directly invariant fashion, i.e. without using a basis or an extra Euclidean structure. To this end, we will define a reduction process which produces a new system out of the original one. The various invariants may then be defined from operators related to the repeated application of the reduction process. We then show the relation between the invariants and the reduced subspace dimensions, and the relation with the regular pencil condition. This is all done using invariant tools only. Making special choices of basis then allows to construct the Kronecker canonical form. In a related manner, we construct the strangeness canonical form associated to weak equivalence.

0 citations0 reviewsNo code linkedOpen access