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...

We reinvestigate the validity of mapping the problem of two onsite interacting particles in a random potential onto an eeective random matrix model. To this end we rst study numerically how the non-interacting basis is coupled by the interaction. Our results indicate that the typical coupling matrix element decreases signiicantly faster with increasing single-particle localization length than i...

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...

In this paper we deal with a single-antenna discrete-time flat-fading channel. The fading process is assumed to be stationary for the duration of a single data block. From block to block the fading process is allowed to be non-stationary. The number of scatterers bounds the rank of the channels covariance matrix. The signal-to-noise ratio (SNR), the user velocity, and the data block-length defi...

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...