Télécharger ddotpv.eso

Retour à la liste

Numérotation des lignes :

ddotpv
  1. C DDOTPV SOURCE PV090527 23/01/09 21:15:04 11549
  2. function ddotpv(lon,val1,val2)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. dimension val1(*),val2(*)
  6. if (lon.eq.0) then
  7. ddotpv=0.d0
  8. return
  9. endif
  10. resu1=0.d0
  11. resu2=0.d0
  12. resu3=0.d0
  13. resu4=0.d0
  14. do 10 i=0,lon-4,4
  15. resu1=resu1+val1(i+1)*val2(i+1)
  16. resu2=resu2+val1(i+2)*val2(i+2)
  17. resu3=resu3+val1(i+3)*val2(i+3)
  18. resu4=resu4+val1(i+4)*val2(i+4)
  19. 10 continue
  20. if (i.le.lon-2) then
  21. resu1=resu1+val1(i+1)*val2(i+1)
  22. resu2=resu2+val1(i+2)*val2(i+2)
  23. i=i+2
  24. endif
  25. if (i.le.lon-1) then
  26. resu1=resu1+val1(i+1)*val2(i+1)
  27. endif
  28. ddotpv=resu1+resu2+resu3+resu4
  29. return
  30. end
  31.  
  32.  
  33.  
  34.  

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