SymPy has some routines to make formulas more palatable. For instance, it can print sympy.Expr objects (expressions) in LATEX:
The second argument of solve() indicates the set of “output” variables. Indeed, we have three equations for twelve variables. Each equation can be used to express one variable as function of the others. Thus, we can pick three variables and express them as functions of the remaining nine.
in IPython qtconsole (which one could embed inside a custom PyQt/PySide application, and keep track of all the defined symbols
T. Cooklev has published over 20 journal papers and over 30 conference papers in the areas of digital signal processing, audio and video compression, fast algorithms, and wireless communications. For one of his papers he received the Best Paper Award at the 1994 IEEE Asia-Pacific Conference on Circuits and Systems, Taipei, Taiwan. He is the inventor on 5 issued and about 12 pending patents in the United States.
Tom Holton, Ph.D.Tom Holton Electrical/Computer Engineering Program Head &Professor of Electrical EngineeringDiscipline: DigitalOffice: SCI 170Office Hours: Click herePhone: (415) 338-1529Fax: (415) 338-0525Email: firstname.lastname@example.org
I know this doesn’t address your problem directly but typical technique is to divide into 2^n samples chunks and process; possibly with overlapping blocks, possibly applying a window function (Google it) depending on desired frequency response. If you are modifying the FFT and applying an inverse you will want overlapping blocks, cross-faded in output, because you will get audible clicks between blocks if the apparent phase (or the 0Hz constant term) changes. BTW, frequency (Hz) is index * sample_rate / block_size