Télécharger ddotpv.eso

Retour à la liste

Numérotation des lignes :

  1. C DDOTPV SOURCE PV 19/07/04 21:15:05 10249
  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. do 20 i=4*(lon/4)+1,lon
  21. resu1=resu1+val1(i)*val2(i)
  22. 20 continue
  23. ddotpv=resu1+resu2+resu3+resu4
  24. return
  25. end
  26.  
  27.  
  28.  

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