Source author record

Hanne Kekkonen

Hanne Kekkonen 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

6works
6topics
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

6 published item(s)

preprint2022arXiv

Exploring Mathematics with Curvagon Tiles

Building blocks and tiles are an excellent way of learning about geometry and mathematics in general. There are several versions of tiles that are either snapped together or connected with magnets that can be used to introduce topics like volume, tessellations, and Platonic solids. However, since these tiles are made of hard plastic, they are not very suitable for creating hyperbolic surfaces or shapes where the tiles need to bend. Curvagons are flexible regular polygon building blocks that allow you to quickly build anything from hyperbolic surfaces and tori to dinosaurs and shoes. They can be used to introduce mathematical concepts from Archimedean solids to Gauss-Bonnet theorem. You can also let your imagination run free and build whatever comes to mind.

preprint2021arXiv

Consistency of Bayesian inference with Gaussian process priors for a parabolic inverse problem

We consider the statistical nonlinear inverse problem of recovering the absorption term $f>0$ in the heat equation $$ \partial_tu-\frac{1}{2}Δu+fu=0 \quad \text{on $\mathcal{O}\times(0,\textbf{T})$}\quad u = g \quad \text{on $\partial\mathcal{O}\times(0,\textbf{T})$}\quad u(\cdot,0)=u_0 \quad \text{on $\mathcal{O}$}, $$ where $\mathcal{O}\in\mathbb{R}^d$ is a bounded domain, $\textbf{T}<\infty$ is a fixed time, and $g,u_0$ are given sufficiently smooth functions describing boundary and initial values respectively. The data consists of $N$ discrete noisy point evaluations of the solution $u_f$ on $\mathcal{O}\times(0,\textbf{T})$. We study the statistical performance of Bayesian nonparametric procedures based on a large class of Gaussian process priors. We show that, as the number of measurements increases, the resulting posterior distributions concentrate around the true parameter generating the data, and derive a convergence rate for the reconstruction error of the associated posterior means. We also consider the optimality of the contraction rates and prove a lower bound for the minimax convergence rate for inferring $f$ from the data, and show that optimal rates can be achieved with truncated Gaussian priors.

preprint2021arXiv

Random tree Besov priors -- Towards fractal imaging

We propose alternatives to Bayesian a priori distributions that are frequently used in the study of inverse problems. Our aim is to construct priors that have similar good edge-preserving properties as total variation or Mumford-Shah priors but correspond to well defined infinite-dimensional random variables, and can be approximated by finite-dimensional random variables. We introduce a new wavelet-based model, where the non zero coefficient are chosen in a systematic way so that prior draws have certain fractal behaviour. We show that realisations of this new prior take values in some Besov spaces and have singularities only on a small set $τ$ that has a certain Hausdorff dimension. We also introduce an efficient algorithm for calculating the MAP estimator, arising from the the new prior, in denoising problem.

preprint2020arXiv

Bernstein-von Mises theorems and uncertainty quantification for linear inverse problems

We consider the statistical inverse problem of recovering an unknown function $f$ from a linear measurement corrupted by additive Gaussian white noise. We employ a nonparametric Bayesian approach with standard Gaussian priors, for which the posterior-based reconstruction of $f$ corresponds to a Tikhonov regulariser $\bar f$ with a reproducing kernel Hilbert space norm penalty. We prove a semiparametric Bernstein-von Mises theorem for a large collection of linear functionals of $f$, implying that semiparametric posterior estimation and uncertainty quantification are valid and optimal from a frequentist point of view. The result is applied to study three concrete examples that cover both the mildly and severely ill-posed cases: specifically, an elliptic inverse problem, an elliptic boundary value problem and the heat equation. For the elliptic boundary value problem, we also obtain a nonparametric version of the theorem that entails the convergence of the posterior distribution to a prior-independent infinite-dimensional Gaussian probability measure with minimal covariance. As a consequence, it follows that the Tikhonov regulariser $\bar f$ is an efficient estimator of $f$, and we derive frequentist guarantees for certain credible balls centred at $\bar{f}$.

preprint2016arXiv

Analysis of regularized inversion of data corrupted by white Gaussian noise

Tikhonov regularization is studied in the case of linear pseudodifferential operator as the forward map and additive white Gaussian noise as the measurement error. The measurement model for an unknown function $u(x)$ is \begin{eqnarray*} m(x) = Au(x) + δ\hspace{.2mm}\varepsilon(x), \end{eqnarray*} where $δ>0$ is the noise magnitude. If $\varepsilon$ was an $L^2$-function, Tikhonov regularization gives an estimate \begin{eqnarray*} T_α(m) = \text{argmin}_{u\in H^r}\big\{\|A u-m\|_{L^2}^2+ α\|u\|_{H^r}^2 \big\}\end{eqnarray*} for $u$ where $α=α(δ)$ is the regularization parameter. Here penalization of the Sobolev norm $ \|u\|_{H^r}$ covers the cases of standard Tikhonov regularization ($r=0$) and first derivative penalty ($r=1$). Realizations of white Gaussian noise are almost never in $L^2$, but do belong to $H^s$ with probability one if $s<0$ is small enough. A modification of Tikhonov regularization theory is presented, covering the case of white Gaussian measurement noise. Furthermore, the convergence of regularized reconstructions to the correct solution as $δ\rightarrow 0$ is proven in appropriate function spaces using microlocal analysis. The convergence of the related finite-dimensional problems to the infinite-dimensional problem is also analysed.

preprint2016arXiv

Posterior consistency and convergence rates for Bayesian inversion with hypoelliptic operators

Bayesian approach to inverse problems is studied in the case where the forward map is a linear hypoelliptic pseudodifferential operator and measurement error is additive white Gaussian noise. The measurement model for an unknown Gaussian random variable $U(x,ω)$ is \begin{eqnarray*} M(y,ω) = A(U(x,ω) )+ δ\hspace{.2mm}\mathcal{E}(y,ω), \end{eqnarray*} where $A$ is a finitely many times smoothing linear hypoelliptic operator and $δ>0$ is the noise magnitude. The covariance operator $C_U$ of $U$ is $2r$ times smoothing, self-adjoint, injective and elliptic pseudodifferential operator. If $\mathcal{E}$ was taking values in $L^2$ then in Gaussian case solving the conditional mean (and maximum a posteriori) estimate is linked to solving the minimisation problem \begin{eqnarray*} T_δ(M) = \text{argmin}_{u\in H^r} \big\{\|A u-m\|_{L^2}^2+ δ^2\|C_U^{-1/2}u\|_{L^2}^2 \big\}. \end{eqnarray*} However, Gaussian white noise does not take values in $L^2$ but in $H^{-s}$ where $s>0$ is big enough. A modification of the above approach to solve the inverse problem is presented, covering the case of white Gaussian measurement noise. Furthermore, the convergence of conditional mean estimate to the correct solution as $δ\rightarrow 0$ is proven in appropriate function spaces using microlocal analysis. Also the contraction of the confidence regions is studied.