Télécharger mamu33.eso

Retour à la liste

Numérotation des lignes :

  1. C MAMU33 SOURCE PV 20/01/14 21:15:13 10498
  2. subroutine mamu33(long,val11,lpl1,val1,lpl,pt)
  3. implicit real*8 (a-h,o-z)
  4. real*8 val11(*),val1(*)
  5. real*8 pt(6,6)
  6. id12=lpl1
  7. id13=lpl1+lpl1+1
  8. id2=lpl
  9. id3=lpl+lpl+1
  10. pt1=0.d0
  11. pt2=0.d0
  12. pt3=0.d0
  13. pt4=0.d0
  14. pt5=0.d0
  15. pt6=0.d0
  16. pt7=0.d0
  17. pt8=0.d0
  18. pt9=0.d0
  19. do 110 i=1,long
  20. xval11=val11(i)
  21. xval12=val11(i+id12)
  22. xval13=val11(i+id13)
  23. xval1=val1(i)
  24. xval2=val1(i+id2)
  25. xval3=val1(i+id3)
  26. pt1=pt1+xval1*xval11
  27. pt2=pt2+xval1*xval12
  28. pt3=pt3+xval1*xval13
  29. pt4=pt4+xval2*xval11
  30. pt5=pt5+xval2*xval12
  31. pt6=pt6+xval2*xval13
  32. pt7=pt7+xval3*xval11
  33. pt8=pt8+xval3*xval12
  34. pt9=pt9+xval3*xval13
  35. 110 continue
  36. pt(1,1)=pt(1,1)+pt1
  37. pt(2,1)=pt(2,1)+pt2
  38. pt(3,1)=pt(3,1)+pt3
  39. pt(1,2)=pt(1,2)+pt4
  40. pt(2,2)=pt(2,2)+pt5
  41. pt(3,2)=pt(3,2)+pt6
  42. pt(1,3)=pt(1,3)+pt7
  43. pt(2,3)=pt(2,3)+pt8
  44. pt(3,3)=pt(3,3)+pt9
  45.  
  46.  
  47. end
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  

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