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    PV        22/10/05    21:15:02     11474          
  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)=XKT
  18.        CALL ERREUR(-259)
  19.       ENDIF
  20.       IF(IIMPI.EQ.1)WRITE(IOIMP,10)
  21.   10  FORMAT(' L''IMPRESSION PRECDENTE EST  AVANT ASSEM1 ')
  22.       CALL ASSEM1(KRIGI,MMATRX,INUINX,ITOPOX,IMINIX,IPOX,INCTRX,IITOPX)
  23.       IF(IERR.NE.0) RETURN
  24.       IF(IIMPI.EQ.1)THEN
  25.        CALL GIBTEM(XKT)
  26.        INTERR(1)=XKT
  27.        CALL ERREUR(-259)
  28.       ENDIF
  29.       IF(IIMPI.EQ.1)WRITE(IOIMP,11)
  30.    11  FORMAT(' L''IMPRESSION PRECEDENTE EST  AVANT ASSEM2')
  31.       CALL ASSEM2(KRIGI,ITOPOX,INUINX,IMINIX,MMATRX,IPOX,INCTRX,IITOPX)
  32.       IF(IERR.NE.0) RETURN
  33.       IF(IIMPI.EQ.1)THEN
  34.        CALL GIBTEM(XKT)
  35.        INTERR(1)=XKT
  36.        CALL ERREUR(-259)
  37.       ENDIF
  38.       IF(IIMPI.EQ.1)WRITE(IOIMP,12)
  39.   12  FORMAT(' L''IMPRESSION PRECEDENTE EST AVANT CHOLE ')
  40.       IF(IERR.NE.0) GO TO 5000
  41. *  sauvegarde de la matrice assemblee sous forme creuse
  42.       call graco11(mmatrx,ilicre,1)
  43.       MRIGID=KRIGI
  44.       SEGACT MRIGID*MOD
  45.       nbnnma=0
  46.       nligra=0
  47.       xmatri=0
  48.       CALL CHOLE(MMATRX,PREC,istab,nbnnma,nligra,xmatri)
  49.       mmatri=mmatrx
  50.       segact mmatri
  51.       IF(IERR.NE.0) RETURN
  52.       IF(IIMPI.EQ.1)THEN
  53.        CALL GIBTEM(XKT)
  54.        INTERR(1)=XKT
  55.        CALL ERREUR(-259)
  56.       ENDIF
  57.       IF(IERR.NE.0) GO TO 5000
  58.       IF(IIMPI.EQ.1) WRITE(IOIMP,13)
  59.    13 FORMAT(' L''IMPRESSION PRECEDENTE EST APRES CHOLE')
  60.       MRIGID=KRIGI
  61.       SEGACT MRIGID*MOD
  62.       ICHOLE=MMATRX
  63.       SEGDES MRIGID
  64.  5000 CONTINUE
  65.       RETURN
  66.       END
  67.  
  68.  
  69.  
  70.  
  71.  
  72.  
  73.  
  74.  
  75.  
  76.  
  77.  
  78.  

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