CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Main CFD Forum (http://www.cfd-online.com/Forums/main/)
-   -   Tridiagonal System (http://www.cfd-online.com/Forums/main/3496-tridiagonal-system.html)

 Zhong Lei May 21, 2001 04:47

Tridiagonal System

Hi, I wanna if there is any method to efficiently solve a 1-d tridiagonal linear system quation on a vector machine. Any idea?

 Ray-zoo May 22, 2001 22:04

Re: Tridiagonal System

Solution from tridiagonal linear equations could be computed with Thoma's algorithm, wich is "pivot de Gauss" method for this particular case.

Do you really need a vector machine for this ?

Best regards,

Christophe

 Zhong Lei May 23, 2001 00:04

Re: Tridiagonal System

The vector rate of Thoma's algorithm is too low when it runs on a parallel machine. I want to improve the vectorization. of the following do-loop.

<DD> do I=1,IMAX <DD> D(I) = D(I) -W(I)*D(I-1) <DD> enddo </DD>

Thanks Zhong Lei

 Ray-zoo May 24, 2001 15:18

Re: Tridiagonal System

Dear Zhong Lei,

Apologises for this late answer due to sparse internet connection.

After a read of J.C. Chien's post about vector machines, I could not guess how to vectorize a loop which requires previous iteration step result.

As I am mostly newbie to this forum, I can not confirm own opinion that meaningless requests may be a reason for none, or very few answers.

If you need running some parallel machines, couldn't you run some domain decomposition method (ddm) ?

I had an internet link to those methods, but could not find it before this post.

Maybe someone else here knows where to find informations about ddm(s).

Yours,

Christophe

 All times are GMT -4. The time now is 16:53.