Nonlocal finite dffierence discretization of a class of renewal equation models for epidemics

In this paper we consider a non-standard discretization to a Volterra integro-dierential system which includes a number of age-of-infection models in the literature. The aim is to provide a general framework to analyze the proposed scheme for the numerical solution of a class of problems whose continuous dynamic is well known in the literature and allow a deeper analysis in cases where the theory lacks

A non-standard numerical method preserving global properties of integro-differential systems

We propose a numerical method for a general integro-differential system of equations which includes a number of age-of-infection epidemic models in the literature [1, 2]. The numerical solution is obtained by a non-standard discretization of the nonlinear terms in the system, and agrees with the analytical solution in many important qualitative aspects. Both the behaviour at finite time and the asymptotic properties of the solution are preserved for any value of the discretization parameter.

Computing integrals with an exponential weight on the real axis in floating point arithmetic

The aim of this work is to propose a fast and reliable algorithm for computing integrals of the type $$\int_{-\infty}^{\infty} f(x) e^{\scriptstyle -x^2 -\frac{\scriptstyle 1}{\scriptstyle x^2}} dx,$$ where $f(x)$ is a sufficiently smooth function, in floating point arithmetic. The algorithm is based on a product integration rule, whose rate of convergence depends only on the regularity of $f$, since the coefficients of the rule are ``exactly'' computed by means of suitable recurrence relations here derived. We prove stability and convergence in the space of locally continuous functions

A Preliminary Investigation of a Single Shock Impact on Italian Mortality Rates Using STMF Data: A Case Study of COVID-19

Mortality shocks, such as pandemics, threaten the consolidated longevity improvements, confirmed in the last decades for the majority of western countries. Indeed, just before the COVID-19 pandemic, mortality was falling for all ages, with a different behavior according to different ages and countries. It is indubitable that the changes in the population longevity induced by shock events, even transitory ones, affecting demographic projections, have financial implications in public spending as well as in pension plans and life insurance.

A multi-GPU aggregation-based AMG preconditioner for iterative linear solvers

We present and release in open source format a sparse linear solver which efficiently exploits heterogeneous parallel computers. The solver can be easily integrated into scientific applications that need to solve large and sparse linear systems on modern parallel computers made of hybrid nodes hosting Nvidia Graphics Processing Unit (GPU) accelerators.

Rational QZ Steps with perfect shifts

In this paper we analyze the stability of the problem of performing a rational QZ$step with a shift that is an eigenvalue of a given regular pencil H-lambda K in unreduced Hessenberg-Hessenberg form. In exact arithmetic, the backward rational QZ step moves the eigenvalue to the top of the pencil, while the rest of the pencil is maintained in Hessenberg-Hessenberg form, which then yields a deflation of the given shift. But in finite-precision the rational QZ step gets ``blurred'' and precludes the deflation of the given shift at the top of the pencil.

The variational approach to s-fractional heat flows and the limit cases s-+0+and s-+1-

This paper deals with the limit cases for s-fractional heat flows in a cylindrical domain, with homogeneous Dirichlet boundary conditions, as s-+ 0+ and s-+ 1-. We describe the fractional heat flows as minimizing move-ments of the corresponding Gagliardo seminorms, with re-spect to the L2 metric. To this end, we first provide a Gamma-convergence analysis for the s-Gagliardo seminorms as s-+ 0+ and s-+ 1-; then, we exploit an abstract stability result for minimizing movements in Hilbert spaces, with respect to a sequence of Gamma-converging uniformly lambda-convex energy function-als.

Seeking critical nodes in digraphs

The Critical Node Detection Problem (CNDP) consists in finding the set of nodes, defined critical, whose removal maximally degrades the graph. In this work we focus on finding the set of critical nodes whose removal minimizes the pairwise connectivity of a direct graph (digraph). Such problem has been proved to be NP-hard, thus we need efficient heuristics to detect critical nodes in real-world applications. We aim at understanding which is the best heuristic we can apply to identify critical nodes in practice, i.e., taking into account time constrains and real-world networks.

RNASeqGUI: a GUI for analysing RNA-Seq data

Summary: We present RNASeqGUI R package, a graphical user interface (GUI) for the identification of differentially expressed genes across multiple biological conditions. This R package includes some wellk-nown RNA-Seq tools, available at www.bioconductor.org. RNASeqGUI package is not just a collection of some known methods and functions, but it is designed to guide the user during the entire analysis process. RNASeqGUI package is mainly addressed to those users who have little experience with command-line software.