Sponsored Content
Full Discussion: Slow FFT in ksh93 and awk.
Operating Systems OS X (Apple) Slow FFT in ksh93 and awk. Post 303023171 by wisecracker on Wednesday 12th of September 2018 04:21:48 PM
Old 09-12-2018
Pure awk FFT.

Hi Corona688...

I just knew you would create a translation and get it working on awk. ;o)
You could remove the SIN_ANGLE and COS_ANGLE and put sin(ANGLE) and cos(ANGLE) directly into the two lines below them, eliminating two lines...

Thanks for the translation, I will see if I can eliminate the Python version in AudioScope.sh and attach your name to it...

Another thing is that there is no need to issue the imaginary values for this basic FFT. As long as the real values are actually, cropped or padded to power(s) of 2 then the imaginary values of 0.0 can be internally generated to the same array length as the real values, eliminating the IMAGSTR initialisation on the command line call.

But nice to see an awk version that I can steal... <wink>

Last edited by wisecracker; 09-14-2018 at 07:45 AM.. Reason: Add title...
 

10 More Discussions You Might Find Interesting

1. UNIX for Dummies Questions & Answers

echo is too slow. HELP with Awk

Hello All, Below is a simple script i worte to find the 208th char in a file. If the char = "C" then I re-direct the line to a file called change.txt. If it is not "C" then I re-direct it to a file called delete.txt. My problem is I have a file 0f 500K lines. this script is very slow. I am... (4 Replies)
Discussion started by: eja
4 Replies

2. Shell Programming and Scripting

ksh93 deprecation...

Any means of running ksh93 in a ksh88-mode? Might sound odd, but I want/need to restrict U/Win-developed scripts to correspond to the ksh88 version on my Solaris environment(s). Thanks. (2 Replies)
Discussion started by: curleb
2 Replies

3. Shell Programming and Scripting

ksh88 or ksh93

Hi all! Does anybody know how can I check if any UNIX installation has implemented ksh88 or ksh93? Thanks in advance. Néstor. (3 Replies)
Discussion started by: Nestor
3 Replies

4. Shell Programming and Scripting

The builtin split function in AWK is too slow

I have a text file that contains 4 million lines, each line contains 2 fields(colon as field separator). as shown: 123:444,555,666,777,888,345 233:5444,555,666,777,888,345 623:454,585,664,773,888,345 ...... Here I have to split the second field(can be up to 40,000 fields) by comma into an... (14 Replies)
Discussion started by: kevintse
14 Replies

5. Shell Programming and Scripting

cut, sed, awk too slow to retrieve line - other options?

Hi, I have a script that, basically, has two input files of this type: file1 key1=value1_1_1 key2=value1_2_1 key4=value1_4_1 ... file2 key2=value2_2_1 key2=value2_2_2 key3=value2_3_1 key4=value2_4_1 ... My files are 10k lines big each (approx). The keys are strings that don't... (7 Replies)
Discussion started by: fzd
7 Replies

6. Shell Programming and Scripting

Making a faster alternative to a slow awk command

Hi, I have a large number of input files with two columns of numbers. For example: 83 1453 99 3255 99 8482 99 7372 83 175 I only wish to retain lines where the numbers fullfil two requirements. E.g: =83 1000<=<=2000 To do this I use the following... (10 Replies)
Discussion started by: s052866
10 Replies

7. UNIX for Advanced & Expert Users

Ksh93 on Linux compatible with ksh93 on AIX

Hi Experts, I have several shell scripts that have been developed on a Linux box for korn ksh93. If we want to run this software on an AIX 6.1 box that runs ksh88 by default can we just change the she-bang line to reference /bin/ksh93 which ships with AIX as its "enhanced shell" to ensure... (6 Replies)
Discussion started by: Keith Turley
6 Replies

8. OS X (Apple)

FFT for the AMIGA through ksh88 shell.

I don't know if anyone is interested but I have been meddling with FFT for the AMIGA. (Sadly we AMIGAns don't have these luxuries through any scripting language. Below is a Python snippet that uses the builtin 'cmath' module to work with the lowly Python 2.0.1 for the AMIGA. It is part of a... (0 Replies)
Discussion started by: wisecracker
0 Replies

9. AIX

Ksh93/AIX compatibility

Hi everyone ! Im trying to know from wich version of AIX KSH93 is available ? Internet tell me 6.x and 7.x AIX are available, bue what about 5.x ? Is KSH93 available on AIX 5.x ? Is it the same way to manipulate variables as KSH93 on 7.x ? Thanks for your support and have a nice day ! (2 Replies)
Discussion started by: majinfrede
2 Replies

10. OS X (Apple)

FFT for Python 2.0.x to 3.7.0.

Hi guys... This is code that was originally designed to work on an upgraded AMIGA A1200 using Python 2.0.x. Unfortunately it broke inside much later versions, NOT because of the print statement/function but other minor subtleties. So this is the final result tested on various machines including... (0 Replies)
Discussion started by: wisecracker
0 Replies
mlib_SignalFFT_4(3MLIB) 				    mediaLib Library Functions					   mlib_SignalFFT_4(3MLIB)

NAME
mlib_SignalFFT_4, mlib_SignalFFT_4_S16_S16, mlib_SignalFFT_4_S16C_S16C, mlib_SignalFFT_4_S16C_S16, mlib_SignalFFT_4_S16, mlib_Sig- nalFFT_4_S16C - signal Fast Fourier Transform (FFT) SYNOPSIS
cc [ flag... ] file... -lmlib [ library... ] #include <mlib.h> mlib_status mlib_SignalFFT_4_S16_S16(mlib_s16 *dstr, mlib_s16 *dsti, const mlib_s16 *srcr, const mlib_s16 *srci, mlib_s32 order, mlib_s32 *scale); mlib_status mlib_SignalFFT_4_S16C_S16C(mlib_s16 *dstc, const mlib_s16 *srcc, mlib_s32 order, mlib_s32 *scale); mlib_status mlib_SignalFFT_4_S16C_S16(mlib_s16 *dstc, const mlib_s16 *srcr, mlib_s32 order, mlib_s32 *scale); mlib_status mlib_SignalFFT_4_S16(mlib_s16 *srcdstr, mlib_s16 *srcdsti, mlib_s32 order, mlib_s32 *scale); mlib_status mlib_SignalFFT_4_S16C(mlib_s16 *srcdstc, mlib_s32 order, mlib_s32 *scale); DESCRIPTION
Each of the functions in this group performs Fast Fourier Transform (FFT). The following equation is used for forward FFT: 1 N-1 dst[k] = ---- SUM {src[n] * exp(-j2*PI*n*k/N)} C1 n=0 and the following equation is used for inverse FFT (IFFT): 1 N-1 dst[n] = ---- SUM {src[k] * exp(j2*PI*n*k/N)} C2 k=0 where k = 0, 1, ..., (N - 1) n = 0, 1, ..., (N - 1) N = 2**order The signal FFT/IFFT functions can be categorized into four groups according to the ScaleMode in the function names in the following form: mlib_Signal[FFT|IFFT]_ScaleMode_OutType_InType_OpMode() mlib_Signal[FFT|IFFT]_ScaleMode_DataType_OpMode() The scaling factors C1 and C2 used in the equations are defined as follows: o For ScaleMode = 1, C1 = 1 and C2 = 2**order. o For ScaleMode = 2, C1 = 2**order and C2 = 1. o For ScaleMode = 3, C1 = C2 = 2**(order/2) when order is even, or C1 = 2**((order+1)/2) and C2 = 2**((order-1)/2) when order is odd. o For ScaleMode = 4, C1 = 2**P and C2 = 2**Q, where P and Q are adaptive scaling factors and are generated by the functions. For functions with only real parts for the source signal, the imaginary parts are assumed to be all zero. For functions with only real parts for the destination signal, the imaginary parts are discarded. The functions with only one data type in their names perform the oper- ation in place. PARAMETERS
Each function takes some of the following arguments: dstr Destination signal array that contains the real parts. dsti Destination signal array that contains the imaginary parts. srcr Source signal array that contains the real parts. srci Source signal array that contains the imaginary parts. dstc Complex destination signal array. dstc[2*i] contains the real parts, and dstc[2*i+1] contains the imaginary parts. srcc Complex source signal array. srcc[2*i] contains the real parts, and srcc[2*i+1] contains the imaginary parts. srcdstr Source and destination signal array that contains the real parts. srcdsti Source and destination signal array that contains the imaginary parts. srcdstc Complex source and destination signal array. srcdstc[2*i] contains the real parts, and srcdstc[2*i+1] contains the imaginary parts. order Order of the transformation. The base-2 logarithm of the number of data samples. scale Adaptive scaling factor. RETURN VALUES
The function returns MLIB_SUCCESS if successful. Otherwise it returns MLIB_FAILURE. ATTRIBUTES
See attributes(5) for descriptions of the following attributes: +-----------------------------+-----------------------------+ | ATTRIBUTE TYPE | ATTRIBUTE VALUE | +-----------------------------+-----------------------------+ |Interface Stability |Committed | +-----------------------------+-----------------------------+ |MT-Level |MT-Safe | +-----------------------------+-----------------------------+ SEE ALSO
mlib_SignalFFT_1(3MLIB), mlib_SignalFFT_2(3MLIB), mlib_SignalFFT_3(3MLIB), mlib_SignalIFFT_1(3MLIB), mlib_SignalIFFT_2(3MLIB), mlib_Sig- nalIFFT_3(3MLIB), mlib_SignalIFFT_4(3MLIB), attributes(5) SunOS 5.11 2 Mar 2007 mlib_SignalFFT_4(3MLIB)
All times are GMT -4. The time now is 10:21 AM.
Unix & Linux Forums Content Copyright 1993-2022. All Rights Reserved.
Privacy Policy