Esope source files

   1 : C SESPA2    SOURCE    CHAT      05/01/13    03:14:36     5004           
   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 : *
  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 SESPA2 ( IPLVEC, IPRIGI, IPRIG1 )
  33 : 
  34 :       IMPLICIT INTEGER(I-N)
  35 :       IMPLICIT REAL*8 (A-H,O-Z)
  36 : 
  37 : -INC  CCOPTIO

  38 : -INC  SMLCHPO

  39 : -INC  SMLMOTS

  40 : -INC  SMCHPOI

  41 : 
  42 :       SEGMENT ,MATRIX
  43 :         REAL*8 A(N,N)
  44 :       ENDSEGMENT
  45 : 
  46 : ******
  47 : *     -- ARGUMENTS --
  48 : ***
  49 :       POINTEUR IPLVEC.MLCHPO
  50 :       POINTEUR IPRIG1.MATRIX
  51 :       INTEGER  IPRIGI
  52 : 
  53 : ******
  54 : *     -- VARIABLES LOCALES --
  55 : ***
  56 :       INTEGER ILDIM, IB100, IPCHP1, IPCHP2
  57 :       REAL*8  XITMXJ
  58 : 
  59 : 
  60 : ******
  61 : *     -- K' = Xt K X --
  62 : ***
  63 :       SEGACT ,IPLVEC
  64 :       ILDIM = IPLVEC.ICHPOI(/1)
  65 :       N = ILDIM
  66 :       SEGINI ,IPRIG1
  67 : 
  68 :       DO 100 IB100 = 1, ILDIM
  69 :          IPCHP1 = IPLVEC.ICHPOI( IB100 )
  69 :          CALL MUCPRI(IPCHP1,IPRIGI,IPCHP3)
  70 :          DO 200 IB200 = 1, ILDIM
  71 :             IPCHP2 = IPLVEC.ICHPOI( IB200 )
  72 :             IF(IB200.EQ.1.AND.IB100.EQ.1)
  72 :      >          CALL CORRSP(iprigi,IPCHP2,IPCHP3,MLMOT1,MLMOT2)
  72 :             CALL XTY1 (IPCHP2,IPCHP3,MLMOT1,MLMOT2,  XITMXJ)
  73 :             IF ( IERR .NE. 0 ) RETURN
  74 :             IPRIG1.A(IB100,IB200) = XITMXJ
  75 :  200     CONTINUE
  75 :         CALL DTCHPO(IPCHP3)
  76 :  100  CONTINUE
  77 :       SEGSUP MLMOT1,MLMOT2
  78 :       SEGDES ,IPLVEC, IPRIG1
  79 : 
  80 :       RETURN
  81 :       END
  82 :  
  83 :