
Geometric Numerical Integration in Ecological Modelling
A major neglected weakness of many ecological models is the numerical method used to solve the governing systems of differential equations. Indeed, the discrete dynamics described by numerical integrators can provide spurious solution of the corresponding continuous model. The approach represented by the geometric numerical integration, by preserving qualitative properties of the solution, leads to improved numerical behaviour expecially in the long-time integration.
Computing the eigenvectors of nonsymmetric tridiagonal matrices
The computation of the eigenvalue decomposition of matrices is
one of the most investigated problems in numerical linear algebra. In particular,
real nonsymmetric tridiagonal eigenvalue problems arise in a variety of
applications. In this paper the problem of computing an eigenvector corresponding
to a known eigenvalue of a real nonsymmetric tridiagonal matrix is
considered, developing an algorithm that combines part of a QR sweep and
part of a QL sweep, both with the shift equal to the known eigenvalue. The
numerical tests show the reliability of the proposed method.
AMG based on compatible weighted matching on GPUs
GPU version of AMG preconditioner
A fast and Robust spectrogram reassignment method
The improvement of the readability of time-frequency transforms is an important topic in the field of fast-oscillating signal processing. The reassignment method is often used due to its adaptivity to different transforms and nice formal properties. However, it strongly depends on the selection of the analysis window and it requires the computation of the same transform using three different but well-defined windows.
On QZ Steps with Perfect Shifts and Computing the Index of a Differential Algebraic Equation
In this paper we revisit the problem of performing a QZ step with a so-called "perfect shift", which is an
"exact" eigenvalue of a given regular pencil lambda B-A in unreduced Hessenberg-Triangular form. In exact
arithmetic, the QZ step moves that eigenvalue to the bottom of the pencil, while the rest of the pencil is
maintained in Hessenberg-Triangular form, which then yields a deflation of the given eigenvalue. But
in finite-precision the QZ step gets "blurred" and precludes the deflation of the given eigenvalue.





