Sources Esope | Cast3M

ift
C IFT       SOURCE    OF166741  26/05/11    21:15:19     12538          
 
*=======================================================================
* IFT: Transformee de Fourier Inverse
*
* Entrée: XF , coefficients de Fourier
*        GAM , opérateur IFT discret
* Sortie: XT , déplacements temporaux
*
*=======================================================================
 
      SUBROUTINE IFT(XF,NT,NDDL,NHBM,NFFT,GAM,XT)
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
 
      REAL*8 XF(NT),XT(NDDL,NFFT),GAM(NFFT,2*NHBM+1)
 
*     Boucle sur le temps
      DO K = 1, NFFT
*       Boucle sur les DDL
        DO J=1, NDDL
*         Harmonique 0
          r_z = XF(J)
*         Boucle sur les harmoniques > 0
          DO I = 2, 2*NHBM+1
            r_z = r_z + GAM(K,I)*XF(NDDL*(I-1)+J)
          ENDDO
          XT(J,K) = r_z
        ENDDO
      ENDDO
 
      RETURN
      END