Télécharger la plaquette Cast3M


fr en

Télécharger sespa0.eso

Retour à la liste

Numérotation des lignes : 

sespa0

  1. C SESPA0    SOURCE    CHAT      05/01/13    03:14:27     5004
  2. C SESPA0    SOURCE    WP        23/08/94
  3. C     SUBROUTINE SESPA0 ( IPLVEC, IPLVAL, IPRIGI, IPMASS )
  4. ************************************************************************
  5. *
  6. *                                SESPA0
  7. *                             -----------
  8. *
  9. * FONCTION:
  10. * ---------
  11. *
  12. *     FAIRE UN PAS D'ITERATION DE SOUS-ESPACE.
  13. *
  14. *
  15. * MODE D'APPEL:
  16. *
  17. *     CALL SESPA1 ( IPLVEC, IPLVAL, IPRIGI, IPMASS )
  18. *
  19. * PARAMETRES:   (E)=ENTREE   (S)=SORTIE
  20. * -----------
  21. *
  22. *     IPLVAL  ENTIER    (S)  POINTEUR DE L'OBJET 'LISTREEL' CONTENANT
  23. *                            LA SUITE DE VALEURS PROPRES CALCULEES
  24. *
  25. *     IPLVEC  ENTIER    (E)  POINTEUR DE L'OBJET 'LISTCHPO' CONTENANT
  26. *                            LA SUITE DE 'CHPOINT' X(i)
  27. *     IPLVEC  ENTIER    (S)  POINTEUR DE L'OBJET 'LISTCHPO' CONTENANT
  28. *                            LA SUITE DE 'CHPOINT' X(i+1)
  29. *
  30. *     IPRIGI  ENTIER    (E)  POINTEUR SUR L'OBJET 'RIGIDITE' K
  31. *     IPMASS  ENTIER    (E)  POINTEUR SUR L'OBJET 'RIGIDITE' M
  32. *
  33. *
  34. * AUTEUR, DATE DE CREATION:
  35. * -------------------------
  36. *
  37. *   A.M. JOLIVALT, W. PASILLAS  12 / 07 / 94.  ( ESOPE )
  38. *
  39. ************************************************************
  40.  
  41.       SUBROUTINE SESPA0 ( IPLVEC, IPLVAL, IPRIGI, IPMASS )
  42.  
  43.       IMPLICIT INTEGER(I-N)
  44.       IMPLICIT REAL*8 (A-H,O-Z)
  45.  
  46. -INC PPARAM
  47. -INC CCOPTIO
  48.  
  49. ******
  50. *     -- ARGUMENTS --
  51. ***
  52.       INTEGER  IPLVEC, IPLVAL, IPRIGI, IPMASS
  53.  
  54. ******
  55. *     -- VARIABLES LOCALES --
  56. ***
  57.       INTEGER IPRIG1, IPMAS1, IPMPHI
  58.  
  59.  
  60. ************************************************************
  61. *            PROCEDURE DE RAYLEIGH-RITZ                    *
  62. ************************************************************
  63.  
  64. ******
  65. *     -- CALCUL DU SOUS-ESPACE SUIVANT: --
  66. ***
  67.       CALL SESPA1 ( IPLVEC, IPRIGI, IPMASS )
  68.       IF ( IERR .NE. 0 ) RETURN
  69.  
  70. ******
  71. *     -- PROJECTION DE K ET DE M SUR CE SOUS-ESPACE --
  72. ***
  73.       CALL SESPA2 ( IPLVEC, IPRIGI, IPRIG1 )
  74.       IF ( IERR .NE. 0 ) RETURN
  75.       CALL SESPA2 ( IPLVEC, IPMASS, IPMAS1 )
  76.       IF ( IERR .NE. 0 ) RETURN
  77.  
  78. ******
  79. *     -- RESOLUTION DU PROBLEME REDUIT PAR LA METHODE DE JACOBI --
  80. ***
  81.       CALL SESPA3 ( IPRIG1, IPMAS1, IPMPHI, IPLVAL )
  82.       IF ( IERR .NE. 0 ) RETURN
  83.  
  84. ******
  85. *     -- CALCUL D'UNE APPROX. DES VECTEURS PROPRES A PARTIR
  86. *     DE LEURS PROJECTIONS SUR LE SOUS-ESPACE --
  87. ***
  88.       CALL SESPA4 ( IPLVEC, IPMPHI )
  89.       IF ( IERR .NE. 0 ) RETURN
  90.  
  91.       RETURN
  92.       END
  93.  
  94.  

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