Télécharger elfree.eso

Retour à la liste

Numérotation des lignes :

  1. C ELFREE SOURCE CHAT 05/01/12 23:32:33 5004
  2. SUBROUTINE ELFREE(V,A,D)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. C
  6. C ====================================================================
  7. C SOUS-PROGRAMME FORTRAN APPELE PAR ELFRES POUR LE CALCUL DU PRODUIT
  8. C DE CONVOLUTION PARTIEL SUR 1 ELEMENT ET 1 PAS DE TEMPS :
  9. C
  10. C VN=AN*DN+BN*CN
  11. C
  12. C V : VECTEUR SECOND MEMBRE VN
  13. C A : TERMES DES MATRICES AN ET BN
  14. C D : VECTEUR DN ET CN
  15. C
  16. C CREATION : 3/11/86
  17. C PROGRAMMEUR : GUILBAUD
  18. C ====================================================================
  19. C
  20. DIMENSION V(*),A(*),D(*)
  21. C
  22. V( 1)=V( 1)+A( 2)*D( 1) +A(16)*D(13)
  23. V( 2)=V( 2)+A( 8)*D( 2)+A( 7)*D( 6)+A(22)*D(14)-A(21)*D(18)
  24. V( 3)=V( 3)+A(13)*D( 3)-A(12)*D( 5)+A(27)*D(15)+A(26)*D(17)
  25. V( 4)=V( 4) +A( 4)*D( 4) +A(18)*D(16)
  26. V( 5)=V( 5)-A(14)*D( 3)+A(13)*D( 5)+A(28)*D(15)+A(27)*D(17)
  27. V( 6)=V( 6)+A( 9)*D( 2)+A( 8)*D( 6)-A(23)*D(14)+A(22)*D(18)
  28. V( 7)=V( 7)+A(16)*D( 1) +A( 2)*D(13)
  29. V( 8)=V( 8)+A(22)*D( 2)+A(21)*D( 6)+A( 8)*D(14)-A( 7)*D(18)
  30. V( 9)=V( 9)+A(27)*D( 3)-A(26)*D( 5)+A(13)*D(15)+A(12)*D(17)
  31. V(10)=V(10) +A(18)*D( 4) +A( 4)*D(16)
  32. V(11)=V(11)-A(28)*D( 3)+A(27)*D( 5)+A(14)*D(15)+A(13)*D(17)
  33. V(12)=V(12)+A(23)*D( 2)+A(22)*D( 6)-A( 9)*D(14)+A( 8)*D(18)
  34. C
  35. V( 1)=V( 1)+A( 1)*D( 7) +A(15)*D(19)
  36. V( 2)=V( 2)+A( 5)*D( 8) +A(19)*D(20)+A(20)*D(24)
  37. V( 3)=V( 3)+A(10)*D( 9) +A(24)*D(21)-A(25)*D(23)
  38. V( 4)=V( 4) +A( 3)*D(10) +A(17)*D(22)
  39. V( 5)=V( 5) -A(12)*D(11)+A(25)*D(21)-A(26)*D(23)
  40. V( 6)=V( 6) -A( 7)*D(12)-A(20)*D(20)-A(21)*D(24)
  41. V( 7)=V( 7)+A(15)*D( 7) +A( 1)*D(19)
  42. V( 8)=V( 8)+A(19)*D( 8)-A(20)*D(12)+A( 5)*D(20)
  43. V( 9)=V( 9)+A(24)*D( 9)+A(25)*D(11)+A(10)*D(21)
  44. V(10)=V(10) +A(17)*D(10) +A( 3)*D(22)
  45. V(11)=V(11)-A(25)*D( 9)-A(26)*D(11) -A(12)*D(23)
  46. V(12)=V(12)+A(20)*D( 8)-A(21)*D(12) -A( 7)*D(24)
  47. C
  48. C WRITE(6,*) ' AN '
  49. C WRITE(6,1001)(A(J),J=1,28)
  50. C WRITE(6,*) ' DN '
  51. C WRITE(6,1001)(D(J),J=1,24)
  52. C WRITE(6,*) ' VN '
  53. C WRITE(6,1001)(V(J),J=1,12)
  54. C1001 FORMAT(1X,12(1PE9.2,1X))
  55. RETURN
  56. END
  57.  
  58.  

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