Télécharger luiw.eso

Retour à la liste

Numérotation des lignes :

  1. C LUIW SOURCE CHAT 05/01/13 01:26:37 5004
  2. CCC
  3. C Producto v=(A)inv*w con A almacenada como LU
  4. CCC
  5. SUBROUTINE LUIW (A,W,V,N)
  6. IMPLICIT INTEGER(I-N)
  7. integer n,i,j
  8. real*8 v(n),w(n),A(n,n)
  9. c forward subst L v = w
  10. v(1)=w(1)/A(1,1)
  11. do i=2,n
  12. v(i)=w(i)
  13. do j=1,i-1
  14. v(i)=v(i)-A(i,j)*v(j)
  15. enddo
  16. v(i)=v(i)/A(i,i)
  17. enddo
  18. c backward subs U v = v
  19. do i=n-1,1,-1
  20. do j=i+1,n
  21. v(i)=v(i)-A(i,j)*v(j)
  22. enddo
  23. enddo
  24. return
  25. end
  26.  
  27.  
  28.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales