Télécharger mamupw.eso

Retour à la liste

Numérotation des lignes :

mamupw
  1. C MAMUPW SOURCE MB234859 26/01/26 21:15:10 12460
  2. SUBROUTINE MAMUPW(IDEB,IFIN,VAL2,IPOSR2,LPL2,NA2,
  3. & VAL1,IPOSR1,LPL1,NA1,PT,NBO)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. -INC CCHOLE
  7. DIMENSION VAL1(*),VAL2(*)
  8. REAL*8 PT(36)
  9. C
  10. DO I=1,MIN(6,NA1)*MIN(6,NA2)
  11. PT(I)=0.D0
  12. ENDDO
  13. C
  14. LON=IFIN-IDEB+1
  15. IF (LON.LE.0) GOTO 999
  16. C
  17. NBO=NBO+LON*(NA1*NA2)
  18. C
  19. IPOS1=IDEB+IPOSR1
  20. IPOS2=IDEB+IPOSR2
  21. IF (NA2.GE.6) THEN
  22. IF (NA1.GE.6) THEN
  23. call mamu66(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  24. ELSEIF (NA1.GE.5) THEN
  25. call mamu65(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  26. ELSEIF (NA1.GE.4) THEN
  27. call mamu64(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  28. ELSEIF (NA1.GE.3) THEN
  29. call mamu63(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  30. ELSEIF (NA1.GE.2) THEN
  31. call mamu62(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  32. ELSEIF (NA1.GE.1) THEN
  33. call mamu61(LON,VAL1(IPOS1), VAL2(IPOS2),LPL2,pt)
  34. ENDIF
  35. ELSEIF (NA2.GE.5) THEN
  36. IF (NA1.GE.6) THEN
  37. call mamu56(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  38. ELSEIF (NA1.GE.5) THEN
  39. call mamu55(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  40. ELSEIF (NA1.GE.4) THEN
  41. call mamu54(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  42. ELSEIF (NA1.GE.3) THEN
  43. call mamu53(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  44. ELSEIF (NA1.GE.2) THEN
  45. call mamu52(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  46. ELSEIF (NA1.GE.1) THEN
  47. call mamu51(LON,VAL1(IPOS1), VAL2(IPOS2),LPL2,pt)
  48. ENDIF
  49. ELSEIF (NA2.GE.4) THEN
  50. IF (NA1.GE.6) THEN
  51. call mamu46(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  52. ELSEIF (NA1.GE.5) THEN
  53. call mamu45(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  54. ELSEIF (NA1.GE.4) THEN
  55. call mamu44(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  56. ELSEIF (NA1.GE.3) THEN
  57. call mamu43(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  58. ELSEIF (NA1.GE.2) THEN
  59. call mamu42(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  60. ELSEIF (NA1.GE.1) THEN
  61. call mamu41(LON,VAL1(IPOS1), VAL2(IPOS2),LPL2,pt)
  62. ENDIF
  63. ELSEIF (NA2.GE.3) THEN
  64. IF (NA1.GE.6) THEN
  65. call mamu36(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  66. ELSEIF (NA1.GE.5) THEN
  67. call mamu35(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  68. ELSEIF (NA1.GE.4) THEN
  69. call mamu34(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  70. ELSEIF (NA1.GE.3) THEN
  71. call mamu33(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  72. ELSEIF (NA1.GE.2) THEN
  73. call mamu32(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  74. ELSEIF (NA1.GE.1) THEN
  75. call mamu31(LON,VAL1(IPOS1), VAL2(IPOS2),LPL2,pt)
  76. ENDIF
  77. ELSEIF (NA2.GE.2) THEN
  78. IF (NA1.GE.6) THEN
  79. call mamu26(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  80. ELSEIF (NA1.GE.5) THEN
  81. call mamu25(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  82. ELSEIF (NA1.GE.4) THEN
  83. call mamu24(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  84. ELSEIF (NA1.GE.3) THEN
  85. call mamu23(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  86. ELSEIF (NA1.GE.2) THEN
  87. call mamu22(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2),LPL2,pt)
  88. ELSEIF (NA1.GE.1) THEN
  89. call mamu21(LON,VAL1(IPOS1), VAL2(IPOS2),LPL2,pt)
  90. ENDIF
  91. ELSEIF (NA2.GE.1) THEN
  92. IF (NA1.GE.6) THEN
  93. call mamu16(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2), pt)
  94. ELSEIF (NA1.GE.5) THEN
  95. call mamu15(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2), pt)
  96. ELSEIF (NA1.GE.4) THEN
  97. call mamu14(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2), pt)
  98. ELSEIF (NA1.GE.3) THEN
  99. call mamu13(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2), pt)
  100. ELSEIF (NA1.GE.2) THEN
  101. call mamu12(LON,VAL1(IPOS1),lpl1,VAL2(IPOS2), pt)
  102. ELSEIF (NA1.GE.1) THEN
  103. pt(1)=pt(1)+ddotpv((LON),VAL1(IPOS1),VAL2(IPOS2))
  104. ENDIF
  105. ENDIF
  106. 999 CONTINUE
  107. END
  108.  
  109.  

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