Effectively dealing with invasive species is a pervasive problem in environmental management. The damages, and associated costs, that stem from invasive species are well known, as is the benefit from their removal. We investigate problems where optimal control theory has been implemented, and we show that these problems can easily become hypersensitive, making their numerical solutions unstable. We show that transforming these problems from state-adjoint systems to state-control systems can provide useful insights into the system dynamics and simplify the numerics. We apply these techniques to two case studies: one of feral cats in Australia, where we use logistic growth; and the other of wild-boars in Italy, where we include an Allee effect. A further development is to optimize the control strategy by taking into account the spatio-temporal features of the invasive species control problems over large and irregular environments. The approach is used in a management scenario where the invasive species to be controlled with an optimal allocation of resources is the deciduous tree Ailanthus Altissima, infesting the Alta Murgia National Park in the south of Italy. This work has been carried out within the H2020 project ECOPOTENTIAL (http://www.ecopotential-project.eu), coordinated by CNR-IGG. The project has received funding from the European Union's Horizon 2020 research and innovation programme (grant agreement No 641762).

Given a population of N elements with their geographical positions and the genetic (or lexical) distances between couples of elements (inferred, for example, from lexical differences between dialects which are spoken in different towns or from genetic differences between animal populations living in different faunal areas) a very interesting problem is to reconstruct the geographical positions of individuals using only genetic/lexical distances. From a technical point of view the program consists in extracting from the genetic/lexical distances a set of reconstructed geographical positions to be compared with the real ones. We show that geographical recovering is successful when the genetic/lexical distances are not a simple consequence of phylogenesis but also of horizontal transfers as, for example, vocabulary borrowings between different languages. Our results go well beyond the simple observation that geographical distances and genetic/lexical distances are correlated. The ascertainment of a correlation, in our perspective, merely is a prerequisite.

We study spreading dynamics of a reaction diffusion process in a special class of heterogeneous graphs with Poissonian degree distribution and composed of both local and long range links. The behavior of the spreading dynamics on such networks are investigated by relating them to the topological features of graphs. We find that the degree of assortativity can give just some indication about the large scale behavior of the spreading dynamics while a detailed description of the process can be addressed by introducing new, more appropriate, topological quantities linked to the distance between nodes. (C) 2016 Published by Elsevier B.V.

The paper presents a method for color quantization (CQ) which uses visual contrast for determining an image-dependent color palette. The proposed method selects image regions in a hierarchical way, according to the visual importance of their colors with respect to the whole image. The method is automatic, image dependent and requires a moderate computational effort. Preliminary results show that the quality of quantized images, measured in terms of Mean Square Error, Color Loss and SSIM, is competitive with some existing CQ approaches.

This paper focuses on an entropy based formalism to speed up the evaluation of the Structural SIMilarity (SSIM) index in images affected by a global distortion. Looking at images as information sources, a visual distortion typical set can be defined for SSIM. This typical set consists of just a subset of information belonging to the original image and the corresponding one in the distorted version. As side effect, some general theoretical criteria for the computation of any full reference quality assessment measure can be given in order to maximize its computational efficiency. Experimental results on various test images show that the proposed approach allows to estimate SSIM with a considerable speed up (about 200 times) and a small relative error (often lower than 5%).

This paper aims at increasing the visual quality of a blurred image according to the contrast sensitivity of a human observer. The main idea is to enhance those image details which can be perceived by a human observer without introducing annoying visible artifacts. To this aim, an adaptive wavelet decomposition is applied to the original blurry image. This decomposition splits the frequency axis into subbands whose central frequency and amplitude width are built according to the contrast sensitivity. The details coefficients of that decomposition are then properly modified according to the just noticeable contrast at each frequency band. Preliminary experimental results show that the proposed method increases the visual quality of the blurred image without introducing visible artifacts. In addition, the contrast sensitivity-based image is a good and recommended initial guess for iterative deblurring methods since it allows them to significantly reduce ringing artifacts and halo effects in the final image.

This paper proposed improving the solve time of the bootstrap AMG proposed previously by the authors. This is achieved by incorporating the information, set of algebraically smooth vectors, generated by the bootstrap algorithm, in a single hierarchy by using sufficiently large aggregates, and these aggregates are compositions of aggregates already built throughout the bootstrap algorithm. The modified AMG method has comparable convergence properties to the original bootstrap one, however with better efficiency. The improvement in solve time with respect to the original bootstrap AMG is illustrated on some difficult linear systems arising from discretization of vector function elliptic Partial Differential Equations (PDEs) in both 2d and 3d.

We study a variational problem for simultaneous video inpainting and motion estimation. We consider a functional proposed by Lauze and Nielsen [25] and we study, by means of the relaxation method of the Calculus of Variations, a slightly modified version of this functional. The domain of the relaxed functional is constituted of functions of bounded variation and we compute a representation formula of the relaxed functional. The representation formula shows the role of discontinuities of the various functions involved in the variational model. The present study clarifies the variational properties of the functional proposed in [25] for motion compensated video inpainting.

Software package in C

Midterm Working Report for EoCoE Project

Release 2.1 of MLD2P4 package

The Optimized Retrieval Code (ORM) was originally designed to be the scientific code used as the reference for the ESA operational Near Real Time Analysis of MIPAS measurements on ENVISAT. After the end of ENVISAT mission, occurred on April 2017, this code has been adapted and improved to perform the reanalysis of the full MIPAS mission.

Position determination of photon emitters and associated strong field parallax effects are investigated using relativistic optics when the photon orbits are confined to the equatorial plane of the Schwarzschild spacetime. We assume the emitter is at a fixed space position and the receiver moves along a circular geodesic orbit. This study requires solving the inverse problem of determining the (spatial) intersection point of two null geodesic initial data problems, serving as a simplified model for applications in relativistic astrometry as well as in radar and satellite communications.

We investigate the scattering of a spinning test particle by a Kerr black hole within the Mathisson-Papapetrou-Dixon model to linear order in spin. The particle's spin and orbital angular momentum are taken to be aligned with the black hole's spin. Both the particle's mass and spin length are assumed to be small in comparison with the characteristic length scale of the background curvature, in order to avoid backreaction effects. We analytically compute the modifications due to the particle's spin to the scattering angle, the periastron shift, and the condition for capture by the black hole, extending previous results valid for the nonrotating Schwarzschild background. Finally, we discuss how to generalize the present analysis beyond the linear approximation in spin, including spin-squared corrections in the case of a black-hole-like quadrupolar structure for the extended test body.

Consider a dynamic general relativistic spacetime in which the proper infinitesimal interval along one spatial coordinate direction decreases monotonically with time, while the corresponding intervals increase along other spatial directions. In a system undergoing such complete anisotropic collapse/expansion, we look for the formation of a cosmic double-jet configuration: free test particles in the ambient medium, relative to the collapsing system, gain energy from the gravitational field and asymptotically line up parallel and antiparallel to the direction of collapse such that their Lorentz factors approach infinity. A strong burst of electromagnetic radiation is expected to accompany this event if some of the free test particles carry electric charge. Previous work in this direction involved mainly Ricci-flat spacetimes; hence, we concentrate here on inhomogeneous perfect fluid spacetimes. We briefly explore the possible connection between these theoretical cosmic jets and astrophysical jets. We also discuss other general relativistic scenarios for the formation of cosmic jets.

We compute the (center-of-mass frame) scattering angle ? of hyperboliclike encounters of two spinning black holes, at the fourth post-Newtonian approximation level for orbital effects, and at the next-to-next-to-leading order for spin-dependent effects. We find it convenient to compute the gauge-invariant scattering angle (expressed as a function of energy, orbital angular momentum and spins) by using the effective-one-body formalism. The contribution to scattering associated with nonlocal, tail effects is computed by generalizing to the case of unbound motions the method of time localization of the action introduced in the case of (small-eccentricity) bound motions by Damour et al. [Phys. Rev. D 91, 084024 (2015)PRVDAQ1550-799810.1103/PhysRevD.91.084024].

In this work we present an analytical gravitational self-force calculation of the spin-orbit precession along an eccentric orbit around a Schwarzschild black hole, following closely the recent prescription of Akcay, Dempsey, and Dolan, giving results to six post-Newtonian orders expanded in small eccentricity through e2. We then transcribe this quantity within the effective-one-body (EOB) formalism, thereby determining several new, linear-in-mass-ratio contributions in the post-Newtonian expansion of the spin-orbit couplings entering the EOB Hamiltonian. Namely, we determine the second gyrogravitomagnetic ratio gS*(r,pr,p?) up to order pr2/r4 included.

The precession angular velocity of a gyroscope moving along a general geodesic in the Kerr spacetime is analyzed using the geometric properties of the spacetime. Natural frames along the gyroscope world line are explicitly constructed by boosting frames adapted to fundamental observers. A novel geometrical description is given to Marck's construction of a parallel propagated orthonormal frame along a general geodesic, identifying and clarifying the special role played by the Carter family of observers in this general context, thus extending previous discussion for the equatorial plane case.