Luka Grubisic

Luka Grubisic, Matko Ljulj, Volker Mehrmann et al.

A new model description for the numerical simulation of elastic stents is proposed. Based on the new formulation an inf-sup inequality for the finite element discretization is proved and the proof of the inf-sup inequality for the continuous problem is simplified. The new formulation also leads to faster simulation times despite an increased number of variables. The techniques also simplify the analysis and numerical solution of the evolution problem describing the movement of the stent under external forces. The results are illustrated via numerical examples.

Luka Grubisic

We give both lower and upper estimates for eigenvalues of unbounded positive definite operators in an arbitrary Hilbert space. We show scaling robust relative eigenvalue estimates for these operators in analogy to such estimates of current interest in Numerical Linear Algebra. Only simple matrix theoretic tools like Schur complements have been used. As prototypes for the strength of our method we discuss a singularly perturbed Schroedinger operator and study convergence estimates for finite element approximations. The estimates can be viewed as a natural quadratic form version of the celebrated Temple--Kato inequality.

Luka Grubisic, Kresimir Veselic

We use a ``weakly formulated'' Sylvester equation $$A^{1/2}TM^{-1/2}-A^{-1/2}TM^{1/2}=F$$ to obtain new bounds for the rotation of spectral subspaces of a nonnegative selfadjoint operator in a Hilbert space. Our bound extends the known results of Davis and Kahan. Another application is a bound for the square root of a positive selfadjoint operator which extends the known rule: ``The relative error in the square root is bounded by the one half of the relative error in the radicand''. Both bounds are illustrated on differential operators which are defined via quadratic forms.

Luka Grubisic

In this article we further develop a perturbation approach to the Rayleigh--Ritz approximations from our earlier work. We both sharpen the estimates and extend the applicability of the theory to nonnegative definite operators . The perturbation argument enables us to solve two problems in one go: We determine which part of the spectrum of the operator is being approximated by the Ritz values and compute the approximation estimates. We also present a Temple--Kato like inequality which --unlike the original Temple--Kato inequality-- applies to any test vectors from the quadratic form domain of the operator.