Esope source files

   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 : *                                                                               
  16 : *     CALL SESPA1 ( IPLVEC, IPLVAL, IPRIGI, IPMASS )
  17 : *                                                                               
  18 : * PARAMETRES:   (E)=ENTREE   (S)=SORTIE                                         
  19 : * -----------                                                                   
  20 : *
  21 : *     IPLVAL  ENTIER    (S)  POINTEUR DE L'OBJET 'LISTREEL' CONTENANT
  22 : *                            LA SUITE DE VALEURS PROPRES CALCULEES
  23 : *                                                                               
  24 : *     IPLVEC  ENTIER    (E)  POINTEUR DE L'OBJET 'LISTCHPO' CONTENANT
  25 : *                            LA SUITE DE 'CHPOINT' X(i)
  26 : *     IPLVEC  ENTIER    (S)  POINTEUR DE L'OBJET 'LISTCHPO' CONTENANT
  27 : *                            LA SUITE DE 'CHPOINT' X(i+1)
  28 : *
  29 : *     IPRIGI  ENTIER    (E)  POINTEUR SUR L'OBJET 'RIGIDITE' K
  30 : *     IPMASS  ENTIER    (E)  POINTEUR SUR L'OBJET 'RIGIDITE' M
  31 : *
  32 : *                                                                               
  33 : * AUTEUR, DATE DE CREATION:                                                     
  34 : * -------------------------                                                     
  35 : *                                                                               
  36 : *   A.M. JOLIVALT, W. PASILLAS  12 / 07 / 94.  ( ESOPE )                                      
  37 : *
  38 : ************************************************************
  39 : 
  40 :       SUBROUTINE SESPA0 ( IPLVEC, IPLVAL, IPRIGI, IPMASS )
  41 :           
  42 :       IMPLICIT INTEGER(I-N)
  43 :       IMPLICIT REAL*8 (A-H,O-Z)
  44 : 
  45 : -INC  CCOPTIO

  46 : 
  47 : ******
  48 : *     -- ARGUMENTS --
  49 : ***
  50 :       INTEGER  IPLVEC, IPLVAL, IPRIGI, IPMASS 
  51 : 
  52 : ******
  53 : *     -- VARIABLES LOCALES --
  54 : ***
  55 :       INTEGER IPRIG1, IPMAS1, IPMPHI
  56 : 
  57 : 
  58 : ************************************************************
  59 : *            PROCEDURE DE RAYLEIGH-RITZ                    *
  60 : ************************************************************
  61 : 
  62 : ******
  63 : *     -- CALCUL DU SOUS-ESPACE SUIVANT: --
  64 : ***
  64 :       CALL SESPA1 ( IPLVEC, IPRIGI, IPMASS )
  65 :       IF ( IERR .NE. 0 ) RETURN
  66 : 
  67 : ******
  68 : *     -- PROJECTION DE K ET DE M SUR CE SOUS-ESPACE --
  69 : ***
  69 :       CALL SESPA2 ( IPLVEC, IPRIGI, IPRIG1 )
  70 :       IF ( IERR .NE. 0 ) RETURN
  70 :       CALL SESPA2 ( IPLVEC, IPMASS, IPMAS1 )
  71 :       IF ( IERR .NE. 0 ) RETURN
  72 : 
  73 : ******
  74 : *     -- RESOLUTION DU PROBLEME REDUIT PAR LA METHODE DE JACOBI --
  75 : ***
  75 :       CALL SESPA3 ( IPRIG1, IPMAS1, IPMPHI, IPLVAL )
  76 :       IF ( IERR .NE. 0 ) RETURN
  77 : 
  78 : ******
  79 : *     -- CALCUL D'UNE APPROX. DES VECTEURS PROPRES A PARTIR
  80 : *     DE LEURS PROJECTIONS SUR LE SOUS-ESPACE --
  81 : ***
  81 :       CALL SESPA4 ( IPLVEC, IPMPHI )
  82 :       IF ( IERR .NE. 0 ) RETURN
  83 : 
  84 :       RETURN
  85 :       END
  86 :