Télécharger inctyp.eso

Retour à la liste

Numérotation des lignes :

  1. C INCTYP SOURCE PV 16/11/17 21:59:40 9180
  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. -INC CCOPTIO
  56. POINTEUR INCX.IZA
  57. -INC SMLENTI
  58. POINTEUR ATYP.MLENTI
  59. C
  60. IRET=0
  61. SEGACT MATRIK
  62. C Initialisation de l'inconnue à 0
  63. JG=KNTTP
  64. SEGINI ATYP
  65. NBVA=KNTTP
  66. MINC=KMINCP
  67. SEGACT MINC
  68. NBI=LISINC(/2)
  69. NPT=NPOS(/1)-1
  70. IDMAT=KIDMAT(1)
  71. SEGACT IDMAT
  72. DO IBI=1,NBI
  73. DO IPT=1,NPT
  74. IF (MPOS(IPT,IBI).NE.0) THEN
  75. ATYP.LECT(NUAN(NPOS(IPT)+
  76. $ MPOS(IPT,IBI)-1))=IBI
  77. ENDIF
  78. ENDDO
  79. ENDDO
  80. SEGDES IDMAT
  81. SEGDES MINC
  82. C
  83. IF (IMPR.GT.5) THEN
  84. WRITE(IOIMP,*) 'inctyp.eso : création du pointeur ATYP=',ATYP
  85. IF (IMPR.GT.7) THEN
  86. WRITE(IOIMP,*) 'ATYP(1..',JG,')= '
  87. WRITE(IOIMP,1002)(ATYP.LECT(II),II=1,JG)
  88. ENDIF
  89. ENDIF
  90. C
  91. SEGDES ATYP
  92. SEGDES MATRIK
  93. *
  94. * Normal termination
  95. *
  96. RETURN
  97. *
  98. * Format handling
  99. *
  100. 1002 FORMAT(10(1X,I7))
  101. *
  102. * Error handling
  103. *
  104. 9999 CONTINUE
  105. WRITE(IOIMP,*) 'An error was detected in inctyp.eso'
  106. RETURN
  107. *
  108. * End of INCTYP
  109. *
  110. END
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  

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