Sources Esope | Cast3M

hbmrww
C HBMRWW    SOURCE    CB215821  26/06/25    21:15:13     12581          
 
      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))
        XJs2=REAL(J/2)
        DO I=1,NDDL
          AA = -2.*XM(I,1)* XJs2**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