Télécharger la plaquette Cast3M


fr en

Télécharger r1f2kf.eso

Retour à la liste

Numérotation des lignes : 

r1f2kf

  1. C R1F2KF    SOURCE    BP208322  18/10/08    21:15:09     9952           
  2. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  3. C
  4. C   FFTPACK 5.1 
  5. C
  6. C   Authors:  Paul N. Swarztrauber and Richard A. Valent
  7. C
  8. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  9.  
  10.       SUBROUTINE R1F2KF (IDO,L1,CC,IN1,CH,IN2,WA1)
  11.       IMPLICIT INTEGER(I-N)
  12.       IMPLICIT REAL*8(A-H,O-Z)
  13.       REAL*8     CH(IN2,IDO,2,L1) ,CC(IN1,IDO,L1,2) , WA1(IDO)
  14. C
  15.       DO 101 K=1,L1
  16.          CH(1,1,1,K) = CC(1,1,K,1)+CC(1,1,K,2)
  17.          CH(1,IDO,2,K) = CC(1,1,K,1)-CC(1,1,K,2)
  18.   101 CONTINUE
  19.       IF (IDO-2) 107,105,102
  20.   102 IDP2 = IDO+2
  21.       DO 104 K=1,L1
  22.          DO 103 I=3,IDO,2
  23.             IC = IDP2-I
  24.             CH(1,I,1,K) = CC(1,I,K,1)+(WA1(I-2)*CC(1,I,K,2)-
  25.      1       WA1(I-1)*CC(1,I-1,K,2))
  26.             CH(1,IC,2,K) = (WA1(I-2)*CC(1,I,K,2)-WA1(I-1)*
  27.      1       CC(1,I-1,K,2))-CC(1,I,K,1)
  28.             CH(1,I-1,1,K) = CC(1,I-1,K,1)+(WA1(I-2)*CC(1,I-1,K,2)+
  29.      1       WA1(I-1)*CC(1,I,K,2))
  30.             CH(1,IC-1,2,K) = CC(1,I-1,K,1)-(WA1(I-2)*CC(1,I-1,K,2)+
  31.      1       WA1(I-1)*CC(1,I,K,2))
  32.   103    CONTINUE
  33.   104 CONTINUE
  34.       IF (MOD(IDO,2) .EQ. 1) RETURN
  35.   105 DO 106 K=1,L1
  36.          CH(1,1,2,K) = -CC(1,IDO,K,2)
  37.          CH(1,IDO,1,K) = CC(1,IDO,K,1)
  38.   106 CONTINUE
  39.   107 RETURN
  40.       END
  41.  
  42.  
  43.  

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