[SciPy-user] phase wrapping algorithm
ryanlists at gmail.com
Tue Oct 10 13:57:09 CDT 2006
For this simple problem, my FORTRAN version is running about 30 times
faster, but neither takes very long:
Python vector time=0.000700950622559
I will probably finish the FORTRAN version since it is almost done.
On 10/10/06, Travis Oliphant <oliphant.travis at ieee.org> wrote:
> Ryan Krauss wrote:
> > Thanks Travis.
> > I looked in scipy.signal and when I didn't see, I just assumed it didn't exist.
> > I was just starting to make good progess on my FORTRAN version too.....
> The FORTRAN version would be faster and not require memory use. The
> NumPy version creates the difference vector, finds discontinuites and
> forms a correction vector that is cumsummed and added to the original.
> > I may need to revisit this when I need more control over messy
> > experimental data. Is there a tricky in-place way to calculate the
> > difference between a vector and itself shift by 1 or more elements?
> c[:-1] -= c[1:] # shift-by-one
> c[:-2] -= c[2:] # shift-by-two
> This will leave un-changed the last elements of the array.
> SciPy-user mailing list
> SciPy-user at scipy.org
More information about the SciPy-user