Télécharger la plaquette Cast3M


fr en

Télécharger sespc2.eso

Retour à la liste

Numérotation des lignes : 

sespc2

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

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