Télécharger la plaquette Cast3M


fr en

Télécharger sespa2.eso

Retour à la liste

Numérotation des lignes : 

sespa2

  1. C SESPA2    SOURCE    CB215821  20/11/25    13:39:41     10792          
  2. C SESPA2    SOURCE    WP        20/08/94
  3. C     SUBROUTINE SESPA2 ( IPLVEC, IPRIGI, IPRIG1 )
  4. ************************************************************************
  5. *
  6. *                                SESPA2
  7. *                             -----------
  8. *
  9. * FONCTION:
  10. * ---------
  11. *
  12. *     CALCULER LA RESTRICTION DE IPRIGI AU LE SOUS-ESPACE IPLVEC.
  13. *
  14. *
  15. * MODE D'APPEL:
  16. *
  17. *     CALL SESPA2 ( IPLVEC, IPRIGI, IPRIG1 )
  18. *
  19. * PARAMETRES:   (E)=ENTREE   (S)=SORTIE
  20. * -----------
  21. *
  22. *     IPLVEC  ENTIER    (E)  POINTEUR DE L'OBJET 'LISTCHPO' CONTENANT
  23. *                            LA SUITE DE 'CHPOINT' X(i)
  24. *
  25. *     IPRIGI  ENTIER    (E)  POINTEUR SUR L'OBJET 'RIGIDITE' K
  26. *     IPRIG1  ENTIER    (E)  POINTEUR SUR L'OBJET 'MATRIX'   K=Xt K X
  27. *
  28. *
  29. *   A.M. JOLIVALT,  W. PASILLAS     12 / 07 / 94.  ( ESOPE )
  30. *
  31. ************************************************************
  32.  
  33.       SUBROUTINE SESPA2 ( IPLVEC, IPRIGI, IPRIG1 )
  34.  
  35.       IMPLICIT INTEGER(I-N)
  36.       IMPLICIT REAL*8 (A-H,O-Z)
  37.  
  38. -INC PPARAM
  39. -INC CCOPTIO
  40. -INC  SMLCHPO
  41. -INC  SMLMOTS
  42. -INC  SMCHPOI
  43.  
  44.       SEGMENT ,MATRIX
  45.         REAL*8 A(N,N)
  46.       ENDSEGMENT
  47.  
  48. ******
  49. *     -- ARGUMENTS --
  50. ***
  51.       POINTEUR IPLVEC.MLCHPO
  52.       POINTEUR IPRIG1.MATRIX
  53.       INTEGER  IPRIGI
  54.  
  55. ******
  56. *     -- VARIABLES LOCALES --
  57. ***
  58.       INTEGER ILDIM, IB100, IPCHP1, IPCHP2
  59.       REAL*8  XITMXJ
  60.  
  61.  
  62. ******
  63. *     -- K' = Xt K X --
  64. ***
  65.       SEGACT ,IPLVEC
  66.       ILDIM = IPLVEC.ICHPOI(/1)
  67.       N = ILDIM
  68.       SEGINI ,IPRIG1
  69.  
  70.       DO 100 IB100 = 1, ILDIM
  71.          IPCHP1 = IPLVEC.ICHPOI( IB100 )
  72.          CALL MUCPRI(IPCHP1,IPRIGI,IPCHP3)
  73.          DO 200 IB200 = 1, ILDIM
  74.             IPCHP2 = IPLVEC.ICHPOI( IB200 )
  75.             IF(IB200.EQ.1.AND.IB100.EQ.1)
  76.      >          CALL CORRSP(iprigi,IPCHP2,IPCHP3,MLMOT1,MLMOT2)
  77.             CALL XTY1 (IPCHP2,IPCHP3,MLMOT1,MLMOT2,  XITMXJ)
  78.             IF ( IERR .NE. 0 ) RETURN
  79.             IPRIG1.A(IB100,IB200) = XITMXJ
  80.  200     CONTINUE
  81.         CALL DTCHPO(IPCHP3)
  82.  100  CONTINUE
  83.       SEGSUP MLMOT1,MLMOT2
  84.       SEGDES ,IPLVEC, IPRIG1
  85.  
  86.       RETURN
  87.       END
  88.  
  89.  
  90.  
  91.  

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