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

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