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.

We study the asymptotic behavior, as the lattice spacing ? tends to zero, of the discrete elastic energy induced by topological singularities in an inhomogeneous ? periodic medium within a two-dimensional model for screw dislocations in the square lattice. We focus on the |log?| regime which, as ?->0 allows the emergence of a finite number of limiting topological singularities.

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.

In this manuscript we propose a numerical method for non-linear integro-differential systems arising in age-of-infection models in a heterogeneously mixed population. The discrete scheme is based on direct quadrature methods and provides an unconditionally positive and bounded solution. Furthermore, we prove the existence of the numerical final size of the epidemic and show that it tends to its continuous equivalent as the discretization steplength vanishes.

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

We consider a discrete model of planar elasticity where the particles, in the reference configuration, sit on a regular triangular lattice and interact through nearest-neighbor pairwise potentials, with bonds modeled as linearized elastic springs. Within this framework, we introduce plastic slip fields, whose discrete circulation around each tri-angle detects the possible presence of an edge dislocation.

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.

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

Epidemic models structured by the age of infection can be formulated in terms of a system of renewal equations and represent a very general mathematical framework for the analysis of infectious diseases ([1, 2]). Here, we propose a formulation of renawal equations that takes into account of the behavioral response of individuals to infection. We use the so called "information index", which is a distributed delay that summarizes the information available on current and past disease trend, and extend some results regarding compartmental behavioral models [3, 4, 5].

In this paper we consider the computation of the modified moments for the system of Laguerre polynomials on the real semiaxis with the Hermite weight. These moments can be used for the computation of integrals with the Hermite weight on the real semiaxis via product rules.