astro-informatics / purify

Next-generation radio interferometric imaging.
https://astro-informatics.github.io/purify
GNU General Public License v2.0
17 stars 13 forks source link

[3.0.1] test_algo_factory fails #283

Open olebole opened 4 years ago

olebole commented 4 years ago

In the same run as #282, there are more failures. The next one is test_algo_factory:

      Start 16: test_algo_factory
16/27 Test #16: test_algo_factory ................***Failed  1125.60 sec
[2020-05-17 11:43:24.354] [sopt] [critical] Performing Proximal ADMM with L1 and L2 operators
[2020-05-17 11:43:24.354] [sopt] [critical] Performing Proximal ADMM

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
test_algo_factory is a Catch v1.12.1 host application.
Run with -? for options

Randomness seeded to: 1589715803

-------------------------------------------------------------------------------
padmm_factory
-------------------------------------------------------------------------------
/build/purify-3.0.1/cpp/tests/algo_factory.cc:22
...............................................................................

/build/purify-3.0.1/cpp/tests/algo_factory.cc:59: FAILED:
  CHECK( diagnostic.niters == 139 )
with expansion:
  1000 (0x3e8) == 139
with message:
  uv_data.vis.head(5) := (-0.0135681,-0.00278705)
   (-0.0107864,0.00278611)
   (0.0195668,-0.00841525)
    (0.00142236,0.0200594)
     (0.0165797,0.0155949)

/build/purify-3.0.1/cpp/tests/algo_factory.cc:64: FAILED:
  CHECK( image.isApprox(solution, 1e-6) )
with expansion:
  false
with messages:
  uv_data.vis.head(5) := (-0.0135681,-0.00278705)
   (-0.0107864,0.00278611)
   (0.0195668,-0.00841525)
    (0.00142236,0.0200594)
     (0.0165797,0.0155949)
  Vector<t_complex>::Map(solution.data(), solution.size()).real().head(10) :=
  0.000134124
  0.000120018
  9.25444e-05
  6.46581e-05
  4.56205e-05
  3.26691e-05
  9.34926e-06
            0
            0
  2.13941e-05
  Vector<t_complex>::Map(image.data(), image.size()).real().head(10) :=       0
        0
  146.661
   340.42
  37.2947
        0
        0
        0
        0
  71.1281
  Vector<t_complex>::Map((image / solution).eval().data(), image.size()).real()
  .head(10) :=           0
            0
  1.58477e+06
  5.26492e+06
       817498
            0
            0
         -nan
         -nan
  3.32467e+06

/build/purify-3.0.1/cpp/tests/algo_factory.cc:71: FAILED:
  CHECK( residual_image.real().isApprox(residual.real(), 1e-6) )
with expansion:
  false
with messages:
  uv_data.vis.head(5) := (-0.0135681,-0.00278705)
   (-0.0107864,0.00278611)
   (0.0195668,-0.00841525)
    (0.00142236,0.0200594)
     (0.0165797,0.0155949)
  Vector<t_complex>::Map(solution.data(), solution.size()).real().head(10) :=
  0.000134124
  0.000120018
  9.25444e-05
  6.46581e-05
  4.56205e-05
  3.26691e-05
  9.34926e-06
            0
            0
  2.13941e-05
  Vector<t_complex>::Map(image.data(), image.size()).real().head(10) :=       0
        0
  146.661
   340.42
  37.2947
        0
        0
        0
        0
  71.1281
  Vector<t_complex>::Map((image / solution).eval().data(), image.size()).real()
  .head(10) :=           0
            0
  1.58477e+06
  5.26492e+06
       817498
            0
            0
         -nan
         -nan
  3.32467e+06
  Vector<t_complex>::Map(residual.data(), residual.size()).real().head(10) :=
  0.00188788
  0.000944006
   0.00412839
   0.00478902
   0.00352017
   0.00123934
   0.00207853
  -0.00163091
  -0.00195286
   0.00415593
  Vector<t_complex>::Map(residuals.data(), residuals.size()).real().head(10) :=
  -0.00220043
   -0.0028312
  0.000602711
   0.00140065
  0.000151452
  -0.00221713
  -0.00154228
  -0.00540908
   -0.0058048
  0.000293272

[2020-05-17 11:48:50.836] [sopt] [critical] Performing Forward Backward with L1 and L2 norms
[2020-05-17 11:48:50.836] [sopt] [critical] Performing Forward Backward Splitting
-------------------------------------------------------------------------------
fb_factory
-------------------------------------------------------------------------------
/build/purify-3.0.1/cpp/tests/algo_factory.cc:74
...............................................................................

/build/purify-3.0.1/cpp/tests/algo_factory.cc:113: FAILED:
  CHECK( diagnostic.niters == 21 )
with expansion:
  974 (0x3ce) == 21
with message:
  uv_data.vis.head(5) := (-0.0135681,-0.00278705)
   (-0.0107864,0.00278611)
   (0.0195668,-0.00841525)
    (0.00142236,0.0200594)
     (0.0165797,0.0155949)

/build/purify-3.0.1/cpp/tests/algo_factory.cc:119: FAILED:
  CHECK( image.isApprox(solution, 1e-6) )
with expansion:
  false
with messages:
  uv_data.vis.head(5) := (-0.0135681,-0.00278705)
   (-0.0107864,0.00278611)
   (0.0195668,-0.00841525)
    (0.00142236,0.0200594)
     (0.0165797,0.0155949)
  Vector<t_complex>::Map(solution.data(), solution.size()).real().head(10) :=
  1.30079e-05
  7.72943e-06
  3.34334e-06
            0
            0
  2.88329e-06
  3.73058e-06
  8.40954e-06
  2.97746e-05
  5.15951e-05
  Vector<t_complex>::Map(image.data(), image.size()).real().head(10) :=  1.
  41211
   1.29521
   1.16625
   1.06191
  0.958439
   0.89535
  0.854436
  0.840068
  0.809206
  0.795173
  Vector<t_complex>::Map((image / solution).eval().data(), image.size()).real()
  .head(10) :=  108558
   167569
   348827
     -nan
     -nan
   310531
   229036
  99894.6
  27177.7
  15411.8

/build/purify-3.0.1/cpp/tests/algo_factory.cc:127: FAILED:
  CHECK( residual_image.real().isApprox(residual.real(), 1e-6) )
with expansion:
  false
with messages:
  uv_data.vis.head(5) := (-0.0135681,-0.00278705)
   (-0.0107864,0.00278611)
   (0.0195668,-0.00841525)
    (0.00142236,0.0200594)
     (0.0165797,0.0155949)
  Vector<t_complex>::Map(solution.data(), solution.size()).real().head(10) :=
  1.30079e-05
  7.72943e-06
  3.34334e-06
            0
            0
  2.88329e-06
  3.73058e-06
  8.40954e-06
  2.97746e-05
  5.15951e-05
  Vector<t_complex>::Map(image.data(), image.size()).real().head(10) :=  1.
  41211
   1.29521
   1.16625
   1.06191
  0.958439
   0.89535
  0.854436
  0.840068
  0.809206
  0.795173
  Vector<t_complex>::Map((image / solution).eval().data(), image.size()).real()
  .head(10) :=  108558
   167569
   348827
     -nan
     -nan
   310531
   229036
  99894.6
  27177.7
  15411.8
  Vector<t_complex>::Map(residual.data(), residual.size()).real().head(10) :=
  0.00242299
    0.00134493
    0.00437652
    0.00495848
    0.00365283
    0.00143164
    0.00250585
  -0.000837817
  -0.000778659
    0.00551395
  Vector<t_complex>::Map(residuals.data(), residuals.size()).real().head(10) :=
  -0.00220068
  -0.00283152
  0.000602778
   0.00140081
  0.000151467
  -0.00221739
  -0.00154246
  -0.00540972
  -0.00580548
  0.000293302

[2020-05-17 11:53:12.644] [sopt] [critical] Performing Joint MAP estimation
[2020-05-17 11:53:12.650] [sopt] [critical] Performing Forward Backward with L1 and L2 norms
[2020-05-17 11:53:12.650] [sopt] [critical] Performing Forward Backward Splitting
[2020-05-17 11:57:42.218] [sopt] [critical] Performing Forward Backward with L1 and L2 norms
[2020-05-17 11:57:42.218] [sopt] [critical] Performing Forward Backward Splitting
[…]
[2020-05-17 12:02:08.318] [sopt] [critical] Performing Forward Backward with L1 and L2 norms
[2020-05-17 12:02:08.319] [sopt] [critical] Performing Forward Backward Splitting
===============================================================================
test cases: 3 | 1 passed | 2 failed
assertions: 9 | 3 passed | 6 failed