Télécharger la plaquette Cast3M


fr en

Télécharger inctyp.eso

Retour à la liste

Numérotation des lignes : 

inctyp

  1. C INCTYP    SOURCE    CB215821  20/11/25    13:30:19     10792          
  2.       SUBROUTINE INCTYP(MATRIK,
  3.      $     ATYP,
  4.      $     IMPR,IRET)
  5.       IMPLICIT INTEGER(I-N)
  6.       IMPLICIT REAL*8 (A-H,O-Z)
  7. C***********************************************************************
  8. C NOM         : INCTYP
  9. C DESCRIPTION :
  10. C     Ce sp initialise le vecteur inconnue primale
  11. C     (i.e. x de Ax=b) avec un chpoint d'inconnues primales.
  12. C     On a en vue une méthode itérative de résolution.
  13. C
  14. C     Ce sp est quasiment identique à mesmbr.eso.
  15. C     C'est une interface à ch2vec.
  16. C
  17. C
  18. C LANGAGE     : ESOPE
  19. C AUTEUR      : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/TTMF)
  20. C               mél : gounand@semt2.smts.cea.fr
  21. C***********************************************************************
  22. C APPELE  : CH2VEC
  23. C***********************************************************************
  24. C ENTREES            : MATRIK, MCHINI, IMPR
  25. C ENTREES/SORTIES    : -
  26. C SORTIES            : INCX, IRET
  27. C CODE RETOUR (IRET) : 0 si ok
  28. C                     <0 si problème
  29. C     MATRIK : pointeur sur segment MATRIK de l'include SMMATRIK
  30. C              on pioche dedans les informations nécessaires
  31. C              (numérotations, inconnues, support géométriques)
  32. C              pour effectuer la conversion MCHINI->INCX.
  33. C     MCHINI : pointeur sur segment MCHPOI de l'include SMCHPOI
  34. C              chpoint d'initialisation du vecteur INCX
  35. C     IMPR   : niveau d'impression
  36. C     INCX   : pointeur sur segment IZA de l'include SMMATRIK
  37. C              vecteur x initial pour la résolution de Ax=b.
  38. C              Il est initialisé à 0 si MCHINI est nul
  39. C              sinon les valeurs non nulles de MCHINI
  40. C              correspondant à des inconnues primales sont
  41. C              chargées dans INCX.
  42. C              Les valeurs non nulles de MCHINI NE correspondant
  43. C              PAS à des inconnues primales donnent lieu à
  44. C              l'émission d'un warning (cf. CH2VEC).
  45. C***********************************************************************
  46. C VERSION    : v1, 01/04/98, version initiale
  47. C HISTORIQUE : v1, 01/04/98, création
  48. C HISTORIQUE :
  49. C HISTORIQUE :
  50. C***********************************************************************
  51. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  52. C en cas de modification de ce sous-programme afin de faciliter
  53. C la maintenance !
  54. C***********************************************************************
  55.  
  56. -INC PPARAM
  57. -INC CCOPTIO
  59.       POINTEUR INCX.IZA
  60. -INC SMLENTI
  61.       POINTEUR ATYP.MLENTI
  62. C
  63.       IRET=0
  64.       SEGACT MATRIK
  65. C Initialisation de l'inconnue à 0
  66.       JG=KNTTP
  67.       SEGINI ATYP
  68.       NBVA=KNTTP
  69.       MINC=KMINCP
  70.       SEGACT MINC
  71.       NBI=LISINC(/2)
  72.       NPT=NPOS(/1)-1
  73.       IDMAT=KIDMAT(1)
  74.       SEGACT IDMAT
  75.       DO IBI=1,NBI
  76.          DO IPT=1,NPT
  77.             IF (MPOS(IPT,IBI).NE.0) THEN
  78.                ATYP.LECT(NUAN(NPOS(IPT)+
  79.      $              MPOS(IPT,IBI)-1))=IBI
  80.             ENDIF
  81.          ENDDO
  82.       ENDDO
  83.       SEGDES IDMAT
  84.       SEGDES MINC
  85. C
  86.       IF (IMPR.GT.5) THEN
  87.          WRITE(IOIMP,*) 'inctyp.eso : création du pointeur ATYP=',ATYP
  88.          IF (IMPR.GT.7) THEN
  89.             WRITE(IOIMP,*) 'ATYP(1..',JG,')= '
  90.             WRITE(IOIMP,1002)(ATYP.LECT(II),II=1,JG)
  91.          ENDIF
  92.       ENDIF
  93. C
  94.       SEGDES ATYP
  95.       SEGDES MATRIK
  96. *
  97. *     Normal termination
  98. *
  99.       RETURN
  100. *
  101. * Format handling
  102. *
  103.  1002 FORMAT(10(1X,I7))
  104. *
  105. * Error handling
  106. *
  107.  9999 CONTINUE
  108.       WRITE(IOIMP,*) 'An error was detected in inctyp.eso'
  109.       RETURN
  110. *
  111. * End of INCTYP
  112. *
  113.       END
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  

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