Researcher profile

Davide Bianchi

Davide Bianchi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2025arXiv

Galaxy-Multiplet Clustering from DESI DR2

We present an efficient estimator for higher-order galaxy clustering using small groups of nearby galaxies, or multiplets. Using the Luminous Red Galaxy (LRG) sample from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2, we identify galaxy multiplets as discrete objects and measure their cross-correlations with the general galaxy field. Our results show that the multiplets exhibit stronger clustering bias as they trace more massive dark matter halos than individual galaxies. When comparing the observed clustering statistics with the mock catalogs generated from the N-body simulation AbacusSummit, we find that the mocks underpredict multiplet clustering despite reproducing the galaxy two-point auto-correlation reasonably well. This discrepancy indicates that the standard Halo Occupation Distribution (HOD) model is insufficient to describe the properties of galaxy multiplets, revealing the greater constraining power of this higher-order statistic on galaxy-halo connection and the possibility that multiplets are specific to additional assembly bias. We demonstrate that incorporating secondary biases into the HOD model improves agreement with the observed multiplet statistics, specifically by allowing galaxies to preferentially occupy halos in denser environments. Our results highlight the potential of utilizing multiplet clustering, beyond traditional two-point correlation measurements, to break degeneracies in models describing the galaxy-dark matter connection.

preprint2022arXiv

The generalized porous medium equation on graphs: existence and uniqueness of solutions with $\ell^1$ data

We study solutions of the generalized porous medium equation on infinite graphs. For nonnegative or nonpositive integrable data, we prove the existence and uniqueness of mild solutions on any graph. For changing sign integrable data, we show existence and uniqueness under extra assumptions such as local finiteness or a uniform lower bound on the node measure.

preprint2021arXiv

Compatibility, embedding and regularization of non-local random walks on graphs

Several variants of the graph Laplacian have been introduced to model non-local diffusion processes, which allow a random walker to {\textquotedblleft jump\textquotedblright} to non-neighborhood nodes, most notably the transformed path graph Laplacians and the fractional graph Laplacian. From a rigorous point of view, this new dynamics is made possible by having replaced the original graph $G$ with a weighted complete graph $G'$ on the same node-set, that depends on $G$ and wherein the presence of new edges allows a direct passage between nodes that were not neighbors in $G$. We show that, in general, the graph $G'$ is not compatible with the dynamics characterizing the original model graph $G$: the random walks on $G'$ subjected to move on the edges of $G$ are not stochastically equivalent, in the wide sense, to the random walks on $G$. From a purely analytical point of view, the incompatibility of $G'$ with $G$ means that the normalized graph $\hat{G}$ can not be embedded into the normalized graph $\hat{G}'$. Eventually, we provide a regularization method to guarantee such compatibility and preserving at the same time all the nice properties granted by $G'$.

preprint2021arXiv

Graph Laplacian for image deblurring

Image deblurring is relevant in many fields of science and engineering. To solve this problem, many different approaches have been proposed and among the various methods, variational ones are extremely popular. These approaches are characterized by substituting the original problem with a minimization one where the functional is composed of two terms, a data fidelity term and a regularization term. In this paper we propose, in the classical $\ell^2-\ell^1$ minimization with the non-negativity constraint of the solution, the use of the graph Laplacian as regularization operator. Firstly, we describe how to construct the graph Laplacian from the observed noisy and blurred image. Once the graph Laplacian has been built, we solve efficiently the proposed minimization problem splitting the convolution operator and the graph Laplacian by the alternating direction method of multipliers (ADMM). Some selected numerical examples show the good performances of the proposed algorithm.

preprint2020arXiv

Analysis of the spectral symbol associated to discretization schemes of linear self-adjoint differential operators

Given a linear self-adjoint differential operator $\mathcal{L}$ along with a discretization scheme (like Finite Differences, Finite Elements, Galerkin Isogeometric Analysis, etc.), in many numerical applications it is crucial to understand how good the (relative) approximation of the whole spectrum of the discretized operator $\mathcal{L}^{(\mathbf{n})}$ is, compared to the spectrum of the continuous operator $\mathcal{L}$. The theory of Generalized Locally Toeplitz sequences allows to compute the spectral symbol function $ω$ associated to the discrete matrix $\mathcal{L}^{(\mathbf{n})}$. We prove that the symbol $ω$ can measure, asymptotically, the maximum spectral relative error $\mathcal{E}\geq 0$. It measures how the scheme is far from a good relative approximation of the whole spectrum of $\mathcal{L}$, and it suggests a suitable (possibly non-uniform) grid such that, if coupled to an increasing refinement of the order of accuracy of the scheme, guarantees $\mathcal{E}=0$.

preprint2020arXiv

Asymptotic spectra of large (grid) graphs with a uniform local structure

We are concerned with sequences of graphs having a grid geometry, with a uniform local structure in a bounded domain $Ω\subset {\mathbb R}^d$, $d\ge 1$. We assume $Ω$ to be Lebesgue measurable with regular boundary and contained, for convenience, in the cube $[0,1]^d$. When $Ω=[0,1]$, such graphs include the standard Toeplitz graphs and, for $Ω=[0,1]^d$, the considered class includes $d$-level Toeplitz graphs. In the general case, the underlying sequence of adjacency matrices has a canonical eigenvalue distribution, in the Weyl sense, and we show that we can associate to it a symbol $f$. The knowledge of the symbol and of its basic analytical features provide many informations on the eigenvalue structure, of localization, spectral gap, clustering, and distribution type. Few generalizations are also considered in connection with the notion of generalized locally Toeplitz sequences and applications are discussed, stemming e.g. from the approximation of differential operators via numerical schemes.

preprint2020arXiv

Generalized Structure Preserving Preconditioners for Frame-Based Image Deblurring

We are interested in fast and stable iterative regularization methods for image deblurring problems with space invariant blur. The associated coefficient matrix has a Block Toeplitz Toeplitz Blocks (BTTB) like structure plus a small rank correction depending on the boundary conditions imposed on the imaging model. In the literature, several strategies have been proposed in the attempt to define proper preconditioner for iterative regularization methods that involve such linear systems. Usually, the preconditioner is chosen to be a Block Circulant with Circulant Blocks (BCCB) matrix because it can be efficiently exploit Fast Fourier Transform (FFT) for any computation, including the (pseudo-)inversion. Nevertheless, for ill-conditioned problems, it is well known that BCCB preconditioners cannot provide a strong clustering of the eigenvalues. Moreover, in order to get an effective preconditioner, it is crucial to preserve the structure of the coefficient matrix. On the other hand, thresholding iterative methods have been recently successfully applied to image deblurring problems, exploiting the sparsity of the image in a proper wavelet domain. Motivated by the results of recent papers, we combine a nonstationary preconditioned iteration with the modified linearized Bregman algorithm (MLBA) and proper regularization operators. Several numerical experiments shows the performances of our methods in terms of quality of the restorations.

preprint2020arXiv

Qualitative properties of bounded subsolutions of nonlinear PDEs

We study decay and compact support properties of positive and bounded solutions of $Δ_{p} u \geq Λ(u)$ on the exterior of a compact set of a complete manifold with rotationally symmetry. In the same setting, we also give a new characterization of stochastic completeness for the $p$-Laplacian in terms of a global $W^{1,p}$-regularity of such solutions. One of the tools we use is a nonlinear version of the Feller property which we investigate on general Riemannian manifolds and which we establish under integral Ricci curvature conditions.

preprint2019arXiv

Confronting missing observations with probability weights: Fourier space and generalised formalism

Due to instrumental limitations, the nature of which vary from case to case, spectroscopic galaxy redshift surveys usually do not collect redshifts for all galaxies in the population of potential targets. Especially problematic is the entanglement between this incompleteness and the true cosmological signal, arising from the fact that the proportion of successful observations is typically lower in regions with higher density of galaxies. The result is a fictitious suppression of the galaxy clustering that, if not properly corrected, can impact severely on cosmological-parameter inference. Recent developments in the field have shown that an unbiased estimate of the 2-point correlation function in the presence of missing observations can be obtained by weighting each pair by its inverse probability of being targeted. In this work we expand on the concept of probability weights by developing a more mature statistical formalism, which provides us with a deeper understanding of their fundamental properties. We take advantage of this novel perspective to handle the problem of estimating the inverse probability, specifically, we discuss how to efficiently determine the weights from a finite set of realisations of the targeting and how to model exactly the resulting sampling effects. This allows us to derive an inverse-probability-based estimator for the power spectrum, which is the main result of this work, but also to improve robustness and computational efficiency of the already existing configuration-space estimator. Finally, we propose a strategy to further extend the concept of inverse probability, providing examples of how traditional missing-observation countermeasures can be included in this more general picture. The effectiveness of the different models and weighting schemes discussed in this work is demonstrated using realisations of an idealised simple survey strategy.

preprint2016arXiv

Laplacian cut-offs, porous and fast diffusion on manifolds and other applications

We construct exhaustion and cut-off functions with controlled gradient and Laplacian on manifolds with Ricci curvature bounded from below by a (possibly unbounded) nonpositive function of the distance from a fixed reference point, without any assumptions on the topology or the injectivity radius. Along the way we prove a generalization of the Li-Yau gradient estimate which is of independent interest. We then apply our cut-offs to the study of the fast and porous media diffusion, of $L^q$-properties of the gradient and of the self-adjointness of Schroedinger-type operators.

preprint2015arXiv

Regularization preconditioners for frame-based image deblurring with reduced boundary artifacts

Thresholding iterative methods are recently successfully applied to image deblurring problems. In this paper, we investigate the modified linearized Bregman algorithm (MLBA) used in image deblurring problems, with a proper treatment of the boundary artifacts. We consider two standard approaches: the imposition of boundary conditions and the use of the rectangular blurring matrix. The fast convergence of the MLBA depends on a regularizing preconditioner that could be computationally expensive and hence it is usually chosen as a block circulant circulant block (BCCB) matrix, diagonalized by discrete Fourier transform. We show that the standard approach based on the BCCB preconditioner may provide low quality restored images and we propose different preconditioning strategies, that improve the quality of the restoration and save some computational cost at the same time. Motivated by a recent nonstationary preconditioned iteration, we propose a new algorithm that combines such method with the MLBA.We prove that it is a regularizing and convergent method. A variant with a stationary preconditioner is also considered. Finally, a large number of numerical experiments shows that our methods provide accurate and fast restorations, when compared with the state of the art.

preprint2014arXiv

Iterated fractional Tikhonov regularization

Fractional Tikhonov regularization methods have been recently proposed to reduce the oversmoothing property of the Tikhonov regularization in standard form, in order to preserve the details of the approximated solution. Their regularization and convergence properties have been previously investigated showing that they are of optimal order. This paper provides saturation and converse results on their convergence rates. Using the same iterative refinement strategy of iterated Tikhonov regularization, new iterated fractional Tikhonov regularization methods are introduced. We show that these iterated methods are of optimal order and overcome the previous saturation results. Furthermore, nonstationary iterated fractional Tikhonov regularization methods are investigated, establishing their convergence rate under general conditions on the iteration parameters. Numerical results confirm the effectiveness of the proposed regularization iterations.