Télécharger prbpol.eso

Retour à la liste

Numérotation des lignes :

  1. C PRBPOL SOURCE GOUNAND 05/12/21 21:34:49 5281
  2. SUBROUTINE PRBPOL(BAPOLY,
  3. $ IMPR,IRET)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. IMPLICIT INTEGER (I-N)
  6. C***********************************************************************
  7. C NOM : PRBPOL
  8. C DESCRIPTION : Imprime un segment décrivant une liste de polynômes.
  9. C
  10. C LANGAGE : ESOPE
  11. C AUTEUR : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  12. C mél : gounand@semt2.smts.cea.fr
  13. C***********************************************************************
  14. C APPELES : -
  15. C APPELES (E/S) : OOOETA
  16. C APPELE PAR : PRLRF
  17. C***********************************************************************
  18. C ENTREES : BAPOLY
  19. C ENTREES/SORTIES : -
  20. C SORTIES : -
  21. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  22. C***********************************************************************
  23. C VERSION : v1, 10/05/2000, version initiale
  24. C HISTORIQUE : v1, 10/05/2000, création
  25. C HISTORIQUE :
  26. C HISTORIQUE :
  27. C***********************************************************************
  28. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  29. C en cas de modification de ce sous-programme afin de faciliter
  30. C la maintenance !
  31. C***********************************************************************
  32. -INC CCOPTIO
  33. CBEGININCLUDE SPOLYNO
  34. SEGMENT POLYNO
  35. REAL*8 COEMON(NBMON)
  36. INTEGER EXPMON(NDIML,NBMON)
  37. ENDSEGMENT
  38. SEGMENT POLYNS
  39. POINTEUR LIPOLY(NBPOLY).POLYNO
  40. ENDSEGMENT
  41. CENDINCLUDE SPOLYNO
  42. POINTEUR BAPOLY.POLYNS
  43. POINTEUR MYPOLY.POLYNO
  44. *
  45. INTEGER IMPR,IRET
  46. *
  47. INTEGER BAPETA,MYPETA
  48. INTEGER NMONO,NPOLY,NDIML
  49. INTEGER IMONO,IPOLY,INDIML
  50. *
  51. * Executable statements
  52. *
  53. IF (IMPR.GT.6) WRITE(IOIMP,*) 'Entrée dans prbpol.eso'
  54. CALL OOOETA(BAPOLY,BAPETA)
  55. IF (BAPETA.NE.1) SEGACT BAPOLY
  56. NPOLY=BAPOLY.LIPOLY(/1)
  57. DO 1 IPOLY=1,NPOLY
  58. WRITE(IOIMP,*) 'Polynôme ',IPOLY
  59. MYPOLY=BAPOLY.LIPOLY(IPOLY)
  60. CALL OOOETA(MYPOLY,MYPETA)
  61. IF (MYPETA.NE.1) SEGACT MYPOLY
  62. NMONO=MYPOLY.COEMON(/1)
  63. NDIML=MYPOLY.EXPMON(/1)
  64. DO 12 IMONO=1,NMONO
  65. WRITE(IOIMP,*) 'Monôme : ',IMONO,
  66. $ ' coeff. et puissances de ksi_i'
  67. WRITE(IOIMP,4004) MYPOLY.COEMON(IMONO)
  68. WRITE(IOIMP,3004)
  69. $ (MYPOLY.EXPMON(INDIML,IMONO),INDIML=1,NDIML)
  70. 12 CONTINUE
  71. IF (MYPETA.NE.1) SEGDES MYPOLY
  72. 1 CONTINUE
  73. IF (BAPETA.NE.1) SEGDES BAPOLY
  74. *
  75. * Normal termination
  76. *
  77. IRET=0
  78. RETURN
  79. *
  80. * Format handling
  81. *
  82. 3004 FORMAT (2X,10(1X,I4))
  83. 4004 FORMAT (2X,6(1X,1PE13.5))
  84. *
  85. * Error handling
  86. *
  87. 9999 CONTINUE
  88. IRET=1
  89. WRITE(IOIMP,*) 'An error was detected in subroutine prbpol'
  90. RETURN
  91. *
  92. * End of subroutine PRBPOL
  93. *
  94. END
  95.  
  96.  
  97.  
  98.  

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