Télécharger la plaquette Cast3M


fr en

Télécharger triang.eso

Retour à la liste

Numérotation des lignes : 

triang

  1. C TRIANG    SOURCE    GOUNAND   24/11/12    21:15:10     12076          
  2.       SUBROUTINE TRIANG(KRIGI,PREC,istab)
  3. C=======================================================================
  4. C   ASSEMBLE LES PETITES MATRICES RIGIDITE ET LES MET SOUS LA FORME
  5. C   Lt.D.L.
  6. C   IL LE POINTEUR DE LA MATRICE RESULTANTE DANS ICHOLE
  7. C=======================================================================
  8. C
  9.       IMPLICIT INTEGER(I-N)
  10.       IMPLICIT REAL*8 (A-H,O-Z)
  11. -INC SMRIGID
  12. -INC PPARAM
  13. -INC CCOPTIO
  14. -INC SMMATRI
  15.       IF(IIMPI.EQ.1)THEN
  16.          CALL GIBTEM(XKT)
  17.          INTERR(1)=INT(XKT)
  18.          CALL ERREUR(-259)
  19.       ENDIF
  20.       IF(IIMPI.EQ.1)WRITE(IOIMP,10)
  21.  10   FORMAT(' L''IMPRESSION PRECEDENTE EST  AVANT ASSEM1 ')
  22. *
  23.       CALL ASSEM1(KRIGI,MMATRX,INUINX,ITOPOX,IMINIX,IPOX,INCTRX,IITOPX)
  24.       IF(IERR.NE.0) RETURN
  25.       IF(IIMPI.EQ.1)THEN
  26.          CALL GIBTEM(XKT)
  27.          INTERR(1)=INT(XKT)
  28.          CALL ERREUR(-259)
  29.       ENDIF
  30.       IF(IIMPI.EQ.1)WRITE(IOIMP,11)
  31.    11  FORMAT(' L''IMPRESSION PRECEDENTE EST  AVANT ASSEM2')
  32.        CALL ASSEM2(KRIGI,ITOPOX,INUINX,IMINIX,MMATRX,IPOX,INCTRX,IITOPX)
  33.       IF(IERR.NE.0) RETURN
  34. *
  35.       IF(IIMPI.EQ.1)THEN
  36.          CALL GIBTEM(XKT)
  37.          INTERR(1)=INT(XKT)
  38.          CALL ERREUR(-259)
  39.       ENDIF
  40.       IF(IIMPI.EQ.1)WRITE(IOIMP,12)
  41.   12  FORMAT(' L''IMPRESSION PRECEDENTE EST AVANT CHOLE ')
  42.       IF(IERR.NE.0) GO TO 5000
  43. *  sauvegarde de la matrice assemblee sous forme creuse
  44.       call graco11(mmatrx,ilicre,1)
  45.       MRIGID=KRIGI
  46.       SEGACT MRIGID*MOD
  47.       nbnnma=0
  48.       nligra=0
  49.       xmatri=0
  50.       CALL CHOLE(MMATRX,PREC,istab,nbnnma,nligra,xmatri)
  51.       mmatri=mmatrx
  52.       segact mmatri
  53.       IF(IERR.NE.0) RETURN
  54.       IF(IIMPI.EQ.1)THEN
  55.          CALL GIBTEM(XKT)
  56.          INTERR(1)=INT(XKT)
  57.          CALL ERREUR(-259)
  58.       ENDIF
  59.       IF(IERR.NE.0) GO TO 5000
  60.       IF(IIMPI.EQ.1) WRITE(IOIMP,13)
  61.    13 FORMAT(' L''IMPRESSION PRECEDENTE EST APRES CHOLE')
  62.       MRIGID=KRIGI
  63.       SEGACT MRIGID*MOD
  64.       ICHOLE=MMATRX
  65.       SEGDES MRIGID
  66.  5000 CONTINUE
  67.       RETURN
  68.       END
  69.  
  70.  

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