Source author record

Vesa Kaarnioja

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

math.NA math.CO Numerical Analysis math.NT

Catalog footprint

Research overview

5works

4topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2021arXiv

Fast approximation by periodic kernel-based lattice-point interpolation with application in uncertainty quantification

This paper deals with the kernel-based approximation of a multivariate periodic function by interpolation at the points of an integration lattice -- a setting that, as pointed out by Zeng, Leung, Hickernell (MCQMC2004, 2006) and Zeng, Kritzer, Hickernell (Constr. Approx., 2009), allows fast evaluation by fast Fourier transform, so avoiding the need for a linear solver. The main contribution of the paper is the application to the approximation problem for uncertainty quantification of elliptic partial differential equations, with the diffusion coefficient given by a random field that is periodic in the stochastic variables, in the model proposed recently by Kaarnioja, Kuo, Sloan (SIAM J. Numer. Anal., 2020). The paper gives a full error analysis, and full details of the construction of lattices needed to ensure a good (but inevitably not optimal) rate of convergence and an error bound independent of dimension. Numerical experiments support the theory.

preprint2020arXiv

Bounds on the spectrum of nonsingular triangular $(0,1)$ -matrices

Let $K_n$ be the set of all nonsingular $n\times n$ lower triangular $(0,1)$ -matrices. Hong and Loewy (2004) introduced the numbers $c_n={\rm min}\{λ\mid λ~\text{is an eigenvalue of}~XX^{\rm T},~X\in K_n\},\quad n\in\mathbb{Z}_+.$ A related family of numbers was considered by Ilmonen, Haukkanen, and Merikoski (2008): $C_n={\rm max}\{λ\mid λ~\text{is an eigenvalue of}~XX^{\rm T},~X\in K_n\},\quad n\in\mathbb{Z}_+.$ These numbers can be used to bound the singular values of matrices belonging to $K_n$ and they appear, e.g., in eigenvalue bounds for power GCD matrices, lattice-theoretic meet and join matrices, and related number-theoretic matrices. In this paper, it is shown that for $n$ odd, one has the lower bound $c_n\geq \frac{1}{\sqrt{\frac{1}{25}φ^{-4n}+\frac{2}{25}φ^{-2n}-\frac{2}{5\sqrt{5}}nφ^{-2n}-\frac{23}{25}+n+\frac{2}{25}φ^{2n}+\frac{2}{5\sqrt{5}}nφ^{2n}+\frac{1}{25}φ^{4n}}},$ and for $n$ even, one has $c_n\geq \frac{1}{\sqrt{\frac{1}{25}φ^{-4n}+\frac{4}{25}φ^{-2n}-\frac{2}{5\sqrt{5}}nφ^{-2n}-\frac{2}{5}+n+\frac{4}{25}φ^{2n}+\frac{2}{5\sqrt{5}}nφ^{2n}+\frac{1}{25}φ^{4n}}},$ where $φ$ denotes the golden ratio. These lower bounds improve the estimates derived previously by Mattila (2015) and Altinişik et al. (2016). The sharpness of these lower bounds is assessed numerically and it is conjectured that $c_n\sim 5φ^{-2n}$ as $n\to\infty$ . In addition, a new closed form expression is derived for the numbers $C_n$ , viz. $C_n=\frac14 \csc^2\bigg(\fracπ{4n+2}\bigg)=\frac{4n^2}{π^2}+\frac{4n}{π^2}+\bigg(\frac{1}{12}+\frac{1}{π^2}\bigg)+\mathcal{O}\bigg(\frac{1}{n^2}\bigg),\quad n\in\mathbb{Z}_+.$

preprint2020arXiv

Positive definite functions on semilattices

We introduce a notion of positive definiteness for functions $f\!:P\to\mathbb{R}$ defined on meet semilattices $(P,\preceq,\wedge)$ and prove several properties for these functions. In addition, we utilize the $LDL^{\rm T}$ decomposition of meet matrices in order to explore the properties of multivariate positive definite arithmetic functions $f\!:\mathbb{Z}_+^d\to\mathbb{R}$ . Finally, we give a series of examples and counterexamples of positive definite functions.

preprint2020arXiv

Uncertainty quantification using periodic random variables

Many studies in uncertainty quantification have been carried out under the assumption of an input random field in which a countable number of independent random variables are each uniformly distributed on an interval, with these random variables entering linearly in the input random field (the so-called affine model). In this paper we consider an alternative model of the random field, in which the random variables have the same uniform distribution on an interval, but the random variables enter the input field as periodic functions. The field is constructed in such a way as to have the same mean and covariance function as the affine random field. Higher moments differ from the affine case, but in general the periodic model seems no less desirable. The new model of the random field is used to compute expected values of a quantity of interest arising from an elliptic PDE with random coefficients. The periodicity is shown to yield a higher order cubature convergence rate of $\mathcal{O}(n^{-1/p})$ independently of the dimension when used in conjunction with rank-1 lattice cubature rules constructed using suitably chosen smoothness-driven product and order dependent weights, where $n$ is the number of lattice points and $p$ is the summability exponent of the fluctuations in the series expansion of the random coefficient. We present numerical examples that assess the performance of our method.

preprint2016arXiv

Polynomial collocation for handling an inaccurately known measurement configuration in electrical impedance tomography

The objective of electrical impedance tomography is to reconstruct the internal conductivity of a physical body based on measurements of current and potential at a finite number of electrodes attached to its boundary. Although the conductivity is the quantity of main interest in impedance tomography, a real-world measurement configuration includes other unknown parameters as well: the information on the contact resistances, electrode positions and body shape is almost always incomplete. In this work, the dependence of the electrode measurements on all aforementioned model properties is parametrized via polynomial collocation. The availability of such a parametrization enables efficient simultaneous reconstruction of the conductivity and other unknowns by a Newton-type output least squares algorithm, which is demonstrated by two-dimensional numerical experiments based on both noisy simulated data and experimental data from two water tanks.

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.

Vesa Kaarnioja

Research overview

Is this your work?

5 publications

Fast approximation by periodic kernel-based lattice-point interpolation with application in uncertainty quantification

Bounds on the spectrum of nonsingular triangular $(0,1)$ -matrices

Positive definite functions on semilattices

Uncertainty quantification using periodic random variables

Polynomial collocation for handling an inaccurately known measurement configuration in electrical impedance tomography

See the researcher in context

Building this map preview

Vesa Kaarnioja

Research overview

Is this your work?

5 publications

Fast approximation by periodic kernel-based lattice-point interpolation with application in uncertainty quantification

Bounds on the spectrum of nonsingular triangular (0,1)(0,1)(0,1)-matrices

Positive definite functions on semilattices

Uncertainty quantification using periodic random variables

Polynomial collocation for handling an inaccurately known measurement configuration in electrical impedance tomography

See the researcher in context

Building this map preview

Bounds on the spectrum of nonsingular triangular $(0,1)$ -matrices