Mercurial > octave-nkf
comparison libcruft/fftpack/passf3.f @ 7789:82be108cc558
First attempt at single precision tyeps
* * *
corrections to qrupdate single precision routines
* * *
prefer demotion to single over promotion to double
* * *
Add single precision support to log2 function
* * *
Trivial PROJECT file update
* * *
Cache optimized hermitian/transpose methods
* * *
Add tests for tranpose/hermitian and ChangeLog entry for new transpose code
author | David Bateman <dbateman@free.fr> |
---|---|
date | Sun, 27 Apr 2008 22:34:17 +0200 |
parents | fccab8e7d35f |
children |
comparison
equal
deleted
inserted
replaced
7788:45f5faba05a2 | 7789:82be108cc558 |
---|---|
1 subroutine passf3 (ido,l1,cc,ch,wa1,wa2) | 1 subroutine passf3 (ido,l1,cc,ch,wa1,wa2) |
2 implicit double precision (a-h,o-z) | |
3 dimension cc(ido,3,l1) ,ch(ido,l1,3) , | 2 dimension cc(ido,3,l1) ,ch(ido,l1,3) , |
4 1 wa1(1) ,wa2(1) | 3 1 wa1(1) ,wa2(1) |
5 data taur,taui /-.5d0,-.866025403784439d0/ | 4 data taur,taui /-.5,-.866025403784439/ |
6 if (ido .ne. 2) go to 102 | 5 if (ido .ne. 2) go to 102 |
7 do 101 k=1,l1 | 6 do 101 k=1,l1 |
8 tr2 = cc(1,2,k)+cc(1,3,k) | 7 tr2 = cc(1,2,k)+cc(1,3,k) |
9 cr2 = cc(1,1,k)+taur*tr2 | 8 cr2 = cc(1,1,k)+taur*tr2 |
10 ch(1,k,1) = cc(1,1,k)+tr2 | 9 ch(1,k,1) = cc(1,1,k)+tr2 |