Sources Esope | Cast3M

hbmrww
C HBMRWW    SOURCE    OF166741  26/05/11    21:15:15     12538          
 
      SUBROUTINE HBMRWW(NT,NDDL,Q1,PHI,XM,RwwPhi)
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
      INTEGER NT,NHBM,NDDL,I,J
*      REAL*8 Q1(NT),PHI(NT),XM(NDDL,NDDL),RwwPhi
      REAL*8 Q1(*),PHI(*),XM(NDDL,*),RwwPhi
 
      REAL*8 AA
 
      NHBM = (NT/NDDL - 1)/2
 
*   Etape intermediaire: calcul de 2*(L^2(o)M)*Q1
*   Produit scalaire avec PHI
      RwwPhi = 0.D0
      DO J=2,2*NHBM,2
        j_z = NDDL*(1+(J-2))
        DO I=1,NDDL
          AA = -2.*XM(I,1)*(J/2)**2
          i_z = j_z + I
          RwwPhi = RwwPhi + AA*Q1(i_z)*PHI(i_z)
          i_z = i_z + NDDL
          RwwPhi = RwwPhi + AA*Q1(i_z)*PHI(i_z)
        ENDDO
      ENDDO
 
      RETURN
      END