This research investigates the implementation mechanism of block-wise ILU(k) preconditioner on GPU. The block-wise ILU(k) algorithm requires both the level k and the block size to be designed as variables. A decoupled ILU(k) algorithm consists of a symbolic phase and a factorization phase. In the symbolic phase, a ILU(k) nonzero pattern is established from the point-wise structure extracted fro...

Run-time support for the CYCLIC(k) redistribution on the SPMD computation model is presently very relevant for the scientific community. This work is focused to the characterization of the sparse matrix redistribution and its associate problematic due to the use of compressed representations. Two main improvements about the buffering and the coordinates calculation modify the original algorithm...

In this paper, the numerical technique based on hybrid Bernoulli and Block-Pulse functions has been developed to approximate the solution of system of linear Volterra integral equations. System of Volterra integral equations arose in many physical problems such as elastodynamic, quasi-static visco-elasticity and magneto-electro-elastic dynamic problems. These functions are formed by the hybridi...

A set of new formulae for the inverse of a block Hankel (or block Toeplitz) matrix is given. The formulae are expressed in terms of certain matrix Pad6 forms, which approximate a matrix power series associated with the block Hankel matrix. By using Frobenius-type identities between certain matrix Pad6 forms, the inversion formulae are shown to generalize the formulae of Gohberg-Heinig and, in t...

In this paper, a so-called auxiliary matrix polynomial Xn(z) and a true right formal orthogonal matrix polynomial (FOMP) A 1 n (z) is connected to each well-conditioned leading principal block submatrix of a given block Toeplitz matrix. From these two matrix polynomi-als, all other right FOMPs of block n of a system of block biorthogonal matrix polynomials with respect to the block Toeplitz mom...

Block coherence of matrices plays an important role in analyzing the performance of block compressed sensing recovery algorithms (Bajwa and Mixon, 2012). In this paper, we characterize two block coherence metrics: worstcase and average block coherence. First, we present lower bounds on worst-case block coherence, in both the general case and also when the matrix is constrained to be a union of ...

lsmr (least squares minimal residual) is an iterative method for the solution of the linear system of equations and leastsquares problems. this paper presents a block version of the lsmr algorithm for solving linear systems with multiple right-hand sides. the new algorithm is based on the block bidiagonalization and derived by minimizing the frobenius norm of the resid ual matrix of normal equa...

Let an n x -matrix A have m < (m ? 2) different eigenvalues ?j of the algebraic multiplicity (j = 1,..., m). It is proved that there are ?j-matrices Aj, each which has a unique eigenvalue ?j, such similar to block-diagonal matrix ?D diag (A1,A2,..., Am). I.e. invertible T, T-1AT ?D. Besides, sharp bound for number kT := ||T||||T-1|| derived. As applications these results we obtain norm estim...

In this paper, block distance matrices are introduced. Suppose F is a square block matrix in which each block is a symmetric matrix of some given order. If F is positive semidefinite, the block distance matrix D is defined as a matrix whose (i, j)-block is given by D ij = F ii +F jj −2F ij. When each block in F is 1 × 1 (i.e., a real number), D is a usual Euclidean distance matrix. Many interes...