Télécharger mamu55.eso

Retour à la liste

Numérotation des lignes :

mamu55
  1. C MAMU55 SOURCE PV 20/02/22 21:15:20 10531
  2. subroutine mamu55(long,val11,lpl1,val1,lpl,pt)
  3. implicit real*8 (a-h,o-z)
  4. real*8 val11(*),val1(*)
  5. real*8 pt(36)
  6. id12=lpl1
  7. id13=lpl1+ lpl1+1
  8. id14=lpl1+ lpl1+1+ lpl1+2
  9. id15=lpl1+ lpl1+1+ lpl1+2+ lpl1+3
  10. id2=lpl
  11. id3=lpl+ lpl+1
  12. id4=lpl+ lpl+1 +lpl+2
  13. id5=lpl+ lpl+1 +lpl+2 +lpl+3
  14. pt1=pt(1)
  15. pt2=pt(2)
  16. pt3=pt(3)
  17. pt4=pt(4)
  18. pt5=pt(5)
  19. pt7=pt(6)
  20. pt8=pt(7)
  21. pt9=pt(8)
  22. pt10=pt(9)
  23. pt11=pt(10)
  24. pt13=pt(11)
  25. pt14=pt(12)
  26. pt15=pt(13)
  27. pt16=pt(14)
  28. pt17=pt(15)
  29. pt19=pt(16)
  30. pt20=pt(17)
  31. pt21=pt(18)
  32. pt22=pt(19)
  33. pt23=pt(20)
  34. pt25=pt(21)
  35. pt26=pt(22)
  36. pt27=pt(23)
  37. pt28=pt(24)
  38. pt29=pt(25)
  39. do 110 i=1,long
  40. xval1=val1(i)
  41. xval2=val1(i+id2)
  42. xval3=val1(i+id3)
  43. xval4=val1(i+id4)
  44. xval5=val1(i+id5)
  45. xval11=val11(i)
  46. xval12=val11(i+id12)
  47. xval13=val11(i+id13)
  48. xval14=val11(i+id14)
  49. xval15=val11(i+id15)
  50. pt1=pt1+xval1*xval11
  51. pt2=pt2+xval1*xval12
  52. pt3=pt3+xval1*xval13
  53. pt4=pt4+xval1*xval14
  54. pt5=pt5+xval1*xval15
  55. pt7=pt7+xval2*xval11
  56. pt8=pt8+xval2*xval12
  57. pt9=pt9+xval2*xval13
  58. pt10=pt10+xval2*xval14
  59. pt11=pt11+xval2*xval15
  60. pt13=pt13+xval3*xval11
  61. pt14=pt14+xval3*xval12
  62. pt15=pt15+xval3*xval13
  63. pt16=pt16+xval3*xval14
  64. pt17=pt17+xval3*xval15
  65. pt19=pt19+xval4*xval11
  66. pt20=pt20+xval4*xval12
  67. pt21=pt21+xval4*xval13
  68. pt22=pt22+xval4*xval14
  69. pt23=pt23+xval4*xval15
  70. pt25=pt25+xval5*xval11
  71. pt26=pt26+xval5*xval12
  72. pt27=pt27+xval5*xval13
  73. pt28=pt28+xval5*xval14
  74. pt29=pt29+xval5*xval15
  75. 110 continue
  76. pt(1)=pt1
  77. pt(2)=pt2
  78. pt(3)=pt3
  79. pt(4)=pt4
  80. pt(5)=pt5
  81. pt(6)=pt7
  82. pt(7)=pt8
  83. pt(8)=pt9
  84. pt(9)=pt10
  85. pt(10)=pt11
  86. pt(11)=pt13
  87. pt(12)=pt14
  88. pt(13)=pt15
  89. pt(14)=pt16
  90. pt(15)=pt17
  91. pt(16)=pt19
  92. pt(17)=pt20
  93. pt(18)=pt21
  94. pt(19)=pt22
  95. pt(20)=pt23
  96. pt(21)=pt25
  97. pt(22)=pt26
  98. pt(23)=pt27
  99. pt(24)=pt28
  100. pt(25)=pt29
  101.  
  102. end
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  

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