Fixation probabilities in weakly compressible fluid flows

Competition between biological species in marine environments is affected by the motion of the surrounding fluid. An effective 2D compressibility can arise, for example, from the convergence and divergence of water masses at the depth at which passively traveling photosynthetic organisms are restricted to live. In this report, we seek to quantitatively study genetics under flow. To this end, we couple an off-lattice agent-based simulation of two populations in 1D to a weakly compressible velocity field--first a sine wave and then a shell model of turbulence.

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.

A free boundary model for the evolution of a geothermal system

The evolution of a geothermal system is studied. A mathematical model is proposed and the corresponding free boundary problem is formulated in a one-dimensional geometry. A situation corresponding to the geothermal field in Larderello, Tuscany (Italy) is considered, showing that the problem has two characteristic time scales, related to the motion of interface and diffusion of vapor.

An inversion method based on random sampling for real-time MEG neuroimaging

The MagnetoEncephaloGraphy (MEG) has gained great interest in neurorehabilitation training due to its high temporal resolution. The challenge is to localize the active regions of the brain in a fast and accurate way. In this paper we use an inversion method based on random spatial sampling to solve the real-time MEG inverse problem. Several numerical tests on synthetic but realistic data show that the method takes just a few hundredths of a second on a laptop to produce an accurate map of the electric activity inside the brain. Moreover, it requires very little memory storage.

Uniform weighted approximation by multivariate filtered polynomials

The paper concerns the weighted uniform approximation of a real function on the d-cube [-1, 1]^d, with d > 1, by means of some multivariate filtered polynomials. These polynomials have been deduced, via tensor product, from certain de la Vallée Poussin type means on [-1, 1], which generalize classical delayed arithmetic means of Fourier partial sums. They are based on arbitrary sequences of filter coefficients, not necessarily connected with a smooth filter function.

A moving-grid approach for fluid-structure interaction problems with hybrid lattice Boltzmann method

In this paper, we propose a hybrid lattice Boltzmann method (HLBM) for solving fluid-structure interaction problems. The proposed numerical approach is applied to model the flow induced by a vibrating thin lamina submerged in a viscous quiescent fluid. The hydrodynamic force exerted by the fluid on the solid body is described by means of a complex hydrodynamic function, whose real and imaginary parts are determined via parametric analysis. Numerical results are validated by comparison with those from other numerical as well as experimental works available in the literature.

Noise Removal from Remote Sensed Images by NonLocal Means with OpenCL Algorithm

We introduce a multi-platform portable implementation of the NonLocal Means methodology aimed at noise removal from remotely sensed images. It is particularly suited for hyperspectral sensors for which real-time applications are not possible with only CPU based algorithms. In the last decades computational devices have usually been a compound of cross-vendor sets of specifications (heterogeneous system architecture) that bring together integrated central processing (CPUs) and graphics processor (GPUs) units.