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.

Exit-time approach to epsilon-entropy

An efficient approach to the calculation of the E-entropy is proposed. The method is based on the idea of looking at the information content of a string nf data hv annalyzing the signal only nt thp instants when the fluctuations are larger than a certain threshold is an element of, i.e., by looking at the exit-time statistics. The practical and theoretical advantages of our method with respect to the usual one are shown by the examples of a deterministic map and a self-affine stochastic process.

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.

Inverse velocity statistics in two-dimensional turbulence

We present a numerical study of two-dimensional turbulent flows in the enstropy cascade regime, with different large-scale energy sinks. In particular, we study the statistics of more-than-differentiable velocity fluctuations by means of two sets of statistical estimators, namely inverse statistics and second-order differences. In this way, we are able to probe statistical fluctuations that are not captured by the usual spectral analysis. We show that a new set of exponents associated to more-than-differentiable fluctuations of the velocity field exists.

A candidate multi-epitope vaccine against SARS-CoV-2

In the past two decades, 7 coronaviruses have infected the human population, with two major outbreaks caused by SARS-CoV and MERS-CoV in the year 2002 and 2012, respectively. Currently, the entire world is facing a pandemic of another coronavirus, SARS-CoV-2, with a high fatality rate. The spike glycoprotein of SARS-CoV-2 mediates entry of virus into the host cell and is one of the most important antigenic determinants, making it a potential candidate for a vaccine. In this study, we have computationally designed a multi-epitope vaccine using spike glycoprotein of SARS-CoV-2.

Less Is Enough: Assessment of the Random Sampling Method for the Analysis of Magnetoencephalography (MEG) Data

Magnetoencephalography (MEG) aims at reconstructing the unknown neuroelectric activity in the brain from non-invasive measurements of the magnetic field induced by neural sources. The solution of this ill-posed, ill-conditioned inverse problem is usually dealt with using regularization techniques that are often time-consuming, and computationally and memory storage demanding. In this paper we analyze how a slimmer procedure, random sampling, affects the estimation of the brain activity generated by both synthetic and real sources.

Relative entropy in diffusive relaxation for a class of discrete velocities BGK models

We provide a framework to extend the relative entropy method to a class of diffusive relaxation systems with discrete velocities. The methodology is detailed in the toy case of the 1D Jin-Xin model under the diffusive scaling, and provides a direct proof of convergence to the limit parabolic equation in any interval of time, in the regime where the solutions are smooth. Recently, the same approach has been successfully used to show the strong convergence of a vector-BGK model to the 2D incompressible Navier-Stokes equations.