I'm writing some polyphase resampling <http://nbviewer.ipython.org/github/jaykickliter/Notebooks/blob/master/Polyphase.ipynb> code. Sometimes there are symmetries that cut down on the number of required multiplies. I'm still learning how to use Julia's metaprogramming facilities, but assuming it's possible to recognize the symmetry and generate some optimized code for a particular instance of a filter object, how would I execute it? The filters will be composite types that store previous state, filter taps, resampling rations, and other information. Is it possible to store the code in the object itself, and execute it efficiently?
This more of a curiosity. If I could get the code to actually compile to vectorized SIMD code, the "optimization" might be slower. Thanks!
