prbpol
C PRBPOL SOURCE GOUNAND 21/06/02 21:17:18 11022 $ IMPR,IRET) IMPLICIT REAL*8 (A-H,O-Z) IMPLICIT INTEGER (I-N) C*********************************************************************** C NOM : PRBPOL C DESCRIPTION : Imprime un segment décrivant une liste de polynômes. C C LANGAGE : ESOPE C AUTEUR : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF) C mél : gounand@semt2.smts.cea.fr C*********************************************************************** C APPELES : - C APPELES (E/S) : OOOETA C APPELE PAR : PRLRF C*********************************************************************** C ENTREES : BAPOLY C ENTREES/SORTIES : - C SORTIES : - C CODE RETOUR (IRET) : = 0 si tout s'est bien passé C*********************************************************************** C VERSION : v1, 10/05/2000, version initiale C HISTORIQUE : v1, 10/05/2000, création C HISTORIQUE : C HISTORIQUE : C*********************************************************************** C Prière de PRENDRE LE TEMPS de compléter les commentaires C en cas de modification de ce sous-programme afin de faciliter C la maintenance ! C*********************************************************************** -INC PPARAM -INC CCOPTIO -INC TNLIN *-INC SPOLYNO POINTEUR BAPOLY.POLYNS POINTEUR MYPOLY.POLYNO * INTEGER IMPR,IRET * INTEGER BAPETA,MYPETA INTEGER NMONO,NPOLY,NDIML INTEGER IMONO,IPOLY,INDIML * * Executable statements * IF (IMPR.GT.6) WRITE(IOIMP,*) 'Entrée dans prbpol.eso' CALL OOOETA(BAPOLY,BAPETA,IMOD) IF (BAPETA.NE.1) SEGACT BAPOLY NPOLY=BAPOLY.LIPOLY(/1) DO 1 IPOLY=1,NPOLY WRITE(IOIMP,*) 'Polynôme ',IPOLY MYPOLY=BAPOLY.LIPOLY(IPOLY) CALL OOOETA(MYPOLY,MYPETA,IMOD) IF (MYPETA.NE.1) SEGACT MYPOLY NMONO=MYPOLY.COEMON(/1) NDIML=MYPOLY.EXPMON(/1) DO 12 IMONO=1,NMONO WRITE(IOIMP,*) 'Monôme : ',IMONO, $ ' coeff. et puissances de ksi_i' WRITE(IOIMP,4004) MYPOLY.COEMON(IMONO) WRITE(IOIMP,3004) $ (MYPOLY.EXPMON(INDIML,IMONO),INDIML=1,NDIML) 12 CONTINUE IF (MYPETA.NE.1) SEGDES MYPOLY 1 CONTINUE IF (BAPETA.NE.1) SEGDES BAPOLY * * Normal termination * IRET=0 RETURN * * Format handling * 3004 FORMAT (2X,10(1X,I4)) 4004 FORMAT (2X,6(1X,1PE13.5)) * * Error handling * 9999 CONTINUE IRET=1 WRITE(IOIMP,*) 'An error was detected in subroutine prbpol' RETURN * * End of subroutine PRBPOL * END
© Cast3M 2003 - Tous droits réservés.
Mentions légales