Skip to content

Instantly share code, notes, and snippets.

@dpo
Created March 20, 2023 18:05
Show Gist options
  • Save dpo/84bacdcffe675ee122a83fd06e03136a to your computer and use it in GitHub Desktop.
Save dpo/84bacdcffe675ee122a83fd06e03136a to your computer and use it in GitHub Desktop.
~/dev/ralna/GALAHAD multi-precision *2 !3 ?10 ❯ ./builddir_gnu/glrtt_single 34s 14:02:27
==== normal exits ======
stopping tolerance = 3.4527E-03, sigma = 1.0000E+03
Iter objective pgnorm lambda gamma it info
0 0.00000000E+00 1.00E+01 - - - -
iteration 0 leftmost eigenvalue = -1.99999995529652E-02
tridiagonal KKT residual = 1.0000E+01
tridiagonal KKT residual = 0.0000E+00
tridiagonal KKT residual = 0.0000E+00
lambda error x
2.0006416365504265E-02 -2.0006416889373213E-05 1.0842021724855044E-19
hard case
iteration 6 hard case: leftmost eigenvector found
Iter objective pgnorm lambda gamma it info
1 -7.79203050E+06 NaN 1.99999995529652E-02 1.40E-01 1 0
it neg l e u piv
1 1 -2.14040828E+00 -2.00009998E-02 -1.99999996E-02 -1.95972695E+04
2 1 -2.14040828E+00 -2.00938344E+00 -2.00009998E-02 -8.77141953E-04
iteration 3 leftmost eigenvalue = -2.01025605201721E+00
tridiagonal KKT residual = 1.0000E+01
tridiagonal KKT residual = 0.0000E+00
tridiagonal KKT residual = 0.0000E+00
lambda error x
2.0102627277374268E+00 -2.0102628041058779E-03 0.0000000000000000E+00
hard case
iteration 6 hard case: leftmost eigenvector found
Iter objective pgnorm lambda gamma it info
2 -3.74613047E+04 NaN 2.01025605201721E+00 1.00E+00 1 0
it neg l e u piv
1 1 -3.14019990E+00 -2.01025701E+00 -2.01025605E+00 -1.04813919E+06
2 1 -3.14019990E+00 -2.01044989E+00 -2.01025701E+00 -5.13163867E+03
3 1 -3.14019990E+00 -2.93637371E+00 -2.01044989E+00 -1.39745712E-01
4 0 -3.14019990E+00 -3.00349927E+00 -2.93637371E+00 1.66893005E-06
leftmost eigenvalue = -3.00349926948547E+00
iteration 4 leftmost eigenvalue = -3.00349926948547E+00
tridiagonal KKT residual = 1.0000E+01
tridiagonal KKT residual = 3.9062E-03
tridiagonal KKT residual = 1.9531E-03
lambda error x
3.0035061836242676E+00 -3.0035062227398157E-03 -1.0842021724855044E-19
hard case
iteration 6 hard case: leftmost eigenvector found
Iter objective pgnorm lambda gamma it info
3 -7.11314111E+03 NaN 3.00349926948547E+00 1.00E+00 1 0
it neg l e u piv
1 1 -3.99999976E+00 -3.00350022E+00 -3.00349927E+00 -2.79497906E+05
2 0 -3.99999976E+00 -3.50174999E+00 -3.00350022E+00 3.01384449E-01
3 0 -3.50174999E+00 -3.41608047E+00 -3.00350022E+00 1.18637085E-03
4 1 -3.41608047E+00 -3.41578460E+00 -3.00350022E+00 -1.19209290E-06
leftmost eigenvalue = -3.41578459739685E+00
iteration 4 leftmost eigenvalue = -3.41578459739685E+00
tridiagonal KKT residual = 1.0000E+01
tridiagonal KKT residual = 1.9531E-03
tridiagonal KKT residual = 1.9531E-03
lambda error x
3.4157915115356445E+00 -3.4157915506511927E-03 -0.0000000000000000E+00
hard case
iteration 6 hard case: leftmost eigenvector found
Iter objective pgnorm lambda gamma it info
4 -3.20910547E+03 Infinity 3.41578459739685E+00 1.00E+00 1 0
it neg l e u piv
1 1 -4.00000000E+00 -3.41578555E+00 -3.41578460E+00 -3.81126344E+05
2 1 -4.00000000E+00 -3.41578770E+00 -3.41578555E+00 -8.73403594E+04
3 1 -4.00000000E+00 -3.41582322E+00 -3.41578770E+00 -6.48837598E+03
4 1 -4.00000000E+00 -3.42160749E+00 -3.41582322E+00 -4.17879448E+01
5 0 -4.00000000E+00 -3.62540913E+00 -3.42160749E+00 4.62012291E-02
6 0 -3.62540913E+00 -3.61887765E+00 -3.42160749E+00 4.07695770E-05
7 0 -3.61887765E+00 -3.61887240E+00 -3.42160749E+00 3.09944153E-06
leftmost eigenvalue = -3.61887240409851E+00
iteration 7 leftmost eigenvalue = -3.61887240409851E+00
tridiagonal KKT residual = 1.0000E+01
tridiagonal KKT residual = 1.4648E-03
tridiagonal KKT residual = 9.7656E-04
lambda error x
3.6188793182373047E+00 -3.6188792437314987E-03 1.0842021724855044E-19
hard case
iteration 6 hard case: leftmost eigenvector found
Iter objective pgnorm lambda gamma it info
5 -1.45280505E+03 NaN 3.61887240409851E+00 1.00E+00 1 0
it neg l e u piv
1 1 -4.00000000E+00 -3.61887336E+00 -3.61887240E+00 -1.02250148E+05
2 0 -4.00000000E+00 -3.80943680E+00 -3.61887336E+00 5.92562199E-01
3 0 -3.80943680E+00 -3.73380184E+00 -3.61887336E+00 1.48906708E-02
4 0 -3.73380184E+00 -3.73255301E+00 -3.61887336E+00 1.45435333E-05
leftmost eigenvalue = -3.73255300521851E+00
iteration 4 leftmost eigenvalue = -3.73255300521851E+00
tridiagonal KKT residual = 1.0000E+01
tridiagonal KKT residual = 4.8828E-04
tridiagonal KKT residual = 9.7656E-04
lambda error x
3.7325599193572998E+00 1.9962674379348755E+00 2.0000000000000000E+00
Iter objective pgnorm lambda gamma it info
6 -7.55384094E+02 3.90E+03 3.73255991935730E+00 1.00E+00 2 1d
it neg l e u piv
1 1 -4.00000000E+00 -3.73255396E+00 -3.73255301E+00 -3.91780117E+04
2 0 -4.00000000E+00 -3.86627698E+00 -3.73255396E+00 6.66471481E-01
3 0 -3.86627698E+00 -3.80357003E+00 -3.73255396E+00 2.39460468E-02
4 1 -3.80357003E+00 -3.80226231E+00 -3.73255396E+00 -4.29153442E-06
leftmost eigenvalue = -3.80226230621338E+00
iteration 4 leftmost eigenvalue = -3.80226230621338E+00
tridiagonal KKT residual = 1.0000E+01
tridiagonal KKT residual = 7.3242E-04
tridiagonal KKT residual = 9.7656E-04
lambda error x
3.8022692203521729E+00 -3.8022692315280437E-03 0.0000000000000000E+00
hard case
iteration 6 hard case: leftmost eigenvector found
Iter objective pgnorm lambda gamma it info
7 -5.66669556E+02 NaN 3.80226230621338E+00 1.00E+00 1 0
it neg l e u piv
1 1 -4.00000000E+00 -3.80226326E+00 -3.80226231E+00 -6.98671562E+04
2 1 -4.00000000E+00 -3.80226707E+00 -3.80226326E+00 -1.17738662E+04
3 1 -4.00000000E+00 -3.80236912E+00 -3.80226707E+00 -5.02606506E+02
4 1 -4.00000000E+00 -3.82354236E+00 -3.80236912E+00 -1.38976836E+00
5 0 -4.00000000E+00 -3.84832478E+00 -3.82354236E+00 9.29498672E-03
6 1 -3.84832478E+00 -3.84798265E+00 -3.82354236E+00 -1.66893005E-06
leftmost eigenvalue = -3.84798264503479E+00
iteration 6 leftmost eigenvalue = -3.84798264503479E+00
tridiagonal KKT residual = 1.0000E+01
tridiagonal KKT residual = 2.4414E-04
tridiagonal KKT residual = 4.8828E-04
lambda error x
3.8479895591735840E+00 -2.0038480758666992E+00 -2.0000000000000000E+00
hard case
iteration 6 hard case: leftmost eigenvector found
Iter objective pgnorm lambda gamma it info
8 -3.74046814E+02 NaN 3.84798264503479E+00 1.00E+00 1 0
it neg l e u piv
1 1 -4.00000000E+00 -3.84798360E+00 -3.84798265E+00 -3.95456562E+04
2 1 -4.00000000E+00 -3.84801292E+00 -3.84798360E+00 -1.20902356E+03
3 1 -4.00000000E+00 -3.85548115E+00 -3.84801292E+00 -3.75981665E+00
4 0 -4.00000000E+00 -3.87998414E+00 -3.85548115E+00 1.66389942E-02
5 0 -3.87998414E+00 -3.87954760E+00 -3.85548115E+00 4.05311584E-05
6 1 -3.87954760E+00 -3.87954569E+00 -3.85548115E+00 -3.19480896E-05
leftmost eigenvalue = -3.87954568862915E+00
iteration 6 leftmost eigenvalue = -3.87954568862915E+00
tridiagonal KKT residual = 1.0000E+01
tridiagonal KKT residual = 4.8828E-04
tridiagonal KKT residual = 7.3242E-04
lambda error x
3.8795526027679443E+00 -3.8795527070760727E-03 -0.0000000000000000E+00
hard case
iteration 6 hard case: leftmost eigenvector found
Iter objective pgnorm lambda gamma it info
9 -2.96618225E+02 NaN 3.87954568862915E+00 1.00E+00 1 0
it neg l e u piv
1 1 -4.00000000E+00 -3.87954664E+00 -3.87954569E+00 -3.81082281E+05
2 1 -4.00000000E+00 -3.87954688E+00 -3.87954664E+00 -7.35406719E+04
3 1 -4.00000000E+00 -3.87954736E+00 -3.87954688E+00 -3.12811484E+04
4 1 -4.00000000E+00 -3.87954926E+00 -3.87954736E+00 -9.54692676E+03
5 1 -4.00000000E+00 -3.87956023E+00 -3.87954926E+00 -1.90614380E+03
6 1 -4.00000000E+00 -3.87976837E+00 -3.87956023E+00 -1.15812027E+02
7 1 -4.00000000E+00 -3.89654565E+00 -3.87976837E+00 -3.94239187E-01
8 0 -4.00000000E+00 -3.90226245E+00 -3.89654565E+00 1.50275230E-03
9 1 -3.90226245E+00 -3.90223336E+00 -3.89654565E+00 -1.93119049E-05
leftmost eigenvalue = -3.90223336219788E+00
iteration 9 leftmost eigenvalue = -3.90223336219788E+00
tridiagonal KKT residual = 1.0000E+01
tridiagonal KKT residual = 2.4414E-04
tridiagonal KKT residual = 2.4414E-04
lambda error x
3.9022402763366699E+00 3.6893488147419103E+19 3.6893488147419103E+19
tridiagonal KKT residual = 1.0000E+01
tridiagonal KKT residual = 1.3752E-23
tridiagonal KKT residual = 1.3752E-23
lambda error x
3.6893488147419103E+22 -3.6893488147419103E+19 2.0000000000000000E+00
real, recurred f = 2.66666674804688E+03 -2.45956572637462E+22
Iter objective pgnorm lambda gamma it info
10 -2.45956573E+22 0.00E+00 3.68934881474191E+22 1.00E+00 2 0d
pass 1 GLRT_solve exit status = 0
stopping tolerance = 2.4414E-03, sigma = 1.0000E+03
Iter objective pgnorm lambda gamma it info
0 0.00000000E+00 7.07E+00 - - - -
iteration 0 leftmost eigenvalue = 1.99000000953674E+00
tridiagonal KKT residual = 7.0711E+00
tridiagonal KKT residual = 0.0000E+00
tridiagonal KKT residual = 0.0000E+00
lambda error x
0.0000000000000000E+00 -2.0000000000000000E+00 -2.0000000000000000E+00
hard case
^C
@dpo
Copy link
Author

dpo commented Mar 20, 2023

NB: there are some NaNs in the output above, along with other suspicious values (1E+22, 1E+39, etc.)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment