On Mar 17, 2013, at 6:01 PM, Mike Jarvis <notifications@github.com> wrote:

> In galsim/lensing.py:
> 
> >  
> > +    def recompute_power(self):
> What is the point of this function? It doesn't seem to ever be called. When we need a new realization of the power (in buildGriddedShears), we just make a new PowerSpectrumRealizer instance.
> 
> Of course, it doesn't hurt to have it, so it's fine to leave this as is. I'm just curious as to why it was put in. (Or why init was broken up into set_size and set_power for that matter, which seems to have been specifically to enable this function.)
> 
> 
Hmm.  Perhaps Barney will set me straight, but I *think* what's going on here is that originally, Joe only modified PowerSpectrumRealizer in order to put in the convergence calculation.  Then he wrote a test script (which later got revamped and put into unit tests) that called PowerSpectrumRealizer directly.  That meant that in the course of the testing, the user who was dealing with PowerSpectrumRealizer might have a reason to use recompute_power.