Télécharger r1f2kf.eso

Retour à la liste

Numérotation des lignes :

  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