Télécharger la plaquette Cast3M


fr en

Télécharger mexini.eso

Retour à la liste

Numérotation des lignes : 

mexini

  1. C MEXINI    SOURCE    CB215821  20/11/25    13:34:26     10792          
  2.       SUBROUTINE MEXINI(MATRIK,MCHINI,
  3.      $     INCX,
  4.      $     IMPR,IRET)
  5.       IMPLICIT INTEGER(I-N)
  6.       IMPLICIT REAL*8 (A-H,O-Z)
  7. C***********************************************************************
  8. C NOM         : MEXINI
  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 SMCHPOI
  61.       POINTEUR MCHINI.MCHPOI
  62. -INC SMELEME
  63.       POINTEUR KISPG.MELEME
  64. C
  65.       IRET=0
  66.       SEGACT MATRIK
  67. C Initialisation de l'inconnue à 0
  68.       NBVA=KNTTP
  69.       SEGINI INCX
  70.       IF(MCHINI.NE.0) THEN
  71. C
  72. C On transforme le chpoint en vecteur
  73. C
  74. C     Entrées
  75.          MCHPOI=MCHINI
  76.          MINC  =KMINCP
  77.          KISPG =KISPGP
  78.          IDMAT =KIDMAT(1)
  79. C     Sortie : INCX
  80. C
  81.          CALL CH2VEC(MCHPOI,MINC,KISPG,IDMAT,INCX,IMPR,IRET)
  82.          IF (IRET.NE.0) GOTO 9999
  83. C
  84.       ENDIF
  85. C
  86.       IF (IMPR.GT.5) THEN
  87.          WRITE(IOIMP,*) 'mexini.eso : création du pointeur INCX=',INCX
  88.          WRITE(IOIMP,*) '            a partir du CHPOINT MCHINI=',MCHINI
  89.          IF (IMPR.GT.7) THEN
  90.             WRITE(IOIMP,*) 'INCX(1..',NBVA,')= '
  91.             WRITE(IOIMP,1002)(INCX.A(II),II=1,NBVA)
  92.          ENDIF
  93.       ENDIF
  94. C
  95.       SEGDES INCX
  96.       SEGDES MATRIK
  97. *
  98. *     Normal termination
  99. *
  100.       RETURN
  101. *
  102. * Format handling
  103. *
  104.  1002 FORMAT(10(1X,1PE11.4))
  105. *
  106. * Error handling
  107. *
  108.  9999 CONTINUE
  109.       WRITE(IOIMP,*) 'An error was detected in mexini.eso'
  110.       RETURN
  111. *
  112. * End of MEXINI
  113. *
  114.       END
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  

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