Télécharger vpmax3.eso

Retour à la liste

Numérotation des lignes :

  1. C VPMAX3 SOURCE CHAT 05/01/13 04:09:49 5004
  2. SUBROUTINE VPMAX3 (SIG,SIGMAX,VECP)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. -INC CCOPTIO
  6. C
  7. DIMENSION SIG(*),VECP(*),A(3,3),D(3),S(3,3)
  8. C
  9. C CETTE ROUTINE SERT A CALCULER LA VALEUR MAX DE SIG ET LE
  10. C VECTEUR PROPRE QUI CORRESPOND
  11. C ENTREES :
  12. C SIG : TENSEUR DE CONTRAINTES
  13. C IIMPI : INDICE D IMPRESSION INTERNE
  14. C SORTIES :
  15. C SIGMAX : CONTRAINTE MAX
  16. C VECP : VECTEUR PROPRE CORRESPONDANT
  17. C
  18. C INITIALISATION
  19. C
  20. A(1,1)=SIG(1)
  21. A(2,2)=SIG(2)
  22. A(3,3)=SIG(3)
  23. A(1,2)=SIG(4)
  24. A(2,1)=SIG(4)
  25. A(2,3)=SIG(5)
  26. A(3,2)=SIG(5)
  27. A(3,1)=SIG(6)
  28. A(1,3)=SIG(6)
  29. C
  30. CALL JACOB3(A,3,D,S)
  31. C
  32. C CALCUL DE LA VALEUR PROPRE MAXI ET SON VECTEUR PROPRE
  33. C
  34. VECP(1)=S(1,1)
  35. VECP(2)=S(2,1)
  36. VECP(3)=S(3,1)
  37. SIGMAX=D(1)
  38. C
  39. DO 1 I=2,3
  40. IF(SIGMAX.GE.D(I)) GO TO 1
  41. SIGMAX=D(I)
  42. DO 2 J=1,3
  43. VECP(J)=S(J,I)
  44. 2 CONTINUE
  45. 1 CONTINUE
  46. C
  47. DO 3 I=1,3
  48. IF(ABS(VECP(I)).LT.1.D-9) VECP(I)=0.D0
  49. 3 CONTINUE
  50. C
  51. IF(IIMPI.EQ.9)
  52. . WRITE(IOIMP,1000) (SIG(I),I=1,6),SIGMAX,(VECP(I),I=1,3),
  53. . ((A(I,J),J=1,3),I=1,3),(D(I),I=1,3),
  54. . ((S(I,J),J=1,3),I=1,3)
  55. 1000 FORMAT(1X,'SIG =',6(3X,D12.5),/,1X,'SIGMAX=',D12.5,/,
  56. .1X,'VECP=',3(3X,D12.5),/,1X,'AAA =',3(3(3X,D12.5),/),
  57. .1X,'DDD =',3(3X,D12.5),/,1X,'SSS =',3(3(3X,D12.5),/))
  58. C
  59. RETURN
  60. END
  61.  
  62.  

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