Re: [Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT)

2008-08-07 Thread Matthias Hillenbrand
Hello, today I found time to do concentrate on the proposed zero padding. My Gaussian beam and the lenses have a diameter of approximately 2^16 array elements while the total array has a size of 2^18. The absolute value of the Gaussian beam multiplied by the lenses converges to zero. Out of this

Re: [Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT)

2008-08-06 Thread Matthias Hillenbrand
Hello, > When you convolve two signals, of lengths N and M, you need to pad the > FFTs to length (N+M-1) before multiplication. > You can take a look at my linear position-invariant filtering code at: > http://mentat.za.net/hg/filter I understand your comments that I need zero padding when doin

Re: [Numpy-discussion] Calculating roots with negative numbers

2008-08-06 Thread Matthias Hillenbrand
Sorry, I have chosen the wrong title for this email. It should be a follow up to the discussion 'Horizontal lines in diffraction image (NumPy FFT)' ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/nu

Re: [Numpy-discussion] Calculating roots with negative numbers

2008-08-06 Thread Matthias Hillenbrand
Hello, > When you convolve two signals, of lengths N and M, you need to pad the > FFTs to length (N+M-1) before multiplication. > You can take a look at my linear position-invariant filtering code at: > http://mentat.za.net/hg/filter I understand your comments that I need zero padding when doi

[Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT)

2008-08-05 Thread Matthias Hillenbrand
Hello, I am trying to calculate the propagation of a Gaussian beam with an angular spectrum propagation algorithm. My program does the following steps: 1. Generation and multiplication of the Gaussian beam, an aperture, and a lens -> u 2. FFT(u) 3. Multiplication of FFT(u) with the angular spectru

Re: [Numpy-discussion] Calculating roots with negative numbers

2008-07-15 Thread Matthias Hillenbrand
Thanks Robert, that is exactely what I was looking for. Matthias ___ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion

[Numpy-discussion] Calculating roots with negative numbers

2008-07-15 Thread Matthias Hillenbrand
Hello, I want to calculate the root of a numpy array with negative numbers. Here is an example: x = linspace(-10,10,100) h = zeros(100) h[where(abs(x) < 2)] = sqrt( abs(x) -2 ) h[where(2 <= abs(x))] = 1j * sqrt( 2 - abs(x) ) Unfortunately I get the following error: Warning: invalid value