On Mar 20, 2013, at 7:12 PM, Barnaby Rowe <notifications@github.com> wrote:

> In galsim/lensing.py:
> 
> > -        kx = self.kx
> > -        ky = self.ky
> > -        TwoPsi=2*np.arctan2(1.0*self.ky, 1.0*self.kx)
> > -        C[kx,ky]=np.cos(TwoPsi)
> > -        S[kx,ky]=np.sin(TwoPsi)
> > -        C[-kx,ky]=-np.cos(TwoPsi)
> > -        S[-kx,ky]=np.sin(TwoPsi)
> > +        # Internal function to generate the cosine and sine spin weightings for the current array
> > +        # set-up.
> > +        C = np.zeros((self.nx, self.ny / 2 + 1))
> > +        S = np.zeros((self.nx, self.ny / 2 + 1))
> > +        i_kx = self.i_kx
> > +        i_ky = self.i_ky
> > +        TwoPsi = 2 * np.arctan2(1.0*self.i_ky, 1.0*self.i_kx)
> > +        C[ i_kx, i_ky] = np.cos(TwoPsi)
> > +        S[ i_kx, i_ky] = np.sin(TwoPsi)
> I agree, I'll push a change along those lines unless anyone objects?
> 
No objection