Télécharger profe1.eso

Retour à la liste

Numérotation des lignes :

  1. C PROFE1 SOURCE BR232186 16/12/06 10:42:34 9242
  2. SUBROUTINE PROFE1(X,Y,II,JJ,KK,LL,z)
  3. C
  4. C====&===1=========2=========3=========4=========5=========6=========7==
  5. C Commentaires : Subroutine permettant de calculer le produit
  6. C de Kronecker entre deux tenseurs X et Y.
  7. C Si X est une matrice 3x2 et B une matrice 2x4
  8. C le produit de Kronecker Z sera de dimension (3*2)x(2*4)
  9. C Exemple :
  10. C [ X(1,1)*Y X(1,2)*Y
  11. C X(2,1)*Y X(2,2)*Y
  12. C X(3,1)*Y X(3,2)*Y ]
  13. C
  14. C
  15. C Argumets : Matrices X(II,JJ) et Y(KK,LL),retourne la matrice Z
  16. C
  17. C AUTEUR : E. KISHTA - CEA/DEN/DANS/DM2S/SEMT/EMSI
  18. C
  19. C====&===1=========2=========3=========4=========5=========6=========7==
  20. C
  21. C----DECLARATION GENERALES----------------------------------------------
  22. C
  23. IMPLICIT INTEGER(I-N)
  24. IMPLICIT REAL*8(A-H,O-Z)
  25.  
  26. REAl*8 x(ii,jj), y(kk,ll), z(ii*kk,jj*ll)
  27.  
  28.  
  29. do i = 1,ii
  30. do j = 1,jj
  31. do k = 1,kk
  32. do l = 1,ll
  33.  
  34. z((kk*(i-1) + k), (ll*(j-1) + l)) = x(i,j)*y(k,l)
  35.  
  36. enddo
  37. enddo
  38. enddo
  39. enddo
  40.  
  41.  
  42.  
  43. RETURN
  44. END
  45.  
  46.  
  47.  

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