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.

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.

Excluded volume effects on active ring polymers

The properties of semiflexible active ring polymers are studied by numerical simulations. The two-dimensionally confined polymer is modeled as a closed bead-spring chain subject to tangential active forces, and the interaction with the fluid is described by the Brownian multiparticle collision dynamics approach. Both phantom polymers and chains with excluded-volume interactions are considered. The size and shape strongly depend on the relative ratio of the persistence length to the ring length as well as on the active force.

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

Screw dislocations in periodic media: Variational coarse graining of the discrete elastic energy

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.

On the probability of (falsely) connecting two distinct components when learning a GGM

In this paper, we extend the result on the probability of (falsely) connecting two distinct components when learning a GGM (Gaussian Graphical Model) by the joint regression based technique. While the classical method of regression based technique learns the neighbours of each node one at a time through a Lasso penalized regression, its joint modification, considered here, learns the neighbours of each node simultaneously through a group Lasso penalized regression.

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.

A renewal equation approach to behavioural epidemic models: analytical and numerical issues

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].