Télécharger hbmrww.eso

Retour à la liste

Numérotation des lignes :

hbmrww
  1. C HBMRWW SOURCE FANDEUR 22/05/02 21:15:23 11359
  2. SUBROUTINE HBMRWW(NT,NDDL,Q1,PHI,XM,RwwPhi)
  3.  
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6.  
  7. INTEGER NT,NHBM,NDDL,I,J
  8. REAL*8 AA,Q1(NT),PHI(NT),AUX(NT),XM(NDDL,NDDL),RwwPhi
  9.  
  10. C Fonctions BLAS/LAPACK
  11. REAL*8 DDOT
  12. EXTERNAL DDOT
  13. *
  14. NHBM = (NT/NDDL - 1)/2
  15. *
  16. * Etape intermediaire: calcul de 2*(L^2(o)M)*Q1
  17. DO I = 1,NDDL
  18. AUX(I)=0.
  19. ENDDO
  20. DO J=2,2*NHBM,2
  21. DO I=1,NDDL
  22. AA = -2.*XM(I,1)*(J/2)**2
  23. AUX(NDDL*(1+(J-2))+I) = AA*Q1(NDDL*(1+(J-2))+I)
  24. AUX(NDDL*(1+(J-1))+I) = AA*Q1(NDDL*(1+(J-1))+I)
  25. ENDDO
  26. ENDDO
  27. * Produit scalaire avec PHI
  28. RwwPhi = DDOT(NT,AUX,1,PHI,1)
  29.  
  30. END
  31.  
  32.  
  33.  

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