[Users] SphericalHarmonicRecon and SphericalHarmonicReconGen

Yosef Zlochower yosef at astro.rit.edu
Thu Jul 13 10:57:50 CDT 2017



On 07/13/2017 11:49 AM, Frank Loeffler wrote:
> On Thu, Jul 13, 2017 at 10:46:33AM -0500, Frank Loeffler wrote:
>> On Thu, Jul 13, 2017 at 11:38:28AM -0400, Yosef Zlochower wrote:
>>> stereo_q should be constant along the second coordinate and stereo_p
>>> along the first.
>>
>> Ok, so they do form a grid. Is there some other restriction? Because 
>> what the code does essentially is check that there is at least one 
>> point closer to 1+i than 1.e-10, but maybe that is just not the case? 
>> The code snippet in the mail just looks for the one with the shortest 
>> distance.
> 
> I take this back. I forgot about the following line within which all the 
> other code is:
> 
> if (minval(abs(zeta - dcmplx(1,1))) < 1.0d-10) then
> 
> This is indeed a little strange. Possibly some round-off or inproper 
> grid setup?
> 

I worry about the latter. Otherwise I would suggest to just remove the 
test. It has all the earmarks of some old test  to debug some 
long-forgotten issue that somehow got committed into the official code.


> Frank
> 


More information about the Users mailing list