On Mar 17, 2013, at 7:31 PM, Mike Jarvis <notifications@github.com> wrote:

> In tests/lensing_reference_data/shearfield_ref.py:
> 
> > @@ -32,14 +32,15 @@
> >  dx = 1.
> >  
> >  # define power spectrum
> > -ps = galsim.lensing.PowerSpectrum(e_power_function="k**2", b_power_function="k")
> > -# get shears
> > -g1, g2 = ps.buildGriddedShears(grid_spacing = dx, ngrid = n, rng=rng)
> > +ps = galsim.lensing.PowerSpectrum(e_power_function="k**0.5", b_power_function="k")
> Another "just curious" question. Why the switch from k**2 to k**0.5?
> 
Purely used as a test/reference case, k**2 would be fine.  But, in #248 one of the things Barney found was a requirement on how the Fourier-space grid samples the power function, such that really steep power functions cannot be adequately represented on our typical grids.  I worried slightly that users might see our use of k**2 in a test as an example of one reasonable use case that the code does well, and since that isn't typically true, I decided to go to something shallower.