Télécharger r1f2kb.eso

Retour à la liste

Numérotation des lignes :

  1. C R1F2KB SOURCE BP208322 18/10/08 21:15:08 9952
  2. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  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 R1F2KB (IDO,L1,CC,IN1,CH,IN2,WA1)
  11. IMPLICIT INTEGER(I-N)
  12. IMPLICIT REAL*8(A-H,O-Z)
  13. REAL*8 CC(IN1,IDO,2,L1), CH(IN2,IDO,L1,2), WA1(IDO)
  14. C
  15. DO 101 K=1,L1
  16. CH(1,1,K,1) = CC(1,1,1,K)+CC(1,IDO,2,K)
  17. CH(1,1,K,2) = CC(1,1,1,K)-CC(1,IDO,2,K)
  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.  
  25. CH(1,I-1,K,1) = CC(1,I-1,1,K)+CC(1,IC-1,2,K)
  26. CH(1,I,K,1) = CC(1,I,1,K)-CC(1,IC,2,K)
  27.  
  28. CH(1,I-1,K,2) = WA1(I-2)*(CC(1,I-1,1,K)-CC(1,IC-1,2,K))
  29. 1 -WA1(I-1)*(CC(1,I,1,K)+CC(1,IC,2,K))
  30. CH(1,I,K,2) = WA1(I-2)*(CC(1,I,1,K)+CC(1,IC,2,K))+WA1(I-1)
  31. 1 *(CC(1,I-1,1,K)-CC(1,IC-1,2,K))
  32.  
  33. 103 CONTINUE
  34. 104 CONTINUE
  35. IF (MOD(IDO,2) .EQ. 1) RETURN
  36. 105 DO 106 K=1,L1
  37. CH(1,IDO,K,1) = CC(1,IDO,1,K)+CC(1,IDO,1,K)
  38. CH(1,IDO,K,2) = -(CC(1,1,2,K)+CC(1,1,2,K))
  39. 106 CONTINUE
  40. 107 RETURN
  41. END
  42.  
  43.  
  44.  

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