rank of F = 10 of 15;  |F x_pinv - g| = 2.8e-16
 CG iterations   |x_cg - x_pinv| / |x_pinv|
             1                     4.22e-01
             2                     2.31e-01
             3                     9.63e-02
             5                     3.50e-03
             8                     2.86e-06
            12                     2.32e-14
            15                     2.32e-14

   delta  KL(old||new)        g.d    true dJ  surrogate
     0.1       0.10452    0.11150    0.11714    0.11868
    0.03       0.03079    0.06107    0.06302    0.06343
    0.01       0.01016    0.03526    0.03595    0.03608
   0.003       0.00303    0.01931    0.01953    0.01956
   0.001       0.00101    0.01115    0.01122    0.01124
