Sources Esope | Cast3M

ldmt1
C LDMT1     SOURCE    PV        20/09/26    21:18:31     10724          
      SUBROUTINE LDMT1(KRIGI,PREC)
C=======================================================================
C   ASSEMBLE LES PETITES MATRICES de RIGIDITE ET LES MET SOUS LA FORME
C        t
C   L.D.M
C   IL LE POINTEUR DE LA MATRICE RESULTANTE DANS ICHOLE (segment MRIGID)
C
C   Cette subroutine est équivalente à TRIANG dans le cas de
C   l'inversion des matrices symétrique
C
C   Appelée par : LDMT
C
C   Auteur : Michel BULIK
C
C   Date : Printemps '95
C
C   Langage : ESOPE + FORTRAN77
C
C=======================================================================
C
      IMPLICIT INTEGER(I-N)
      REAL*8 XKT,PREC
-INC SMRIGID
-INC SMMATRI
 
-INC PPARAM
-INC CCOPTIO
 
C ... Ces variables ont pour but, de diriger le comportement de LDMT2 ...
C       TRSUP  - TRiangle SUPérieur
C       MENAGE - évident
C       LDIAG  - initialisation et remplissage de MDIAG et MDNOR demandés
 
      LOGICAL TRSUP,MENAGE,LDIAG
 
 
      IF(IIMPI.EQ.1)THEN
         CALL GIBTEM(XKT)
         INTERR(1)=XKT
         CALL ERREUR(-259)
      ENDIF
      IF(IIMPI.EQ.1)WRITE(IOIMP,10)
  10  FORMAT(' L''IMPRESSION PRECDENTE EST  AVANT ASSEM1 ')
 
C ... MMATRI est initialisé dans ASSEM1 et renvoyé en tant que résultat
C     dans la variable MMATRX, il est désactivé à la sortie ...
      CALL ASNS1(KRIGI,MMATRX,INUINX,ITOPOX,IMINIX,IPOX,INCTRX,INCTRZ,
     &           IITOPX,ITOPOD,IITOPD,IPODD)
      IF(IERR.NE.0) RETURN
 
      IF(IIMPI.EQ.1) THEN
         CALL GIBTEM(XKT)
         INTERR(1)=XKT
         CALL ERREUR(-259)
      ENDIF
      IF(IIMPI.EQ.1)WRITE(IOIMP,11)
   11 FORMAT(' L''IMPRESSION PRECEDENTE EST  AVANT LDMT2')
 
C ... On initialise IJMAX ici et non dans LDMT2, car celui-ci est
C     appelé deux fois ...
      MMATRI=MMATRX
      SEGACT,MMATRI*MOD
      IJMAX=0
      SEGDES,MMATRI
 
      TRSUP =.FALSE.
      MENAGE=.FALSE.
      LDIAG =.TRUE.
      njtot=0
*      write(6,*) ' premier appel'
      CALL LDMT2(KRIGI,ITOPOX,INUINX,IMINIX,MMATRX,IPOX,INCTRX,INCTRZ,
     &           IITOPX,TRSUP,MENAGE,LDIAG,IITOPD,ITOPOD,IPODD,njtot,1)
      IF(IERR.NE.0) RETURN
      TRSUP =.TRUE.
      MENAGE=.TRUE.
      LDIAG =.FALSE.
*      write(6,*) ' deucxieme appel'
      CALL LDMT2(KRIGI,ITOPOX,INUINX,IMINIX,MMATRX,IPOX,INCTRX,INCTRZ,
     &           IITOPX,TRSUP,MENAGE,LDIAG,IITOPD,ITOPOD,IPODD,njtot,1)
      IF(IERR.NE.0) RETURN
 
      IF(IIMPI.EQ.1)THEN
         CALL GIBTEM(XKT)
         INTERR(1)=XKT
         CALL ERREUR(-259)
      ENDIF
      IF(IIMPI.EQ.1) WRITE(IOIMP,12)
  12  FORMAT(' L''IMPRESSION PRECEDENTE EST AVANT LDMT3')
 
      IF(IERR.NE.0) GO TO 5000
 
      CALL LDMT3(MMATRX,PREC)
      IF(IERR.NE.0) RETURN
 
      IF(IIMPI.EQ.1)THEN
         CALL GIBTEM(XKT)
         INTERR(1)=XKT
         CALL ERREUR(-259)
      ENDIF
      IF(IERR.NE.0) GO TO 5000
 
      IF(IIMPI.EQ.1) WRITE(IOIMP,13)
   13 FORMAT(' L''IMPRESSION PRECEDENTE EST APRES LDMT3')
 
      MRIGID=KRIGI
      SEGACT MRIGID*MOD
      ICHOLE=MMATRX
      SEGDES MRIGID
 
 5000 CONTINUE
      RETURN
      END