Next: Matrix multiplication with reduced
Up: Array Sections
Previous: Two-dimensional Fourier transform
If is a symmetric positive definite matrix, associated
linear equations are often using Choleski decomposition:
where is a lower triangular matrix. In practise this is
followed by forward and back substitutions:
to complete the solution of .
A pseudocode algorithm for Cholesky decomposition is
A parallel version, assuming the main array is stored by columns
with the rows cyclically distributed, is given in figure 4.4.
The array is accumulated in the lower part of the input array a. Note that the array b has a replicated distribution, so the
remap operation is a broadcast of the relevant part of column