Algoritma Thomas
Algoritma Thomas
Algoritma Thomas
Tu=q
Dimana T adalah matriks tridiagonal, u adalah vektor variabel, dan q adalah vektor
konstanta.
SPL tersebut dapat diselesaikan dengan menggunakan Algoritma Thomas yang
merupakan penyederhanaan dari Eliminasi Gauss.
Algoritma Thomas adalah sebagai berikut.
1. Memfaktorkan matriks T menjadi matriks segitiga bawah B dan matriks
segitiga atas A, dimana B dan A adalah matriks bidiagonal. Langkah ini
dikenal dengan LU decomposition.
t 11 t 12
t 21 t 22
0 0
t 23 0
0 0 t nn1
][
1 0 0 0
b
1 0 0
= 2
0 0 bn 1
t nn
][
d 1 a1
0 d2
0 0 0 dn
T
Dimana
0 0
a2 0
d 1=t 11 , a1=t 12
Untuk
i=2, 3, , n1
bi=t i ,i 1 /d i1
d i=t i ,ibi ai1
ai=t i ,i +1
bn =t n , n1 /d n1
d n=t n ,n bn an 1
2. Substitusi maju dengan melibatkan matriks B;
1 0 0 0
b2 1 0 0
0 0 bn 1
Dimana
][ ] [ ]
r1
q1
r2
q
= 2
rn
qn
r 1=q 1
r i=qib i r i1 ,
untuk
i=2, 3, , n
d 1 a1
0 d2
0 0
a2 0
Br =q
Au=r
][ ] [ ]
u1
r1
u2
r
= 2
u
r
0 0 0 dn n
n
Dimana
un=r n /d n
ui=(r iai ui +1)/d i ,
untuk
FALSE
Reff: http://www.cfm.brown.edu/people/gk/chap6/node13.html
Start
LU decomposition
j=0
j=j+1
j <= nt-1