/ Usenet 1994 January / usenetsourcesnewsgroupsinfomagicjanuary1994.iso / sources / misc / volume1 / 8707 / 48

< prev

next >

view JSON data

view as text

open on a Mac

open on a PC

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 46 72 6f 6d 3a 20 61 6c | 65 65 40 74 79 62 61 6c |From: al|ee@tybal|
|00000010| 74 2e 63 61 6c 74 65 63 | 68 2e 65 64 75 20 28 41 |t.caltec|h.edu (A|
|00000020| 6e 64 72 65 77 20 4c 65 | 65 29 0a 4e 65 77 73 67 |ndrew Le|e).Newsg|
|00000030| 72 6f 75 70 73 3a 20 63 | 6f 6d 70 2e 73 6f 75 72 |roups: c|omp.sour|
|00000040| 63 65 73 2e 6d 69 73 63 | 0a 53 75 62 6a 65 63 74 |ces.misc|.Subject|
|00000050| 3a 20 4e 2d 64 69 6d 65 | 6e 73 69 6f 6e 61 6c 2c |: N-dime|nsional,|
|00000060| 20 52 61 64 69 78 20 32 | 20 46 46 54 20 52 6f 75 | Radix 2| FFT Rou|
|00000070| 74 69 6e 65 0a 4d 65 73 | 73 61 67 65 2d 49 44 3a |tine.Mes|sage-ID:|
|00000080| 20 3c 32 39 32 37 40 6e | 63 6f 61 73 74 2e 55 55 | <2927@n|coast.UU|
|00000090| 43 50 3e 0a 44 61 74 65 | 3a 20 31 37 20 4a 75 6c |CP>.Date|: 17 Jul|
|000000a0| 20 38 37 20 32 32 3a 32 | 36 3a 32 39 20 47 4d 54 | 87 22:2|6:29 GMT|
|000000b0| 0a 53 65 6e 64 65 72 3a | 20 61 6c 6c 62 65 72 79 |.Sender:| allbery|
|000000c0| 40 6e 63 6f 61 73 74 2e | 55 55 43 50 0a 4f 72 67 |@ncoast.|UUCP.Org|
|000000d0| 61 6e 69 7a 61 74 69 6f | 6e 3a 20 43 61 6c 66 6f |anizatio|n: Calfo|
|000000e0| 72 6e 69 61 20 49 6e 73 | 74 69 74 75 74 65 20 6f |rnia Ins|titute o|
|000000f0| 66 20 54 65 63 68 6e 6f | 6c 6f 67 79 0a 4c 69 6e |f Techno|logy.Lin|
|00000100| 65 73 3a 20 32 38 31 0a | 41 70 70 72 6f 76 65 64 |es: 281.|Approved|
|00000110| 3a 20 61 6c 6c 62 65 72 | 79 40 6e 63 6f 61 73 74 |: allber|y@ncoast|
|00000120| 2e 55 55 43 50 0a 58 2d | 41 72 63 68 69 76 65 3a |.UUCP.X-|Archive:|
|00000130| 20 63 6f 6d 70 2e 73 6f | 75 72 63 65 73 2e 6d 69 | comp.so|urces.mi|
|00000140| 73 63 2f 38 37 30 37 2f | 34 38 0a 0a 09 54 68 69 |sc/8707/|48...Thi|
|00000150| 73 20 69 73 20 61 6e 20 | 6e 2d 64 69 6d 65 6e 73 |s is an |n-dimens|
|00000160| 69 6f 6e 61 6c 20 46 46 | 54 20 72 6f 75 74 69 6e |ional FF|T routin|
|00000170| 65 20 49 20 77 72 6f 74 | 65 2e 20 20 54 68 65 72 |e I wrot|e.  Ther|
|00000180| 65 20 61 72 65 20 74 77 | 6f 0a 63 6f 6d 70 69 6c |e are tw|o.compil|
|00000190| 61 74 69 6f 6e 20 6f 70 | 74 69 6f 6e 73 2c 20 53 |ation op|tions, S|
|000001a0| 50 45 45 44 20 61 6e 64 | 20 4d 55 4c 54 34 2c 20 |PEED and| MULT4, |
|000001b0| 77 68 69 63 68 20 61 72 | 65 20 64 65 73 63 72 69 |which ar|e descri|
|000001c0| 62 65 64 20 69 6e 20 74 | 68 65 0a 64 6f 63 75 6d |bed in t|he.docum|
|000001d0| 65 6e 74 61 74 69 6f 6e | 20 66 69 6c 65 2c 20 66 |entation| file, f|
|000001e0| 66 74 6e 2e 64 6f 63 2e | 0a 0a 2d 2d 2d 2d 2d 2d |ftn.doc.|..------|
|000001f0| 2d 2d 2d 2d 2d 2d 2d 2d | 20 43 75 74 20 48 65 72 |--------| Cut Her|
|00000200| 65 20 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 0a |e ------|-------.|
|00000210| 23 21 20 2f 62 69 6e 2f | 73 68 0a 23 20 54 68 69 |#! /bin/|sh.# Thi|
|00000220| 73 20 69 73 20 61 20 73 | 68 65 6c 6c 20 61 72 63 |s is a s|hell arc|
|00000230| 68 69 76 65 2c 20 6d 65 | 61 6e 69 6e 67 3a 0a 23 |hive, me|aning:.#|
|00000240| 20 31 2e 20 52 65 6d 6f | 76 65 20 65 76 65 72 79 | 1. Remo|ve every|
|00000250| 74 68 69 6e 67 20 61 62 | 6f 76 65 20 74 68 65 20 |thing ab|ove the |
|00000260| 23 21 20 2f 62 69 6e 2f | 73 68 20 6c 69 6e 65 2e |#! /bin/|sh line.|
|00000270| 0a 23 20 32 2e 20 53 61 | 76 65 20 74 68 65 20 72 |.# 2. Sa|ve the r|
|00000280| 65 73 75 6c 74 69 6e 67 | 20 74 65 78 74 20 69 6e |esulting| text in|
|00000290| 20 61 20 66 69 6c 65 2e | 0a 23 20 33 2e 20 45 78 | a file.|.# 3. Ex|
|000002a0| 65 63 75 74 65 20 74 68 | 65 20 66 69 6c 65 20 77 |ecute th|e file w|
|000002b0| 69 74 68 20 2f 62 69 6e | 2f 73 68 20 28 6e 6f 74 |ith /bin|/sh (not|
|000002c0| 20 63 73 68 29 20 74 6f | 20 63 72 65 61 74 65 20 | csh) to| create |
|000002d0| 74 68 65 20 66 69 6c 65 | 73 3a 0a 23 09 66 66 74 |the file|s:.#.fft|
|000002e0| 6e 2e 64 6f 63 0a 23 09 | 63 6f 6d 70 6c 65 78 2e |n.doc.#.|complex.|
|000002f0| 68 0a 23 09 66 66 74 6e | 2e 63 0a 23 09 64 65 6d |h.#.fftn|.c.#.dem|
|00000300| 6f 2e 63 0a 23 20 54 68 | 69 73 20 61 72 63 68 69 |o.c.# Th|is archi|
|00000310| 76 65 20 63 72 65 61 74 | 65 64 3a 20 57 65 64 20 |ve creat|ed: Wed |
|00000320| 4a 75 6c 20 31 35 20 32 | 33 3a 31 34 3a 34 36 20 |Jul 15 2|3:14:46 |
|00000330| 31 39 38 37 0a 65 78 70 | 6f 72 74 20 50 41 54 48 |1987.exp|ort PATH|
|00000340| 3b 20 50 41 54 48 3d 2f | 62 69 6e 3a 24 50 41 54 |; PATH=/|bin:$PAT|
|00000350| 48 0a 69 66 20 74 65 73 | 74 20 2d 66 20 27 66 66 |H.if tes|t -f 'ff|
|00000360| 74 6e 2e 64 6f 63 27 0a | 74 68 65 6e 0a 09 65 63 |tn.doc'.|then..ec|
|00000370| 68 6f 20 73 68 61 72 3a | 20 77 69 6c 6c 20 6e 6f |ho shar:| will no|
|00000380| 74 20 6f 76 65 72 2d 77 | 72 69 74 65 20 65 78 69 |t over-w|rite exi|
|00000390| 73 74 69 6e 67 20 66 69 | 6c 65 20 22 27 66 66 74 |sting fi|le "'fft|
|000003a0| 6e 2e 64 6f 63 27 22 0a | 65 6c 73 65 0a 63 61 74 |n.doc'".|else.cat|
|000003b0| 20 3c 3c 20 5c 53 48 41 | 52 5f 45 4f 46 20 3e 20 | << \SHA|R_EOF > |
|000003c0| 27 66 66 74 6e 2e 64 6f | 63 27 0a 09 48 65 72 65 |'fftn.do|c'..Here|
|000003d0| 20 69 73 20 61 6e 20 6e | 2d 64 69 6d 65 6e 73 69 | is an n|-dimensi|
|000003e0| 6f 6e 61 6c 20 46 46 54 | 20 72 6f 75 74 69 6e 65 |onal FFT| routine|
|000003f0| 20 49 20 77 72 6f 74 65 | 20 72 65 63 65 6e 74 6c | I wrote| recentl|
|00000400| 79 2e 20 20 49 74 27 73 | 20 62 61 73 65 64 0a 6f |y.  It's| based.o|
|00000410| 6e 20 74 68 65 20 46 6f | 72 74 72 61 73 68 20 61 |n the Fo|rtrash a|
|00000420| 6e 64 20 50 61 73 63 61 | 6c 20 72 6f 75 74 69 6e |nd Pasca|l routin|
|00000430| 65 73 20 69 6e 20 74 68 | 65 20 62 6f 6f 6b 20 22 |es in th|e book "|
|00000440| 4e 75 6d 65 72 69 63 61 | 6c 20 52 65 63 69 70 65 |Numerica|l Recipe|
|00000450| 73 22 0a 28 49 20 66 6f | 72 67 6f 74 20 74 68 65 |s".(I fo|rgot the|
|00000460| 20 61 75 74 68 6f 72 2c | 20 61 6e 64 20 74 68 65 | author,| and the|
|00000470| 20 62 6f 6f 6b 20 69 73 | 20 61 74 20 77 6f 72 6b | book is| at work|
|00000480| 2e 29 2c 20 62 75 74 20 | 6d 75 63 68 20 6f 66 20 |.), but |much of |
|00000490| 69 74 20 68 61 73 0a 62 | 65 65 6e 20 63 68 61 6e |it has.b|een chan|
|000004a0| 67 65 64 20 28 74 6f 20 | 74 61 6b 65 20 61 64 76 |ged (to |take adv|
|000004b0| 61 6e 74 61 67 65 20 6f | 66 20 43 27 73 20 6c 6f |antage o|f C's lo|
|000004c0| 67 69 63 61 6c 20 6f 70 | 65 72 61 74 69 6f 6e 73 |gical op|erations|
|000004d0| 20 66 6f 72 20 74 68 65 | 20 62 69 74 0a 72 65 76 | for the| bit.rev|
|000004e0| 65 72 73 61 6c 2c 20 74 | 6f 20 75 73 65 20 43 27 |ersal, t|o use C'|
|000004f0| 73 20 61 72 72 61 79 20 | 69 6e 64 65 78 20 6f 72 |s array |index or|
|00000500| 64 65 72 69 6e 67 20 61 | 6e 64 20 73 74 61 72 74 |dering a|nd start|
|00000510| 69 6e 67 20 69 6e 64 65 | 78 20 6f 66 20 30 2c 0a |ing inde|x of 0,.|
|00000520| 65 74 63 2e 29 0a 09 49 | 66 20 74 68 65 20 70 72 |etc.)..I|f the pr|
|00000530| 65 70 72 6f 63 65 73 73 | 6f 72 20 76 61 72 69 61 |eprocess|or varia|
|00000540| 62 6c 65 20 53 50 45 45 | 44 20 69 73 20 64 65 66 |ble SPEE|D is def|
|00000550| 69 6e 65 64 2c 20 64 61 | 74 61 5b 5d 20 6d 75 73 |ined, da|ta[] mus|
|00000560| 74 20 62 65 20 61 6e 0a | 61 72 72 61 79 20 6f 66 |t be an.|array of|
|00000570| 20 66 63 6f 6d 70 6c 65 | 78 2c 20 6e 6f 74 20 63 | fcomple|x, not c|
|00000580| 6f 6d 70 6c 65 78 2c 20 | 77 68 69 63 68 20 75 73 |omplex, |which us|
|00000590| 65 73 20 66 6c 6f 61 74 | 27 73 20 69 6e 73 74 65 |es float|'s inste|
|000005a0| 61 64 20 6f 66 20 64 6f | 75 62 6c 65 27 73 2e 0a |ad of do|uble's..|
|000005b0| 49 20 75 73 65 64 20 74 | 68 65 20 66 6c 6f 61 74 |I used t|he float|
|000005c0| 27 73 20 69 6e 20 74 68 | 65 20 61 70 70 6c 69 63 |'s in th|e applic|
|000005d0| 61 74 69 6f 6e 20 49 20 | 77 72 6f 74 65 20 74 68 |ation I |wrote th|
|000005e0| 69 73 20 66 6f 72 20 74 | 6f 20 72 65 64 75 63 65 |is for t|o reduce|
|000005f0| 20 74 68 65 20 73 69 7a | 65 0a 6f 66 20 74 68 65 | the siz|e.of the|
|00000600| 20 61 72 72 61 79 2c 20 | 61 6e 64 20 74 6f 20 72 | array, |and to r|
|00000610| 65 64 75 63 65 20 74 68 | 65 20 74 69 6d 65 20 74 |educe th|e time t|
|00000620| 61 6b 65 6e 20 62 79 20 | 6d 65 6d 6f 72 79 20 61 |aken by |memory a|
|00000630| 63 63 65 73 73 65 73 2c | 20 73 69 6e 63 65 20 68 |ccesses,| since h|
|00000640| 69 67 68 0a 61 63 63 75 | 72 61 63 79 20 77 61 73 |igh.accu|racy was|
|00000650| 6e 27 74 20 6e 65 65 64 | 65 64 2e 0a 09 54 68 65 |n't need|ed...The|
|00000660| 20 72 6f 75 74 69 6e 65 | 20 68 61 73 20 62 65 65 | routine| has bee|
|00000670| 6e 20 6f 70 74 69 6d 69 | 7a 65 64 20 66 6f 72 20 |n optimi|zed for |
|00000680| 73 70 65 65 64 20 66 6f | 72 20 4d 53 43 20 34 2e |speed fo|r MSC 4.|
|00000690| 30 2c 20 63 6f 6d 70 69 | 6c 69 6e 67 0a 66 6f 72 |0, compi|ling.for|
|000006a0| 20 61 6e 20 38 30 32 38 | 36 2c 20 75 73 69 6e 67 | an 8028|6, using|
|000006b0| 20 69 6e 6c 69 6e 65 20 | 38 30 38 37 2f 32 38 37 | inline |8087/287|
|000006c0| 20 63 6f 64 65 20 28 72 | 65 71 75 69 72 65 73 20 | code (r|equires |
|000006d0| 74 68 65 20 38 37 2f 32 | 38 37 29 2c 20 61 6e 64 |the 87/2|87), and|
|000006e0| 20 69 74 20 69 73 0a 70 | 72 6f 62 61 62 6c 79 20 | it is.p|robably |
|000006f0| 6e 6f 74 20 6f 70 74 69 | 6d 61 6c 20 66 6f 72 20 |not opti|mal for |
|00000700| 6f 74 68 65 72 20 73 79 | 73 74 65 6d 73 2e 0a 09 |other sy|stems...|
|00000710| 49 20 68 61 76 65 20 69 | 6e 63 6c 75 64 65 64 20 |I have i|ncluded |
|00000720| 32 20 6f 70 74 69 6f 6e | 73 20 66 6f 72 20 74 68 |2 option|s for th|
|00000730| 65 20 63 6f 6d 70 6c 65 | 78 20 6d 75 6c 74 69 70 |e comple|x multip|
|00000740| 6c 69 63 61 74 69 6f 6e | 2c 20 61 20 63 6f 6d 70 |lication|, a comp|
|00000750| 6c 65 78 0a 6d 75 6c 74 | 69 70 6c 79 20 75 73 69 |lex.mult|iply usi|
|00000760| 6e 67 20 34 20 72 65 61 | 6c 20 6d 75 6c 74 69 70 |ng 4 rea|l multip|
|00000770| 6c 69 65 73 2c 20 61 6e | 64 20 6f 6e 65 20 75 73 |lies, an|d one us|
|00000780| 69 6e 67 20 33 20 72 65 | 61 6c 20 6d 75 6c 74 69 |ing 3 re|al multi|
|00000790| 70 6c 69 65 73 2e 20 0a | 54 68 65 20 33 20 6d 75 |plies. .|The 3 mu|
|000007a0| 6c 74 69 70 6c 79 20 6d | 65 74 68 6f 64 20 69 73 |ltiply m|ethod is|
|000007b0| 20 75 73 65 64 20 75 6e | 6c 65 73 73 20 74 68 65 | used un|less the|
|000007c0| 20 70 72 65 70 72 6f 63 | 65 73 73 6f 72 20 76 61 | preproc|essor va|
|000007d0| 72 69 61 62 6c 65 20 4d | 55 4c 54 34 20 69 73 0a |riable M|ULT4 is.|
|000007e0| 64 65 66 69 6e 65 64 20 | 28 77 68 69 63 68 20 69 |defined |(which i|
|000007f0| 73 20 61 75 74 6f 6d 61 | 74 69 63 61 6c 6c 79 20 |s automa|tically |
|00000800| 64 6f 6e 65 20 66 6f 72 | 20 4d 53 43 29 2e 20 20 |done for| MSC).  |
|00000810| 55 73 75 61 6c 6c 79 2c | 20 74 68 65 20 33 20 6d |Usually,| the 3 m|
|00000820| 75 6c 74 69 70 6c 79 0a | 6d 65 74 68 6f 64 20 73 |ultiply.|method s|
|00000830| 68 6f 75 6c 64 20 62 65 | 20 66 61 73 74 65 72 2c |hould be| faster,|
|00000840| 20 62 75 74 20 75 6e 66 | 6f 72 74 75 6e 61 74 65 | but unf|ortunate|
|00000850| 6c 79 2c 20 4d 53 43 20 | 34 2e 30 20 64 6f 65 73 |ly, MSC |4.0 does|
|00000860| 20 61 6c 6d 6f 73 74 20 | 6e 6f 0a 6f 70 74 69 6d | almost |no.optim|
|00000870| 69 7a 61 74 69 6f 6e 20 | 66 6f 72 20 74 68 65 20 |ization |for the |
|00000880| 38 30 32 38 37 2c 20 73 | 6f 20 69 74 20 77 61 73 |80287, s|o it was|
|00000890| 20 73 74 6f 72 69 6e 67 | 20 74 31 20 61 6e 64 20 | storing| t1 and |
|000008a0| 74 32 20 62 61 63 6b 20 | 74 6f 20 6d 65 6d 6f 72 |t2 back |to memor|
|000008b0| 79 0a 61 6e 64 20 72 65 | 61 64 69 6e 67 20 74 68 |y.and re|ading th|
|000008c0| 65 6d 20 62 61 63 6b 2c | 20 77 68 69 63 68 20 61 |em back,| which a|
|000008d0| 6c 6c 20 74 6f 6f 6b 20 | 6d 6f 72 65 20 74 69 6d |ll took |more tim|
|000008e0| 65 20 74 68 61 6e 20 64 | 6f 69 6e 67 20 74 68 65 |e than d|oing the|
|000008f0| 20 65 78 74 72 61 20 6d | 75 6c 74 69 70 6c 79 21 | extra m|ultiply!|
|00000900| 0a 0a 09 4e 6f 77 20 66 | 6f 72 20 74 68 65 20 75 |...Now f|or the u|
|00000910| 73 61 67 65 20 28 66 69 | 6e 61 6c 6c 79 29 3a 0a |sage (fi|nally):.|
|00000920| 64 61 74 61 5b 5d 20 69 | 73 20 74 68 65 20 61 72 |data[] i|s the ar|
|00000930| 72 61 79 20 6f 66 20 63 | 6f 6d 70 6c 65 78 20 6e |ray of c|omplex n|
|00000940| 75 6d 62 65 72 73 20 74 | 6f 20 62 65 20 74 72 61 |umbers t|o be tra|
|00000950| 6e 73 66 6f 72 6d 65 64 | 2c 0a 6e 6e 5b 5d 20 69 |nsformed|,.nn[] i|
|00000960| 73 20 74 68 65 20 61 72 | 72 61 79 20 67 69 76 69 |s the ar|ray givi|
|00000970| 6e 67 20 74 68 65 20 64 | 69 6d 65 6e 73 69 6f 6e |ng the d|imension|
|00000980| 73 20 28 49 20 6d 65 61 | 6e 20 73 69 7a 65 29 20 |s (I mea|n size) |
|00000990| 6f 66 20 74 68 65 20 61 | 72 72 61 79 2c 0a 6e 64 |of the a|rray,.nd|
|000009a0| 69 6d 20 69 73 20 74 68 | 65 20 6e 75 6d 62 65 72 |im is th|e number|
|000009b0| 20 6f 66 20 64 69 6d 65 | 6e 73 69 6f 6e 73 20 6f | of dime|nsions o|
|000009c0| 66 20 74 68 65 20 61 72 | 72 61 79 2c 20 61 6e 64 |f the ar|ray, and|
|000009d0| 0a 69 73 69 67 6e 20 69 | 73 20 2b 31 20 66 6f 72 |.isign i|s +1 for|
|000009e0| 20 61 20 66 6f 72 77 61 | 72 64 20 74 72 61 6e 73 | a forwa|rd trans|
|000009f0| 66 6f 72 6d 2c 20 61 6e | 64 20 2d 31 20 66 6f 72 |form, an|d -1 for|
|00000a00| 20 61 6e 20 69 6e 76 65 | 72 73 65 20 74 72 61 6e | an inve|rse tran|
|00000a10| 73 66 6f 72 6d 2e 0a 0a | 64 61 74 61 5b 5d 20 61 |sform...|data[] a|
|00000a20| 6e 64 20 6e 6e 5b 5d 20 | 61 72 65 20 73 74 6f 72 |nd nn[] |are stor|
|00000a30| 65 64 20 69 6e 20 74 68 | 65 20 22 6e 61 74 75 72 |ed in th|e "natur|
|00000a40| 61 6c 22 20 6f 72 64 65 | 72 20 66 6f 72 20 43 3a |al" orde|r for C:|
|00000a50| 0a 6e 6e 5b 30 5d 20 67 | 69 76 65 73 20 74 68 65 |.nn[0] g|ives the|
|00000a60| 20 6e 75 6d 62 65 72 20 | 6f 66 20 65 6c 65 6d 65 | number |of eleme|
|00000a70| 6e 74 73 20 61 6c 6f 6e | 67 20 74 68 65 20 6c 65 |nts alon|g the le|
|00000a80| 66 74 6d 6f 73 74 20 69 | 6e 64 65 78 2c 0a 6e 6e |ftmost i|ndex,.nn|
|00000a90| 5b 6e 64 69 6d 20 2d 20 | 31 5d 20 67 69 76 65 73 |[ndim - |1] gives|
|00000aa0| 20 74 68 65 20 6e 75 6d | 62 65 72 20 6f 66 20 65 | the num|ber of e|
|00000ab0| 6c 65 6d 65 6e 74 73 20 | 61 6c 6f 6e 67 20 74 68 |lements |along th|
|00000ac0| 65 20 72 69 67 68 74 6d | 6f 73 74 20 69 6e 64 65 |e rightm|ost inde|
|00000ad0| 78 2c 20 61 6e 64 0a 64 | 61 74 61 20 73 68 6f 75 |x, and.d|ata shou|
|00000ae0| 6c 64 20 62 65 20 64 65 | 63 6c 61 72 65 64 20 61 |ld be de|clared a|
|00000af0| 6c 6f 6e 67 20 74 68 65 | 20 6c 69 6e 65 73 20 6f |long the| lines o|
|00000b00| 66 0a 09 73 74 72 75 63 | 74 20 28 66 29 63 6f 6d |f..struc|t (f)com|
|00000b10| 70 6c 65 78 20 64 61 74 | 61 5b 6e 6e 5b 30 5d 2c |plex dat|a[nn[0],|
|00000b20| 20 6e 6e 5b 31 5d 2c 20 | 2e 2e 2e 2c 20 6e 6e 5b | nn[1], |..., nn[|
|00000b30| 6e 64 69 6d 20 2d 20 31 | 5d 5d 0a 0a 09 41 64 64 |ndim - 1|]]...Add|
|00000b40| 69 74 69 6f 6e 61 6c 20 | 6e 6f 74 65 73 3a 20 54 |itional |notes: T|
|00000b50| 68 65 20 72 6f 75 74 69 | 6e 65 20 64 6f 65 73 20 |he routi|ne does |
|00000b60| 4e 4f 20 4e 4f 52 4d 41 | 4c 49 5a 41 54 49 4f 4e |NO NORMA|LIZATION|
|00000b70| 2c 20 73 6f 20 69 66 20 | 79 6f 75 20 64 6f 20 61 |, so if |you do a|
|00000b80| 0a 66 6f 72 77 61 72 64 | 2c 20 61 6e 64 20 74 68 |.forward|, and th|
|00000b90| 65 6e 20 61 6e 20 69 6e | 76 65 72 73 65 20 74 72 |en an in|verse tr|
|00000ba0| 61 6e 73 66 6f 72 6d 20 | 6f 6e 20 61 6e 20 61 72 |ansform |on an ar|
|00000bb0| 72 61 79 2c 20 74 68 65 | 20 72 65 73 75 6c 74 20 |ray, the| result |
|00000bc0| 77 69 6c 6c 0a 62 65 20 | 69 64 65 6e 74 69 63 61 |will.be |identica|
|00000bd0| 6c 20 74 6f 20 74 68 65 | 20 6f 72 69 67 69 6e 61 |l to the| origina|
|00000be0| 6c 20 61 72 72 61 79 20 | 4d 55 4c 54 49 50 4c 49 |l array |MULTIPLI|
|00000bf0| 45 44 20 42 59 20 54 48 | 45 20 4e 55 4d 42 45 52 |ED BY TH|E NUMBER|
|00000c00| 20 4f 46 0a 45 4c 45 4d | 45 4e 54 53 20 49 4e 20 | OF.ELEM|ENTS IN |
|00000c10| 54 48 45 20 41 52 52 41 | 59 2e 20 20 41 6c 73 6f |THE ARRA|Y.  Also|
|00000c20| 2c 20 6f 66 20 63 6f 75 | 72 73 65 2c 20 74 68 65 |, of cou|rse, the|
|00000c30| 20 64 69 6d 65 6e 73 69 | 6f 6e 73 20 6f 66 20 64 | dimensi|ons of d|
|00000c40| 61 74 61 5b 5d 0a 6d 75 | 73 74 20 61 6c 6c 20 62 |ata[].mu|st all b|
|00000c50| 65 20 70 6f 77 65 72 73 | 20 6f 66 20 32 2e 20 20 |e powers| of 2.  |
|00000c60| 49 20 68 61 76 65 20 61 | 6c 73 6f 20 65 6e 63 6c |I have a|lso encl|
|00000c70| 6f 73 65 64 20 61 20 76 | 65 72 73 69 6f 6e 20 6f |osed a v|ersion o|
|00000c80| 66 20 74 68 65 20 64 65 | 6d 6f 0a 70 72 6f 67 72 |f the de|mo.progr|
|00000c90| 61 6d 20 66 72 6f 6d 20 | 74 68 65 20 65 78 61 6d |am from |the exam|
|00000ca0| 70 6c 65 73 20 62 6f 6f | 6b 20 66 6f 72 20 22 4e |ples boo|k for "N|
|00000cb0| 75 6d 65 72 69 63 61 6c | 20 52 65 63 69 70 65 73 |umerical| Recipes|
|00000cc0| 22 2e 0a 53 48 41 52 5f | 45 4f 46 0a 66 69 20 23 |"..SHAR_|EOF.fi #|
|00000cd0| 20 65 6e 64 20 6f 66 20 | 6f 76 65 72 77 72 69 74 | end of |overwrit|
|00000ce0| 69 6e 67 20 63 68 65 63 | 6b 0a 69 66 20 74 65 73 |ing chec|k.if tes|
|00000cf0| 74 20 2d 66 20 27 63 6f | 6d 70 6c 65 78 2e 68 27 |t -f 'co|mplex.h'|
|00000d00| 0a 74 68 65 6e 0a 09 65 | 63 68 6f 20 73 68 61 72 |.then..e|cho shar|
|00000d10| 3a 20 77 69 6c 6c 20 6e | 6f 74 20 6f 76 65 72 2d |: will n|ot over-|
|00000d20| 77 72 69 74 65 20 65 78 | 69 73 74 69 6e 67 20 66 |write ex|isting f|
|00000d30| 69 6c 65 20 22 27 63 6f | 6d 70 6c 65 78 2e 68 27 |ile "'co|mplex.h'|
|00000d40| 22 0a 65 6c 73 65 0a 63 | 61 74 20 3c 3c 20 5c 53 |".else.c|at << \S|
|00000d50| 48 41 52 5f 45 4f 46 20 | 3e 20 27 63 6f 6d 70 6c |HAR_EOF |> 'compl|
|00000d60| 65 78 2e 68 27 0a 23 69 | 66 6e 64 65 66 20 44 4f |ex.h'.#i|fndef DO|
|00000d70| 53 09 2f 2a 09 44 65 66 | 69 6e 69 74 69 6f 6e 20 |S./*.Def|inition |
|00000d80| 6f 66 20 63 6f 6d 70 6c | 65 78 20 6e 6f 74 20 6e |of compl|ex not n|
|00000d90| 65 65 64 65 64 20 66 6f | 72 20 4d 53 43 09 2a 2f |eeded fo|r MSC.*/|
|00000da0| 0a 23 69 66 6e 64 65 66 | 20 4d 53 44 4f 53 09 2f |.#ifndef| MSDOS./|
|00000db0| 2a 09 41 6c 72 65 61 64 | 79 20 69 6e 20 3c 6d 61 |*.Alread|y in <ma|
|00000dc0| 74 68 2e 68 3e 09 2a 2f | 0a 0a 73 74 72 75 63 74 |th.h>.*/|..struct|
|00000dd0| 20 63 6f 6d 70 6c 65 78 | 20 7b 0a 09 64 6f 75 62 | complex| {..doub|
|00000de0| 6c 65 20 78 2c 20 79 3b | 0a 7d 3b 0a 0a 23 65 6e |le x, y;|.};..#en|
|00000df0| 64 69 66 0a 0a 73 74 72 | 75 63 74 20 66 63 6f 6d |dif..str|uct fcom|
|00000e00| 70 6c 65 78 20 7b 0a 09 | 66 6c 6f 61 74 20 78 2c |plex {..|float x,|
|00000e10| 20 79 3b 0a 7d 3b 0a 53 | 48 41 52 5f 45 4f 46 0a | y;.};.S|HAR_EOF.|
|00000e20| 66 69 20 23 20 65 6e 64 | 20 6f 66 20 6f 76 65 72 |fi # end| of over|
|00000e30| 77 72 69 74 69 6e 67 20 | 63 68 65 63 6b 0a 69 66 |writing |check.if|
|00000e40| 20 74 65 73 74 20 2d 66 | 20 27 66 66 74 6e 2e 63 | test -f| 'fftn.c|
|00000e50| 27 0a 74 68 65 6e 0a 09 | 65 63 68 6f 20 73 68 61 |'.then..|echo sha|
|00000e60| 72 3a 20 77 69 6c 6c 20 | 6e 6f 74 20 6f 76 65 72 |r: will |not over|
|00000e70| 2d 77 72 69 74 65 20 65 | 78 69 73 74 69 6e 67 20 |-write e|xisting |
|00000e80| 66 69 6c 65 20 22 27 66 | 66 74 6e 2e 63 27 22 0a |file "'f|ftn.c'".|
|00000e90| 65 6c 73 65 0a 63 61 74 | 20 3c 3c 20 5c 53 48 41 |else.cat| << \SHA|
|00000ea0| 52 5f 45 4f 46 20 3e 20 | 27 66 66 74 6e 2e 63 27 |R_EOF > |'fftn.c'|
|00000eb0| 0a 23 69 6e 63 6c 75 64 | 65 20 3c 6d 61 74 68 2e |.#includ|e <math.|
|00000ec0| 68 3e 0a 23 69 6e 63 6c | 75 64 65 20 22 63 6f 6d |h>.#incl|ude "com|
|00000ed0| 70 6c 65 78 2e 68 22 0a | 0a 2f 2a 09 55 73 65 20 |plex.h".|./*.Use |
|00000ee0| 34 2d 6d 75 6c 74 69 70 | 6c 79 20 63 6f 6d 70 6c |4-multip|ly compl|
|00000ef0| 65 78 20 6d 75 6c 74 69 | 70 6c 79 20 66 6f 72 20 |ex multi|ply for |
|00000f00| 4d 53 43 09 2a 2f 0a 23 | 69 66 64 65 66 20 44 4f |MSC.*/.#|ifdef DO|
|00000f10| 53 09 2f 2a 09 49 42 4d | 20 43 09 2a 2f 0a 23 64 |S./*.IBM| C.*/.#d|
|00000f20| 65 66 69 6e 65 20 4d 55 | 4c 54 34 0a 23 65 6e 64 |efine MU|LT4.#end|
|00000f30| 69 66 0a 23 69 66 64 65 | 66 20 4d 53 44 4f 53 09 |if.#ifde|f MSDOS.|
|00000f40| 2f 2a 09 4d 69 63 72 6f | 73 6f 66 74 20 43 09 2a |/*.Micro|soft C.*|
|00000f50| 2f 0a 23 64 65 66 69 6e | 65 20 4d 55 4c 54 34 0a |/.#defin|e MULT4.|
|00000f60| 23 65 6e 64 69 66 0a 0a | 23 64 65 66 69 6e 65 20 |#endif..|#define |
|00000f70| 50 49 20 33 2e 31 34 31 | 35 39 32 36 35 33 35 39 |PI 3.141|59265359|
|00000f80| 37 39 0a 0a 66 66 74 6e | 28 64 61 74 61 2c 20 6e |79..fftn|(data, n|
|00000f90| 6e 2c 20 6e 64 69 6d 2c | 20 69 73 69 67 6e 29 0a |n, ndim,| isign).|
|00000fa0| 0a 23 69 66 64 65 66 20 | 53 50 45 45 44 0a 73 74 |.#ifdef |SPEED.st|
|00000fb0| 72 75 63 74 20 66 63 6f | 6d 70 6c 65 78 20 64 61 |ruct fco|mplex da|
|00000fc0| 74 61 5b 5d 3b 0a 23 65 | 6c 73 65 0a 73 74 72 75 |ta[];.#e|lse.stru|
|00000fd0| 63 74 20 63 6f 6d 70 6c | 65 78 20 64 61 74 61 5b |ct compl|ex data[|
|00000fe0| 5d 3b 0a 23 65 6e 64 69 | 66 0a 0a 75 6e 73 69 67 |];.#endi|f..unsig|
|00000ff0| 6e 65 64 20 6e 6e 5b 5d | 3b 0a 69 6e 74 20 6e 64 |ned nn[]|;.int nd|
|00001000| 69 6d 2c 20 69 73 69 67 | 6e 3b 0a 7b 0a 09 69 6e |im, isig|n;.{..in|
|00001010| 74 20 69 64 69 6d 3b 0a | 09 75 6e 73 69 67 6e 65 |t idim;.|.unsigne|
|00001020| 64 20 69 31 2c 20 69 32 | 72 65 76 2c 20 69 33 72 |d i1, i2|rev, i3r|
|00001030| 65 76 2c 20 69 62 69 74 | 3b 0a 09 75 6e 73 69 67 |ev, ibit|;..unsig|
|00001040| 6e 65 64 20 69 70 32 2c | 20 69 66 70 31 2c 20 69 |ned ip2,| ifp1, i|
|00001050| 66 70 32 2c 20 6b 32 2c | 20 6e 3b 0a 09 75 6e 73 |fp2, k2,| n;..uns|
|00001060| 69 67 6e 65 64 20 6e 70 | 72 65 76 20 3d 20 31 2c |igned np|rev = 1,|
|00001070| 20 6e 72 65 6d 2c 20 6e | 74 6f 74 20 3d 20 31 3b | nrem, n|tot = 1;|
|00001080| 0a 09 72 65 67 69 73 74 | 65 72 20 75 6e 73 69 67 |..regist|er unsig|
|00001090| 6e 65 64 20 69 32 2c 20 | 69 33 3b 0a 09 64 6f 75 |ned i2, |i3;..dou|
|000010a0| 62 6c 65 20 74 68 65 74 | 61 3b 0a 09 73 74 72 75 |ble thet|a;..stru|
|000010b0| 63 74 20 63 6f 6d 70 6c | 65 78 20 77 2c 20 77 70 |ct compl|ex w, wp|
|000010c0| 3b 0a 0a 23 69 66 64 65 | 66 20 53 50 45 45 44 0a |;..#ifde|f SPEED.|
|000010d0| 09 66 6c 6f 61 74 20 77 | 74 65 6d 70 3b 0a 09 73 |.float w|temp;..s|
|000010e0| 74 72 75 63 74 20 66 63 | 6f 6d 70 6c 65 78 20 74 |truct fc|omplex t|
|000010f0| 65 6d 70 2c 20 77 74 3b | 0a 23 65 6c 73 65 0a 09 |emp, wt;|.#else..|
|00001100| 64 6f 75 62 6c 65 20 77 | 74 65 6d 70 3b 0a 09 73 |double w|temp;..s|
|00001110| 74 72 75 63 74 20 63 6f | 6d 70 6c 65 78 20 74 65 |truct co|mplex te|
|00001120| 6d 70 2c 20 77 74 3b 0a | 23 65 6e 64 69 66 0a 0a |mp, wt;.|#endif..|
|00001130| 23 69 66 6e 64 65 66 20 | 4d 55 4c 54 34 09 2f 2a |#ifndef |MULT4./*|
|00001140| 09 54 65 6d 70 6f 72 61 | 72 79 20 76 61 72 69 61 |.Tempora|ry varia|
|00001150| 62 6c 65 73 20 6e 65 65 | 64 65 64 20 66 6f 72 20 |bles nee|ded for |
|00001160| 33 2d 6d 75 6c 74 69 70 | 6c 79 20 63 6f 6d 70 6c |3-multip|ly compl|
|00001170| 65 78 20 6d 75 6c 74 2e | 09 2a 2f 0a 09 64 6f 75 |ex mult.|.*/..dou|
|00001180| 62 6c 65 20 74 31 2c 20 | 74 32 3b 0a 23 65 6e 64 |ble t1, |t2;.#end|
|00001190| 69 66 0a 0a 09 2f 2a 20 | 20 20 20 20 20 43 6f 6d |if.../* |     Com|
|000011a0| 70 75 74 65 20 74 6f 74 | 61 6c 20 6e 75 6d 62 65 |pute tot|al numbe|
|000011b0| 72 20 6f 66 20 63 6f 6d | 70 6c 65 78 20 76 61 6c |r of com|plex val|
|000011c0| 75 65 73 20 20 2a 2f 0a | 09 66 6f 72 20 28 69 64 |ues  */.|.for (id|
|000011d0| 69 6d 20 3d 20 30 3b 20 | 69 64 69 6d 20 3c 20 6e |im = 0; |idim < n|
|000011e0| 64 69 6d 3b 20 2b 2b 69 | 64 69 6d 29 0a 09 09 6e |dim; ++i|dim)...n|
|000011f0| 74 6f 74 20 2a 3d 20 6e | 6e 5b 69 64 69 6d 5d 3b |tot *= n|n[idim];|
|00001200| 0a 0a 09 66 6f 72 20 28 | 69 64 69 6d 20 3d 20 6e |...for (|idim = n|
|00001210| 64 69 6d 20 2d 20 31 3b | 20 69 64 69 6d 20 3e 3d |dim - 1;| idim >=|
|00001220| 20 30 3b 20 2d 2d 69 64 | 69 6d 29 20 7b 0a 09 09 | 0; --id|im) {...|
|00001230| 6e 20 3d 20 6e 6e 5b 69 | 64 69 6d 5d 3b 0a 0a 09 |n = nn[i|dim];...|
|00001240| 09 6e 72 65 6d 20 3d 20 | 6e 74 6f 74 20 2f 20 28 |.nrem = |ntot / (|
|00001250| 6e 20 2a 20 6e 70 72 65 | 76 29 3b 0a 09 09 69 70 |n * npre|v);...ip|
|00001260| 32 20 3d 20 6e 70 72 65 | 76 20 2a 20 6e 3b 20 20 |2 = npre|v * n;  |
|00001270| 20 20 20 20 20 20 2f 2a | 20 20 20 20 20 20 55 6e |      /*|      Un|
|00001280| 69 74 20 73 74 65 70 20 | 66 6f 72 20 6e 65 78 74 |it step |for next|
|00001290| 20 64 69 6d 65 6e 73 69 | 6f 6e 20 2a 2f 0a 09 09 | dimensi|on */...|
|000012a0| 69 32 72 65 76 20 3d 20 | 30 3b 20 20 20 20 20 20 |i2rev = |0;      |
|000012b0| 20 20 20 20 20 20 20 20 | 2f 2a 20 20 20 20 20 20 |        |/*      |
|000012c0| 42 69 74 20 72 65 76 65 | 72 73 65 64 20 69 32 20 |Bit reve|rsed i2 |
|000012d0| 2a 2f 0a 0a 09 09 2f 2a | 20 20 20 20 20 20 54 68 |*/..../*|      Th|
|000012e0| 69 73 20 69 73 20 74 68 | 65 20 62 69 74 20 72 65 |is is th|e bit re|
|000012f0| 76 65 72 73 61 6c 20 73 | 65 63 74 69 6f 6e 20 6f |versal s|ection o|
|00001300| 66 20 74 68 65 20 72 6f | 75 74 69 6e 65 20 2a 2f |f the ro|utine */|
|00001310| 0a 09 09 2f 2a 20 20 20 | 20 20 20 4c 6f 6f 70 20 |.../*   |   Loop |
|00001320| 6f 76 65 72 20 63 75 72 | 72 65 6e 74 20 64 69 6d |over cur|rent dim|
|00001330| 65 6e 73 69 6f 6e 20 20 | 20 20 20 2a 2f 0a 09 09 |ension  |   */...|
|00001340| 66 6f 72 20 28 69 32 20 | 3d 20 30 3b 20 69 32 20 |for (i2 |= 0; i2 |
|00001350| 3c 20 69 70 32 3b 20 69 | 32 20 2b 3d 20 6e 70 72 |< ip2; i|2 += npr|
|00001360| 65 76 29 20 7b 0a 09 09 | 09 69 66 20 28 69 32 20 |ev) {...|.if (i2 |
|00001370| 3c 20 69 32 72 65 76 29 | 0a 09 09 09 09 2f 2a 20 |< i2rev)|...../* |
|00001380| 20 20 20 20 20 4c 6f 6f | 70 20 6f 76 65 72 20 6c |     Loo|p over l|
|00001390| 6f 77 65 72 20 64 69 6d | 65 6e 73 69 6f 6e 73 20 |ower dim|ensions |
|000013a0| 20 20 20 20 20 2a 2f 0a | 09 09 09 09 66 6f 72 20 |     */.|....for |
|000013b0| 28 69 31 20 3d 20 69 32 | 3b 20 69 31 20 3c 20 69 |(i1 = i2|; i1 < i|
|000013c0| 32 20 2b 20 6e 70 72 65 | 76 3b 20 2b 2b 69 31 29 |2 + npre|v; ++i1)|
|000013d0| 0a 09 09 09 09 09 2f 2a | 20 20 20 20 20 20 4c 6f |....../*|      Lo|
|000013e0| 6f 70 20 6f 76 65 72 20 | 68 69 67 68 65 72 20 64 |op over |higher d|
|000013f0| 69 6d 65 6e 73 69 6f 6e | 73 20 20 2a 2f 0a 09 09 |imension|s  */...|
|00001400| 09 09 09 66 6f 72 20 28 | 69 33 20 3d 20 69 31 3b |...for (|i3 = i1;|
|00001410| 20 69 33 20 3c 20 6e 74 | 6f 74 3b 20 69 33 20 2b | i3 < nt|ot; i3 +|
|00001420| 3d 20 69 70 32 29 20 7b | 0a 09 09 09 09 09 09 69 |= ip2) {|.......i|
|00001430| 33 72 65 76 20 3d 20 69 | 33 20 2b 20 69 32 72 65 |3rev = i|3 + i2re|
|00001440| 76 20 2d 20 69 32 3b 0a | 09 09 09 09 09 09 74 65 |v - i2;.|......te|
|00001450| 6d 70 20 3d 20 64 61 74 | 61 5b 69 33 5d 3b 0a 09 |mp = dat|a[i3];..|
|00001460| 09 09 09 09 09 64 61 74 | 61 5b 69 33 5d 20 3d 20 |.....dat|a[i3] = |
|00001470| 64 61 74 61 5b 69 33 72 | 65 76 5d 3b 0a 09 09 09 |data[i3r|ev];....|
|00001480| 09 09 09 64 61 74 61 5b | 69 33 72 65 76 5d 20 3d |...data[|i3rev] =|
|00001490| 20 74 65 6d 70 3b 0a 09 | 09 09 09 09 7d 0a 0a 09 | temp;..|....}...|
|000014a0| 09 09 69 62 69 74 20 3d | 20 69 70 32 3b 0a 09 09 |..ibit =| ip2;...|
|000014b0| 09 2f 2a 20 20 20 20 20 | 20 49 6e 63 72 65 6d 65 |./*     | Increme|
|000014c0| 6e 74 20 66 72 6f 6d 20 | 68 69 67 68 20 65 6e 64 |nt from |high end|
|000014d0| 20 6f 66 20 69 32 72 65 | 76 20 74 6f 20 6c 6f 77 | of i2re|v to low|
|000014e0| 20 2a 2f 0a 09 09 09 64 | 6f 20 7b 0a 09 09 09 09 | */....d|o {.....|
|000014f0| 69 62 69 74 20 3e 3e 3d | 20 31 3b 0a 09 09 09 09 |ibit >>=| 1;.....|
|00001500| 69 32 72 65 76 20 5e 3d | 20 69 62 69 74 3b 0a 09 |i2rev ^=| ibit;..|
|00001510| 09 09 7d 20 77 68 69 6c | 65 20 28 69 62 69 74 20 |..} whil|e (ibit |
|00001520| 3e 3d 20 6e 70 72 65 76 | 20 26 26 20 21 28 69 62 |>= nprev| && !(ib|
|00001530| 69 74 20 26 20 69 32 72 | 65 76 29 29 3b 0a 09 09 |it & i2r|ev));...|
|00001540| 7d 0a 0a 09 09 2f 2a 20 | 20 20 20 20 20 48 65 72 |}..../* |     Her|
|00001550| 65 20 62 65 67 69 6e 73 | 20 74 68 65 20 44 61 6e |e begins| the Dan|
|00001560| 69 65 6c 73 6f 6e 2d 4c | 61 6e 63 7a 6f 73 20 73 |ielson-L|anczos s|
|00001570| 65 63 74 69 6f 6e 20 6f | 66 20 74 68 65 20 72 6f |ection o|f the ro|
|00001580| 75 74 69 6e 65 20 2a 2f | 0a 09 09 2f 2a 20 20 20 |utine */|.../*   |
|00001590| 20 20 20 4c 6f 6f 70 20 | 6f 76 65 72 20 73 74 65 |   Loop |over ste|
|000015a0| 70 20 73 69 7a 65 73 20 | 20 20 20 2a 2f 0a 09 09 |p sizes |   */...|
|000015b0| 66 6f 72 20 28 69 66 70 | 31 20 3d 20 6e 70 72 65 |for (ifp|1 = npre|
|000015c0| 76 3b 20 69 66 70 31 20 | 3c 20 69 70 32 3b 20 69 |v; ifp1 |< ip2; i|
|000015d0| 66 70 31 20 3c 3c 3d 20 | 31 29 20 7b 0a 09 09 09 |fp1 <<= |1) {....|
|000015e0| 69 66 70 32 20 3d 20 69 | 66 70 31 20 3c 3c 20 31 |ifp2 = i|fp1 << 1|
|000015f0| 3b 0a 09 09 09 2f 2a 20 | 20 49 6e 69 74 69 61 6c |;..../* | Initial|
|00001600| 69 7a 65 20 66 6f 72 20 | 74 68 65 20 74 72 69 67 |ize for |the trig|
|00001610| 2e 20 72 65 63 75 72 72 | 65 6e 63 65 20 2a 2f 0a |. recurr|ence */.|
|00001620| 09 09 09 74 68 65 74 61 | 20 3d 20 69 73 69 67 6e |...theta| = isign|
|00001630| 20 2a 20 32 2e 30 20 2a | 20 50 49 20 2f 20 28 69 | * 2.0 *| PI / (i|
|00001640| 66 70 32 20 2f 20 6e 70 | 72 65 76 29 3b 0a 09 09 |fp2 / np|rev);...|
|00001650| 09 77 70 2e 78 20 3d 20 | 73 69 6e 28 30 2e 35 20 |.wp.x = |sin(0.5 |
|00001660| 2a 20 74 68 65 74 61 29 | 3b 0a 09 09 09 77 70 2e |* theta)|;....wp.|
|00001670| 78 20 2a 3d 20 2d 32 2e | 30 20 2a 20 77 70 2e 78 |x *= -2.|0 * wp.x|
|00001680| 3b 0a 09 09 09 77 70 2e | 79 20 3d 20 73 69 6e 28 |;....wp.|y = sin(|
|00001690| 74 68 65 74 61 29 3b 0a | 09 09 09 77 2e 78 20 3d |theta);.|...w.x =|
|000016a0| 20 31 2e 30 3b 0a 09 09 | 09 77 2e 79 20 3d 20 30 | 1.0;...|.w.y = 0|
|000016b0| 2e 30 3b 0a 0a 09 09 09 | 2f 2a 20 20 4c 6f 6f 70 |.0;.....|/*  Loop|
|000016c0| 20 62 79 20 75 6e 69 74 | 20 73 74 65 70 20 69 6e | by unit| step in|
|000016d0| 20 63 75 72 72 65 6e 74 | 20 64 69 6d 65 6e 73 69 | current| dimensi|
|000016e0| 6f 6e 20 20 2a 2f 0a 09 | 09 09 66 6f 72 20 28 69 |on  */..|..for (i|
|000016f0| 33 20 3d 20 30 3b 20 69 | 33 20 3c 20 69 66 70 31 |3 = 0; i|3 < ifp1|
|00001700| 3b 20 69 33 20 2b 3d 20 | 6e 70 72 65 76 29 20 7b |; i3 += |nprev) {|
|00001710| 0a 09 09 09 09 2f 2a 20 | 20 20 20 20 20 4c 6f 6f |...../* |     Loo|
|00001720| 70 20 6f 76 65 72 20 6c | 6f 77 65 72 20 64 69 6d |p over l|ower dim|
|00001730| 65 6e 73 69 6f 6e 73 20 | 20 20 20 20 20 2a 2f 0a |ensions |     */.|
|00001740| 09 09 09 09 66 6f 72 20 | 28 69 31 20 3d 20 69 33 |....for |(i1 = i3|
|00001750| 3b 20 69 31 20 3c 20 69 | 33 20 2b 20 6e 70 72 65 |; i1 < i|3 + npre|
|00001760| 76 3b 20 2b 2b 69 31 29 | 0a 09 09 09 09 09 2f 2a |v; ++i1)|....../*|
|00001770| 20 20 4c 6f 6f 70 20 6f | 76 65 72 20 68 69 67 68 |  Loop o|ver high|
|00001780| 65 72 20 64 69 6d 65 6e | 73 69 6f 6e 73 20 2a 2f |er dimen|sions */|
|00001790| 0a 09 09 09 09 09 66 6f | 72 20 28 69 32 20 3d 20 |......fo|r (i2 = |
|000017a0| 69 31 3b 20 69 32 20 3c | 20 6e 74 6f 74 3b 20 69 |i1; i2 <| ntot; i|
|000017b0| 32 20 2b 3d 20 69 66 70 | 32 29 20 7b 0a 09 09 09 |2 += ifp|2) {....|
|000017c0| 09 09 09 2f 2a 20 20 20 | 20 20 20 44 61 6e 69 65 |.../*   |   Danie|
|000017d0| 6c 73 6f 6e 2d 4c 61 6e | 63 7a 6f 73 20 66 6f 72 |lson-Lan|czos for|
|000017e0| 6d 75 6c 61 20 2a 2f 0a | 09 09 09 09 09 09 6b 32 |mula */.|......k2|
|000017f0| 20 3d 20 69 32 20 2b 20 | 69 66 70 31 3b 0a 09 09 | = i2 + |ifp1;...|
|00001800| 09 09 09 09 77 74 20 3d | 20 64 61 74 61 5b 6b 32 |....wt =| data[k2|
|00001810| 5d 3b 0a 0a 23 69 66 64 | 65 66 20 4d 55 4c 54 34 |];..#ifd|ef MULT4|
|00001820| 0a 2f 2a 09 43 6f 6d 70 | 6c 65 78 20 6d 75 6c 74 |./*.Comp|lex mult|
|00001830| 69 70 6c 79 20 75 73 69 | 6e 67 20 34 20 72 65 61 |iply usi|ng 4 rea|
|00001840| 6c 20 6d 75 6c 74 69 70 | 6c 69 65 73 2e 20 20 46 |l multip|lies.  F|
|00001850| 61 73 74 65 72 20 69 6e | 20 4d 53 43 09 2a 2f 0a |aster in| MSC.*/.|
|00001860| 09 09 09 09 09 09 64 61 | 74 61 5b 6b 32 5d 2e 78 |......da|ta[k2].x|
|00001870| 20 3d 20 64 61 74 61 5b | 69 32 5d 2e 78 20 2d 20 | = data[|i2].x - |
|00001880| 28 74 65 6d 70 2e 78 20 | 3d 20 77 2e 78 20 2a 20 |(temp.x |= w.x * |
|00001890| 77 74 2e 78 20 2d 20 77 | 2e 79 20 2a 20 77 74 2e |wt.x - w|.y * wt.|
|000018a0| 79 29 3b 0a 09 09 09 09 | 09 09 64 61 74 61 5b 6b |y);.....|..data[k|
|000018b0| 32 5d 2e 79 20 3d 20 64 | 61 74 61 5b 69 32 5d 2e |2].y = d|ata[i2].|
|000018c0| 79 20 2d 20 28 74 65 6d | 70 2e 79 20 3d 20 77 2e |y - (tem|p.y = w.|
|000018d0| 78 20 2a 20 77 74 2e 79 | 20 2b 20 77 2e 79 20 2a |x * wt.y| + w.y *|
|000018e0| 20 77 74 2e 78 29 3b 0a | 23 65 6c 73 65 0a 2f 2a | wt.x);.|#else./*|
|000018f0| 09 43 6f 6d 70 6c 65 78 | 20 6d 75 6c 74 69 70 6c |.Complex| multipl|
|00001900| 79 20 75 73 69 6e 67 20 | 33 20 72 65 61 6c 20 6d |y using |3 real m|
|00001910| 75 6c 74 69 70 6c 69 65 | 73 2e 20 20 53 68 6f 75 |ultiplie|s.  Shou|
|00001920| 6c 64 20 75 73 75 61 6c | 6c 79 20 62 65 20 66 61 |ld usual|ly be fa|
|00001930| 73 74 65 72 2e 09 2a 2f | 0a 09 09 09 09 09 09 64 |ster..*/|.......d|
|00001940| 61 74 61 5b 6b 32 5d 2e | 78 20 3d 20 64 61 74 61 |ata[k2].|x = data|
|00001950| 5b 69 32 5d 2e 78 20 2d | 20 28 74 65 6d 70 2e 78 |[i2].x -| (temp.x|
|00001960| 20 3d 0a 09 09 09 09 09 | 09 09 28 74 31 20 3d 20 | =......|..(t1 = |
|00001970| 77 2e 78 20 2a 20 77 74 | 2e 78 29 20 2d 20 28 74 |w.x * wt|.x) - (t|
|00001980| 32 20 3d 20 77 2e 79 20 | 2a 20 77 74 2e 79 29 29 |2 = w.y |* wt.y))|
|00001990| 3b 0a 09 09 09 09 09 09 | 64 61 74 61 5b 6b 32 5d |;.......|data[k2]|
|000019a0| 2e 79 20 3d 20 64 61 74 | 61 5b 69 32 5d 2e 79 20 |.y = dat|a[i2].y |
|000019b0| 2d 20 28 74 65 6d 70 2e | 79 20 3d 0a 09 09 09 09 |- (temp.|y =.....|
|000019c0| 09 09 09 28 77 2e 78 20 | 2b 20 77 2e 79 29 20 2a |...(w.x |+ w.y) *|
|000019d0| 20 28 77 74 2e 78 20 2b | 20 77 74 2e 79 29 20 2d | (wt.x +| wt.y) -|
|000019e0| 20 74 31 20 2d 20 74 32 | 29 3b 0a 23 65 6e 64 69 | t1 - t2|);.#endi|
|000019f0| 66 0a 09 09 09 09 09 09 | 64 61 74 61 5b 69 32 5d |f.......|data[i2]|
|00001a00| 2e 78 20 2b 3d 20 74 65 | 6d 70 2e 78 3b 0a 09 09 |.x += te|mp.x;...|
|00001a10| 09 09 09 09 64 61 74 61 | 5b 69 32 5d 2e 79 20 2b |....data|[i2].y +|
|00001a20| 3d 20 74 65 6d 70 2e 79 | 3b 0a 09 09 09 09 09 7d |= temp.y|;......}|
|00001a30| 0a 09 09 09 09 2f 2a 20 | 20 20 20 20 20 54 72 69 |...../* |     Tri|
|00001a40| 67 6f 6e 6f 6d 65 74 72 | 69 63 20 72 65 63 75 72 |gonometr|ic recur|
|00001a50| 72 65 6e 63 65 20 20 20 | 20 20 20 20 20 2a 2f 0a |rence   |     */.|
|00001a60| 09 09 09 09 77 74 65 6d | 70 20 3d 20 77 2e 78 3b |....wtem|p = w.x;|
|00001a70| 0a 23 69 66 64 65 66 20 | 4d 55 4c 54 34 0a 2f 2a |.#ifdef |MULT4./*|
|00001a80| 09 43 6f 6d 70 6c 65 78 | 20 6d 75 6c 74 69 70 6c |.Complex| multipl|
|00001a90| 79 20 75 73 69 6e 67 20 | 34 20 72 65 61 6c 20 6d |y using |4 real m|
|00001aa0| 75 6c 74 69 70 6c 69 65 | 73 2e 09 2a 2f 0a 09 09 |ultiplie|s..*/...|
|00001ab0| 09 09 77 2e 78 20 2b 3d | 20 77 2e 78 20 2a 20 77 |..w.x +=| w.x * w|
|00001ac0| 70 2e 78 20 2d 20 77 2e | 79 20 2a 20 77 70 2e 79 |p.x - w.|y * wp.y|
|00001ad0| 3b 0a 09 09 09 09 77 2e | 79 20 2b 3d 20 77 74 65 |;.....w.|y += wte|
|00001ae0| 6d 70 20 2a 20 77 70 2e | 79 20 2b 20 77 2e 79 20 |mp * wp.|y + w.y |
|00001af0| 2a 20 77 70 2e 78 3b 0a | 23 65 6c 73 65 0a 2f 2a |* wp.x;.|#else./*|
|00001b00| 09 43 6f 6d 70 6c 65 78 | 20 6d 75 6c 74 69 70 6c |.Complex| multipl|
|00001b10| 79 20 75 73 69 6e 67 20 | 33 20 72 65 61 6c 20 6d |y using |3 real m|
|00001b20| 75 6c 74 69 70 6c 69 65 | 73 2e 09 2a 2f 0a 09 09 |ultiplie|s..*/...|
|00001b30| 09 09 77 2e 78 20 2b 3d | 20 28 74 31 20 3d 20 77 |..w.x +=| (t1 = w|
|00001b40| 2e 78 20 2a 20 77 70 2e | 78 29 20 2d 20 28 74 32 |.x * wp.|x) - (t2|
|00001b50| 20 3d 20 77 2e 79 20 2a | 20 77 70 2e 79 29 3b 0a | = w.y *| wp.y);.|
|00001b60| 09 09 09 09 77 2e 79 20 | 2b 3d 20 28 77 74 65 6d |....w.y |+= (wtem|
|00001b70| 70 20 2b 20 77 2e 79 29 | 20 2a 20 28 77 70 2e 78 |p + w.y)| * (wp.x|
|00001b80| 20 2b 20 77 70 2e 79 29 | 20 2d 20 74 31 20 2d 20 | + wp.y)| - t1 - |
|00001b90| 74 32 3b 0a 23 65 6e 64 | 69 66 0a 09 09 09 7d 0a |t2;.#end|if....}.|
|00001ba0| 09 09 7d 0a 09 6e 70 72 | 65 76 20 2a 3d 20 6e 3b |..}..npr|ev *= n;|
|00001bb0| 0a 09 7d 0a 7d 0a 53 48 | 41 52 5f 45 4f 46 0a 66 |..}.}.SH|AR_EOF.f|
|00001bc0| 69 20 23 20 65 6e 64 20 | 6f 66 20 6f 76 65 72 77 |i # end |of overw|
|00001bd0| 72 69 74 69 6e 67 20 63 | 68 65 63 6b 0a 69 66 20 |riting c|heck.if |
|00001be0| 74 65 73 74 20 2d 66 20 | 27 64 65 6d 6f 2e 63 27 |test -f |'demo.c'|
|00001bf0| 0a 74 68 65 6e 0a 09 65 | 63 68 6f 20 73 68 61 72 |.then..e|cho shar|
|00001c00| 3a 20 77 69 6c 6c 20 6e | 6f 74 20 6f 76 65 72 2d |: will n|ot over-|
|00001c10| 77 72 69 74 65 20 65 78 | 69 73 74 69 6e 67 20 66 |write ex|isting f|
|00001c20| 69 6c 65 20 22 27 64 65 | 6d 6f 2e 63 27 22 0a 65 |ile "'de|mo.c'".e|
|00001c30| 6c 73 65 0a 63 61 74 20 | 3c 3c 20 5c 53 48 41 52 |lse.cat |<< \SHAR|
|00001c40| 5f 45 4f 46 20 3e 20 27 | 64 65 6d 6f 2e 63 27 0a |_EOF > '|demo.c'.|
|00001c50| 23 69 6e 63 6c 75 64 65 | 20 3c 73 74 64 69 6f 2e |#include| <stdio.|
|00001c60| 68 3e 0a 23 69 6e 63 6c | 75 64 65 20 3c 6d 61 74 |h>.#incl|ude <mat|
|00001c70| 68 2e 68 3e 0a 23 69 6e | 63 6c 75 64 65 20 22 63 |h.h>.#in|clude "c|
|00001c80| 6f 6d 70 6c 65 78 2e 68 | 22 0a 0a 23 64 65 66 69 |omplex.h|"..#defi|
|00001c90| 6e 65 20 6e 64 69 6d 20 | 33 0a 23 64 65 66 69 6e |ne ndim |3.#defin|
|00001ca0| 65 20 6e 64 61 74 20 35 | 31 32 0a 0a 6d 61 69 6e |e ndat 5|12..main|
|00001cb0| 28 29 0a 7b 0a 09 69 6e | 74 20 69 2c 20 6a 2c 20 |().{..in|t i, j, |
|00001cc0| 6b 2c 20 6c 2c 20 6c 6c | 2c 20 69 73 69 67 6e 20 |k, l, ll|, isign |
|00001cd0| 3d 20 31 3b 0a 09 69 6e | 74 20 6e 6e 5b 6e 64 69 |= 1;..in|t nn[ndi|
|00001ce0| 6d 5d 3b 0a 0a 23 69 66 | 64 65 66 20 53 50 45 45 |m];..#if|def SPEE|
|00001cf0| 44 0a 09 73 74 72 75 63 | 74 20 66 63 6f 6d 70 6c |D..struc|t fcompl|
|00001d00| 65 78 20 64 61 74 61 5b | 6e 64 61 74 5d 3b 0a 23 |ex data[|ndat];.#|
|00001d10| 65 6c 73 65 0a 09 73 74 | 72 75 63 74 20 63 6f 6d |else..st|ruct com|
|00001d20| 70 6c 65 78 20 64 61 74 | 61 5b 6e 64 61 74 5d 3b |plex dat|a[ndat];|
|00001d30| 0a 23 65 6e 64 69 66 0a | 0a 09 66 6f 72 20 28 69 |.#endif.|..for (i|
|00001d40| 20 3d 20 30 3b 20 69 20 | 3c 20 6e 64 69 6d 3b 20 | = 0; i |< ndim; |
|00001d50| 2b 2b 69 29 0a 09 6e 6e | 5b 69 5d 20 3d 20 34 20 |++i)..nn|[i] = 4 |
|00001d60| 3c 3c 20 69 3b 0a 0a 09 | 66 6f 72 20 28 69 20 3d |<< i;...|for (i =|
|00001d70| 20 30 3b 20 69 20 3c 20 | 6e 6e 5b 30 5d 3b 20 2b | 0; i < |nn[0]; +|
|00001d80| 2b 69 29 0a 09 09 66 6f | 72 20 28 6a 20 3d 20 30 |+i)...fo|r (j = 0|
|00001d90| 3b 20 6a 20 3c 20 6e 6e | 5b 31 5d 3b 20 2b 2b 6a |; j < nn|[1]; ++j|
|00001da0| 29 0a 09 09 09 66 6f 72 | 20 28 6b 20 3d 20 30 3b |)....for| (k = 0;|
|00001db0| 20 6b 20 3c 20 6e 6e 5b | 32 5d 3b 20 2b 2b 6b 29 | k < nn[|2]; ++k)|
|00001dc0| 20 7b 0a 09 09 09 09 6c | 20 3d 20 28 69 20 2a 20 | {.....l| = (i * |
|00001dd0| 6e 6e 5b 31 5d 20 2b 20 | 6a 29 20 2a 20 6e 6e 5b |nn[1] + |j) * nn[|
|00001de0| 32 5d 20 2b 20 6b 3b 0a | 09 09 09 09 6c 6c 20 3d |2] + k;.|....ll =|
|00001df0| 20 28 6c 20 3c 3c 20 31 | 29 20 2b 20 31 3b 0a 09 | (l << 1|) + 1;..|
|00001e00| 09 09 09 64 61 74 61 5b | 6c 5d 2e 78 20 3d 20 6c |...data[|l].x = l|
|00001e10| 6c 3b 0a 09 09 09 09 64 | 61 74 61 5b 6c 5d 2e 79 |l;.....d|ata[l].y|
|00001e20| 20 3d 20 6c 6c 20 2b 20 | 31 3b 0a 09 09 09 7d 0a | = ll + |1;....}.|
|00001e30| 09 66 66 74 6e 28 64 61 | 74 61 2c 20 6e 6e 2c 20 |.fftn(da|ta, nn, |
|00001e40| 6e 64 69 6d 2c 20 69 73 | 69 67 6e 29 3b 0a 09 69 |ndim, is|ign);..i|
|00001e50| 73 69 67 6e 20 3d 20 2d | 31 3b 0a 0a 09 70 72 69 |sign = -|1;...pri|
|00001e60| 6e 74 66 28 22 5c 6e 44 | 6f 75 62 6c 65 20 33 2d |ntf("\nD|ouble 3-|
|00001e70| 64 69 6d 65 6e 73 69 6f | 6e 61 6c 20 54 72 61 6e |dimensio|nal Tran|
|00001e80| 73 66 6f 72 6d 22 29 3b | 0a 09 70 72 69 6e 74 66 |sform");|..printf|
|00001e90| 28 22 5c 6e 44 6f 75 62 | 6c 65 20 54 72 61 6e 73 |("\nDoub|le Trans|
|00001ea0| 66 2e 20 20 4f 72 69 67 | 69 6e 61 6c 20 44 61 74 |f.  Orig|inal Dat|
|00001eb0| 61 20 20 52 61 74 69 6f | 22 29 3b 0a 09 70 72 69 |a  Ratio|");..pri|
|00001ec0| 6e 74 66 28 22 5c 6e 52 | 65 61 6c 20 20 20 49 6d |ntf("\nR|eal   Im|
|00001ed0| 61 67 2e 09 52 65 61 6c | 20 20 20 49 6d 61 67 2e |ag..Real|   Imag.|
|00001ee0| 20 20 20 52 65 61 6c 20 | 20 20 49 6d 61 67 2e 22 |   Real |  Imag."|
|00001ef0| 29 3b 0a 20 0a 09 66 66 | 74 6e 28 64 61 74 61 2c |);. ..ff|tn(data,|
|00001f00| 20 6e 6e 2c 20 6e 64 69 | 6d 2c 20 69 73 69 67 6e | nn, ndi|m, isign|
|00001f10| 29 3b 0a 0a 09 66 6f 72 | 20 28 69 20 3d 20 30 3b |);...for| (i = 0;|
|00001f20| 20 69 20 3c 20 34 3b 20 | 2b 2b 69 29 20 7b 0a 09 | i < 4; |++i) {..|
|00001f30| 09 6a 20 3d 20 69 20 3c | 3c 20 31 3b 0a 09 09 6b |.j = i <|< 1;...k|
|00001f40| 20 3d 20 6a 20 3c 3c 20 | 31 3b 0a 09 09 6c 20 3d | = j << |1;...l =|
|00001f50| 20 28 69 20 2a 20 6e 6e | 5b 31 5d 20 2b 20 6a 29 | (i * nn|[1] + j)|
|00001f60| 20 2a 20 6e 6e 5b 32 5d | 20 2b 20 6b 3b 0a 09 09 | * nn[2]| + k;...|
|00001f70| 6c 6c 20 3d 20 28 6c 20 | 3c 3c 20 31 29 20 2b 20 |ll = (l |<< 1) + |
|00001f80| 31 3b 0a 09 09 70 72 69 | 6e 74 66 28 22 5c 6e 25 |1;...pri|ntf("\n%|
|00001f90| 2e 32 6c 66 20 25 2e 32 | 6c 66 20 25 64 20 25 64 |.2lf %.2|lf %d %d|
|00001fa0| 20 25 2e 32 6c 66 20 25 | 2e 32 6c 66 22 2c 20 64 | %.2lf %|.2lf", d|
|00001fb0| 61 74 61 5b 6c 5d 2e 78 | 2c 20 64 61 74 61 5b 6c |ata[l].x|, data[l|
|00001fc0| 5d 2e 79 2c 20 6c 6c 2c | 20 6c 6c 2b 31 2c 20 64 |].y, ll,| ll+1, d|
|00001fd0| 61 74 61 5b 6c 5d 2e 78 | 20 2f 20 6c 6c 2c 20 64 |ata[l].x| / ll, d|
|00001fe0| 61 74 61 5b 6c 5d 2e 79 | 20 2f 20 28 6c 6c 2b 31 |ata[l].y| / (ll+1|
|00001ff0| 29 29 3b 0a 09 7d 0a 0a | 09 70 72 69 6e 74 66 28 |));..}..|.printf(|
|00002000| 22 5c 6e 54 68 65 20 70 | 72 6f 64 75 63 74 20 6f |"\nThe p|roduct o|
|00002010| 66 20 74 72 61 6e 73 66 | 6f 72 6d 20 6c 65 6e 67 |f transf|orm leng|
|00002020| 74 68 73 20 69 73 3a 20 | 25 64 5c 6e 22 2c 20 6e |ths is: |%d\n", n|
|00002030| 6e 5b 30 5d 20 2a 20 6e | 6e 5b 31 5d 20 2a 20 6e |n[0] * n|n[1] * n|
|00002040| 6e 5b 32 5d 29 3b 0a 7d | 0a 53 48 41 52 5f 45 4f |n[2]);.}|.SHAR_EO|
|00002050| 46 0a 66 69 20 23 20 65 | 6e 64 20 6f 66 20 6f 76 |F.fi # e|nd of ov|
|00002060| 65 72 77 72 69 74 69 6e | 67 20 63 68 65 63 6b 0a |erwritin|g check.|
|00002070| 23 09 45 6e 64 20 6f 66 | 20 73 68 65 6c 6c 20 61 |#.End of| shell a|
|00002080| 72 63 68 69 76 65 0a 65 | 78 69 74 20 30 0a 0a 09 |rchive.e|xit 0...|
|00002090| 41 6e 64 72 65 77 20 4c | 65 65 0a 09 61 6c 65 65 |Andrew L|ee..alee|
|000020a0| 40 63 61 6c 74 65 63 68 | 2e 62 69 74 6e 65 74 2c |@caltech|.bitnet,|
|000020b0| 20 61 6c 65 65 40 74 79 | 62 61 6c 74 2e 63 61 6c | alee@ty|balt.cal|
|000020c0| 74 65 63 68 2e 65 64 75 | 2c 0a 09 2e 2e 2e 21 73 |tech.edu|,.....!s|
|000020d0| 65 69 73 6d 6f 21 63 69 | 74 2d 76 61 78 21 74 79 |eismo!ci|t-vax!ty|
|000020e0| 62 61 6c 74 21 61 6c 65 | 65 0a                   |balt!ale|e.      |
+--------+-------------------------+-------------------------+--------+--------+