WFS update, and added function for rectangular setups in the main

Increased the readability of mainWFS.py.
Normalized the intensity distribution and the lensletCentx/y.
Added support length for FFT in the parameters for the sensors.
Added function in main to handle rectangular lenslet configurations.
 
But isn't the range  independent from the size of the pixels.
If you have a screen of 2 meters and only 2 pixels, then the length of
every pixel is 1 meter.  So I can't understand why you write dxo =
lx/(Nx-1.0). Or do dx and dxo representing something different than the
sample length?


On 15 June 2014 12:36, Nick Kant <notifications@github.com> wrote:

> In src/WFS/mainWFS.py:
>
> >  	k = 2*pi/lam # Wavenumber
> >  	
> > -	# Create grid in the pupil plane
> > -	dxo = lx/(Nx - 1.0) # Sample length on x-axis [m]
> > -	dyo = ly/(Ny - 1.0) # Sample length on y-axis [m]
> > -	xo = arange(0.0,lx + dxo,dxo) # Sample positions on x-axis [m]
> > -	yo = arange(0.0,ly + dyo,dyo) # Sample positions on y-axis [m]
> > -	Xo, Yo = meshgrid(xo,yo) # Create spatial grid
> > +	# Create grid for in the focal plane
> > +	dx = lx/(Nx - 1.0) # Sample length on x-axis [m]
>
> It is necessary to overspan te range of 0 to and with lx with Nx samples.
> So not doing it would make the range from dx to and with lx or 0 to lx.
>
> —
> Reply to this email directly or view it on GitHub
> <https://github.com/faitdivers/pyao/pull/78/files#r13785392>.
>