5 December 2010 9:00:08.404 PM EXPOKIT: FORTRAN77 version Tests for the EXPOKIT library. 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.80E-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.74E-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.99388628512518640E-002 0.13849879082076080 0.13038023839275395 0.10822981419271124 0.17750548647807074 With DSCHBV: exp(t*A)e_1 = 9.99388628512555138E-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.53582622006686407E-002 1.20255209080964744E-002 8.77613542815507304E-003 3.49148113952841253E-003 6.46066721625751278E-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.10E+00 ||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.95E-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.88E-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.21814038901391416E-002 1.35176293168635907E-002 1.26164540290726857E-003 1.77127556565795617E-002 1.65319052794742571E-003 1.83453540728859780E-003 1.71223304680269184E-004 1.14179972943451464E-002 1.06567974747221405E-003 final report---------------------------------------- runtime = 0.54E-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.21814038901391416E-002 1.35176293168635907E-002 1.26164540290726857E-003 1.77127556565795617E-002 1.65319052794742571E-003 1.83453540728859780E-003 1.71223304680269184E-004 1.14179972943451464E-002 1.06567974747221405E-003 final report---------------------------------------- runtime = 0.54E-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.97E-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.98E-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.12061291039622999E-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.97E-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.75378622767 , -113.89587171936590 ) ( 6.74919204175804873E-008, 6.86578673704936263E-008) ( -3.0966858825504842 , 4.6959293376943450 ) (-1.67768922280969565E-006,-7.88929404446561570E-006) ( 9.16756563521883951E-006,-2.68608029386347638E-006) ( 6.89110050472356225E-003,-8.28642323619984272E-002) ( 2.63648479299239106E-005, 8.23661994586515414E-004) (-0.20614403650421770 ,-0.20907747166749444 ) ( 5.9467621461748950 ,-0.95593374832015676 ) ( 0.14657967066943786 , 0.23806647697786873 ) final report---------------------------------------- runtime = 0.21E+00 ||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.75378622825 , -113.89587171923773 ) ( 6.74919204175696480E-008, 6.86578673704930969E-008) ( -3.0966858825504757 , 4.6959293376943254 ) (-1.67768922280969544E-006,-7.88929404446575123E-006) ( 9.16756563521874295E-006,-2.68608029386345901E-006) ( 6.89110050472352235E-003,-8.28642323619975946E-002) ( 2.63648479299214915E-005, 8.23661994586511077E-004) (-0.20614403650421689 ,-0.20907747166749352 ) ( 5.9467621461748674 ,-0.95593374832014988 ) ( 0.14657967066943731 , 0.23806647697786812 ) final report---------------------------------------- runtime = 0.72E-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: Normal end of execution. 5 December 2010 9:00:09.770 PM