Télécharger asksgl.eso

Retour à la liste

Numérotation des lignes :

  1. C ASKSGL SOURCE CHAT 07/10/22 21:15:07 5921
  2. SUBROUTINE ASKSGL(XK,R)
  3. IMPLICIT REAL*8(A-H,O-Z)
  4. IMPLICIT INTEGER (I-N)
  5.  
  6. * SUBROUTINE ASSE_KSTAB_GL(K,R)
  7. C ANTOINE LEGAY LE 15/02/2001
  8. C TOURNE KSTAB VERS LE REPERE GLOBAL
  9. * REAL *8 K(24,24),R(3,3),KRT(24,24),TMP(24,24),KR(24,24)
  10. * INTEGER I,J,KK
  11. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  12.  
  13. DIMENSION XK(24,24),R(3,3),XKRT(24,24),TMP(24,24),XKR(24,24)
  14. CALL ZDANUL(XKRT,576)
  15. CALL ZDANUL(XKR,576)
  16. DO KK=1,8
  17. DO J=1,3
  18. DO I=1,3
  19. CCCCCCC ATTENTION< INVERSE LE 3 NOVEMBRE 2001
  20. XKR((KK-1)*3+I,(KK-1)*3+J) = R(I,J)
  21. XKRT((KK-1)*3+I,(KK-1)*3+J) = R(J,I)
  22. ENDDO
  23. ENDDO
  24. ENDDO
  25. C WRITE(6,*)'****************** KR ***************'
  26. C DO I=1,24
  27. C WRITE(6,*)(KR(I,J),J=1,24)
  28. C ENDDO
  29. C WRITE(6,*)'*********************************'
  30. * CALL MULMAT(24,24,24,XK,XKR,TMP)
  31. * CALL MULMAT(24,24,24,XKRT,TMP,XKR)
  32. call mulmat ( tmp,xk,xkr,24,24,24)
  33. call mulmat(xkr,xkrt,tmp,24,24,24)
  34. CALL SHIFTD(XKR,XK,576)
  35.  
  36. C
  37. RETURN
  38. END
  39.  
  40.  
  41.  

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