Télécharger prbpol.eso

Retour à la liste

Numérotation des lignes :

prbpol
  1. C PRBPOL SOURCE GOUNAND 21/06/02 21:17:18 11022
  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 PPARAM
  33. -INC CCOPTIO
  34. -INC TNLIN
  35. *-INC SPOLYNO
  36. POINTEUR BAPOLY.POLYNS
  37. POINTEUR MYPOLY.POLYNO
  38. *
  39. INTEGER IMPR,IRET
  40. *
  41. INTEGER BAPETA,MYPETA
  42. INTEGER NMONO,NPOLY,NDIML
  43. INTEGER IMONO,IPOLY,INDIML
  44. *
  45. * Executable statements
  46. *
  47. IF (IMPR.GT.6) WRITE(IOIMP,*) 'Entrée dans prbpol.eso'
  48. CALL OOOETA(BAPOLY,BAPETA,IMOD)
  49. IF (BAPETA.NE.1) SEGACT BAPOLY
  50. NPOLY=BAPOLY.LIPOLY(/1)
  51. DO 1 IPOLY=1,NPOLY
  52. WRITE(IOIMP,*) 'Polynôme ',IPOLY
  53. MYPOLY=BAPOLY.LIPOLY(IPOLY)
  54. CALL OOOETA(MYPOLY,MYPETA,IMOD)
  55. IF (MYPETA.NE.1) SEGACT MYPOLY
  56. NMONO=MYPOLY.COEMON(/1)
  57. NDIML=MYPOLY.EXPMON(/1)
  58. DO 12 IMONO=1,NMONO
  59. WRITE(IOIMP,*) 'Monôme : ',IMONO,
  60. $ ' coeff. et puissances de ksi_i'
  61. WRITE(IOIMP,4004) MYPOLY.COEMON(IMONO)
  62. WRITE(IOIMP,3004)
  63. $ (MYPOLY.EXPMON(INDIML,IMONO),INDIML=1,NDIML)
  64. 12 CONTINUE
  65. IF (MYPETA.NE.1) SEGDES MYPOLY
  66. 1 CONTINUE
  67. IF (BAPETA.NE.1) SEGDES BAPOLY
  68. *
  69. * Normal termination
  70. *
  71. IRET=0
  72. RETURN
  73. *
  74. * Format handling
  75. *
  76. 3004 FORMAT (2X,10(1X,I4))
  77. 4004 FORMAT (2X,6(1X,1PE13.5))
  78. *
  79. * Error handling
  80. *
  81. 9999 CONTINUE
  82. IRET=1
  83. WRITE(IOIMP,*) 'An error was detected in subroutine prbpol'
  84. RETURN
  85. *
  86. * End of subroutine PRBPOL
  87. *
  88. END
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  

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