Block LU Decomposition: explained

1
Block LU Decomposition explained

• Keiran OHaire

2
Introduction

• Method for performing LU decomposition on large
  matrices
• Requires more calculations
• Faster due to memory locality

3
The Matrix
4
Notes

• Each of the boxes represent sub matrices of the
  original matrix
• The sub matrices size are determined by the
  block size
• Example sub matrix with block size 4

5
Matrix Dimensions
b
n-b
b
n-b
6
Step 1

• LU Decomposition of left-most column

7
Step 2

• Solve the rest of the top row using the topmost L

8
Step 2 Breakdown
Example L
Example U
9
Column Version

• Use one column at a time from the matrix you wish
  to solve

10
Column Version Contd
The result!
11
Column Version Contd

• Now insert the next column
• Repeat until all columns are solved

12
Row Version

• Same matrices as before, but this time were doing
  rows
• The version is dependent on the language you are
  using, such as C or FORTRAN due to array
  orientation

13
Row Version Contd

• This time, solve each of these (one cell at a
  time) in a row.
• i.e. solve 1 0 3
• Clearly, for the first, x13, and for the second,
  x11

14
Row Version Contd

• These calculated xs are now used in the next
  row.
• i.e., 3(x13) x2 2, x2-7
• Then, 3(x11) x2 5, x22

15
Continued

• So now we have calculated this portion. Now..

16
The Next Step

• Continue recursively

17
The Completed Decomposition

• Eventually you end up with the solution

18
Questions?
Questions?