Télécharger la plaquette Cast3M


fr en

Télécharger deviat.eso

Retour à la liste

Numérotation des lignes : 

deviat

  1. C DEVIAT    SOURCE    BR232186  12/08/01    21:15:22     7456
  2. C DEVIAT    SOURCE    BR      06/01/10    21:15:29     5284
  3.       SUBROUTINE DEVIAT(A,B)
  4. C======================================================================
  5. C Routine permettant de calculer le deviateur de A
  6. C
  7. C     OBJET
  8. C     -----
  9. C     CALCUL DE B = DEVIAT(A) = A - (1.0D0/3.0D0) * TRACE(A) * UNIT
  10. C
  11. C     ENTREES
  12. C     -------
  13. C     A(3,3) = MATRICE SYMETRIQUE
  14. C
  15. C     SORTIES
  16. C     -------
  17. C     B      = DEVIAT(A)
  18. C
  19. C===============================================================
  20.       IMPLICIT INTEGER(I-N)
  21.       IMPLICIT REAL*8 (A-H,O-Z)
  22.  
  23.       REAL*8 A(3,3),UNIT(3,3),B(3,3)
  24.  
  25.       TRA = A(1,1)+A(2,2)+A(3,3)
  26.  
  27.       DO I=1,3
  28.          DO J=1,3
  29.  
  30.             UNIT(I,J) = 0.0D0
  31.  
  32.             IF (I.EQ.J) THEN
  33.                UNIT(I,J) = 1.0D0
  34.             ENDIF
  35.  
  36.          ENDDO
  37.       ENDDO
  38.  
  39.       DO I=1,3
  40.          DO J=1,3
  41.  
  42.             B(I,J) = 0.0D0
  43.             B(I,J) = A(I,J) - (1.0D0/3.0D0)*TRA*UNIT(I,J)
  44.  
  45.          ENDDO
  46.       ENDDO
  47.  
  48.  
  49.       RETURN
  50.       END
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  

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