Télécharger la plaquette Cast3M


fr en

Télécharger ecimat.eso

Retour à la liste

Numérotation des lignes : 

ecimat

  1. C ECIMAT    SOURCE    PV        20/09/26    21:16:32     10724          
  2.       SUBROUTINE ECIMAT(IMATRI,NIVIMP)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8 (A-H,O-Z)
  5. C***********************************************************************
  6. C NOM         : ECIMAT
  7. C DESCRIPTION : Impression d'un segment de type IMATRI.
  8. C               Ce segment contient des groupes de
  9. C               matrices élémentaires non assemblées.
  10. C               cf. include SMMATRIK
  11. C
  12. C LANGAGE     : ESOPE
  13. C AUTEUR      : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/TTMF)
  14. C               mél : gounand@semt2.smts.cea.fr
  15. C***********************************************************************
  16. C APPELES : -
  17. C***********************************************************************
  18. C ENTREES            : IMATRI, NIVIMP
  19. C ENTREES/SORTIES    : -
  20. C SORTIES            : -
  21. C CODE RETOUR (IRET) : -
  22. C    IMATRI : segment de type IMATRI (include SMMATRIK)
  23. C    NIVIMP : niveau d'impression. Suivant sa valeur, on obtient :
  24. C  Convention (probablement non totalement respectée) :
  25. C  ----------   0 : presque rien (numéro de pointeur)
  26. C               1 : affichage du chapeau IMATRI
  27. C               2 : affichage des données concernant les objets
  28. C                   pointés par IMATRI (s'il y en a)
  29. C               3 : affichage du contenu des objets vectoriels
  30. C               4 : affichage du contenu des objets matriciels
  31. C On ne change pas l'état (actif ou inactif) du segment IMATRI.
  32. C***********************************************************************
  33. C VERSION    : v1, 01/04/98, version initiale
  34. C HISTORIQUE : v1, 01/04/98, création
  35. C HISTORIQUE : 29/10/98, modif. l'état du segment reste inchangé
  36. C                        en sortie
  37. C HISTORIQUE :
  38. C***********************************************************************
  39. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  40. C en cas de modification de ce sous-programme afin de faciliter
  41. C la maintenance !
  42. C***********************************************************************
  43. -INC PPARAM
  44. -INC CCOPTIO
  46. * Variable d'état du segment IMATRK
  47.       INTEGER IMAETA
  48. *
  49. * Executable statements
  50. *
  51.       IF (IMATRI.EQ.0) THEN
  52.          WRITE(IOIMP,*) 'Nil pointer transmitted to ecimat'
  53.          GOTO 9999
  54.       ENDIF
  55.       CALL OOOETA(IMATRI,IMAETA,IMOD)
  56.       IF (IMAETA.NE.1) SEGACT IMATRI
  57.       WRITE(IOIMP,2001) 'Segment IMATRI de pointeur ',IMATRI
  58. *
  59.       NBME=LIZAFM(/2)
  60.       NBSOUS=LIZAFM(/1)
  61.       WRITE(IOIMP,2003) 'NBME  =',NBME,'(nb groupe mat. élém.)'
  62.       WRITE(IOIMP,2003) 'NBSOUS=',NBSOUS,'(nb sous-obj.)'
  63.       WRITE(IOIMP,2004) 'KSPGP=',KSPGP,'KSPGD=',KSPGD
  64.       WRITE(IOIMP,1999) 'Liste inconnues primales et duales'
  65.       WRITE(IOIMP,2999)  'LISPRI(NBME)=',(LISPRI(N),N=1,NBME)
  66.       WRITE(IOIMP,2999)  'LISDUA(NBME)=',(LISDUA(N),N=1,NBME)
  67.       WRITE(IOIMP,1999) 'Matrices élémentaires non assemblées'
  68.       WRITE(IOIMP,1999) 'Tableau LIZAFM(NBSOUS,NBME)'
  69.       DO 4 J=1,NBME
  70.          DO 5 K=1,NBSOUS
  71.             IZAFM=LIZAFM(K,J)
  72.             WRITE(IOIMP,2007) 'LIZAFM(',K,',',J,')=',IZAFM,
  73.      $           'pointeur IZAFM'
  74.             IF ((NIVIMP.GT.2).AND.(IZAFM.NE.0)) THEN
  75. * Affichage des IZAFM
  76.                SEGACT IZAFM
  77.                NBEL=AM(/1)
  78.                NP=AM(/2)
  79.                MP=AM(/3)
  80.                WRITE(IOIMP,2008) 'NBEL=',NBEL,'(nb élém)'
  81.                WRITE(IOIMP,2009) 'NP  =',NP,'MP  =',MP,
  82.      $              '(nb points)'
  83.                DO 6 L=1,NBEL
  84.                   WRITE(IOIMP,2010) 'Matrice pour l''élément',L
  85.                   DO 7 N=1,NP
  86. * Affichage des AM
  87.                      WRITE(IOIMP,2011) (AM(L,N,M),M=1,MP)
  88.  7                CONTINUE
  89.  6             CONTINUE
  90.                SEGDES IZAFM
  91.             ENDIF
  92.             WRITE(IOIMP,1998) 'End segment IZAFM'
  93.             WRITE(IOIMP,1998) '---'
  94.  5       CONTINUE
  95.  4    CONTINUE
  96.       IF (IMAETA.NE.1) SEGDES IMATRI
  97. *
  98. * Normal termination
  99. *
  100.       RETURN
  101. *
  102. * Format handling
  103. *
  104.  1998 FORMAT (2X,A)
  105.  1999 FORMAT (A)
  106.  2001 FORMAT (A,1X,I6)
  107.  2003 FORMAT (2X,A,I2,1X,A)
  108.  2004 FORMAT (2X,A,I6,4X,A,I6)
  109.  2007 FORMAT (2X,A,I2,A,I2,A,I6,2X,A)
  110.  2008 FORMAT (4X,A,I6,1X,A)
  111.  2009 FORMAT (4X,A,I6,2X,A,I6,1X,A)
  112.  2010 FORMAT (6X,A,1X,I6)
  113.  2011 FORMAT (4(1X,1PE18.10))
  114.  2999 FORMAT (2X,A,100(1X,A8))
  115. *
  116. * Error handling
  117. *
  118.  9999 CONTINUE
  119.       WRITE(IOIMP,*) 'An error was detected in subroutine ecimat'
  120.       RETURN
  121. *
  122. * End of subroutine ECIMAT
  123. *
  124.       END
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  

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