作者cl36260 (小小小宝贝)
看板Fortran
标题Re: [问题] TDMA
时间Sun Jan 17 11:56:13 2010
我回应在板上一次好了QQ 有人寄信来跟我要
= = = = = =
| a b | | x | | d |
| a b c | | x | | d |
| a b c | x | x | = | d |
| . . . | | x | | d |
| . . . | | x | | d |
| b c | | x | | d |
= = = = = =
subroutine TDMA(a,b,c,d,x,n)
real :: a(19),b(19),c(19),d(19),x(19)
integer::N, I
c(2)=c(2)/b(2)
d(2)=d(2)/b(2)
do i=3,N
c(i)=c(i)/(b(i)-a(i)*c(i-1))
d(i)=(d(i)-a(i)*d(i-1))/(b(i)-a(i)*c(i-1))
end do
x(N)=d(N)
do i=N-1,2,-1
x(i)=d(i)-x(i+1)*c(i)
end do
return
end subroutine
这样应该看的懂吧 囧 不是很会画图..
※ 引述《libra1005 (PTT乡民)》之铭言：
: 请问一下，现在我正在修计算流体力学这门课，有用到TDMA这东西
: 想请问一下这要如何写出来??
: 我目前对他所认识只是知道他是三带宽矩阵而已，却不知如何写成
: 程式码？请板上高手解惑
: 感激不尽~~
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.134.115
※ 编辑: cl36260 来自: 140.113.134.115 (01/18 19:14)
1F：推 imys:谢谢分享喔～希望这个版能够越来越多人 02/03 16:06