5 October 2023 7:40:10.576 AM expokit_test(): FORTRAN77 version Test expokit(). SAMPLE_B 1SYMMETRIC MATRIX FROM NINE POINT START ON A 30 X 30 GRID. type :RSA size : 900 900 order : 900 number of nonzero : 4322 Harwell-Boeing matrix loaded ||A||_inf= 0.16E+02 ---------------------------------------------------- DSEXPV (Forward) has completed: ---------------------------------------------------- w(1:10) = 3456.5698306801160 7.3427169843682139 4094.7323184930651 1275.0417533588898 2939.0163458164616 2130.1604660263356 2471.7750081012391 2345.0013653129258 2387.0442087025081 2374.4520162188674 final report---------------------------------------- runtime = 0.42E-02 ||A||_inf = 0.16E+02 nz = 7744 n = 900 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 2 nreject = 0 nmult = 62 nexph = 2 nscale = 8 tol = 0.10E-09 t = 0.10E+01 tbrkdwn = 0.00E+00 step_min = 0.30E+00 step_max = 0.70E+00 max_round = 0.00E+00 sum_round = 0.00E+00 max_error = 0.36E-14 sum_error = 0.71E-14 hump = 0.21E+04 scale-norm= 0.21E+04 ---------------------------------------------------- DSEXPV (Backward) has completed: ---------------------------------------------------- w(1:10) = 1.0000000000001383 1.0000000000002518 1.0000000000003229 1.0000000000003466 1.0000000000003335 1.0000000000003058 1.0000000000002760 1.0000000000002542 1.0000000000002425 1.0000000000002411 final report---------------------------------------- runtime = 0.38E-02 ||A||_inf = 0.16E+02 nz = 7744 n = 900 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 2 nreject = 0 nmult = 62 nexph = 2 nscale = 8 tol = 0.10E-09 t = -.10E+01 tbrkdwn = 0.00E+00 step_min = 0.24E+00 step_max = 0.76E+00 max_round = 0.00E+00 sum_round = 0.00E+00 max_error = 0.36E-14 sum_error = 0.71E-14 hump = 0.10E+01 scale-norm= 0.48E-03 SAMPLE_D t = 2.0000000000000000 REAL SYMMETRIC CASE ******************************* A = -0.2265D+01 0.3401D+00 0.1605D+00 0.3114D+00 0.6003D+00 0.3401D+00 -0.2390D+01 0.8465D+00 0.3564D+00 0.8367D+00 0.1605D+00 0.8465D+00 -0.2290D+01 0.5922D+00 0.5924D+00 0.3114D+00 0.3564D+00 0.5922D+00 -0.1982D+01 0.3205D+00 0.6003D+00 0.8367D+00 0.5924D+00 0.3205D+00 -0.1648D+01 With DGPADM: exp(t*A) = 0.9994D-01 0.1385D+00 0.1304D+00 0.1082D+00 0.1775D+00 0.1385D+00 0.2125D+00 0.2014D+00 0.1626D+00 0.2648D+00 0.1304D+00 0.2014D+00 0.1953D+00 0.1582D+00 0.2502D+00 0.1082D+00 0.1626D+00 0.1582D+00 0.1388D+00 0.2020D+00 0.1775D+00 0.2648D+00 0.2502D+00 0.2020D+00 0.3404D+00 With DSPADM: exp(t*A) = 0.9994D-01 0.1385D+00 0.1304D+00 0.1082D+00 0.1775D+00 0.1385D+00 0.2125D+00 0.2014D+00 0.1626D+00 0.2648D+00 0.1304D+00 0.2014D+00 0.1953D+00 0.1582D+00 0.2502D+00 0.1082D+00 0.1626D+00 0.1582D+00 0.1388D+00 0.2020D+00 0.1775D+00 0.2648D+00 0.2502D+00 0.2020D+00 0.3404D+00 With DGCHBV: exp(t*A)e_1 = 9.9938862851251864E-002 0.13849879082076080 0.13038023839275395 0.10822981419271124 0.17750548647807074 With DSCHBV: exp(t*A)e_1 = 9.9938862851255514E-002 0.13849879082076086 0.13038023839275442 0.10822981419271194 0.17750548647807241 REAL UPPER HESSENBERG CASE ************************ H = -0.2265D+01 0.3401D+00 0.1605D+00 0.3114D+00 0.6003D+00 0.3401D+00 -0.2390D+01 0.8465D+00 0.3564D+00 0.8367D+00 0.0000D+00 0.8465D+00 -0.2290D+01 0.5922D+00 0.5924D+00 0.0000D+00 0.0000D+00 0.5922D+00 -0.1982D+01 0.3205D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.3205D+00 -0.1648D+01 With DGPADM: exp(t*H) = 0.1536D-01 0.2051D-01 0.2913D-01 0.4193D-01 0.6832D-01 0.1203D-01 0.3804D-01 0.5222D-01 0.6524D-01 0.1043D+00 0.8776D-02 0.3649D-01 0.5509D-01 0.6574D-01 0.9729D-01 0.3491D-02 0.1879D-01 0.3486D-01 0.5122D-01 0.5677D-01 0.6461D-03 0.4463D-02 0.1077D-01 0.2400D-01 0.5080D-01 With DNCHBV: exp(t*A)e_1 = 1.5358262200668641E-002 1.2025520908096474E-002 8.7761354281550730E-003 3.4914811395284125E-003 6.4606672162575128E-004 COMPLEX HERMITIAN CASE **************************** Re(H) = -0.1697D+01 0.1930D+00 0.8407D-01 0.2124D+00 0.2147D+00 0.1930D+00 -0.1643D+01 -0.2380D-01 -0.1988D+00 0.4879D+00 0.8407D-01 -0.2380D-01 -0.1876D+01 -0.4748D+00 0.4947D+00 0.2124D+00 -0.1988D+00 -0.4748D+00 -0.1295D+01 -0.3451D-01 0.2147D+00 0.4879D+00 0.4947D+00 -0.3451D-01 -0.1287D+01 Im(H) = 0.0000D+00 0.3918D+00 0.1505D+00 -0.2909D+00 0.3722D-01 -0.3918D+00 0.0000D+00 0.1402D+00 -0.3691D+00 -0.2446D+00 -0.1505D+00 -0.1402D+00 0.0000D+00 0.1096D+00 -0.3141D+00 0.2909D+00 0.3691D+00 -0.1096D+00 0.0000D+00 -0.3874D+00 -0.3722D-01 0.2446D+00 0.3141D+00 0.3874D+00 0.0000D+00 With ZGPADM: exp(t*H)e_1 = (9.06590329711321613E-002,3.88957790507077711E-018) (4.95706631370926992E-002,-7.57284642532656538E-002) (-3.96759983975782278E-003,-5.15579943320459630E-002) (3.09066819242459039E-002,6.97046558676174682E-002) (5.50148381156331523E-002,-3.60400966614133275E-002) With ZHPADM: exp(t*H)_e_1 = (9.06590329711320780E-002,-1.21480386307847623E-018) (4.95706631370926576E-002,-7.57284642532655983E-002) (-3.96759983975782625E-003,-5.15579943320459561E-002) (3.09066819242458830E-002,6.97046558676174405E-002) (5.50148381156331384E-002,-3.60400966614133206E-002) With ZGCHBV: exp(t*H)e_1 = (9.06590329708769765E-002,-8.45995331893913821E-016) (4.95706631370486164E-002,-7.57284642531964314E-002) (-3.96759983975812029E-003,-5.15579943319686290E-002) (3.09066819242636744E-002,6.97046558675448458E-002) (5.50148381155987423E-002,-3.60400966613236770E-002) SAMPLE_G 1UNSYMMETRIC MATRIX ORANI678 ST ERTAR 24OCT84 type :RUA size : 2529 2529 order : 2529 number of nonzero : 90158 There is a second hand Harwell-Boeing matrix loaded ||A||_inf= 0.95E+01 ---------------------------------------------------- DGEXPV (COO) has completed: ---------------------------------------------------- w(1:10) = 6464.4009480839586 4828.8899856184116 5450.5845968253079 5592.7251537899101 4575.6820139573229 4199.7768387207225 5087.4782505485646 3746.9681136961131 4678.1202261774997 3797.7801341101217 final report---------------------------------------- runtime = 0.50E-01 ||A||_inf = 0.95E+01 nz = 90158 n = 2529 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 3 nreject = 0 nmult = 93 nexph = 3 nscale = 18 tol = 0.15E-07 t = 0.10E+02 tbrkdwn = 0.00E+00 step_min = 0.61E+00 step_max = 0.49E+01 max_round = 0.00E+00 sum_round = 0.00E+00 max_error = 0.31E-08 sum_error = 0.37E-08 hump = 0.25E+05 scale-norm= 0.25E+05 ---------------------------------------------------- DGEXPV (CCS) has completed: ---------------------------------------------------- w(1:10) = 6464.4009480839586 4828.8899856184116 5450.5845968253079 5592.7251537899101 4575.6820139573229 4199.7768387207225 5087.4782505485646 3746.9681136961131 4678.1202261774997 3797.7801341101217 final report---------------------------------------- runtime = 0.47E-01 ||A||_inf = 0.95E+01 nz = 90158 n = 2529 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 3 nreject = 0 nmult = 93 nexph = 3 nscale = 18 tol = 0.15E-07 t = 0.10E+02 tbrkdwn = 0.00E+00 step_min = 0.61E+00 step_max = 0.49E+01 max_round = 0.00E+00 sum_round = 0.00E+00 max_error = 0.31E-08 sum_error = 0.37E-08 hump = 0.25E+05 scale-norm= 0.25E+05 ---------------------------------------------------- DGEXPV (CRS) has completed: ---------------------------------------------------- w(1:10) = 6464.4009480839586 4828.8899856184116 5450.5845968253079 5592.7251537899101 4575.6820139573229 4199.7768387207225 5087.4782505485646 3746.9681136961131 4678.1202261774997 3797.7801341101217 final report---------------------------------------- runtime = 0.49E-01 ||A||_inf = 0.95E+01 nz = 90158 n = 2529 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 3 nreject = 0 nmult = 93 nexph = 3 nscale = 18 tol = 0.15E-07 t = 0.10E+02 tbrkdwn = 0.00E+00 step_min = 0.61E+00 step_max = 0.49E+01 max_round = 0.00E+00 sum_round = 0.00E+00 max_error = 0.31E-08 sum_error = 0.37E-08 hump = 0.25E+05 scale-norm= 0.25E+05 SAMPLE_M ||A||_inf= 0.32E+02 ---------------------------------------------------- DMEXPV has completed: ---------------------------------------------------- w(1:10) = 0.13051504168006217 1.2181403890139142E-002 1.3517629316863591E-002 1.2616454029072686E-003 1.7712755656579562E-002 1.6531905279474257E-003 1.8345354072885978E-003 1.7122330468026918E-004 1.1417997294345146E-002 1.0656797474722140E-003 final report---------------------------------------- runtime = 0.27E-01 ||A||_inf = 0.32E+02 nz = 11264 n = 1024 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 6 nreject = 0 nmult = 186 nexph = 6 nscale = 44 tol = 0.10E-09 t = 0.10E+02 tbrkdwn = 0.00E+00 step_min = 0.17E+00 step_max = 0.32E+01 max_round = 0.32E-16 sum_round = 0.81E-16 max_error = 0.72E-14 sum_error = 0.43E-13 hump = 0.10E+01 scale-norm= 0.24E+00 ---------------------------------------------------- DGEXPV has completed: ---------------------------------------------------- w(1:10) = 0.13051504168006217 1.2181403890139142E-002 1.3517629316863591E-002 1.2616454029072686E-003 1.7712755656579562E-002 1.6531905279474257E-003 1.8345354072885978E-003 1.7122330468026918E-004 1.1417997294345146E-002 1.0656797474722140E-003 final report---------------------------------------- runtime = 0.27E-01 ||A||_inf = 0.32E+02 nz = 11264 n = 1024 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 6 nreject = 0 nmult = 186 nexph = 6 nscale = 44 tol = 0.10E-09 t = 0.10E+02 tbrkdwn = 0.00E+00 step_min = 0.17E+00 step_max = 0.32E+01 max_round = 0.00E+00 sum_round = 0.00E+00 max_error = 0.72E-14 sum_error = 0.43E-13 hump = 0.10E+01 scale-norm= 0.24E+00 SAMPLE_P 1UNSYMMETRIC MATRIX ORANI678 ST ERTAR 24OCT84 type :RUA size : 2529 2529 order : 2529 number of nonzero : 90158 There is a second hand Harwell-Boeing matrix loaded ||A||_inf= 0.95E+01 ---------------------------------------------------- DGPHIV (CCS) has completed: ---------------------------------------------------- w(1:10) = 6464.4009480834493 4828.8899856180487 5450.5845968248868 5592.7251537894754 4575.6820139569709 4199.7768387204142 5087.4782505481780 3746.9681136958407 4678.1202261771368 3797.7801341098088 final report---------------------------------------- runtime = 0.49E-01 ||A||_inf = 0.95E+01 nz = 90158 n = 2529 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 3 nreject = 0 nmult = 96 nexph = 3 nscale = 17 tol = 0.15E-07 t = 0.10E+02 tbrkdwn = 0.00E+00 step_min = 0.62E+00 step_max = 0.52E+01 max_round = 0.00E+00 sum_round = 0.00E+00 max_error = 0.49E-08 sum_error = 0.49E-08 ---------------------------------------------------- DGPHIV (CCS) has completed: ---------------------------------------------------- w(1:10) = 4962.2511930988721 3658.2951908942173 4152.9628069217488 4279.0497069392668 3442.7780495195943 3155.7521413695863 3861.9898164862652 2781.6080219244691 3525.0002206756062 2933.5009201612670 final report---------------------------------------- runtime = 0.48E-01 ||A||_inf = 0.95E+01 nz = 90158 n = 2529 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 3 nreject = 0 nmult = 96 nexph = 3 nscale = 17 tol = 0.15E-07 t = 0.10E+02 tbrkdwn = 0.00E+00 step_min = 0.61E+00 step_max = 0.52E+01 max_round = 0.00E+00 sum_round = 0.00E+00 max_error = 0.32E-08 sum_error = 0.32E-08 relative difference (phi vs. exp) = 5.1206129103962300E-014 ---------------------------------------------------- DGPHIV (CCS) has completed: ---------------------------------------------------- w(1:10) = 11426.652141182287 8487.1851765123174 9603.5474037466338 9871.7748607287613 8018.4600634765529 7355.5289800900528 8949.4680670344824 6528.5761356203893 8203.1204468527321 6731.2810542710522 final report---------------------------------------- runtime = 0.49E-01 ||A||_inf = 0.95E+01 nz = 90158 n = 2529 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 3 nreject = 0 nmult = 96 nexph = 3 nscale = 17 tol = 0.15E-07 t = 0.10E+02 tbrkdwn = 0.00E+00 step_min = 0.61E+00 step_max = 0.52E+01 max_round = 0.00E+00 sum_round = 0.00E+00 max_error = 0.56E-08 sum_error = 0.56E-08 SAMPLE_Z 1SYMMETRIC STRUCTURE OF REPRESENTATION OF ENTIRE U.S. POWER NETWORK type :PSA size : 5300 5300 order : 5300 number of nonzero : 13571 Harwell-Boeing pattern loaded ||A||_inf= 0.61E+02 ---------------------------------------------------- ZGEXPV has completed: ---------------------------------------------------- w(1:10) = (305250.71372538898,-113.89588308939845) (6.74919321355144169E-008,6.86578771621437937E-008) (-3.0966853417659617,4.6959293837327314) (-1.67768882951225342E-006,-7.88929583291198018E-006) (9.16756637825495513E-006,-2.68608052331716616E-006) (6.89109212036109146E-003,-8.28641975910175271E-002) (2.63648739146995889E-005,8.23662197501563568E-004) (-0.20614403044848950,-0.20907746417618769) (5.9467627640614289,-0.95593305646046267) (0.14657965977112086,0.23806649030555826) final report---------------------------------------- runtime = 0.96E-01 ||A||_inf = 0.61E+02 nz = 21842 n = 5300 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 3 nreject = 0 nmult = 93 nexph = 3 nscale = 11 tol = 0.10E-04 t = 0.10E+01 tbrkdwn = 0.00E+00 step_min = 0.80E-01 step_max = 0.80E+00 max_round = 0.00E+00 sum_round = 0.00E+00 max_error = 0.40E-06 sum_error = 0.40E-06 hump = 0.15E+07 scale-norm= 0.15E+07 ---------------------------------------------------- ZHEXPV has completed: ---------------------------------------------------- w(1:10) = (305250.71372538852,-113.89588308936860) (6.74919321355037364E-008,6.86578771621655519E-008) (-3.0966853417659750,4.6959293837327429) (-1.67768882951258799E-006,-7.88929583291203439E-006) (9.16756637825358633E-006,-2.68608052331702725E-006) (6.89109212036058925E-003,-8.28641975910162087E-002) (2.63648739146868088E-005,8.23662197501586445E-004) (-0.20614403044848995,-0.20907746417618778) (5.9467627640614831,-0.95593305646044713) (0.14657965977112170,0.23806649030555979) final report---------------------------------------- runtime = 0.33E-01 ||A||_inf = 0.61E+02 nz = 21842 n = 5300 m = 30 itrace = 0 iflag = 0 ibrkflag = 0 mbrkdwn = 30 nstep = 3 nreject = 0 nmult = 93 nexph = 3 nscale = 11 tol = 0.10E-04 t = 0.10E+01 tbrkdwn = 0.00E+00 step_min = 0.80E-01 step_max = 0.80E+00 max_round = 0.00E+00 sum_round = 0.00E+00 max_error = 0.40E-06 sum_error = 0.40E-06 hump = 0.15E+07 scale-norm= 0.15E+07 expokit_test(): Normal end of execution. 5 October 2023 7:40:11.172 AM