10 April 2012 9:37:50.877 AM TEST_MAT_PRB FORTRAN77 version Test the TEST_MAT library. TEST_COND Compute the condition number of an example of each test matrix Matrix title N COND AEGERTER 5 24.0000 BAB 5 8.42986 BODEWIG 4 10.4366 COMBIN 3 18.4545 CONEX3 5 80.0000 RUTIS5 4 62608.0 SUMMATION 5 10.0000 TRI_UPPER 5 524.600 WILK03 3 0.260000E+12 WILSON 4 4488.00 TEST_DETERMINANT Compute the determinants of an example of each test matrix; compare with the determinant routine, if available. Print the matrix Frobenius norm for an estimate of magnitude. Matrix title N Determ Determ ||A|| AEGERTER 5 -25.0000 -25.0000 9.43398 ANTICIRCULANT 3 -231.800 -231.800 10.8499 ANTICIRCULANT 4 1397.76 1397.76 12.5857 ANTICIRCULANT 5 7207.92 7207.92 14.1845 ANTIHADAMARD 5 1.00000 1.00000 3.31662 ANTISYMM_RANDOM 5 0.00000 2.87369 ANTISYMM_RANDOM 6 0.973530E-01 3.33451 BAB 5 -2075.15 -2075.15 14.4388 BIMARKOV_RANDOM 5 -0.862803E-04 1.38793 BIS 5 51.5363 51.5363 19.8202 BODEWIG 4 568.000 568.000 12.7279 BOOTHROYD 5 1.00000 1.00000 886.710 BORDERBAND 5 -0.328125 -0.328125 2.76699 CARRY 5 0.165382E-07 0.165382E-07 1.41391 CAUCHY 5 0.689186E-09 0.689186E-09 8.73726 CHEBY_DIFF1 5 -0.213163E-13 13.4722 CHEBY_DIFF1 6 -0.511591E-12 20.7702 CHEBY_T 5 64.0000 64.0000 12.6886 CHEBY_U 5 1024.00 1024.00 22.4277 CHEBY_VAN1 5 0.188744E+08 4896.26 CHEBY_VAN2 2 -2.00000 -2.00000 2.00000 CHEBY_VAN2 3 -1.41421 -1.41421 2.00000 CHEBY_VAN2 4 1.00000 1.00000 2.08167 CHEBY_VAN2 5 0.707107 0.707107 2.17945 CHEBY_VAN2 6 -0.500000 -0.500000 2.28035 CHEBY_VAN2 7 -0.353553 -0.353553 2.38048 CHEBY_VAN2 8 0.250000 0.250000 2.47848 CHEBY_VAN2 9 0.176777 0.176777 2.57391 CHEBY_VAN2 10 -0.125000 -0.125000 2.66667 CHEBY_VAN3 5 13.9754 13.9754 3.87298 CHOW 5 201639. 201639. 69.9110 CIRCULANT 5 7207.92 7207.92 14.1845 CIRCULANT2 3 18.0000 18.0000 6.48074 CIRCULANT2 4 -160.000 -160.000 10.9545 CIRCULANT2 5 1875.00 1875.00 16.5831 CLEMENT1 5 0.00000 0.00000 7.74597 CLEMENT1 6 -225.000 -225.000 10.4881 CLEMENT2 5 0.00000 0.00000 6.32456 CLEMENT2 6 -225.000 -225.000 8.36660 CLEMENT3 5 0.00000 0.00000 9.00000 CLEMENT3 6 -144.000 -144.000 10.3441 COMBIN 5 1175.97 1175.97 50.3925 COMPANION 5 -3.00000 -3.00000 7.00000 COMPLEX_I 2 1.00000 1.00000 1.41421 CONEX1 4 2.20000 2.20000 7.63937 CONEX2 3 0.454545 0.454545 2.65381 CONEX3 5 -1.00000 -1.00000 3.87298 CONFERENCE 6 -125.000 -125.000 5.47723 CREATION 5 0.00000 0.00000 5.47723 DAUB2 4 1.00000 1.00000 2.00000 DAUB4 8 -1.00000 -1.00000 2.82843 DAUB6 12 1.00000 1.00000 3.46410 DAUB8 16 -1.00000 -1.00000 4.00000 DAUB10 20 1.00000 1.00000 4.47214 DAUB12 24 -1.00000 -1.00000 4.89898 DIAGONAL 5 45.0000 45.0000 6.70820 DIF1 5 0.00000 0.00000 2.82843 DIF1CYCLIC 5 0.00000 0.00000 3.16228 DIF2 5 6.00000 6.00000 5.29150 DIF2CYCLIC 5 0.00000 0.00000 5.47723 DORR 5 0.190140E+11 422.499 DOWNSHIFT 5 1.00000 1.00000 2.23607 EBERLEIN 5 0.00000 0.184173E-10 38.0316 EULERIAN 5 1.00000 1.00000 77.2981 EXCHANGE 5 1.00000 1.00000 2.23607 FIBONACCI1 5 0.00000 -0.155847E-40 369.525 FIBONACCI2 5 -1.00000 -1.00000 3.00000 FIBONACCI3 5 8.00000 8.00000 3.60555 FIEDLER 7 1502.87 1502.87 30.1171 FORSYTHE 5 81539.5 81539.5 21.6712 FOURIER_COSINE 5 1.00000 1.00000 2.23607 FOURIER_SINE 5 1.00000 1.00000 2.23607 FRANK 5 1.00000 1.00000 11.6190 GEAR 4 -0.244929E-15 -0.00000 2.82843 GEAR 5 2.00000 2.00000 3.16228 GEAR 6 -4.00000 -4.00000 3.46410 GEAR 7 2.00000 2.00000 3.74166 GEAR 8 0.489859E-15 0.00000 4.00000 GFPP 5 104.858 104.858 7.57628 GIVENS 5 16.0000 16.0000 20.6155 GK316 5 -25.0000 -25.0000 9.43398 GK323 5 32.0000 32.0000 10.0000 GK324 5 12.0384 12.0384 11.4315 GRCAR 5 8.00000 3.60555 HADAMARD 5 0.00000 4.00000 HANKEL 5 -3057.94 15.2774 HANOWA 6 714.498 714.498 7.55248 HARMAN 8 0.954779E-03 0.954779E-03 5.05359 HARTLEY 5 55.9017 55.9017 5.00000 HARTLEY 6 -216.000 -216.000 6.00000 HARTLEY 7 -907.493 -907.493 7.00000 HARTLEY 8 -4096.00 -4096.00 8.00000 HELMERT 5 1.00000 1.00000 2.23607 HELMERT2 5 1.00000 1.00000 2.23607 HERMITE 5 1024.00 1024.00 54.1941 HERNDON 5 -0.400000E-01 -0.400000E-01 1.77133 HILBERT 5 0.374930E-11 0.374930E-11 1.58091 HOUSEHOLDER 5 -1.00000 -1.00000 2.23607 IDEM_RANDOM 5 0.00000 0.00000 1.00000 IDENTITY 5 1.00000 1.00000 2.23607 IJFACT1 5 0.716636E+10 0.716636E+10 0.366559E+07 IJFACT2 5 0.149480E-20 0.149480E-20 0.557720 ILL3 3 6.00000 6.00000 817.763 INTEGRATION 6 1.00000 1.00000 3.61716 INVOL 5 -1.00000 -1.00000 1942.46 INVOL_RANDOM 5 -1.00000 2.23607 JACOBI 5 0.00000 0.00000 1.49071 JACOBI 6 -0.216450E-01 -0.216450E-01 1.65145 JORDAN 6 113.380 113.380 5.83438 KAHAN 5 0.412260E-01 0.412260E-01 1.80785 KERSHAW 4 1.00000 1.00000 8.24621 KERSHAWTRI 5 1.00000 538.766 8.73155 KMS 5 217.433 217.433 41.6571 LAGUERRE 5 0.347222E-02 0.347222E-02 6.85376 LEHMER 5 0.656250E-01 3.28041 LESLIE 4 0.605244 0.605244 1.78414 LESP 5 -42300.0 -42300.0 22.3487 LIETZKE 5 48.0000 48.0000 18.0278 LIGHTS_OUT 25 -0.325405E-29 10.2470 LINE_ADJ 5 0.00000 0.00000 2.82843 LINE_LOOP_ADJ 5 0.00000 -0.00000 3.60555 LOEWNER 5 -29.0825 20.5227 LOTKIN 5 0.187465E-10 0.187465E-10 2.45676 MARKOV_RANDOM 5 0.488558E-02 1.33584 MAXIJ 5 5.00000 5.00000 19.8746 MILNES 5 12.0384 12.0384 11.4315 MINIJ 5 1.00000 1.00000 12.4499 MOLER1 5 1.00000 1.00000 47.5848 MOLER2 5 0.00000 0.102538E-06 101035. MOLER3 5 1.00000 1.00000 8.66025 NEUMANN 25 0.00000 0.124631E-02 23.2379 ONE 5 0.00000 -0.00000 5.00000 ORTEGA 5 -20.5920 -20.5920 248.221 ORTH_RANDOM 5 1.00000 1.00000 2.23607 ORTH_SYMM 5 1.00000 1.00000 2.23607 OTO 5 6.00000 6.00000 5.29150 PARTER 5 131.917 131.917 6.34077 PASCAL1 5 1.00000 1.00000 9.94987 PASCAL2 5 1.00000 1.00000 92.4608 PASCAL3 5 1.00000 1.00000 122.629 PDS_RANDOM 5 0.404187E-01 0.404187E-01 1.46230 PEI 5 168.664 168.664 8.43801 PERMUTATION_RANDOM 5 1.00000 1.00000 2.23607 PLU 5 0.193261E+08 0.193261E+08 152.462 POISSON 25 0.325655E+14 0.325655E+14 21.9089 PROLATE 5 1644.22 9.84461 RECTANGLE_ADJ 25 0.00000 0.00000 8.94427 REDHEFFER 5 -2.00000 -2.00000 3.74166 REF_RANDOM 5 0.00000 0.00000 2.63560 REF_RANDOM 5 1.00000 1.00000 2.81894 RIEMANN 5 96.0000 8.83176 RING_ADJ 1 1.00000 1.00000 1.00000 RING_ADJ 2 -1.00000 -1.00000 1.41421 RING_ADJ 3 2.00000 2.00000 2.44949 RING_ADJ 4 0.00000 0.00000 2.82843 RING_ADJ 5 2.00000 2.00000 3.16228 RING_ADJ 6 -4.00000 -4.00000 3.46410 RING_ADJ 7 2.00000 2.00000 3.74166 RING_ADJ 8 0.00000 0.00000 4.00000 RIS 5 4.12239 4.12239 3.17039 RODMAN 5 -2126.84 -2126.84 12.7201 ROSSER1 8 0.00000 -9480.58 2482.26 ROUTH 5 75.7680 75.7680 7.74855 RUTIS1 4 -375.000 -375.000 16.6132 RUTIS2 4 100.000 100.000 11.4018 RUTIS3 4 624.000 624.000 14.1421 RUTIS4 5 216.000 216.000 59.1270 RUTIS5 4 1.00000 1.00000 23.7697 SCHUR_BLOCK 5 572.032 572.032 8.38093 SKEW_CIRCULANT 5 -10057.0 -10057.0 14.1845 SPLINE 5 -2444.28 -2444.28 20.6417 STIRLING 5 1.00000 1.00000 67.9191 STRIPE 5 2112.00 14.8324 SUMMATION 5 1.00000 1.00000 3.87298 SWEET1 6 -0.204682E+08 70.1997 SWEET2 6 9562.52 30.1433 SWEET3 6 -0.540561E+08 73.4234 SWEET4 13 -0.646348E+17 119.704 SYLVESTER 5 -219.525 12.5092 SYMM_RANDOM 5 19.7837 19.7837 6.34350 TOEPLITZ 5 -3057.94 15.2774 TOEPLITZ_5DIAG 5 -672.905 12.7562 TOEPLITZ_5S 25 -0.141078E+18 40.2651 TOEPLITZ_PDS 5 0.424591E-03 2.82557 TOURNAMENT_RANDOM 5 0.00000 0.00000 4.47214 TRANSITION_RANDOM 5 0.486764E-02 1.32331 TRENCH 5 -35.0147 6.99866 TRI_UPPER 5 1.00000 1.00000 9.13236 TRIS 5 6656.04 6656.04 13.3462 TRIV 5 -712.622 -712.622 11.1893 TRIW 5 1.00000 1.00000 9.46784 UPSHIFT 5 1.00000 1.00000 2.23607 VAND1 5 132934. 132934. 476.382 VAND2 5 132934. 132934. 476.382 WATHEN 96 0.116235+164 1324.87 WILK03 3 0.900000E-20 0.900000E-20 1.39284 WILK04 4 0.442923E-16 0.442923E-16 1.89545 WILK05 5 0.379950E-14 0.379947E-14 1.51485 WILK12 12 1.00000 1.00000 53.5910 WILK20 20 0.146801E+26 102.361 WILK21 21 -0.415825E+13 -0.415825E+13 28.4605 WILSON 4 1.00000 1.00000 30.5450 ZERO 5 0.00000 0.00000 0.00000 ZIELKE 5 446.457 13.7099 TEST_EIGEN Compute the Frobenius norm of the eigenvalue error: A * X - X * LAMBDA given a set of K eigenvectors X and eigenvalues LAMBDA. Matrix title N K ||A|| ||(A-Lambda*I)*X|| BODEWIG 4 4 12.7279 0.917346E-14 CARRY 5 5 1.41391 0.117642E-14 CHOW 5 5 69.9110 0.241480E-12 COMBIN 5 5 50.3925 0.753644E-14 DIF2 5 5 5.29150 0.119584E-14 EXCHANGE 5 5 2.23607 0.00000 IDEM_RANDOM 5 5 1.73205 0.136010E-14 IDENTITY 5 5 2.23607 0.00000 ILL3 3 3 817.763 0.162356E-10 KMS 5 5 2.32288 0.320550E-07 ONE 5 5 5.00000 0.00000 ORTEGA 5 5 248.221 0.632798E-12 OTO 5 5 5.29150 0.119584E-14 PDS_RANDOM 5 5 1.46230 0.172673E-14 PEI 5 5 8.43801 0.00000 RODMAN 5 5 10.0896 0.00000 ROSSER1 8 8 2482.26 0.261994E-10 RUTIS1 4 4 16.6132 0.00000 RUTIS2 4 4 11.4018 0.00000 RUTIS5 4 4 23.7697 0.146286E-13 WILK12 12 12 53.5910 0.101528E-06 WILSON 4 4 30.5450 0.248731E-13 ZERO 5 5 0.00000 0.00000 TEST_INVERSE A = a test matrix of order N; B = inverse as computed by a routine. C = inverse as computed by R8MAT_INVERSE. ||A|| = Frobenius norm of A. ||I-AB|| = Frobenius norm of I-A*B. ||I-AC|| = Frobenius norm of I-A*C. Matrix title N ||A|| ||I-AB|| ||I-AC|| AEGERTER 5 9.43398 0.710890E-15 0.705075E-15 BAB 5 14.4388 0.677600E-15 0.536803E-15 BERNSTEIN 5 25.2784 0.00000 0.00000 BIS 5 19.8202 0.00000 0.00000 BODEWIG 4 12.7279 0.708784E-15 0.866075E-15 BOOTHROYD 5 886.710 0.00000 0.511553E-10 BORDERBAND 5 2.76699 0.00000 0.00000 CARRY 5 1.41391 0.117630E-12 0.240612E-12 CAUCHY 5 8.73726 0.367556E-09 0.158004E-09 CHEBY_T 5 12.6886 0.00000 0.00000 CHEBY_U 5 22.4277 0.00000 0.00000 CHEBY_VAN2 5 2.17945 0.591396E-15 0.474151E-15 CHEBY_VAN3 5 3.87298 0.954577E-15 0.114785E-14 CHOW 5 69.9110 0.215951E-14 0.166914E-14 CIRCULANT --- NOT READY CIRCULANT2 5 16.5831 0.180915E-14 0.114595E-14 CLEMENT1 6 10.4881 0.00000 0.00000 CLEMENT2 6 8.36660 0.00000 0.647646E-15 CLEMENT3 6 10.2215 0.114947E-14 0.127415E-14 COMBIN 5 50.3925 0.287803E-14 0.615447E-14 COMPANION 5 7.00000 0.222045E-15 0.00000 COMPLEX_I 2 1.41421 0.00000 0.00000 CONEX1 4 7.63937 0.00000 0.00000 CONEX2 3 2.65381 0.00000 0.00000 CONEX3 5 3.87298 0.00000 0.00000 CONFERENCE 6 5.47723 0.00000 0.742465E-15 DAUB2 4 2.00000 0.888178E-15 0.00000 DAUB4 8 2.82843 0.210011E-14 0.396675E-15 DAUB6 12 3.46410 0.139744E-14 0.113170E-14 DAUB8 16 4.00000 0.464930E-14 0.176073E-14 DAUB10 20 4.47214 0.869739E-14 0.165358E-14 DAUB12 24 4.89898 0.195470E-13 0.224395E-14 DIAGONAL 5 6.34350 0.00000 0.00000 DIF2 5 5.29150 0.686635E-15 0.113951E-14 DOWNSHIFT 5 2.23607 0.00000 0.00000 EULERIAN 5 77.2981 0.00000 0.250443E-12 EXCHANGE 5 2.23607 0.00000 0.00000 FIBONACCI2 5 3.00000 0.00000 0.00000 FIBONACCI3 5 3.60555 0.00000 0.157009E-15 FIEDLER 7 30.1171 0.636935E-14 0.182766E-13 FORSYTHE 5 21.6712 0.918189E-17 0.226462E-15 FOURIER_COSINE 5 2.23607 0.117732E-14 0.988355E-15 FOURIER_SINE 5 2.23607 0.198448E-14 0.871782E-15 FRANK 5 11.6190 0.00000 0.347485E-13 GFPP 5 3.07848 0.472623E-15 0.432359E-15 GIVENS 5 20.6155 0.00000 0.00000 GK316 5 9.43398 0.710890E-15 0.705075E-15 GK323 5 10.0000 0.00000 0.00000 GK324 5 11.4315 0.328589E-14 0.174132E-14 HANOWA 8 7.77056 0.628037E-15 0.315243E-15 HARMAN 8 5.05359 0.954191E-06 0.709414E-14 HARTLEY 5 5.00000 0.269465E-14 0.697215E-15 HELMERT 5 2.23607 0.741186E-15 0.509651E-15 HELMERT2 5 2.23607 0.319070E-15 0.735070E-15 HERMITE 5 54.1941 0.00000 0.00000 HERNDON 5 1.77133 0.710890E-15 0.178578E-14 HILBERT 5 1.58091 0.727596E-11 0.201079E-10 HOUSEHOLDER 5 2.23607 0.701160E-15 0.110659E-14 IDENTITY 5 2.23607 0.00000 0.00000 ILL3 3 817.763 0.471661E-10 0.163698E-10 INTEGRATION 6 3.61716 0.522116E-15 0.00000 INVOL 5 1942.46 0.727596E-11 0.840355E-10 JORDAN 6 5.83438 0.196262E-16 0.208167E-16 KAHAN 5 1.80785 0.379646E-15 0.362693E-15 KERSHAW 4 8.24621 0.00000 0.444089E-14 KERSHAWTRI 5 8.73155 0.772264E-15 0.753655E-15 KMS 5 41.6571 0.693334E-14 0.130038E-13 LAGUERRE 5 6.85376 0.00000 0.173921E-13 LEGENDRE 5 6.80762 0.268032E-15 0.248253E-15 LEHMER 5 3.28041 0.141744E-14 0.220341E-14 LIETZKE 5 18.0278 0.695553E-15 0.489906E-14 LOTKIN 5 2.45676 0.00000 0.288028E-10 MAXIJ 5 19.8746 0.00000 0.137775E-14 MILNES 5 11.4315 0.327303E-14 0.174132E-14 MINIJ 5 12.4499 0.00000 0.00000 MOLER1 5 47.5848 0.156584E-12 0.292922E-12 MOLER3 5 8.66025 0.00000 0.00000 ORTEGA 5 248.221 0.194889E-11 0.131369E-11 ORTH_SYMM 5 2.23607 0.196898E-14 0.843169E-15 OTO 5 5.29150 0.686635E-15 0.113951E-14 PARTER 5 6.34077 0.696142E-16 0.748746E-15 PASCAL1 5 9.94987 0.00000 0.00000 PASCAL2 5 92.4608 0.00000 0.00000 PASCAL3 5 122.629 0.901578E-13 0.287134E-12 PDS_RANDOM 5 1.46230 0.960755E-14 0.718843E-15 PEI 5 8.43801 0.108211E-14 0.715962E-15 PERMUTATION_RANDOM 5 2.23607 0.00000 0.00000 PLU 5 152.462 0.129789E-14 0.125372E-14 RIS 5 3.17039 0.837173E-16 0.735818E-15 RODMAN 5 12.7201 0.112127E-14 0.574229E-15 RUTIS1 4 16.6132 0.105471E-14 0.151741E-14 RUTIS2 4 11.4018 0.683824E-15 0.599455E-15 RUTIS3 4 14.1421 0.601660E-15 0.860331E-15 RUTIS4 5 59.1270 0.288694E-12 0.933214E-12 RUTIS5 4 23.7697 0.00000 0.480611E-11 SCHUR_BLOCK 5 8.38093 0.785046E-16 0.961481E-16 SPLINE 5 20.6417 0.967194E-15 0.498071E-15 STIRLING 5 67.9191 0.00000 0.376704E-13 SUMMATION 5 3.87298 0.00000 0.00000 TRI_UPPER 5 9.13236 0.923022E-14 0.923022E-14 TRIS 5 13.3462 0.641687E-15 0.466457E-15 TRIV 5 11.1893 0.128151E-14 0.154380E-14 TRIW 5 9.46784 0.00000 0.00000 UPSHIFT 5 2.23607 0.00000 0.00000 VAND1 5 476.382 0.112190E-13 0.793837E-14 VAND2 5 476.382 0.112190E-13 0.456693E-13 WILK03 3 1.39284 0.674350E-06 0.674350E-06 WILK04 4 1.89545 10.7174 0.482362E-04 WILK05 5 1.51485 0.122740E-08 0.742989E-09 WILK21 21 28.4605 0.382082E-14 0.155131E-14 WILSON 4 30.5450 0.00000 0.636941E-12 TEST_NULL A = a test matrix of order M by N x = an N vector, candidate for a null vector. ||A|| = Frobenius norm of A. ||x|| = L2 norm of x. ||A*x||/||x|| = L2 norm of A*x over L2 norm of x. Matrix title M N ||A|| ||x|| ||A*x||/||x|| ARCHIMEDES 7 8 93.3970 0.187697E+08 0.00000 CHEBY_DIFF1 5 5 13.4722 2.23607 0.649741E-15 CREATION 5 5 5.47723 1.00000 0.00000 DIF1 5 5 2.82843 1.73205 0.00000 DIF1CYCLIC 5 5 3.16228 2.23607 0.00000 DIF2CYCLIC 5 5 5.47723 2.23607 0.00000 EBERLEIN (left) 5 5 18.0665 2.23607 0.561733E-15 FIBONACCI1 5 5 96.9208 1.73205 0.00000 LAUCHLI (left) 6 5 6.64831 3.58329 0.00000 LINE_ADJ 7 7 3.46410 2.00000 0.00000 MOLER2 5 5 101035. 1016.30 0.00000 NEUMANN 25 25 23.2379 5.00000 0.00000 ONE 5 5 5.00000 1.41421 0.00000 RING_ADJ 12 12 4.89898 3.46410 0.00000 ROSSER1 8 8 2482.26 22.3607 0.00000 ZERO 5 5 0.00000 2.23607 0.00000 TEST_PLU A = a test matrix of order M by N P, L, U are the PLU factors. ||A|| = Frobenius norm of A. ||A-PLU|| = Frobenius norm of A-P*L*U. Matrix title M N ||A|| ||A-PLU|| BODEWIG 4 4 12.7279 0.412430E-14 BORDERBAND 5 5 2.76699 0.00000 DIF2 5 5 5.29150 0.00000 GFPP 5 5 7.57628 0.158882E-13 GIVENS 5 5 20.6155 0.00000 KMS 5 5 41.6571 0.131968E-12 MAXIJ 5 5 19.8746 0.00000 MINIJ 5 5 12.4499 0.00000 MOLER1 5 5 47.5848 0.00000 MOLER3 5 5 8.66025 0.00000 OTO 5 5 5.29150 0.00000 PASCAL2 5 5 92.4608 0.00000 WILSON 4 4 30.5450 0.732411E-14 TEST_SOLUTION Compute the Frobenius norm of the solution error: A * X - B given MxN matrix A, NxK solution X, MxK right hand side B. Matrix title M N K ||A|| ||A*X-B|| BODEWIG 4 4 1 12.7279 0.00000 DIF2 10 10 2 7.61577 0.00000 FRANK 10 10 2 38.6652 0.00000 POISSON 20 20 1 19.5448 0.00000 WILK03 3 3 1 1.39284 0.674350E-06 WILK04 4 4 1 1.89545 0.129463E+17 WILSON 4 4 1 30.5450 0.00000 TEST_MAT_PRB Normal end of execution. 10 April 2012 9:37:50.888 AM