| """ | |
| Shows how to do a cross join (i.e. cartesian product) between two pandas DataFrames using an example on | |
| calculating the distances between origin and destination cities. | |
| Tested with pandas 0.17.1 and 0.18 on Python 3.4 and Python 3.5 | |
| Best run this with Spyder (see https://github.com/spyder-ide/spyder) | |
| Author: Markus Konrad <post@mkonrad.net> | |
| April 2016 |
| import numpy as np | |
| import fftw3 | |
| import pyfftw | |
| import multiprocessing | |
| import matplotlib | |
| import matplotlib.pyplot as pl | |
| import time | |
| def fft_comparison_tests(size=2048, dtype=np.complex128, byte_align=False): |
| function x = istft(d, ftsize, w, h) | |
| % X = istft(D, F, W, H) Inverse short-time Fourier transform. | |
| % Performs overlap-add resynthesis from the short-time Fourier transform | |
| % data in D. Each column of D is taken as the result of an F-point | |
| % fft; each successive frame was offset by H points (default | |
| % W/2, or F/2 if W==0). Data is hann-windowed at W pts, or | |
| % W = 0 gives a rectangular window (default); | |
| % W as a vector uses that as window. | |
| % This version scales the output so the loop gain is 1.0 for | |
| % either hann-win an-syn with 25% overlap, or hann-win on |
Empirically testing whether FWT or FFT is faster. Disclaimer: I don't understand wavelets.
Speed of FFT implementation doesn't care what the data is:
In [23]: a = rand(2**20)
In [24]: timeit fft(a)
10 loops, best of 3: 128 ms per loop
In [25]: a = rand(2**20)
| # -*- coding: utf-8 -*- | |
| """ | |
| Created on Wed Feb 03 2016 | |
| Show how the time derivative of phase spectrogram is constant, shows the | |
| frequency offset within a bin | |
| """ | |
| from numpy import sin, linspace, pi, abs, angle, diff, unwrap | |
| from numpy.random import randn |
| #!/bin/bash | |
| # Example Script to forward ssh-agent socket into an unprivileged lxc | |
| # This allows access to the ssh-agent by commands executed with lxc-attach | |
| set -m | |
| CONTAINER=$1 # e.g., test-container | |
| REMOTE=$2 # e.g, user@example.com |
| import subprocess | |
| import tempfile | |
| import time | |
| from IPython.core.magic import register_line_cell_magic | |
| def snakeviz(line, cell=None): | |
| """ | |
| Profile code and display the profile in Snakeviz. |