[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