Télécharger ltl.eso

Retour à la liste

Numérotation des lignes :

ltl
  1. C LTL SOURCE CHAT 05/01/13 01:25:53 5004
  2. SUBROUTINE LTL
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. C======================================================================
  6. C = CALCUL DU PRODUIT SCALAIRE DE DEUX LISTREEL =
  7. C = =
  8. C = SYNTAXE : =
  9. C = =
  10. C = FLOT = LTL LIST1 (LIST2 FLOT1 FLOT2) ; =
  11. C = =
  12. C = LIST1, LIST2 : LISTREEL =
  13. C = =
  14. C = FLOT, FLOT1, FLOT2 : FLOTTANT =
  15. C = =
  16. C = FLOT = (FLOT1*LIST1) SCALAIRE (FLOT2*LIST2) =
  17. C = =
  18. C = SI LIST2 N'EST PAS DONNE, LIST2=LIST1 =
  19. C = SI FLOT2 N'EST PAS DONNE, FLOT2=1. =
  20. C = ST FLOT1 ET FLOT2 NE SONT PAS DONNE ,FLOT1=FLOT2=1. =
  21. C = =
  22. C = CREATION : 24/09/90 =
  23. C = PROGRAMMEUR : P.P. =
  24. C======================================================================
  25. C
  26. -INC PPARAM
  27. -INC CCOPTIO
  28. -INC SMLREEL
  29. C
  30. FLOT1=1.D0
  31. C
  32. CALL LIROBJ('LISTREEL',MLREE1,1,IRET)
  33. IF(IRET.EQ.0) GOTO 666
  34. C
  35. CALL LIROBJ('LISTREEL',MLREE2,0,IRET)
  36. IF(IRET.EQ.0) MLREE2=0
  37. IF(MLREE1.EQ.MLREE2) MLREE2=0
  38. C
  39. CALL LIRREE(FLOT,0,IRET)
  40. IF(IRET.NE.0)THEN
  41. FLOT1=FLOT
  42. CALL LIRREE(FLOT,0,IRET)
  43. IF(IRET.NE.0)FLOT1=FLOT*FLOT1
  44. ENDIF
  45. C
  46. SEGACT, MLREE1
  47. NPT=MLREE1.PROG(/1)
  48. IF (MLREE2.NE.0)THEN
  49. SEGACT, MLREE2
  50. NPTT=MLREE2.PROG(/1)
  51. IF(NPT.NE.NPTT)THEN
  52. CALL ERREUR(577)
  53. SEGDES, MLREE1, MLREE2
  54. GOTO 666
  55. ENDIF
  56. ENDIF
  57. C
  58. FLOT=0.D0
  59. IF (MLREE2.EQ.0)THEN
  60. DO 10 IE1=1,NPT
  61. FLOT=FLOT+FLOT1*MLREE1.PROG(IE1)**2
  62. 10 CONTINUE
  63. SEGDES, MLREE1
  64. ELSE
  65. DO 11 IE1=1,NPT
  66. FLOT=FLOT+FLOT1*MLREE1.PROG(IE1)*MLREE2.PROG(IE1)
  67. 11 CONTINUE
  68. SEGDES, MLREE1, MLREE2
  69. ENDIF
  70. C
  71. CALL ECRREE(FLOT)
  72. C
  73. 666 RETURN
  74. END
  75.  
  76.  
  77.  
  78.  
  79.  

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