Télécharger la plaquette Cast3M


fr en

Télécharger extipd.eso

Retour à la liste

Numérotation des lignes : 

extipd

  1. C EXTIPD    SOURCE    PV        20/09/26    21:16:51     10724          
  2.       SUBROUTINE EXTIPD
  3.       IMPLICIT REAL*8 (A-H,O-Z)
  4.       IMPLICIT INTEGER (I-N)
  5. C***********************************************************************
  6. C NOM         : EXTIPD
  7. C DESCRIPTION : Extrait les noms d'inconnues primales ou duales
  8. C               d'un MATRIK, on réduit à CH*4 pour des raisons
  9. C               de compatibilité
  10. C
  11. C
  12. C LANGAGE     : ESOPE
  13. C AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  14. C               mél : gounand@semt2.smts.cea.fr
  15. C***********************************************************************
  16. C APPELES          :
  17. C APPELES (E/S)    :
  18. C APPELES (BLAS)   :
  19. C APPELES (CALCUL) :
  20. C APPELE PAR       :
  21. C***********************************************************************
  22. C SYNTAXE GIBIANE    :
  23. C    LMOT1  = 'KOPS'  'EXTRINPR'  MATRIK1  ;
  24. C    LMOT1  = 'KOPS'  'EXTRINDU'  MATRIK1  ;
  25. C ENTREES            :
  26. C ENTREES/SORTIES    :
  27. C SORTIES            :
  28. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  29. C***********************************************************************
  30. C VERSION    : v1, 10/05/2006, version initiale
  31. C HISTORIQUE : v1, 10/05/2006, création
  32. C HISTORIQUE :
  33. C HISTORIQUE :
  34. C***********************************************************************
  35. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  36. C en cas de modification de ce sous-programme afin de faciliter
  37. C la maintenance !
  38. C***********************************************************************
  39.  
  40. -INC PPARAM
  41. -INC CCOPTIO
  43. -INC SMLMOTS
  44.       POINTEUR LINC.MLMOTS
  45.       POINTEUR LINC2.MLMOTS
  46. *
  47.       INTEGER IMPR,IRET
  48. *
  49. * Executable statements
  50. *
  51.       IMPR=0
  52.       IF (IMPR.GT.2) WRITE(IOIMP,*) 'Entrée dans extipd.eso'
  53.       CALL LIROBJ('MATRIK  ',MATRIK,1,IRETOU)
  54.       IF (IERR.NE.0) RETURN
  55. * D'abord les primales puis les duales
  56. * On renverse l'ordre car ecrobj écrit sur une pile
  57.       DO IPRIDU=2,1,-1
  58. *
  59. * Dim max de LINC
  60. *
  61.          JGM=0
  62.          SEGACT MATRIK
  63.          NMAT=IRIGEL(/2)
  64.          DO IMAT=1,NMAT
  65.             IMATRI=IRIGEL(4,IMAT)
  66.             SEGACT IMATRI
  67.             JGM=JGM+LIZAFM(/2)
  68.          ENDDO
  69. *
  70. * Remplissage de LINC
  71. *
  72.          JGN=4
  73.          IGM=0
  74.          SEGINI LINC
  75.          DO IMAT=1,NMAT
  76.             IMATRI=IRIGEL(4,IMAT)
  77.             NBME=LIZAFM(/2)
  78.             DO IBME=1,NBME
  79.                IGM=IGM+1
  80.                IF (IPRIDU.EQ.1) THEN
  81.                   LINC.MOTS(IGM)=LISPRI(IBME)
  82.                ELSEIF (IPRIDU.EQ.2) THEN
  83.                   LINC.MOTS(IGM)=LISDUA(IBME)
  84.                ELSE
  85.                   GOTO 9999
  86.                ENDIF
  87.             ENDDO
  88.             SEGDES IMATRI
  89.          ENDDO
  90.          SEGDES MATRIK
  91. *
  92. * Enlever les doublons dans LINC
  93. *
  94.          SEGINI,LINC2=LINC
  95.          CALL CUNIQ(LINC2.MOTS,LINC2.MOTS(/1),LINC2.MOTS(/2),
  96.      $        LINC.MOTS,NIUNIQ,
  97.      $        IMPR,IRET)
  98.          IF (IRET.NE.0) GOTO 9999
  99.          JGN=4
  100.          JGM=NIUNIQ
  101.          SEGADJ,LINC
  102.          SEGSUP LINC2
  103.          SEGDES LINC
  104.          CALL ECROBJ('LISTMOTS',LINC)
  105.       ENDDO
  106. *
  107. * Normal termination
  108. *
  109. *      IRET=0
  110.       RETURN
  111. *
  112. * Format handling
  113. *
  114. *
  115. * Error handling
  116. *
  117.  9999 CONTINUE
  118. *      IRET=1
  119.       WRITE(IOIMP,*) 'An error was detected in subroutine extipd'
  120.       CALL ERREUR(5)
  121.       RETURN
  122. *
  123. * End of subroutine EXTIPD
  124. *
  125.       END
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  

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