mò rß@c@sQdZdZdZyeWnd\ZZnXy(dkZdkZdkZdkZWn@e j o4yei i deƒWn eGHnXe ‚nXei di eiƒdZei ƒdZead„Zd „Zd kTdkZdkZd Zd „Zd „Zd„Zd„Zd„Zd„Zedjo eƒndS(sopackage of variants on fourier transform, including power spectrum see descriptions of individual functionss% power - variants on fourier transform Copyright (C) 2004 John Comeau This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. s5Not all needed libraries found, upgrade or check pathiiNs%s iÿÿÿÿcGs)totiidt|ƒƒndS(Ns%s (t debuggingtsyststderrtwritetreprtwhatever(R((t/home/jcomeau/src/power.pyt DebugPrint0scGsdGHt|ƒdS(s6check if lowercase functions actually work as programss=this is a test. it is only a test. passed arguments follow...N(Rtargs(R((Rttest5s(t*s:this string has to be here, see template.py for reason whycCsËt|ƒ}x‘t|dƒtijp1t|dƒtijpt|dƒtjoEt|dƒ}t|ƒdjot|dƒ}n|d}qWt |ƒ}t t |ƒ}t |ƒ|fS(sreturn a numeric array and size... whether it was passed AS an array, or another sequence type, or even as a list of strings like .12, .22, .65, etc. if the first element is an array, then the second element will be taken as the size parameter to be passed to fft()iiN(tlenRtsizettypettypestListTypet TupleTypet ArrayTypetinttlisttdatatmaptfloattarray(RRR ((RtgetarrayAs K cGs t|ƒS(sfor command-line testingN(tfftR(R((RtFftVscGs"t|ƒ\}}ti||ƒS(s2pydoc FFT.fft, this is just an interface into thatN(RRRR tFFTR(RRR ((RRZscGs t|ƒS(sfor command-line testingN(t powerspectrumR(R((Rt PowerSpectrumascGsIt|ƒ\}}t||ƒ}t|t|ƒƒt|ƒd}|S(srPower spectrum code from University of Cape Town website http://www.phy.uct.ac.za/courses/phy400w/cp/pyfrags.htmiN( RRRR Rt transformtabst conjugateR tspectrum(RRR!RR ((RRes $cCstiditiƒd}|idƒd}y­td|ƒti jo‚t i dƒi |ƒd jotd|ƒqÜyttd|ƒƒGHWqÜtiid|ttidƒfƒ‚qÜXntd |‚Wn¬tj o }y|id}Wn d }nX|d |jo t|‚nttiƒdjotiidƒtƒqŒtiid ttidƒdfƒnXd S( sÀmain routine, only used for command-line testing let's say this program is called 'template.py' don't, by the way, name the program with the name of an included function if you're using this system... now if the program is invoked as './template.py', it will simply output a usage message. if invoked as './template.py MyFunction 1 2 3' we will attempt to call MyFunction(1, 2, 3) and show whatever is found if invoked as './template.py myfunction 1 2 3' we attempt to call myfunction(1, 2, 3) as a program and let it output whatever it decides likewise if you symlink like so: 'ln -s template.py myfunction' then invoke as './myfunction 1 2 3'iiÿÿÿÿt.stype(%s)s ^[a-z0-9]+$s%s(sys.argv[1:])s&Error occurred testing routine %s(%s) is%s not a functionsunknown NameErrors_Usage: %s [OPTIONS] INPUT[...] or for more detailed help, read source or: pydoc %s N(RtargvtsplittostseptprogramtfunctiontevalRt FunctionTypetretcompiletmatchtNoneRRRt NameErrortinstanceRtmessageR tpoptmaintself(R(R0R'R1((RR3ns0'    t__main__(ii(t__doc__t Copyrightt errormessagetTruetFalseRR%RR+t exceptionRRR#R$R&R4t originalselfRRR tNumericRtpowerRRRRRR3t__name__(R8R7RRRRR4R+R R3R9R>R<RRRR:RRR%((Rt?s<(         2