Sources Esope | Cast3M

hbmdvec
C HBMDVEC   SOURCE    OF166741  26/05/11    21:15:10     12538          
 
      SUBROUTINE HBMDVEC(NT,NHBM,NDDL,Q1,OMEG,V1)
 
c-----------------------------------------------------------------------
c       Derivation dans le domaine frequentiel
c                         V1 = w*kron(L,In)*Q1
c       e  NT,NHBM,NDDL : nombre de coefficients, harmoniques, ddl
c       e  Q1           : vecteur de coefficients à deriver
c       e  OMEG         : pulsation du cycle
c       s  V1           : derivee du vecteur Q1 en frequentiel
c-----------------------------------------------------------------------
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
      INTEGER NT,NHBM,NDDL
      REAL*8 OMEG
c*      REAL*8 Q1(NT),V1(NT)
      REAL*8 Q1(*),V1(*)
 
c   Composante statique nulle
      DO I=1,NDDL
        V1(I) = 0.D0
      ENDDO
 
c   Boucle sur les harmoniques > 0
      DO J=2,2*NHBM,2
        r_z = OMEG*(J/2)
        j_z1 = (J-1)*NDDL
        j_z2 = j_z1 + NDDL
        DO I=1,NDDL
          V1(j_z1+I) =  r_z*Q1(j_z2+I)
          V1(j_z2+I) = -r_z*Q1(j_z1+I)
        ENDDO
      ENDDO
 
      RETURN
      END