Riemannian Manifold Learning via Shearlet Approximation

Analyzing Geometry by Applied Harmonic Analysis Methodologies

Applied harmonic analysis provides powerful methodologies to approximate geometric objects, which might be given as a Riemannian manifold itself or as an approximating point cloud. The main tools are specifically designed representation systems such as shearlets. These systems are of a multiscale type, thus an approximation process provides different resolution levels. One might ask: "Which resolution level allows detection of which geometric properties, such as curvature or torsion?" Project A10 aims to analyze such relations between approximations and learning of geometrical properties.


Many processes in geometry and dynamics can be considered as time dependent systems. For instance, dynamic geometry studies elementary geometric constructions in which the base elements of a construction may be moved and the dependent elements move accordingly and consistently. Such a movement may be considered as a process in time. In physics there are also time-dependent processes. One may think of a planet orbiting around a sun or an arrangement of atoms moving under mutual forces. Both systems - moving geometric constructions and particle/force systems - suffer from the fundamental problem that singular situations, which occur at a specific moment in time, may prevent the system from being continued beyond this moment. For example, the orbit of a planet, directly heading into a sun, is stopped at the moment of the collision. The philosophy in research project B01 expands the concept of an ordinary linear time flow to a concept in two dimensions, where the time is represented by a variable over the complex numbers. Complex numbers are inherently two-dimensional objects, but in many aspects they behave like ordinary numbers. In particular all formulas governing geometric constructions or force/particle systems allow direct generalizations to complex numbers. Now, figuratively speaking, the two-dimensional nature of complex numbers allows us to bypass critical singular situations by so-called complex detours. At first, this seems to be a purely esoteric approach to geometry and physics. However, this approach, which is a central element in the research in B01, creates possibilities for new, more stable and consistent methods in dynamic geometry as well as for physics simulations.

Scientific Details+

This project deals with the study of complexified time in ordinary differential equations and other systems that continuously propagate in time. Structural aspects (underlying Riemann surfaces, monodromy, singularities) as well as algorithmic and simulation aspects (complexified integrators, problem of path jumping, discretization, path following, complexity estimates) will be addressed. One of the main aims of the project, is to apply complex detour approaches that have proved fruitful in the context of dynamic geometry to broader contexts including general straight-line and branching programs, numerical integration and the topology of Riemann surfaces. In particular we want to deal with singular situations (or nearly singular situations) by circumventing them in the complex plane rather than applying standard methods like decreasing the stepsize. Applied to numerical integration this kind of complex integrators also allows for a deeper study of the solution structure of a differential equation by providing information on the underlying Riemann surface. In synergy with the project C1 the project also aims for a general algorithmic framework that applies to the treatment of related problems.

One starting point of the project is the experience with complex detour techniques in the context of situations in dynamic geometry. Dynamic geometry deals with geometric configurations (defined by a collection of constraints or a geometric construction) and the modeling of consistent and continuous (i.e. constraint-stable) movements in their configuration space. A typical problem in this context arises if the configuration approaches a singular point of the configuration space, there an intrinsic ambiguity may cause unnatural behavior in the algorithmic simulation of such a configuration. The problem in this context can be completely resolved by embedding the original configuration into a ambient complex space and avoid the singularities by performing complex detours. Besides interesting theoretical insights in the structure of geometric configuration spaces these techniques allow for a satisfactory algorithmic treatment of such problems

One long term aim of the presented research project is to transfer as much as possible from this approach to the situation of ordinary differential equations. In the situation of ODEs one usually assumes a continuous and real propagation of time. Singular or nearly singular situations may make the numerical simulation of an ODE hard or even impossible (like a simulated particle that is moving directly towards the center of an attracting singular force field). Considering the time as a complex variable introduces new freedom that allows for a resolution of the singularity. Singular situations may by this be understood in a broader framework and for nearly singular situations more stable discretization schemes may be found by this approach. Although this approach seems natural and well motivated, there seems to be only few literature on this topic. Several fundamental difficulties arise that at first sight form fundamental obstacles in such a kind of approach. Resolving these difficulties will be a major research goal in this project. Here we list a few corner stones:

  • Avoidance of path jumping: Typical numerical integrates for ODEs (like for instance Runge Kutta methods) involve many evaluations of the right side of the ODE y’(t) = f (g(t)). If the function f is allowed to branch then also the simulation method has to be adapted to this potential branching situation, otherwise path jumping of the solution may occur.
  • Understanding monodromy: The complex interpretation of a singular situation may introduce new monodromy into the problem. This monodromy must be understood as part of the structure inherent to a specific ODE and the numerical simulation has to be adapted to the existence of this monodromy.
  • Analyzing examples; algorithmic complexity: As soon as reliable algorithmic methods in the above sense are available we will investigate specific interesting examples. We also aim for intrinsic algorithmic bounds for these examples. The extraction of crucial examples that can serve as benchmarks for the developed methods is also of relevance.
  • Structure preserving discretization: In particular (and in relation with project A1) we aim for the development for structure preserving or symplectic discretization techniques that still work well for ambient complex spaces of an ODE.


Shearlet approximation of functions with discontinuous derivatives

Author: Petersen, Philipp
Note: preprint
Date: Aug 2015
Download: arXiv

Anisotropic multiscale systems on bounded domains

Authors: Grohs, P. and Kutyniok, G. and Ma, J. and Petersen, P.
Note: preprint
Date: 2015
Download: arXiv

Classification of edges using compactly supported shearlets

Authors: Kutyniok, Gitta and Petersen, Philipp
Journal: Applied and Computational Harmonic Analysis
Date: 2015
DOI: 10.1016/j.acha.2015.08.006
Download: external arXiv

Linear independence of compactly supported separable shearlet systems

Authors: Ma, Jackie and Petersen, Philipp
Journal: Journal of Mathematical Analysis and Applications, 428:238-257
Date: 2015
DOI: 10.1016/j.jmaa.2015.03.001
Download: external arXiv


Authors: Grohs, Philipp and Keiper, Sandra and Kutyniok, Gitta and Schäfer, Martin
Note: preprint
Date: Jul 2014
Download: arXiv

Regularization and Numerical Solution of the Inverse Scattering Problem using Shearlet Frames

Authors: Kutyniok, Gitta and Mehrmann, Volker and Petersen, Philipp
Note: preprint
Date: Jul 2014
Download: arXiv

Dualizable Shearlet Frames and Sparse Approximation

Authors: Kutyniok, Gitta and Lim, Wang-Q
Note: preprint
Date: 2014
Download: internal

ShearLab 3D: Faithful Digital Shearlet Transforms based on Compactly Supported Shearlets

Authors: Kutyniok, Gitta and Lim, Wang-Q and Reisenhofer, Rafael
Note: Preprint
Date: Jan 2014
Download: internal arXiv

Gabor Shearlets

Authors: Bodmann, Bernhard G. and Kutyniok, Gitta and Zhuang, Xiaosheng
Journal: Appl. Comput. Harmon. Anal.
Note: submitted
Date: Mar 2013
Download: external arXiv


Prof. Dr. Gitta Kutyniok   +

Jackie Ma   +

Projects: A10, C03
University: TU Berlin
E-Mail: ma[at]math.tu-berlin.de

Dr. Philipp Petersen   +

Projects: A10, C03
University: TU Berlin
E-Mail: petersen[at]math.tu-berlin.de