Research
My research activities are divided into different strands and projects.
The publications are reported by subject area and project in descending
temporal order.
The same information can be found on either my
Google Scholar Profile, my
ResearchGate, or my
Orcid profile.
Solving Linear Algebra problems is a core task in pretty much every Scientific
challenge, from the numerical solution of partial differential equations to the
new problems arising in the analysis of larger and larger sets of data. Moreover,
we are fastly moving towards the use of exascaleenabled HPC platforms, thus we
need to be capable of pushing our algorithms (and codes) to work efficiently in
this environment. This first research line deals with this type of issue, and
more specifically with the task of solving linear systems of equations,
$$ A \mathbf{x} = \mathbf{b}, \qquad A \in \mathbb{R}^{N \times N}, \; \mathbf{x},\mathbf{b} \in \mathbb{R}^{N}, $$
in which $A$ is a
large, $N \geq 10^9$, and
sparse,
number on nonzero entries proportional to $N$, matrix. The key instruments are
here
Krylov Subspace Methods,
Geometric and
Algebraic Multigrid
preconditioners, and the application of update strategies to deal with more than
one linear system.
Preconditioning for large and sparse linear systems 

Owen, H.; Lehmkuhl, O.;
D'Ambra, P.; Durastante, F.;
Filippone, S.; Alya towards Exascale: Algorithmic Scalability
using PSCToolkit. J. Supercomput (2024).
Open access

Bertaccini, D.; D'Ambra, P.;
Durastante, F.; Filippone, S.;
Why diffusionbased preconditioning of Richards equation works: spectral
analysis and computational experiments at very large scale. (2023).
Numer. Linear Algebra Appl.; e2523.
Open access

D'Ambra, P.; Durastante, F.;
Filippone, S.; Zikatanov, L..
Automatic coarsening in Algebraic Multigrid utilizing quality measures
for matchingbased aggregations. Comput. Math. Appl. 144 (2023),
290305.
MR4615364
Open access

D'Ambra, P.; Durastante, F.;
Ferdous, S.M.; Filippone, S.;
Halappanavar, M.; Pothen A..
AMG Preconditioners based on Parallel Hybrid Coarsening and Multiobjective
Graph Matching. (2023). 31st Euromicro International Conference on Parallel,
Distributed and NetworkBased Processing (PDP), Naples, Italy, 2023, pp. 5967.

P. D'Ambra; F. Durastante;
S. Filippone; Parallel Sparse Computation Toolkit.
Software Impacts 15 (2023), no. 100463.
Open Access.

D'Ambra, P.; Durastante, F.;
Filippone, S..
AMG Preconditioners for Linear Solvers towards Extreme Scale.
SIAM J. Sci. Comput. 43 (2021),
no. 5, S679S703.
MR4331965

Bertaccini, D.; Durastante, F..
Iterative methods and preconditioning for large and sparse linear systems
with applications. Monographs and Research Notes in Mathematics.
CRC Press, Boca Raton, FL, 2018. xviii+353 pp.
ISBN: 9781498764162
MR3793630

Bertaccini, D.; Donatelli, M.;
Durastante, F.; SerraCapizzano, S..
Optimizing a multigrid RungeKutta smoother for variablecoefficient
convectiondiffusion equations. Linear Algebra Appl. 533 (2017),
507535. MR3695922

Bertaccini, D.; Durastante, F..
Interpolating preconditioners for the solution of sequence of linear systems.
Comput. Math. Appl. 72 (2016), no. 4, 11181130.
MR3529065


Fractional calculus is that subset of Analysis that studies the possibilities of
defining differential operators of real or complex order. After the definition
of such objects, the need to numerically approximate the solution of differential
equations defined in terms of them immediately arises. The following papers deal
with the efficient numerical solution of different formulations of such problems.

Fractional Differential Equations: Modeling, Discretization, and Numerical Solvers.
Springer INdAM Series (SINDAMS, volume 50). Editors: A. Cardone, M. Donatelli,
F. Durastante, R. Garrappa, M. Mazza,
M. Popolizio. (2023). Pages XII, 146. ISBN: 9789811977152. ISSN: 2281518X.

Aceto, L.; Durastante, F.. Efficient
computation of the Wright function and its applications to fractional
diffusionwave equations. ESAIM Math. Model. Numer. Anal.
56 (2022), no. 6, 21812196.
MR4516169

Durastante, F.. Efficient solution of timefractional
differential equations with a new adaptive multiterm discretization of
the generalized CaputoDzherbashyan derivative.
Calcolo 56 (2019), no. 4, Paper No. 36, 24 pp.
MR4015149

Bertaccini, D.; Durastante, F..
Efficient preconditioner updates for semilinear spacetime fractional
reactiondiffusion equations. Structured matrices in numerical linear algebra,
285302, Springer INdAM Ser., 30, Springer, Cham, 2019.
MR3931580

Bertaccini, D.; Durastante, F.
Block structured preconditioners in tensor form for the allatonce
solution of a finite volume fractional diffusion equation.
Appl. Math. Lett. 95 (2019), 9297.
MR3936774

Bertaccini, D.; Durastante, F..
Limited memory block preconditioners for fast solution of fractional partial
differential equations. J. Sci. Comput. 77
(2018), no. 2, 950970.
MR3860196

Bertaccini, D.; Durastante, F..
Solving mixed classical and fractional partial differential equations
using shortmemory principle and approximate inverses.
Numer. Algorithms 74 (2017), no. 4, 10611082.
MR3626328

Numerical treatment of Fractional Differential Equations 
Numerical optimal control theory is a part of numerical optimization that deals
with the problem of finding a
control for either a boundary value problem
or a dynamical system over a period of time in such a way that an objective function
is minimized/maximized. The prototypical problem of this type has the form
$$ (y^*,u^*) = \arg\min_{u \in \mathcal{U}} J(y,u) = \frac{1}{2}\y  y_z\_2^2 + \frac{\lambda}{2} \ u \_2^2, \text{ s.t. } e(y,u) = 0,$$
where $\mathcal{U}$ is a function space containing the admissible controls, $y_z$
a target state, and $e(y,u)=0$ a boundary value problem that the quantities $y$
and $u$ have to satisfy. The numerical solution of such a problem is challenging
in many ways, it mixes together optimization, differential, and linear algebra
problems in a nontrivial way.
Numerical optimal control of Partial Differential Equations 

Durastante, F.; Furci, I..
Spectral analysis of saddlepoint matrices from optimization problems
with elliptic PDE constraints. Electron. J. Linear Algebra 36
(2020), 773798.
MR4188673
Open Access.

Durastante, F.; Cipolla, S..
Fractional PDE constrained optimization: box and sparse constrained problems.
Numerical methods for optimal control problems,
111135, Springer INdAM Ser., 29, Springer, Cham, 2018.
MR3889795

Cipolla, S.; Durastante, F..
Fractional PDE constrained optimization: an optimizethendiscretize
approach with LBFGS and approximate inverse preconditioning.
Appl. Numer. Math. 123 (2018), 4357.
MR3711990

A matrix function is simply a function which maps a matrix to another matrix,
indeed to maintain interesting properties one has to extend the concept of
scalar function to a matrix in a clever way. If, e.g., $f$ is a function
admitting a Taylor expansion with radius of convergence $r$, and $A$ is a
matrix such $\A\ < r$ for any compatible norm $\\cdot\$, then $f(A)$ can
be defined as
$$ f(A) = \sum_{j=0}^{+\infty} c_k A^k, $$
for $c_k$ the Taylor coefficients of $f$. Equivalently, if $A$ is diagonalizable,
and $f$ is defined on the eigenvalues of $A$, then
$$ f(A) = X f(\Lambda) X^{1}, \qquad A = X \Lambda X^{1}, \; f(\Lambda)_{i,i} = f(\lambda_i(A)), $$
and in several other ways for more general $A$ and $f$. In any case, one of the
principal numerical task that one can face is the computation of MatrixFunctionvector
products, i.e., $\mathbf{y} = f(A)\mathbf{x}$, for $A$ a
large and sparse
matrix. The following line of research deals with these kind of problems by
different techniques, and for different purposes (solution of Fractional
differential equations, network analysis,
etc.).

Aceto, L.; Durastante, F..
Efficient computation of the sinc matrix function for the integration
of secondorder differential equations. Submitted.

Aceto, L.; Durastante, F..
Theoretical error estimates for computing the matrix logarithm by
Padétype approximants. Submitted.

Durastante, F.; Meini, B..
Stochastic $p$th root approximation of a stochastic matrix: A Riemannian
optimization approach. SIAM J. Matrix Anal. Appl. Accepted.

Bertaccini, D.; Durastante, F..
Computing function of large matrices by a preconditioned rational Krylov
method.
Numerical mathematics and advanced applications ENUMATH 2019,
343351, Lect. Notes Comput. Sci. Eng., 139, Springer, Cham,
[2021], 2021.
MR4266513

Bertaccini, D.; Durastante, F.
Computing functions of very large matrices with small TT/QTT ranks
by quadrature formulas. J. Comput. Appl. Math. 370 (2020),
112663, 15 pp.
MR4046619

Aceto, L.; Bertaccini, D.;
Durastante, F.; Novati, P.
Rational Krylov methods for functions of matrices with applications
to fractional partial differential equations.
J. Comput. Phys. 396 (2019), 470482.
MR3989621

Bertaccini, D.; Popolizio, M.;
Durastante, F. Efficient approximation of functions
of some large matrices by partial fraction expansions.
Int. J. Comput. Math. 96 (2019), no. 9,
17991817.
MR3960343

MatrixFunctions 
Network science is the study of complex networks that are nothing more than
graphs that are both not random and not structured, e.g., one usually does not
consider Cayley graphs encoding the structure of an abstract group.
Examples are telecommunication networks, computer networks, biological
networks, social networks, and many others. In each case, every actor is
represented by a vertex of the graph and their connections as an edge.
My interest in this field lies both in the modeling and the computational
aspects, specifically in the task of extracting information on the networked
phenomena by looking at the topology of the graph.
Graphs and Network Science: ComplexNetworks 

Bini, D. A.; Durastante, F.;
Kim, S.; Meini, B.
On Kemeny's constant and stochastic complement. Submitted.

Bertaccini, D.; Durastante, F.
Nonlocal diffusion of variable order on complex networks.
Int. J. Comput. Math. Comput. Syst. Theory 7
(2022), no. 3, 172191.
MR4483942

Bianchi, D.; Donatelli, M.;
Durastante, F.; Mazza, M..
Compatibility, embedding and regularization of nonlocal random
walks on graphs. J. Math. Anal. Appl. 511
(2022), no. 1, Paper No. 126020, 30 pp.
MR4379318

Arrigo, F.; Durastante, F..
Mittag–Leffler Functions and their Applications in Network Science.
SIAM J. Matrix Anal. Appl. 42 (2021),
no. 4, 15811601. MR4340667

Cipolla, S.; Durastante, F.;
Tudisco, F.. Nonlocal pagerank.
ESAIM Math. Model. Numer. Anal. 55 (2021),
no. 1, 7797.
MR4216832

Benzi, M.; Bertaccini, D.;
Durastante, F.; Simunec, I..
Nonlocal network dynamics via fractional graph Laplacians.
J. Complex Netw. 8 (2020), no. 3, cnaa017,
29 pp. MR4130854

Inverse problems are involved in the task of inferring the parameters of the
model from a set of (possibly noisy) measurements requires. These types of
problems can be difficult to solve for several reasons, e.g., one could be
dealing with a lack of unicity, that is: different values of parameters are
consistent with the measurement, and discovering the parameters may be
challenging from the computational point of view. I have done some work in using
linear algebra techniques for a particular type of inverse problem that is
image denoise and deblurring, that is the process of removing blurring
artifacts and noise from images (e.g., somebody shook the camera while we were
taking the photo).

Cipolla, S.; Donatelli, M.;
Durastante, F.. Regularization of Inverse
Problems by an Approximate Matrix–Function Technique.
Numer. Algorithms 88 (2021), 12751308.
DOI. 10.1007/s1107502101076y.

Cipolla, S.; Di Fiore, C.;
Durastante, F.; Zellini, P..
Regularizing properties of a class of matrices including the optimal and
the superoptimal preconditioners. Numer. Linear Algebra Appl.
26 (2019), no. 2, e2225, 17 pp.
MR3911025

Regularization and Inverse Problems

Mathematical methods play also a role in unexpected settings like archaeology.
They enable the analysis and interpretation of data, making predictions, and
gaining insights into quantitatives aspects. Statistical techniques are commonly
employed to analyze large datasets, such as excavation records, artifact distributions,
and survey data. These methods help archaeologists identify patterns, trends,
and correlations, allowing them to make inferences about past human behaviors
and societal structures. Mathematical modeling is also utilized to reconstruct
ancient landscapes, simulate settlement patterns, and explore the dynamics
of cultural change over time. The integration of mathematics into archaeological
investigations enhances the accuracy, precision, and efficiency of research,
contributing to a quantitative understanding of some aspects of history.

Mauro, C.; Durastante, F..
Nocturnal Seafaring: the Reduction of Visibility at Night and its Impact on Ancient Mediterranean Seafaring.
A Study Based on 84th Centuries BC Evidence Journal of Maritime Archaeology, (2024), 119 pp.
Open access

Mauro, C.; Durastante, F..
Evaluating visibility at sea: Instrumental data and historical nautical
records. Mount Etna from the Calabrian Ionian coast (Italy).
The Journal of Island and Coastal Archaeology, 17.1 (2022),
2142 pp.

Applications, intersections, contaminations 
Projects

Fundings for € 4000.00 for the "Progetti Speciali per la Didattica
Università di Pisa 2023/24". Project: HighPerformance Mathematics.

Funding for 4950.00 € for the 2024 GNCS INdAM project "Metodi di riduzione di
modello ed approssimazioni di rango basso per problemi altodimensionali".
P.I. C. Pagliantini. Participants: L. Aceto, F. Durastante, B, Meini, F.G.
Poloni, L, Robol, M, Strazzullo, A. Bucci, A. A. Casulli, L. Saluzzi, I. Simunec.

Collaborator as research associate to the Istituto per le Applicazioni del
Calcolo "Mauro Picone" of the National Research Council of Italy to the
20242026 Energy Oriented Center of Excellence (EoCoE III): Fostering the European
Energy Transition with Exascale A Center of Excellence in Computing Applications.
EuroHPC Project, Horizon Europe Program for Research and Innovation.

EuroHPC Project DD23167, Sparse Matrix Computations at Extreme Scales for
the access to the Karolina supercomputer at the eINFRA CZ facilities (Czech Republic),
from 5/12/2023.
 Funding for 5600.00 €
for the 2023 GNCS INdAM project "Metodi basati su matrici e tensori strutturati
per problemi di algebra lineare di grandi dimensioni". P.I. S. Massei.
Participants: F. Durastante, L. Robol, B. Meini, F. Poloni,
M. Benzi, B. Iannazzo, L. Aceto, I. Simunec, C. Faccio, M. Rinelli, A.
Bucci, M. Viviani.
 Fundings for
€ 5943,65 for the "Progetti Speciali per la Didattica
Università di Pisa 2022/23". Project: Calcolo Parallelo dall'Infrastruttura alla Matematica.

Leonardo Early Access Program (LEAP). Awarderd computing time on the
Leonardo HPC System for the project
"PSCToolkit for Sparse Matrix Computations at Extreme Scales". Principal
Investigator S. Filippone and CoPrincipal Investigator P. D'Ambra.
Participants: S. Filippone, P. D'Ambra, F. Durastante, V. Cardellini, A. Celestini.
 Funding for 2700.00 €
for the 2022 GNCS INdAM project "Metodi numerici per l’analisi di modelli
innovativi di reti complesse" with Principal Investigator: Fabio Durastante
and participants: B. Meini, P. Boito, S. Massei, G. Del Corso, M. Benzi,
D. Bertaccini, P. D'Ambra, B. Iannazzo, A. Bucci, C. Faccio, M. Rinelli,
I. Simunec.
 Funding for 11656.03 €
for the project "¿Ver o no ver? Una aproximación interdisciplinar a los
estudios de visibilidad en el contexto de la navegación antigua (siglos IXIV a.C.)."
(PR27/21018) cofinanced by the Comunidad de Madrid and the Universidad
Complutense. Principal Investigator C. M. Mauro.
 Member of
the WorkPackage 6 (WP6) "Applications and Use cases" in the European
TEXTAROSSA project.
 Funding for
15000€ for the organization of the INdAM workshop on "Fractional Differential Equations:
Modeling, Discretization, and Numerical Solvers" (
https://fractionalworkshop.github.io/). Organizing committee:
Angelamaria Cardone, Marco Donatelli, Fabio Durastante, Roberto Garrappa,
Mariarosa Mazza, Marina Popolizio.
 Member of
the WorkPackage 3 (WP3) "Scalable Solvers" in the European
EoCoEII project.
 Writer for the Deliverable: ``D3.2
Preliminary results and performance evaluation of
Linear Algebra solvers'',
 Contributor for the Deliverable: ``D3.3
Updated results and new releases of LA solvers.''

Funding for 1200€ for the 2020 GNCSINDAM Project "Nonlocal models for
the analysis of complex networks" with Principal Investigator F. Durastante

Funding for 3400€ for the 2019 GNCSINDAM Project "Tecniche innovative e
parallele per sistemi lineari e non lineari di grandi dimensioni, funzioni ed equazioni
matriciali ed applicazioni" with Principal Investigator D. Bertaccini together
with L. Bergamaschi, P. D’Ambra, M. Ferronato, C. Janna, F. Marcuzzi, A.
Martinez Calomardo, V. Simoncini, M. Porcelli, D. Brandoni

Funding for 4500€ for the 2018 GNCSINDAM project "Tecniche innovative
per problemi di algebra lineare" with Principal Investigator D. Bertaccini together
with D. Bini, G. Del Corso, D. Fasino, L. Gemignani, B. Iannazzo, N.
Mastronardi, B. Meini, F. Poloni, A. Fayyaz, P. Boito, M. Fasi, I. Furci, S.
Massei and L. Robol.
Membership in scientific societies

Istituto Nazionale di Alta Matematica (INdAM)  Gruppo Nazionale di Calcolo Scientifico (GNCS).

Società Italiana di Matematica Applicata e Industriale (SIMAI).

Society of Industrial and Applied Mathematics (SIAM)  Activity group in Linear Algebra.

International Linear Algebra Society (ILAS).

Unione Matematica Italiana (UMI).

ACM: Association for Computing Machinery.
Collaboration Network
Most of the work described here has been done with several collaborators. I am
always working on expanding my network of collaborators by looking for people to
do interesting things with.
Bibliometric indexes
Not everything that counts can be counted, and not everything that can be counted counts.


(Scopus) HIndex: 8 
(Scopus) Subject Areas 
Editorial Work
I serve as editor for:
I have been part of the
 AD/AE Appendices Process & Badges committee for the I Am HPC International Conference for High Performance
Computing, Networking, Storage, and Analysis; Denver, Colorado November 1217, 2023.
I have served as reviewer for manuscripts submitted to the following journals.