[Users] SphericalHarmonicRecon and SphericalHarmonicReconGen

[Yosef Zlochower]

Hi Ian,
   Bela found a bug in the null code and I committed his fix on June 2. 
That commit included an updated testsuite for SphericalHarmonicRecon, 
but not for SphericalHarmonicReconGen. I wonder if the 
SphericalHarmonicReconGen testsuite is out of date.

The failure for SphericalHarmonicRecon is quite weird.
The failing code is

  if (minval(abs(zeta - dcmplx(1,1))) < 1.0d-10) then
         Tarr = minloc(abs(stereo_q(:,1)-1.))
         loc_q = Tarr(1)
         Tarr = minloc(abs(stereo_p(1,:)-1.))
         loc_p = Tarr(1)
         if (abs(zeta(loc_q,loc_p) - dcmplx(1.,1.)) .gt. 1d-10) then
            call CCTK_WARN(0, " Error ")
         endif

      endif
I can't figure out why the test is there. However, I also can't figure 
out how it can fail.

Basically the code uses either two real angular coordinates (q and p) 
or one complex one (zeta). With zeta = q + i * p. The code is checking 
that if zeta == 1+i anywhere, that it is equal to 1 + i at the point 
where q=1 and p=1.

Can you try to output
loc_q, loc_p, stereo_p(1,loc_p), stereo_p(loc_q,1), and zeta(loc_q, loq_p)?




On 07/13/2017 04:31 AM, Ian Hinder wrote:
> Hi,
> 
> We are down to three failing tests in Jenkins:
> 
> SphericalHarmonicRecon.regression_test/2procs
> SphericalHarmonicReconGen.SpEC-dat-test/2procs
> SphericalHarmonicReconGen.SpEC-h5-test/2procs
> 
> These tests all pass on one process but fail on two processes in 
> Jenkins, which uses the ubuntu.cfg optionlist.  They all seem to pass on 
> multiple processes on all other machines 
> (http://einsteintoolkit.org/testsuite_results/index.php), including my 
> laptop with gcc.
> 
> The first, SphericalHarmonicRecon.regression_test, fails like this:
> 
> WARNING level 0 from host 7ce14e5707a0 process 0
>    while executing schedule bin NullEvol_Initial, routine 
> NullEvolve::NullEvol_InitialSlice
>    in thorn NullEvolve, file NullEvol_InitialSlice.F90:42:
>    ->  Error
> 
> The second, SphericalHarmonicReconGen.SpEC-dat-test, fails like this:
> 
> NewsB_scri.L02Mm01.asc: substantial differences
>        significant differences on 1 (out of 2) lines
>        maximum absolute difference in column 1 is 963
>        maximum absolute difference in column 2 is 0.000185770963653907
>        maximum absolute difference in column 3 is 0.000142466608463344
>        maximum relative difference in column 1 is 1
>        maximum relative difference in column 2 is 1
>        maximum relative difference in column 3 is 1
>        ...
> 
> The third, SphericalHarmonicReconGen.SpEC-h5-test, fails like this:
> 
> NewsB_scri.L02Mm01.asc: substantial differences
>        significant differences on 1 (out of 2) lines
>        maximum absolute difference in column 1 is 963
>        maximum absolute difference in column 2 is 0.000185770963653907
>        maximum absolute difference in column 3 is 0.000142466608463344
>        maximum relative difference in column 1 is 1
>        maximum relative difference in column 2 is 1
>        maximum relative difference in column 3 is 1
>        ...
> 
> I suspect the second and third failures have the same cause.  Do we have 
> any idea why these tests fail?  They don't seem to fail on any other 
> machines (http://einsteintoolkit.org/testsuite_results/index.php).
> 
> -- 
> Ian Hinder
> http://members.aei.mpg.de/ianhin
> 
> 
> 
> _______________________________________________
> Users mailing list
> Users at einsteintoolkit.org
> http://lists.einsteintoolkit.org/mailman/listinfo/users
>