Journal of Applied Mathematics

Volume 2012 (2012), Article ID 872901, 14 pages

http://dx.doi.org/10.1155/2012/872901

## Parallel Rayleigh Quotient Optimization with FSAI-Based Preconditioning

Department of Mathematical Methods and Models for Scientific Applications, University of Padova, Via Trieste 63, 35121 Padova, Italy

Received 2 November 2011; Revised 1 February 2012; Accepted 3 February 2012

Academic Editor: Massimiliano Ferronato

Copyright © 2012 Luca Bergamaschi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The present paper describes a parallel preconditioned algorithm for the solution of partial eigenvalue problems for large sparse symmetric matrices, on parallel computers. Namely, we consider the Deflation-Accelerated Conjugate Gradient (DACG) algorithm accelerated by factorized-sparse-approximate-inverse- (FSAI-) type preconditioners. We present an enhanced parallel implementation of the FSAI preconditioner and make use of the recently developed Block FSAI-IC preconditioner, which combines the FSAI and the Block Jacobi-IC preconditioners. Results onto matrices of large size arising from finite element discretization of geomechanical models reveal that DACG accelerated by these type of preconditioners is competitive with respect to the available public parallel *hypre* package, especially in the computation of a few of the leftmost eigenpairs. The parallel DACG code accelerated by FSAI is written in MPI-Fortran 90 language and exhibits good scalability up to one thousand processors.

#### 1. Introduction

The computation by iterative methods of the partial eigenspectrum of the generalized eigenproblem: where are large sparse symmetric positive definite (SPD) matrices, is an important and difficult task in many applications. It has become increasingly widespread owing to the development in the last twenty years of robust and computationally efficient schemes and corresponding software packages. Among the most well-known approaches for the important class of symmetric positive definite (SPD) matrices are the implicitly restarted Arnoldi method (equivalent to the Lanczos technique for this type of matrices) [1, 2], the Jacobi-Davidson (JD) algorithm [3], and schemes based on preconditioned conjugate gradient minimization of the Rayleigh quotient [4, 5].

The basic idea of the latter is to minimize the Rayleigh Quotient
in a subspace which is orthogonal to the previously computed eigenvectors via a preconditioned CG-like procedure. Among the different variants of this technique we chose to use the Deflation-Accelerated Conjugate Gradient (DACG) scheme [4, 6] which has been shown to be competitive with the Jacobi Davidson method and with the PARPACK package [7]. As in any other approach, for our DACG method, the choice of the preconditioning technique is a key factor to accelerate and, in some cases even to allow for, convergence. To accelerate DACG in a parallel environment we selected the Factorized Sparse Approximate inverse (FSAI) preconditioner introduced in [8]. We have developed a parallel implementation of this algorithm which has displayed excellent performances on both the setup phase and the application phase within a Krylov subspace solver [9–11]. The effectiveness of the FSAI preconditioner in the acceleration of DACG is compared to that of the Block FSAI-IC preconditioner, recently developed in [12], which combines the FSAI and the Block Jacobi-IC preconditioners obtaining good results on a small number of processors for the solution of SPD linear systems and for the solution of large eigenproblems [13]. We used the resulting parallel codes to compute a few of the leftmost eigenpairs of a set of test matrices of large size arising from Finite Element discretization of geomechanical models. The reported results show that DACG preconditioned with either FSAI or BFSAI is a scalable and robust algorithm for the partial solution of SPD eigenproblems. The parallel performance of DACG is also compared to that of the publicly available parallel package *hypre* [14] which implements a number of preconditioners which can be used in combination with the Locally Optimal Block PCG (LOBPCG) iterative eigensolver [15]. The results presented in this paper show that the parallel DACG code accelerated by FSAI exhibits good scalability up to one thousand processors and displays comparable performance with respect to *hypre*, specially when a low number of eigenpairs is sought.

The outline of the paper is as follows: in Section 2 we describe the DACG Algorithm; in Sections 3 and 4 we recall the definition and properties of the FSAI and BFSAI preconditioners, respectively. Section 5 contains the numerical results obtained with the proposed algorithm in the eigensolution of very large SPD matrices of size up to almost 7 million unknowns and nonzeros. A comparison with the *hypre* eigensolver code is also included. Section 6 ends the paper with some conclusions.

#### 2. The DACG Iterative Eigensolver and Implementation

The DACG algorithm sequentially computes the eigenpairs, starting from the leftmost one . To evaluate the th eigenpair, , DACG minimizes the Rayleigh Quotient (RQ) in a subspace orthogonal to the eigenvectors previously computed. More precisely, DACG minimizes the Rayleigh Quotient: where The first eigenpair is obtained by minimization of (2.1) with . Indicating with the preconditioning matrix, that is, , the leftmost eigenpairs are computed by the conjugate gradient procedure [6] described in Algorithm 1.

The schemes relying on the Rayleigh quotient optimization are quite attractive for parallel computations; however preconditioning is an essential feature to ensure practical convergence. When seeking for an eigenpair it can be proved that the number of iterations is proportional to the square root of the condition number of the Hessian of the Rayleigh quotient in the stationary point [4]. It turns out that is similar to which is not SPD. However, operates on the orthogonal space spanned by the previous eigenvectors, so that the only important eigenvalues are the positive ones. In the non-preconditioned case (i.e., ) we would have where in the ideal case , we have Therefore, even though is not the optimal preconditioner for , however, if is a good preconditioner of then the condition number will approach .

#### 3. The FSAI Preconditioner

The FSAI preconditioner, initially proposed in [8, 16], has been later developed and implemented in parallel by Bergamaschi and Martínez in [9]. Here, we only shortly recall the main features of this preconditioner. Given an SPD matrix the FSAI preconditioner approximately factorizes its inverse as a product of two sparse triangular matrices as
The choice of nonzeros in is based on a sparsity pattern which in our work may be the same as where is the result of *prefiltration* [10] of , that is, dropping of all elements below of a threshold parameter . The entries of are computed by minimizing the Frobenius norm of , where is the exact Cholesky factor of , without forming explicitly the matrix . The computed is then sparsified by dropping all the elements which are below a second tolerance parameter (. The final FSAI preconditioner is therefore related to the following three parameters: , prefiltration threshold; , power of generating the sparsity pattern (we allow in our experiments); , postfiltration threshold.

##### 3.1. Parallel Implementation of FSAI-DACG

We have developed a parallel code written in FORTRAN 90 and which exploits the MPI library for exchanging data among the processors. We used a block row distribution of all matrices (, , and ), that is, with complete rows assigned to different processors. All these matrices are stored in static data structures in CSR format.

Regarding the preconditioner computation, we stress that any row of matrix of FSAI preconditioner is computed independently of each other, by solving a small SPD dense linear system of size equal to the number of nonzeros allowed in row of . Some of the rows which contribute to form this linear system may be nonlocal to processor and should be received from other processors. To this aim we implemented a routine called *get_extra_rows* which carries out all the row exchanges among the processors, before starting the computation of , which proceed afterwards entirely in parallel. Since the number of nonlocal rows needed by each processor is relatively small we chose to temporarily replicate these rows on auxiliary data structures. Once is obtained a parallel transposition routine provides every processor with its part of .

The DACG iterative solver is essentially based on scalar and matrix-vector products. We made use of an optimized parallel matrix-vector product which has been developed in [17] showing its effectiveness up to 1024 processors.

#### 4. Block FSAI-IC Preconditioning

The Block FSAI-IC preconditioner, BFSAI-IC in the following, is a recent development for the parallel solution to Symmetric Positive Definite (SPD) linear systems. Assume that is an arbitrary nonsingular block diagonal matrix consisting of equal size blocks.

Let and be a sparse lower triangular and a dense block diagonal nonzero pattern, respectively, for an matrix. Even though not strictly necessary, for the sake of simplicity assume that consists of diagonal blocks with equal size and let be an arbitrary full-rank matrix with nonzero pattern .

Consider the set of lower block triangular matrices with a prescribed nonzero pattern and minimize over the Frobenius norm: where is the exact lower Cholesky factor of an SPD matrix . A matrix satisfying the minimality condition (4.1) for a given is the lower block triangular factor of BFSAI-IC. Recalling the definition of the classical FSAI preconditioner, it can be noticed that BFSAI-IC is a block generalization of the FSAI concept.

The differentiation of (4.1) with respect to the unknown entries , , yields the solution to independent dense subsystems which, in the standard FSAI case, do not require the explicit knowledge of . The effect of applying to is to concentrate the largest entries of the preconditioned matrix into diagonal blocks. However, as is arbitrary, it is still not ensured that is better than in an iterative method, so it is necessary to precondition again. As resembles a block diagonal matrix, an efficient technique relies on using a block diagonal matrix which collects an approximation of the inverse of each diagonal block of .

It is easy to show that is guaranteed to exist with SPD matrices and is SPD, too [12]. Using an IC decomposition with partial fill-in for each block and collecting in the lower IC factors, the resulting preconditioned matrix reads with the final preconditioner in (4.3) is the BFSAI-IC preconditioner of .

For its computation BFSAI-IC needs the selection of and . The basic requirement for the number of blocks is to be larger than or equal to the number of computing cores . From a practical viewpoint, however, the most efficient choice in terms of both wall clock time and iteration count is to keep the blocks as large as possible, thus implying . Hence, is by default set equal to . By distinction, the choice of is theoretically more challenging and still not completely clear. A widely accepted option for other approximate inverses, such as FSAI or SPAI, is to select the nonzero pattern of for small values of on the basis of the Neumann series expansion of . Using a similar approach, in the BFSAI construction we select as the lower block triangular pattern of . As the nonzeros located in the diagonal blocks are not used for the computation of a larger value of , say 3 or 4, can still be used.

Though theoretically not necessary, three additional user-specified parameters are worth introducing in order to better control the memory occupation and the BFSAI-IC density:(1) is a postfiltration parameter that allows for dropping the smallest entries of . In particular, is neglected if , where is the th row of ;(2) is a parameter that controls the fill-in of and determines the maximum allowable number of nonzeros for each row of in addition to the corresponding entries of . Quite obviously, the largest entries only are retained;(3) is a parameter that controls the fill-in of each IC factor denoting the maximum allowable number of nonzeros for each row of in addition to the corresponding entries of .

An OpenMP implementation of the algorithms above is available in [18].

#### 5. Numerical Results

In this section we examine the performance of the parallel DACG preconditioned by both FSAI and BFSAI in the partial solution of four large-size sparse eigenproblems. The test cases, which we briefly describe below, are taken from different real engineering mechanical applications. In detail, they are as follows.(i)FAULT-639 is obtained from a structural problem discretizing a faulted gas reservoir with tetrahedral finite elements and triangular interface elements [19]. The interface elements are used with a penalty formulation to simulate the faults behavior. The problem arises from a 3D discretization with three displacement unknowns associated to each node of the grid.(ii)PO-878 arises in the simulation of the consolidation of a real gas reservoir of the Po Valley, Italy, used for underground gas storage purposes (for details, see [20]).(iii)GEO-1438 is obtained from a geomechanical problem discretizing a region of the earth crust subject to underground deformation. The computational domain is a box with an areal extent of 50 × 50 km and 10 km deep consisting of regularly shaped tetrahedral finite elements. The problem arises from a 3D discretization with three displacement unknowns associated to each node of the grid [21].(iv)CUBE-6091 arises from the equilibrium of a concrete cube discretized by a regular unstructured tetrahedral grid.

Matrices FAULT-639 and GEO-1438 are publicly available in the University of Florida Sparse Matrix Collection at http://www.cise.ufl.edu/research/sparse/matrices/.

In Table 1 we report sizes and nonzeros of the four matrices together with three of the most significant eigenvalues for each problem.

The computational performance of FSAI is compared to the one obtained by using BFSAI as implemented in [12]. The comparison is done evaluating the number of iterations to converge at the same tolerance, the wall clock time in seconds and for the preconditioner computation, and the eigensolver to converge, respectively, with the total time . All tests are performed on the IBM SP6/5376 cluster at the CINECA Centre for High Performance Computing, equipped with IBM Power6 processors at 4.7 GHz with 168 nodes, 5376 computing cores, and 21 Tbyte of internal network RAM. The FSAI-DACG code is written in Fortran 90 and compiled with -O4 -q64 -qarch=pwr6 -qtune=pwr6 -qnoipa -qstrict -bmaxdata:0x70000000 options. For the BFSAI-IC code only an OpenMP implementation is presently available.

To study parallel performance we will use a strong scaling measure to see how the CPU times vary with the number of processors for a fixed total problem size. Denote with the total CPU elapsed times expressed in seconds on processors. We introduce a relative measure of the parallel efficiency achieved by the code, , which is the pseudo speedup computed with respect to the smallest number of processors used to solve a given problem. Accordingly, we will denote by the corresponding efficiency:

##### 5.1. FSAI-DACG Results

In this section we report the results of our FSAI-DACG implementation in the computation of the 10 leftmost eigenpairs of the 4 test problems. We used the exit test described in the DACG algorithm (see Algorithm 1) with . The results are summarized in Table 2. As the FSAI parameters, we choose , , and for all the test matrices. This combination of parameters produces, on the average, the best (or close to the best) performance of the iterative procedure. Note that the number of iterations does not change with the number of processors, for a fixed problem. The scalability of the code is very satisfactory in both the setup stage (preconditioner computation) and the iterative phase.

##### 5.2. BFSAI-IC-DACG Results

We present in this section the results of DACG accelerated by the BFSAI-IC preconditioner for the approximation of the leftmost eigenpairs of the matrices described above.

Table 3 provides iteration count and total CPU time for BFSAI-DACG with different combinations of the parameters needed to construct the BFSAI-IC preconditioner for matrix PO-878 and using from 2 to 8 processors. It can be seen from Table 3 that the assessment of the optimal parameters, , , and , is not an easy task, since the number of iterations may highly vary depending on the number of processors. We chose in this case the combination of parameters producing the second smallest total time with processors. After intensive testing for all the test problems, we selected similarly the “optimal” values which are used in the numerical experiments reported in Table 4:(i)FAULT-639: , , , ,(ii)PO-878 , , , ,(iii)GEO-1438: , , , ,(iv)CUBE-6091: , , , .

The user-specified parameters for BFSAI-IC given above provide evidence that it is important to build a dense preconditioner based on the lower nonzero pattern of (except for CUBE-6091, which is built on a regular discretization) with the aim at decreasing the number of DACG iterations. Anyway, the cost for computing such a dense preconditioner appears to be almost negligible with respect to the wall clock time needed to iterate to convergence.

We recall that, presently, the code BFSAI-IC is implemented in OpenMP, and the results in terms of CPU time are significant only for . For this reason the number of iterations reported in Table 4 is obtained with increasing number of blocks and with processors. This iteration number accounts for a potential implementation of BFSAI-DACG under the MPI (or hybrid OpenMP-MPI) environment as the number of iterations depends only on the number of blocks, irrespective of the number of processors.

The only meaningful comparison between FSAI-DACG and BFSAI-DACG can be carried out in terms of iteration numbers which are smaller for BFSAI-DACG for a small number of processors. The gap between FSAI and BFSAI iterations reduces when the number of processors increases.

##### 5.3. Comparison with the LOBPCG Eigensolver Provided by *hypre*

In order to validate the effectiveness of our preconditioning in the proposed DACG algorithm with respect to already available public parallel eigensolvers, the results given in Tables 2 and 4 are compared with those obtained by the schemes implemented in the *hypre* software package [14]. The Locally Optimal Block Preconditioned Conjugate Gradient method (LOBPCG) [15] is experimented with, using the different preconditioners developed in the *hypre* project, that is, algebraic multigrid (AMG), diagonal scaling (DS), approximate inverse (ParaSails), additive Schwarz (Schwarz), and incomplete LU (Euclid). The *hypre* preconditioned CG is used for the inner iterations within LOBPCG. For details on the implementation of the LOBPCG algorithm, see, for instance, [22]. The selected preconditioner, ParaSails, is on its turn based on the FSAI preconditioner, so that the different FSAI-DACG and ParaSails-LOBPCG performances should be ascribed mainly to the different eigensolvers rather than to the preconditioners.

We first carried out a preliminary set of runs with the aim of assessing the optimal value of the block size bl parameter, that is, the size of the subspace where to seek for the eigenvectors. Obviously it must be bl ≥ *s* = 10. We fixed to 16 the number of processors and obtained the results summarized in Table 5 with different values of bl . We found that, only in problem CUBE-6091, a value of bl larger than 10, namely, bl = 12, yields an improvement in the CPU time. Note that we also made this comparison with different number of processors, and we obtained analogous results.

Table 6 presents the number of iterations and timings using the LOBPCG algorithm in the *hypre* package. The LOBPCG wall clock time is obtained with the preconditioner allowing for the best performance in the specific problem at hand, that is, ParaSails for all the problems. Using AMG as the preconditioner did not allow for convergence in three cases out of four, with the only exception of the FAULT-639 problem, in which the CPU timings were however very much larger than using ParaSails.

All matrices have to be preliminarily scaled by their maximum coefficient in order to allow for convergence. To make the comparison meaningful, the outer iterations of the different methods are stopped when the average relative error measure of the computed leftmost eigenpairs gets smaller than , in order to obtain a comparable accuracy as in the other codes. We also report in Table 6 the number of inner preconditioned CG iterations (pcgitr).

To better compare our FSAI DACG with the LOBPCG method, we depict in Figure 1 the total CPU time versus the number of processor for the two codes. FSAI-DACG and LOBPCG provide very similar scalability, being the latter code a little bit more performing on the average. On the FAULT-639 problem, DACG reveals faster than LOBPCG, irrespective of the number of processors employed.

Finally, we have carried out a comparison of the two eigensolvers in the computation of only the leftmost eigenpair. Differently from LOBPCG, which performs a simultaneous approximation of all the selected eigenpairs, DACG proceeds in the computation of the selected eigenpairs in a sequential way. For this reason, DACG should be the better choice, at least in principle, when just one eigenpair is sought. We investigate this feature, and the results are summarized in Table 7. We include the total CPU time and iteration count needed by LOBPCG and FSAI-DACG to compute the leftmost eigenpair with 16 processors. For the LOBPCG code we report only the number of outer iterations.

The parameters used to construct the FSAI preconditioner for these experiments are as follows: (1)FAULT-639. , , ,(2) PO-878. , , ,(3) GEO-1438. , , ,(4) CUBE-6091. , , .

These parameters differ from those employed to compute the FSAI preconditioner in the assessment of the 10 leftmost eigenpairs and have been selected in order to produce a preconditioner relatively cheap to compute. This is so because otherwise the setup time would prevail over the iteration time. Similarly, to compute just one eigenpair with LOBPCG we need to setup a different value for pcgitr, the number of inner iterations. As it can be seen from Table 7, in the majority of the test cases, LOBPCG takes less time to compute 2 eigenpairs than just only 1. FSAI-DACG reveals more efficient than the best LOBPCG on problems PO-878 and GEO-1438. On the remaining two problems the slow convergence exhibited by DACG is probably due to the small relative separation between and .

#### 6. Conclusions

We have presented the parallel DACG algorithm for the partial eigensolution of large and sparse SPD matrices. The scalability of DACG, accelerated with FSAI-type preconditioners, has been studied on a set of test matrices of very large size arising from real engineering mechanical applications. Our FSAI-DACG code has shown comparable performances with the LOBPCG eigensolver within the well-known public domain package, *hypre*. Numerical results reveal that not only the scalability achieved by our code is roughly identical to that of *hypre* but also, in some instances, FSAI-DACG proves more efficient in terms of absolute CPU time. In particular, for the computation of the leftmost eigenpair, FSAI-DACG is more convenient in 2 problems out of 4.

#### Acknowledgment

The authors acknowledge the CINECA Iscra Award SCALPREC (2011) for the availability of HPC resources and support.

#### References

- P. Arbenz, U. Hetmaniuk, R. Lehoucq, and R. Tuminaro, “A comparison of eigensolvers for large-scale 3D modal analysis using AMG-preconditioned
iterative methods,”
*International Journal for Numerical Methods in Engineering*, vol. 64, no. 2, pp. 204–236, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - R. B. Lehoucq and D. C. Sorensen, “Deflation techniques for an implicitly restarted Arnoldi iteration,”
*SIAM Journal on Matrix Analysis and Applications*, vol. 17, no. 4, pp. 789–821, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - G. L. G. Sleijpen and H. A. Van der Vorst, “A Jacobi-Davidson iteration method for linear eigenvalue problems,”
*SIAM Journal on Matrix Analysis and Applications*, vol. 17, no. 2, pp. 401–425, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - L. Bergamaschi, G. Gambolati, and G. Pini, “Asymptotic convergence of conjugate gradient methods for the partial symmetric eigenproblem,”
*Numerical Linear Algebra with Applications*, vol. 4, no. 2, pp. 69–84, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. V. Knyazev and A. L. Skorokhodov, “Preconditioned gradient-type iterative methods in a subspace for partial generalized symmetric eigenvalue problems,”
*SIAM Journal on Numerical Analysis*, vol. 31, no. 4, pp. 1226–1239, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - L. Bergamaschi, G. Pini, and F. Sartoretto, “Approximate inverse preconditioning in the parallel solution of sparse eigenproblems,”
*Numerical Linear Algebra with Applications*, vol. 7, no. 3, pp. 99–116, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - L. Bergamaschi and M. Putti, “Numerical comparison of iterative eigensolvers for large sparse symmetric positive definite matrices,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 191, no. 45, pp. 5233–5247, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - L. Yu. Kolotilina and A. Yu. Yeremin, “Factorized sparse approximate inverse preconditionings. I. Theory,”
*SIAM Journal on Matrix Analysis and Applications*, vol. 14, no. 1, pp. 45–58, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - L. Bergamaschi and A. Martínez, “Parallel acceleration of Krylov solvers by factorized approximate inverse preconditioners,” in
*VECPAR 2004*, M. Dayd et al., Ed., vol. 3402 of*Lecture Notes in Computer Sciences*, pp. 623–636, Springer, Heidelberg, Germany, 2005. View at Google Scholar · View at Zentralblatt MATH - L. Bergamaschi, A. Martínez, and G. Pini, “An efficient parallel MLPG method for poroelastic models,”
*Computer Modeling in Engineering & Sciences*, vol. 49, no. 3, pp. 191–215, 2009. View at Google Scholar - L. Bergamaschi and A. Martínez, “Parallel inexact constraint preconditioners for saddle point problems,” in
*Proceedings of the 17th International Conference on Parallel Processing (Euro-Par'11)*, R. N. E. Jeannot and J. Roman, Eds., vol. 6853, part 2, pp. 78–89, Springer, Bordeaux, France, 2011, Lecture Notes in Computer Sciences. - C. Janna, M. Ferronato, and G. Gambolati, “A block FSAI-ILU parallel preconditioner for symmetric positive definite linear systems,”
*SIAM Journal on Scientific Computing*, vol. 32, no. 5, pp. 2468–2484, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. Ferronato, C. Janna, and G. Pini, “Efficient parallel solution to large-size sparse eigenproblems with block FSAI preconditioning,”
*Numerical Linear Algebra with Applications*. In press. View at Publisher · View at Google Scholar - Lawrence Livermore National Laboratory, hypre user manual, software version 1.6.0. Center for Applied Scientific Computing (CASC), University of California, 2001.
- A. V. Knyazev, “Toward the optimal preconditioned eigensolver: locally optimal block preconditioned conjugate gradient method,”
*SIAM Journal on Scientific Computing*, vol. 23, no. 2, pp. 517–541, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - L. Yu. Kolotilina, A. A. Nikishin, and A. Yu. Yeremin, “Factorized sparse approximate inverse preconditionings. IV. Simple approaches to rising efficiency,”
*Numerical Linear Algebra with Applications*, vol. 6, no. 7, pp. 515–531, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - A. Martínez, L. Bergamaschi, M. Caliari, and M. Vianello, “A massively parallel exponential integrator for advection-diffusion models,”
*Journal of Computational and Applied Mathematics*, vol. 231, no. 1, pp. 82–91, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - C. Janna, M. Ferronato, and N. Castelletto, “BFSAI-IC OpenMP implementation,” Release V1.0, January 2011, http://www.dmsa.unipd.it/~ferronat/software.html.
- M. Ferronato, C. Janna, and G. Gambolati, “Mixed constraint preconditioning in computational contact mechanics,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 197, no. 45–48, pp. 3922–3931, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - N. Castelletto, M. Ferronato, G. Gambolati et al., “3D geomechanics in UGS projects: a comprehensive study in northern Italy,” in
*Proceedings of the 44th US Rock Mechanics Symposium*, Salt Lake City, Utah, USA, 2010. View at Publisher · View at Google Scholar - P. Teatini, M. Ferronato, G. Gambolati, D. Bau, and M. Putti, “An-thropogenic Venice uplift by seawater pumping into a heterogeneous aquifer system,”
*Water Resources Research*, vol. 46, Article ID W11547, 16 pages, 2010. View at Publisher · View at Google Scholar - A. V. Knyazev and M. E. Argentati, “Implementation of a preconditioned eigensolver using hypre,” 2005, http://math.ucdenver.edu/~rargenta/index_files/rep220.pdf.