Source author record

Antti Hannukainen

Antti Hannukainen 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

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

7 published item(s)

preprint2026arXiv

Hybrid Nitsche method for distributed computing

We extend a distributed finite element method built upon model order reduction to arbitrary polynomial degree using a hybrid Nitsche scheme. The new method considerably simplifies the transformation of the finite element system to the reduced basis for large problems. We prove that the error of the reduced Nitsche solution converges optimally with respect to the approximation order of the finite element spaces and linearly with respect to the dimension reduction parameter. Numerical tests with nontrivial tetrahedral meshes using second-degree polynomial bases support the theoretical results.

preprint2020arXiv

Efficient solution of symmetric eigenvalue problems from families of coupled systems

Efficient solution of the lowest eigenmodes is studied for a family of related eigenvalue problems with common $2\times 2$ block structure. It is assumed that the upper diagonal block varies between different versions while the lower diagonal block and the range of the coupling blocks remains unchanged. Such block structure naturally arises when studying the effect of a subsystem to the eigenmodes of the full system. The proposed method is based on interpolation of the resolvent function after some of its singularities have been removed by a spectral projection. Singular value decomposition can be used to further reduce the dimension of the computational problem. Error analysis of the method indicates exponential convergence with respect to the number of interpolation points. Theoretical results are illustrated by two numerical examples related to finite element discretisation of the Laplace operator.

preprint2016arXiv

Efficient inclusion of total variation type priors in quantitative photoacoustic tomography

Quantitative photoacoustic tomography is an emerging imaging technique aimed at estimating the distribution of optical parameters inside tissues from photoacoustic images, which are formed by combining optical information and ultrasonic propagation. This optical parameter estimation problem is ill-posed and needs to be approached within the framework of inverse problems. Photoacoustic images are three-dimensional and high-resolution. Furthermore, high-resolution reconstructions of the optical parameters are targeted. Therefore, in order to provide a practical method for quantitative photoacoustic tomography, the inversion algorithm needs to be able to perform successfully with problems of prominent size. In this work, an efficient approach for the inverse problem of quantitative photoacoustic tomography is proposed, assuming an edge-preferring prior for the optical parameters. The method is based on iteratively combining priorconditioned LSQR with a lagged diffusivity step and a linearisation of the measurement model, with the needed multiplications by Jacobians performed in a matrix-free manner. The algorithm is tested with three-dimensional numerical simulations. The results show that the approach can be used to produce accurate and high quality estimates of absorption and diffusion in complex three-dimensional geometries with moderate computation time and cost.

preprint2015arXiv

Convergence analysis of GMRES for the Helmholtz equation via pseudospectrum

Most finite element methods for solving time-harmonic wave-propagation problems lead to a linear system with a non-normal coefficient matrix. The non-normality is due to boundary conditions and losses. One way to solve these systems is to use a preconditioned iterative method. Detailed mathematical analysis of the convergence properties of these methods is important for developing new and understanding old preconditioners. Due to non-normality, there is currently very little existing literature in this direction. In this paper, we study the convergence of GMRES for such systems by deriving inclusion and exclusion regions for the pseudospectrum of the coefficient matrix. All analysis is done a priori by relating the properties of the weak problem to the coefficient matrix. The inclusion is derived from the stability properties of the problem and the exclusion is established via field of values and boundedness of the weak form. The derived tools are applied to estimate the pseudospectrum of time-harmonic Helmholtz equation with first-order absorbing boundary conditions, with and without a shifted-Laplace preconditioner.

preprint2015arXiv

Edge-promoting reconstruction of absorption and diffusivity in optical tomography

In optical tomography a physical body is illuminated with near-infrared light and the resulting outward photon flux is measured at the object boundary. The goal is to reconstruct internal optical properties of the body, such as absorption and diffusivity. In this work, it is assumed that the imaged object is composed of an approximately homogeneous background with clearly distinguishable embedded inhomogeneities. An algorithm for finding the maximum a posteriori estimate for the absorption and diffusion coefficients is introduced assuming an edge-preferring prior and an additive Gaussian measurement noise model. The method is based on iteratively combining a lagged diffusivity step and a linearization of the measurement model of diffuse optical tomography with priorconditioned LSQR. The performance of the reconstruction technique is tested via three-dimensional numerical experiments with simulated measurement data.

preprint2015arXiv

Spectral Study of the Vocal Tract in Vowel Synthesis: A Comparison between 1D and 3D Acoustic Analysis

A state-of-the-art 1D acoustic synthesizer has been previously developed, and coupled to speaker-specific biomechanical models of oropharynx in ArtiSynth. As expected, the formant frequencies of the synthesized vowel sounds were shown to be different from those of the recorded audio. Such discrepancy was hypothesized to be due to the simplified geometry of the vocal tract model as well as the one dimensional implementation of Navier-Stokes equations. In this paper, we calculate Helmholtz resonances of our vocal tract geometries using 3D finite element method (FEM), and compare them with the formant frequencies obtained from the 1D method and audio. We hope such comparison helps with clarifying the limitations of our current models and/or speech synthesizer.

preprint2011arXiv

Field of values analysis of preconditioners for the Helmholtz equation in lossy media

In this paper, we analyze the convergence of the preconditioned GMRES method for the first order finite element discretizations of the Helmholtz equation in media with losses. We consider a Laplace preconditioner, an inexact Laplace preconditioner and a two-level preconditioner. Our analysis is based on bounding the field of values of the preconditioned system matrix in the complex plane. The analysis takes the non-normal nature of the linear system naturally into account and allows us to easily consider certain type of inexact Laplace preconditioners via a perturbation argument. For the two-level preconditioner, our convergence analysis takes into account a media, which has not been considered in previous works.