subroutine tri dag

Version imprimable

subroutine tri dag

salut tous le monde

A quoi sert la subroutine tri dag (j'ai le code en fortran):

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
subroutine tridag(if,l,a,b,c,d,rr)
      implicit real*8(A-h,o-z)
      dimension A(1),b(1),c(1),d(1),rr(1),beta(20000),gamma(20000)
      beta(if)=b(if)
      gamma(if)=d(if)/beta(if)
      ifp1=if+1
      do 1 i=ifp1,l
      beta(i)=b(i)-a(i)*c(i-1)/beta(i-1)
      gamma(i)=(d(i)-a(i)*gamma(i-1))/beta(i)
   1  continue
      rr(l)=gamma(l)
      last=l-if
      do 2 k=1,last
      i=l-k
      rr(i)=gamma(i)-c(i)*rr(i+1)/beta(i)
   2  continue
      return
      end

merci d'avance

24/03/2011, 00h58
FR119492

Salut!

En plus du fait que tu ne sais pas à quoi sert cette routine, tu peux directement la mettre à la poubelle, car elle est écrite d'une manière totalement obsolète.

Par exemple:
Es-tu certain que son nom est tridag et non tridiag?

Code:

implicit real*8(A-h,o-z)

Cette instruction revient à répéter ce qui était de toute façon l'option par défaut. En revanche, il vaudrait mieux la remplacer par

Code:

implicit none

afin que les variables non déclarées soient signalées par des messages d'erreur à la compilation.

Code:

dimension A(1),b(1),c(1),d(1),rr(1) ...

Cette manière d'écrire les dimensions des tableaux n'est plus acceptée par certains compilateurs. Il faudrait écrire

Code:

dimension A(*),b(*),c(*),d(*),rr(*) ...

Code:

beta(if)=b(if)

C'est introduire un risque inutile de confusion que de donner à une variable le nom d'une instruction.

Il y a peut-être d'autres horreurs, mais je crois que ça suffit pour le moment.

Jean-Marc Blanc
25/03/2011, 22h51
Aleph69

Bonsoir,

a priori, ton code correspond à ceci:
http://www.aoe.vt.edu/~mason/Mason_f/CAtxtAppH1.pdf