Lecture 22 - Forward Substitution [A]cfw_x = cfw_b [L]...

Info iconThis preview shows pages 1–8. Sign up to view the full content.

View Full Document Right Arrow Icon
Forward Substitution [ ] { } = + + + + + + = = 4 3 2 1 4 44 3 43 2 42 1 41 3 33 2 32 1 31 2 22 1 21 1 11 4 3 2 1 44 43 42 41 33 32 31 22 21 11 b b b b d l d l d l d l d l d l d l d l d l d l d d d d l l l l 0 l l l 0 0 l l 0 0 0 l d L - - - = - - = - = = 44 3 43 2 42 1 41 4 4 33 2 32 1 31 3 3 22 1 21 2 2 11 1 1 l d l d l d l b d l d l d l b d l d l b d l b d / ) ( / ) ( / ) ( / Very efficient for large matrices ! Once [ L ] is formed, we can use forward substitution instead of forward elimination for different { b }’s [L] [U]{x} = {b} [L] {d} = {b} [A] {x} = {b}
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Back Substitution [ ] { } = + + + + + + = = 4 3 2 1 4 44 4 34 3 33 4 24 3 23 2 22 4 14 3 13 2 12 1 11 4 3 2 1 44 34 33 24 23 22 14 13 12 11 d d d d x u x u x u x u x u x u x u x u x u x u x x x x u 0 0 0 u u 0 0 u u u 0 u u u u x U - - - = - - = - = = 11 4 14 3 13 2 12 1 1 22 4 24 3 23 2 2 33 4 34 3 3 44 4 4 u x u x u x u d x u x u x u d x u x u d x u d x / ) ( / ) ( / ) ( / Identical to Gauss elimination
Background image of page 2
Forward Substitution 33 14(2) 6 1 14d d 6d 1 d 2 (1/2)d 2 d 0 1 1 d 1 d 1 d 3 2 1 4 2 3 1 2 1 - = - - = - + - = = - = = + - = + - = = [ ] { } { } b 1 2 1 1 d d d d 1 14 1 6 0 1 1/2 0 0 0 1 1 0 0 0 1 d L 4 3 2 1 = - = - = { } - = 33 2 0 1 d Example:
Background image of page 3

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Back-Substitution - = - - = = - = - = - = = - - = 13/70 x 3 x 2 1 x 8/35 x 2 x 4/35 2 x 4 x 33/70 70 33/ x 4 3 1 3 2 4 3 4 - - = 33/70 4/35 8/35 13/70 x } { [ ] { } 33 2 0 1 x x x x 70 0 0 0 4 1 0 0 0 4 2 0 3 2 0 1 x U 4 3 2 1 - = - - = {d}
Background image of page 4
Forward and Back Substitutions Forward-substitution Back-substitution (identical to Gauss elimination) - = = - = l i 1 j j ij i i n , 1,2, i for d l b d 1 2 3 2 n 1 n i for u x u d x a d x ii n 1 i j j ij i i nn n n , , , , , / - - = - = = + =
Background image of page 5

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Forward and Back Substitutions
Background image of page 6
Still need pivoting in LU decomposition Messes up order of [ L ] What to do? Need to pivot both [
Background image of page 7

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Image of page 8
This is the end of the preview. Sign up to access the rest of the document.

This note was uploaded on 11/15/2010 for the course EGM 5403 taught by Professor Mei during the Spring '10 term at University of Florida.

Page1 / 29

Lecture 22 - Forward Substitution [A]cfw_x = cfw_b [L]...

This preview shows document pages 1 - 8. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online