Manuel Graef

Martin Ehler, Manuel Graef, Chris. J. Oates

The worst case integration error in reproducing kernel Hilbert spaces of standard Monte Carlo methods with n random points decays as $n^{-1/2}$. However, re-weighting of random points can sometimes be used to improve the convergence order. This paper contributes general theoretical results for Sobolev spaces on closed Riemannian manifolds, where we verify that such re-weighting yields optimal approximation rates up to a logarithmic factor. We also provide numerical experiments matching the theoretical results for some Sobolev spaces on the unit sphere and on the Grassmannian manifold. Our theoretical findings also cover function spaces on more general sets such as the unit ball, the cube, and the simplex.

Anna Breger, Martin Ehler, Manuel Graef

Numerical integration and function approximation on compact Riemannian manifolds based on eigenfunctions of the Laplace-Beltrami operator have been widely studied in the recent literature. The standard example in numerical experiments is the Euclidean sphere. Here, we derive numerically feasible expressions for the approximation schemes on the Grassmannian manifold, and we present the associated numerical experiments on the Grassmannian. Indeed, our experiments illustrate and match the corresponding theoretical results in the literature.

Anna Breger, Martin Ehler, Manuel Graef et al.

We briefly explain the use of cubature points on Grassmannians in numerical analysis.

Bernhard G. Bodmann, Martin Ehler, Manuel Graef

We aim to compute the first few moments of a high-dimensional random vector from the first few moments of a number of its low-dimensional projections. To this end, we identify algebraic conditions on the set of low-dimensional projectors that yield explicit reconstruction formulas. We also provide a computational framework, with which suitable projectors can be derived by solving an optimization problem. Finally, we show that randomized projections permit approximate recovery.

Anna Breger, Martin Ehler, Manuel Graef

Given a finite set of points on the Euclidean sphere, the worst case quadrature error in Sobolev spaces has recently been shown to provide upper bounds on the covering radius of the point set. Moreover, quasi-Monte Carlo integration points on the sphere achieve the asymptotically optimal covering radius. Here, we extend these results to points on compact smooth Riemannian manifolds and provide numerical experiments illustrating our findings for the Grassmannian manifold.

Martin Ehler, Manuel Graef, Franz J. Kiraly

As a generalization of the standard phase retrieval problem, we seek to reconstruct symmetric rank-1 matrices from inner products with subclasses of positive semidefinite matrices. For such subclasses, we introduce random cubatures for spaces of multivariate polynomials based on moment conditions. The inner products with samples from sufficiently strong random cubatures allow the reconstruction of symmetric rank-1 matrices with a decent probability by solving the feasibility problem of a semidefinite program.