22 January 2020 8:44:21.337 AM EISPACK_TEST FORTRAN77 version. Test the EISPACK library. BALANC_TEST BALANC balances a real general matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 110.000 12.0000 13.0000 0.00000 15.0000 2 0.00000 22.0000 0.00000 0.00000 0.00000 3 310.000 32.0000 33.0000 0.00000 35.0000 4 0.00000 0.00000 43.0000 44.0000 0.00000 5 510.000 0.00000 53.0000 0.00000 55.0000 LOW = 2 IGH = 4 Scaling vector SCALE: 1 4.0000000 2 1.0000000 3 1.0000000 4 0.62500000E-01 5 2.0000000 The balanced matrix A: Col 1 2 3 4 5 Row 1 44.0000 0.00000 43.0000 0.00000 0.00000 2 0.00000 55.0000 53.0000 31.8750 0.00000 3 0.00000 35.0000 33.0000 19.3750 32.0000 4 0.00000 240.000 208.000 110.000 192.000 5 0.00000 0.00000 0.00000 0.00000 22.0000 BANDV_TEST BANDV computes the eigenvectors of a real symmetric band matrix. Matrix order = 5 Half bandwidth + 1 = 2 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 The eigenvector matrix X: Col 1 2 3 4 5 Row 1 -0.288675 -0.500000 0.577350 0.500000 -0.288675 2 -0.500000 -0.500000 -0.00000 -0.500000 0.500000 3 -0.577350 -0.218463E-14 -0.577350 -0.178742E-14 -0.577350 4 -0.500000 0.500000 -0.128198E-14 0.500000 0.500000 5 -0.288675 0.500000 0.577350 -0.500000 -0.288675 The residual (A-Lambda*I)*X: Col 1 2 3 4 5 Row 1 0.596745E-15 0.111022E-14 0.00000 0.111022E-14 0.666134E-15 2 0.388578E-15 0.111022E-14 0.122125E-14 0.333067E-14 0.111022E-14 3 -0.152656E-14 0.107441E-14 0.133227E-14 0.114343E-14 0.133227E-14 4 0.693889E-15 0.116573E-14 0.123168E-14 0.888178E-15 0.444089E-15 5 0.610623E-15 0.105471E-14 0.133227E-14 0.00000 0.666134E-15 BISECT_TEST BISECT computes some eigenvalues of a real symmetric tridiagonal matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 BQR_TEST BQR computes some eigenvalues of a real symmetric band matrix. Matrix order = 5 Half bandwidth+1 = 2 The compressed matrix A: Col 1 2 Row 1 0.00000 2.00000 2 -1.00000 2.00000 3 -1.00000 2.00000 4 -1.00000 2.00000 5 -1.00000 2.00000 Eigenvalues: 1 0.26794919243112286 2 1.0000000000000009 3 2.0000000000000013 4 3.0000000000000004 5 3.7320508075688794 CBAL_TEST CBAL balances a complex general matrix. Matrix order = 5 The matrix AR: Col 1 2 3 4 5 Row 1 110.000 12.0000 13.0000 0.00000 15.0000 2 0.00000 22.0000 0.00000 0.00000 0.00000 3 310.000 32.0000 33.0000 0.00000 35.0000 4 0.00000 0.00000 43.0000 44.0000 0.00000 5 510.000 0.00000 53.0000 0.00000 55.0000 The matrix AI: Col 1 2 3 4 5 Row 1 110.500 12.5000 13.5000 0.00000 15.5000 2 0.00000 22.5000 0.00000 0.00000 0.00000 3 310.500 32.5000 33.5000 0.00000 35.5000 4 0.00000 0.00000 43.5000 44.5000 0.00000 5 510.500 0.00000 53.5000 0.00000 55.5000 LOW = 2 IGH = 4 Scaling vector SCALE: 1 4.0000000 2 1.0000000 3 1.0000000 4 0.62500000E-01 5 2.0000000 The balanced matrix AR: Col 1 2 3 4 5 Row 1 44.0000 0.00000 43.0000 0.00000 0.00000 2 0.00000 55.0000 53.0000 31.8750 0.00000 3 0.00000 35.0000 33.0000 19.3750 32.0000 4 0.00000 240.000 208.000 110.000 192.000 5 0.00000 0.00000 0.00000 0.00000 22.0000 The balanced matrix AI: Col 1 2 3 4 5 Row 1 44.5000 0.00000 43.5000 0.00000 0.00000 2 0.00000 55.5000 53.5000 31.9062 0.00000 3 0.00000 35.5000 33.5000 19.4062 32.5000 4 0.00000 248.000 216.000 110.500 200.000 5 0.00000 0.00000 0.00000 0.00000 22.5000 CG_LR_TEST CG_LR computes the eigenvalues and eigenvectors of a complex general matrix using elementary transformations. Matrix order = 4 Real and imaginary parts of eigenvalues: 1 4.82843 0.00000 2 4.00000 0.00000 3 -0.135688E-15 0.00000 4 -0.828427 0.00000 CG_LR_TEST CG_LR computes the eigenvalues and eigenvectors of a complex general matrix using elementary transformations. Matrix order = 4 Real and imaginary parts of eigenvalues: 1 4.82843 0.00000 2 4.00000 0.00000 3 -0.135688E-15 0.00000 4 -0.828427 0.00000 Eigenvector 1 4.82843 0.00000 -0.436130 0.00000 0.00000 -0.744309 0.00000 0.136143 Eigenvector 2 4.82843 0.00000 0.436130 0.00000 0.00000 0.744309 0.00000 0.136143 Eigenvector 3 0.00000 2.00000 0.00000 0.436130 0.744309 0.00000 0.328679 0.00000 Eigenvector 4 0.00000 -2.00000 0.00000 0.436130 0.744309 0.00000 -0.328679 0.00000 CG_QR_TEST CG_QR computes the eigenvalues and eigenvectors of a complex general matrix using unitary transformations. Matrix order = 4 Real and imaginary parts of eigenvalues: 1 4.82843 0.00000 2 4.00000 0.00000 3 -0.162724E-17 0.00000 4 -0.828427 0.00000 CG_QR_TEST CG_QR computes the eigenvalues and eigenvectors of a complex general matrix using unitary transformations. Matrix order = 4 Real and imaginary parts of eigenvalues: 1 4.82843 0.00000 2 4.00000 0.00000 3 -0.162724E-17 0.00000 4 -0.828427 0.00000 Eigenvector 1 0.653281 0.00000 0.500000 0.00000 -0.500000 0.00000 0.270598 0.00000 Eigenvector 2 0.653281 0.00000 -0.500000 0.00000 0.500000 0.00000 0.270598 0.00000 Eigenvector 3 0.00000 0.270598 0.00000 -0.500000 0.00000 -0.500000 0.00000 -0.653281 Eigenvector 4 0.00000 -0.270598 0.00000 -0.500000 0.00000 -0.500000 0.00000 0.653281 CH_TEST CH computes the eigenvalues and eigenvectors of a complex hermitian matrix. Matrix order = 4 The eigenvalues Lambda: 1 -0.82842712 2 0.12894578E-16 3 4.0000000 4 4.8284271 CH_TEST CH computes the eigenvalues and eigenvectors of a complex hermitian matrix. Matrix order = 4 The eigenvalues Lambda: 1 -0.82842712 2 0.13322676E-14 3 4.0000000 4 4.8284271 Eigenvector 1 0.00000 0.270598 0.00000 0.500000 0.00000 0.500000 0.00000 0.653281 Eigenvector 2 0.00000 0.270598 0.00000 -0.500000 0.00000 -0.500000 0.00000 0.653281 Eigenvector 3 0.653281 0.00000 -0.500000 0.00000 0.500000 0.00000 -0.270598 0.00000 Eigenvector 4 -0.653281 0.00000 -0.500000 0.00000 0.500000 -0.00000 0.270598 -0.00000 CH3_TEST CH3 computes the eigenvalues and eigenvectors of a complex hermitian matrix. Matrix order = 4 The eigenvalues Lambda: 1 -0.82842712 2 0.12894578E-16 3 4.0000000 4 4.8284271 CH3_TEST CH3 computes the eigenvalues and eigenvectors of a complex hermitian matrix. Matrix order = 4 The eigenvalues Lambda: 1 -0.82842712 2 0.13322676E-14 3 4.0000000 4 4.8284271 Eigenvector 1 0.00000 -0.270598 0.00000 -0.500000 0.00000 -0.500000 0.00000 -0.653281 Eigenvector 2 0.00000 -0.270598 0.00000 0.500000 0.00000 0.500000 0.00000 -0.653281 Eigenvector 3 0.653281 0.00000 -0.500000 0.00000 0.500000 0.00000 -0.270598 0.00000 Eigenvector 4 -0.653281 0.00000 -0.500000 0.00000 0.500000 -0.00000 0.270598 -0.00000 CINVIT_TEST CINVIT computes the eigenvectors of a complex Hessenberg matrix. Matrix order = 4 Real and imaginary parts of eigenvalues: 1 20.2603 16.6458 2 0.568541 6.82620 3 3.32443 -2.74203 4 -5.15323 -8.72994 Matrix Real Part Ar: Col 1 2 3 4 Row 1 4.00000 7.00000 7.00000 7.00000 2 9.00000 8.00000 8.00000 2.00000 3 0.00000 8.00000 7.00000 7.00000 4 0.00000 0.00000 4.00000 0.00000 Matrix Imag Part Ai: Col 1 2 3 4 Row 1 8.00000 6.00000 10.0000 10.0000 2 9.00000 1.00000 10.0000 5.00000 3 0.00000 3.00000 2.00000 8.00000 4 0.00000 0.00000 10.0000 1.00000 Real and imaginary parts of eigenvalues: 1 20.2603 16.6458 2 0.568541 6.82620 3 3.32443 -2.74203 4 -5.15323 -8.72994 Eigenvector 1 0.952157 0.250710 1.00000 0.00000 -0.330070 -0.222298 0.335491 -0.680030E-01 Eigenvector 2 -1.00000 -0.444087E-17 -0.503206 0.824242 -1.00000 -0.444087E-17 0.768179 -0.105140E-01 Eigenvector 3 -0.501509 0.172182 0.215017 -0.275480 -0.257418 -0.309094 -1.00000 0.00000 Eigenvector 4 -0.224707 -0.460156E-01 0.244894 -0.613996 0.866321 -0.203311 1.00000 0.106965 IMTQLV_TEST IMTQLV computes the eigenvalues of a real symmetric tridiagonal matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 INVIT_TEST INVIT computes the eigenvectors of a real Hessenberg matrix. Matrix order = 5 Matrix A: Col 1 2 3 4 5 Row 1 9.00000 4.00000 1.00000 3.00000 2.00000 2 4.00000 3.00000 1.00000 7.00000 1.00000 3 0.00000 3.00000 1.00000 2.00000 4.00000 4 0.00000 0.00000 5.00000 5.00000 1.00000 5 0.00000 0.00000 0.00000 1.00000 2.00000 Real and imaginary parts of eigenvalues: 1 11.9421 0.00000 2 -0.484651 3.05040 3 -0.484651 -3.05040 4 7.23209 0.00000 5 1.79507 0.00000 Eigenvector: 1 0.576743 0.00000 2 1.00000 0.00000 3 -0.108856 0.00000 4 -0.795549E-01 0.00000 5 0.800179E-02 0.00000 Eigenvector: 1 -0.929264E-01 -0.660226E-01 2 1.00000 0.00000 3 0.237405E-01 0.323168 4 0.118671 -0.231110 5 0.645954E-01 -0.137115E-01 Eigenvector: 1 -0.929264E-01 0.660226E-01 2 1.00000 -0.00000 3 0.237405E-01 -0.323168 4 0.118671 0.231110 5 0.645954E-01 0.137115E-01 Eigenvector: 1 1.00000 0.00000 2 0.111558 0.00000 3 0.206404 0.00000 4 0.505654 0.00000 5 -0.966447E-01 0.00000 Eigenvector: 1 -0.134922 0.00000 2 1.00000 0.00000 3 -0.188338E-01 0.00000 4 -0.562261E-01 0.00000 5 -0.274369 0.00000 Residuals (A*x-Lambda*x) for eigenvalue 1 1.97163 0.00000 -7.29291 0.00000 4.06401 0.00000 0.160036E-01 0.00000 -0.159110 0.00000 Residuals (A*x-Lambda*x) for eigenvalue 2 3.42618 -0.740324 4.03198 -4.62280 4.51676 -0.109692 0.129191 -0.274230E-01 0.237343 -0.462220 Residuals (A*x-Lambda*x) for eigenvalue 3 3.42618 0.740324 4.03198 4.62280 4.51676 0.109692 0.129191 0.274230E-01 0.237343 0.462220 Residuals (A*x-Lambda*x) for eigenvalue 4 3.74422 0.00000 7.17721 0.00000 -0.326924 0.00000 -0.193289 0.00000 1.01131 0.00000 Residuals (A*x-Lambda*x) for eigenvalue 5 2.29164 0.00000 -0.215472E-01 0.00000 1.80504 0.00000 -0.548739 0.00000 -0.112452 0.00000 MINFIT_TEST MINFIT solves an overdetermined linear system using least squares methods. Matrix rows = 5 Matrix columns = 2 The matrix A: Col 1 2 Row 1 1.00000 1.00000 2 2.05000 -1.00000 3 3.06000 1.00000 4 -1.02000 2.00000 5 4.08000 -1.00000 The right hand side B: Col 1 Row 1 1.98000 2 0.950000 3 3.98000 4 0.920000 5 2.90000 MINFIT error code IERR = 0 The singular values: 1 5.7385075 2 2.7059992 The least squares solution X: 1 0.96310140 2 0.98854334 The residual A * X - B: 1 -0.28355256E-01 2 0.35814526E-01 3 -0.44366372E-01 4 0.74723261E-01 5 0.40910368E-01 RG_ELM_TEST RG_ELM computes the eigenvalues and eigenvectors of a real general matrix, using elementary transforms. Matrix order = 3 The matrix A: Col 1 2 3 Row 1 33.0000 16.0000 72.0000 2 -24.0000 -10.0000 -57.0000 3 -8.00000 -4.00000 -17.0000 Real and imaginary parts of eigenvalues: 1 3.00000 0.00000 2 1.00000 0.00000 3 2.00000 0.00000 RG_ELM_TEST RG_ELM computes the eigenvalues and eigenvectors of a real general matrix, using elementary transforms. Matrix order = 3 The matrix A: Col 1 2 3 Row 1 33.0000 16.0000 72.0000 2 -24.0000 -10.0000 -57.0000 3 -8.00000 -4.00000 -17.0000 Real and imaginary parts of eigenvalues: 1 3.00000 0.00000 2 1.00000 0.00000 3 2.00000 0.00000 Eigenvector: 1 1.00000 0.00000 2 -0.750000 0.00000 3 -0.250000 0.00000 Eigenvector: 1 -1.00000 0.00000 2 0.800000 0.00000 3 0.266667 0.00000 Eigenvector: 1 1.00000 0.00000 2 -0.812500 0.00000 3 -0.250000 0.00000 Residuals (A*x-Lambda*x) for eigenvalue 1 0.00000 0.00000 0.444089E-15 0.00000 0.00000 0.00000 Residuals (A*x-Lambda*x) for eigenvalue 2 0.00000 0.00000 0.666134E-15 0.00000 -0.555112E-16 0.00000 Residuals (A*x-Lambda*x) for eigenvalue 3 -0.355271E-14 0.00000 0.355271E-14 0.00000 0.999201E-15 0.00000 RG_ORT_TEST RG_ORT computes the eigenvalues and eigenvectors of a real general matrix, using orthogonal transforms. Matrix order = 3 The matrix A: Col 1 2 3 Row 1 33.0000 16.0000 72.0000 2 -24.0000 -10.0000 -57.0000 3 -8.00000 -4.00000 -17.0000 Real and imaginary parts of eigenvalues: 1 3.00000 0.00000 2 1.00000 0.00000 3 2.00000 0.00000 RG_ORT_TEST RG_ORT computes the eigenvalues and eigenvectors of a real general matrix, using orthogonal transforms. Matrix order = 3 The matrix A: Col 1 2 3 Row 1 33.0000 16.0000 72.0000 2 -24.0000 -10.0000 -57.0000 3 -8.00000 -4.00000 -17.0000 Real and imaginary parts of eigenvalues: 1 3.00000 0.00000 2 1.00000 0.00000 3 2.00000 0.00000 Eigenvector: 1 1.00000 0.00000 2 -0.750000 0.00000 3 -0.250000 0.00000 Eigenvector: 1 1.00000 0.00000 2 -0.800000 0.00000 3 -0.266667 0.00000 Eigenvector: 1 1.00000 0.00000 2 -0.812500 0.00000 3 -0.250000 0.00000 Residuals (A*x-Lambda*x) for eigenvalue 1 0.355271E-14 0.00000 -0.111022E-13 0.00000 -0.421885E-14 0.00000 Residuals (A*x-Lambda*x) for eigenvalue 2 -0.355271E-14 0.00000 -0.965894E-14 0.00000 -0.233147E-14 0.00000 Residuals (A*x-Lambda*x) for eigenvalue 3 0.532907E-14 0.00000 -0.117684E-13 0.00000 -0.366374E-14 0.00000 RGG_TEST: RGG for real generalized problem. Find scalars LAMBDA and vectors x so that A*x = LAMBDA * B * x Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 -7.00000 7.00000 6.00000 6.00000 2 -10.0000 8.00000 10.0000 8.00000 3 -8.00000 3.00000 10.0000 11.0000 4 -4.00000 0.00000 4.00000 12.0000 The matrix B: Col 1 2 3 4 Row 1 2.00000 1.00000 0.00000 0.00000 2 1.00000 2.00000 1.00000 0.00000 3 0.00000 1.00000 2.00000 1.00000 4 0.00000 0.00000 1.00000 2.00000 Real and imaginary parts of eigenvalues: 1 2.00000 0.00000 2 1.00000 0.00000 3 3.00000 0.00000 4 4.00000 0.00000 RGG_TEST: RGG for real generalized problem. Find scalars LAMBDA and vectors x so that A*x = LAMBDA * B * x Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 -7.00000 7.00000 6.00000 6.00000 2 -10.0000 8.00000 10.0000 8.00000 3 -8.00000 3.00000 10.0000 11.0000 4 -4.00000 0.00000 4.00000 12.0000 The matrix B: Col 1 2 3 4 Row 1 2.00000 1.00000 0.00000 0.00000 2 1.00000 2.00000 1.00000 0.00000 3 0.00000 1.00000 2.00000 1.00000 4 0.00000 0.00000 1.00000 2.00000 Real and imaginary parts of eigenvalues: 1 2.00000 0.00000 2 1.00000 0.00000 3 3.00000 0.00000 4 4.00000 0.00000 The eigenvectors are: Eigenvector 1 1.00000 1.00000 -1.00000 -1.00000 Eigenvector 2 1.00000 0.750000 -1.00000 -1.00000 Eigenvector 3 0.666667 0.500000 -1.00000 -1.00000 Eigenvector 4 0.333333 0.250000 -0.500000 -1.00000 Residuals (A*x-(Alfr+Alfi*I)*B*x) for eigenvalue 1 0.355271E-14 -0.355271E-14 -0.632827E-14 -0.244249E-14 Residuals (A*x-(Alfr+Alfi*I)*B*x) for eigenvalue 2 0.00000 -0.482947E-14 -0.696665E-14 -0.277556E-14 Residuals (A*x-(Alfr+Alfi*I)*B*x) for eigenvalue 3 -0.399680E-14 0.444089E-14 0.488498E-14 0.888178E-15 Residuals (A*x-(Alfr+Alfi*I)*B*x) for eigenvalue 4 0.00000 0.666134E-14 0.710543E-14 0.532907E-14 RS_TEST RS computes the eigenvalues and eigenvectors of a real symmetric matrix. Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 5.00000 4.00000 1.00000 1.00000 2 4.00000 5.00000 1.00000 1.00000 3 1.00000 1.00000 4.00000 2.00000 4 1.00000 1.00000 2.00000 4.00000 The eigenvalues Lambda: 1 1.0000000 2 2.0000000 3 5.0000000 4 10.000000 RS_TEST RS computes the eigenvalues and eigenvectors of a real symmetric matrix. Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 5.00000 4.00000 1.00000 1.00000 2 4.00000 5.00000 1.00000 1.00000 3 1.00000 1.00000 4.00000 2.00000 4 1.00000 1.00000 2.00000 4.00000 The eigenvalues Lambda: 1 1.0000000 2 2.0000000 3 5.0000000 4 10.000000 The eigenvector matrix: Col 1 2 3 4 Row 1 0.707107 -0.971445E-16 0.316228 0.632456 2 -0.707107 -0.277556E-16 0.316228 0.632456 3 0.00000 0.707107 -0.632456 0.316228 4 0.00000 -0.707107 -0.632456 0.316228 The residual (A-Lambda*I)*X: Col 1 2 3 4 Row 1 0.00000 -0.360822E-15 -0.888178E-15 0.00000 2 0.00000 -0.499600E-15 0.00000 0.00000 3 0.00000 0.888178E-15 -0.266454E-14 -0.133227E-14 4 0.00000 -0.111022E-14 0.00000 0.00000 RSB_TEST RSB computes the eigenvalues and eigenvectors of a real symmetric band matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 RSB_TEST RSB computes the eigenvalues and eigenvectors of a real symmetric band matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 The eigenvector matrix X: Col 1 2 3 4 5 Row 1 -0.288675 -0.500000 -0.577350 -0.500000 0.288675 2 -0.500000 -0.500000 -0.414673E-16 0.500000 -0.500000 3 -0.577350 0.00000 0.577350 0.222045E-15 0.577350 4 -0.500000 0.500000 -0.553180E-16 -0.500000 -0.500000 5 -0.288675 0.500000 -0.577350 0.500000 0.288675 The residual (A-Lambda*I)*X: Col 1 2 3 4 5 Row 1 0.416334E-16 0.00000 -0.222045E-15 0.00000 -0.222045E-15 2 0.113798E-14 0.00000 0.829346E-16 -0.222045E-15 0.444089E-15 3 0.277556E-16 -0.111022E-15 0.222045E-15 -0.555112E-16 0.00000 4 0.105471E-14 0.00000 -0.111409E-15 -0.444089E-15 0.444089E-15 5 0.277556E-15 0.00000 -0.222045E-15 0.00000 0.222045E-15 RSG_TEST: RSG for real symmetric generalized problem. Find scalars LAMBDA and vectors X so that A*X = LAMBDA * B * X Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 0.00000 1.00000 2.00000 3.00000 2 1.00000 0.00000 1.00000 2.00000 3 2.00000 1.00000 0.00000 1.00000 4 3.00000 2.00000 1.00000 0.00000 The matrix B: Col 1 2 3 4 Row 1 2.00000 -1.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 4 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 -2.4357817 2 -0.52079729 3 -0.16421833 4 11.520797 RSG_TEST: RSG for real symmetric generalized problem. Find scalars LAMBDA and vectors X so that A*X = LAMBDA * B * X Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 0.00000 1.00000 2.00000 3.00000 2 1.00000 0.00000 1.00000 2.00000 3 2.00000 1.00000 0.00000 1.00000 4 3.00000 2.00000 1.00000 0.00000 The matrix B: Col 1 2 3 4 Row 1 2.00000 -1.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 4 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 -2.4357817 2 -0.52079729 3 -0.16421833 4 11.520797 The eigenvector matrix X: Col 1 2 3 4 Row 1 0.526940 0.251292 0.149448 0.660948 2 0.287038 -0.409656 -0.342942 0.912240 3 -0.287038 -0.409656 0.342942 0.912240 4 -0.526940 0.251292 -0.149448 0.660948 Residuals (A*x-(w*I)*B*x) for eigenvalue 1 -0.222045E-15 0.111022E-15 0.133227E-14 -0.133227E-14 Residuals (A*x-(w*I)*B*x) for eigenvalue 2 0.222045E-15 -0.721645E-15 -0.693889E-15 0.183187E-14 Residuals (A*x-(w*I)*B*x) for eigenvalue 3 -0.929812E-15 0.256739E-15 -0.534295E-15 -0.107553E-14 Residuals (A*x-(w*I)*B*x) for eigenvalue 4 0.355271E-14 -0.177636E-14 -0.177636E-14 0.355271E-14 RSGAB_TEST: RSGAB for real symmetric generalized problem. Find scalars LAMBDA and vectors X so that A*B*X = LAMBDA * X Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 0.00000 1.00000 2.00000 3.00000 2 1.00000 0.00000 1.00000 2.00000 3 2.00000 1.00000 0.00000 1.00000 4 3.00000 2.00000 1.00000 0.00000 The matrix B: Col 1 2 3 4 Row 1 2.00000 -1.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 4 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 -5.0000000 2 -2.0000000 3 -2.0000000 4 3.0000000 RSGAB_TEST: RSGAB for real symmetric generalized problem. Find scalars LAMBDA and vectors X so that A*B*X = LAMBDA * X Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 0.00000 1.00000 2.00000 3.00000 2 1.00000 0.00000 1.00000 2.00000 3 2.00000 1.00000 0.00000 1.00000 4 3.00000 2.00000 1.00000 0.00000 The matrix B: Col 1 2 3 4 Row 1 2.00000 -1.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 4 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 -5.0000000 2 -2.0000000 3 -2.0000000 4 3.0000000 The eigenvector matrix X: Col 1 2 3 4 Row 1 0.547723 -0.171729E-16 -0.314018E-15 0.707107 2 0.182574 0.325082E-01 -0.815849 0.707107 3 -0.182574 -0.690292 -0.436078 0.707107 4 -0.547723 -0.138588E-15 -0.644493E-16 0.707107 The residual matrix (A*B-Lambda*I)*X: Col 1 2 3 4 Row 1 0.355271E-14 -0.343458E-16 -0.405992E-15 -0.444089E-15 2 0.199840E-14 -0.222045E-15 -0.222045E-15 -0.444089E-15 3 0.122125E-14 -0.111022E-14 -0.888178E-15 0.00000 4 -0.444089E-15 -0.551323E-16 -0.128769E-14 0.444089E-15 RSGBA_TEST: RSGBA for real symmetric generalized problem. Find scalars LAMBDA and vectors X so that B*A*X = LAMBDA * X Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 0.00000 1.00000 2.00000 3.00000 2 1.00000 0.00000 1.00000 2.00000 3 2.00000 1.00000 0.00000 1.00000 4 3.00000 2.00000 1.00000 0.00000 The matrix B: Col 1 2 3 4 Row 1 2.00000 -1.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 4 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 -5.0000000 2 -2.0000000 3 -2.0000000 4 3.0000000 RSGBA_TEST: RSGBA for real symmetric generalized problem. Find scalars LAMBDA and vectors X so that B*A*X = LAMBDA * X Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 0.00000 1.00000 2.00000 3.00000 2 1.00000 0.00000 1.00000 2.00000 3 2.00000 1.00000 0.00000 1.00000 4 3.00000 2.00000 1.00000 0.00000 The matrix B: Col 1 2 3 4 Row 1 2.00000 -1.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 4 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 -5.0000000 2 -2.0000000 3 -2.0000000 4 3.0000000 The eigenvector matrix X: Col 1 2 3 4 Row 1 0.912871 -0.325082E-01 0.815849 0.707107 2 -0.222045E-15 0.755308 -1.19562 -0.166533E-15 3 0.388578E-15 -1.41309 -0.563058E-01 -0.277556E-16 4 -0.912871 0.690292 0.436078 0.707107 The residual matrix (B*A-Lambda*I)*X: Col 1 2 3 4 Row 1 0.444089E-14 -0.291434E-15 0.222045E-15 -0.177636E-14 2 -0.444089E-15 0.444089E-15 0.444089E-15 0.183187E-14 3 0.105471E-14 -0.133227E-14 -0.166533E-15 -0.360822E-15 4 -0.266454E-14 0.444089E-15 -0.222045E-15 0.00000 RSM_TEST RSM computes some eigenvalues and eigenvectors of a real symmetric matrix. Matrix order = 4 Number of eigenvectors desired = 4 The matrix A: Col 1 2 3 4 Row 1 5.00000 4.00000 1.00000 1.00000 2 4.00000 5.00000 1.00000 1.00000 3 1.00000 1.00000 4.00000 2.00000 4 1.00000 1.00000 2.00000 4.00000 The eigenvalues Lambda: 1 1.0000000 2 2.0000000 3 5.0000000 4 10.000000 The eigenvector matrix X: Col 1 2 3 4 Row 1 0.707107 -0.222045E-15 0.316228 0.632456 2 -0.707107 -0.555112E-16 0.316228 0.632456 3 0.00000 0.707107 -0.632456 0.316228 4 0.00000 -0.707107 -0.632456 0.316228 The residual (A-Lambda*I)*X: Col 1 2 3 4 Row 1 0.00000 -0.333067E-15 -0.266454E-14 -0.266454E-14 2 0.00000 -0.444089E-15 -0.155431E-14 -0.266454E-14 3 0.00000 0.888178E-15 -0.310862E-14 0.222045E-14 4 0.00000 0.666134E-15 0.888178E-15 0.177636E-14 RSP_TEST RSP computes the eigenvalues and eigenvectors of a real symmetric packed matrix. Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 5.00000 4.00000 1.00000 1.00000 2 4.00000 5.00000 1.00000 1.00000 3 1.00000 1.00000 4.00000 2.00000 4 1.00000 1.00000 2.00000 4.00000 The eigenvalues Lambda: 1 1.0000000 2 2.0000000 3 5.0000000 4 10.000000 RSP_TEST RSP computes the eigenvalues and eigenvectors of a real symmetric packed matrix. Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 5.00000 4.00000 1.00000 1.00000 2 4.00000 5.00000 1.00000 1.00000 3 1.00000 1.00000 4.00000 2.00000 4 1.00000 1.00000 2.00000 4.00000 The eigenvalues Lambda: 1 1.0000000 2 2.0000000 3 5.0000000 4 10.000000 The eigenvector matrix X: Col 1 2 3 4 Row 1 0.707107 0.555112E-16 0.316228 0.632456 2 -0.707107 0.111022E-15 0.316228 0.632456 3 0.00000 0.707107 -0.632456 0.316228 4 0.00000 -0.707107 -0.632456 0.316228 The residual matrix (A-Lambda*I)*X: Col 1 2 3 4 Row 1 0.00000 0.777156E-15 -0.177636E-14 -0.888178E-15 2 0.00000 0.777156E-15 -0.666134E-15 0.00000 3 0.00000 0.133227E-14 -0.222045E-14 -0.888178E-15 4 0.00000 -0.222045E-15 -0.888178E-15 0.00000 RSPP_TEST RSPP finds some eigenvalues and eigenvectors of a real symmetric packed matrix. Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 5.00000 4.00000 1.00000 1.00000 2 4.00000 5.00000 1.00000 1.00000 3 1.00000 1.00000 4.00000 2.00000 4 1.00000 1.00000 2.00000 4.00000 The eigenvalues Lambda: 1 1.0000000 2 2.0000000 3 5.0000000 4 10.000000 RSPP_TEST RSPP finds some eigenvalues and eigenvectors of a real symmetric packed matrix. Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 5.00000 4.00000 1.00000 1.00000 2 4.00000 5.00000 1.00000 1.00000 3 1.00000 1.00000 4.00000 2.00000 4 1.00000 1.00000 2.00000 4.00000 The eigenvalues Lambda: 1 1.0000000 2 2.0000000 3 5.0000000 4 10.000000 The eigenvector matrix X: Col 1 2 3 4 Row 1 0.707107 -0.555112E-16 0.316228 -0.632456 2 -0.707107 0.111022E-15 0.316228 -0.632456 3 0.00000 0.707107 -0.632456 -0.316228 4 0.00000 -0.707107 -0.632456 -0.316228 The residual matrix (A-Lambda*I)*X: Col 1 2 3 4 Row 1 0.00000 0.555112E-15 -0.333067E-14 0.888178E-15 2 0.00000 0.444089E-15 -0.222045E-14 0.00000 3 0.00000 0.133227E-14 -0.310862E-14 0.888178E-15 4 0.00000 0.00000 0.222045E-14 0.00000 RST_TEST RST computes the eigenvalues and eigenvectors of a real symmetric tridiagonal matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 RST_TEST RST computes the eigenvalues and eigenvectors of a real symmetric tridiagonal matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 The eigenvector matrix X: Col 1 2 3 4 5 Row 1 -0.288675 0.500000 -0.577350 -0.500000 -0.288675 2 -0.500000 0.500000 0.346252E-15 0.500000 0.500000 3 -0.577350 -0.226647E-15 0.577350 -0.489055E-15 -0.577350 4 -0.500000 -0.500000 -0.150268E-15 -0.500000 0.500000 5 -0.288675 -0.500000 -0.577350 0.500000 -0.288675 The residual matrix (A-Lambda*I)*X: Col 1 2 3 4 5 Row 1 -0.555112E-15 0.00000 -0.444089E-15 -0.222045E-15 0.00000 2 0.00000 -0.222045E-15 -0.359436E-15 0.00000 0.00000 3 -0.249800E-15 -0.509091E-16 -0.222045E-15 -0.871483E-16 0.00000 4 -0.832667E-16 0.166533E-15 -0.325319E-16 0.00000 0.00000 5 -0.180411E-15 0.444089E-15 0.222045E-15 -0.222045E-15 0.444089E-15 RT_TEST RT computes the eigenvalues and eigenvectors of a real sign-symmetric tridiagonal matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 RT_TEST RT computes the eigenvalues and eigenvectors of a real sign-symmetric tridiagonal matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 The eigenvector matrix X: Col 1 2 3 4 5 Row 1 -0.288675 0.500000 -0.577350 -0.500000 -0.288675 2 -0.500000 0.500000 0.346252E-15 0.500000 0.500000 3 -0.577350 -0.226647E-15 0.577350 -0.489055E-15 -0.577350 4 -0.500000 -0.500000 -0.150268E-15 -0.500000 0.500000 5 -0.288675 -0.500000 -0.577350 0.500000 -0.288675 The residual matrix (A-Lambda*I)*X: Col 1 2 3 4 5 Row 1 -0.555112E-15 0.00000 -0.444089E-15 -0.222045E-15 0.00000 2 0.00000 -0.222045E-15 -0.359436E-15 0.00000 0.00000 3 -0.249800E-15 -0.509091E-16 -0.222045E-15 -0.871483E-16 0.00000 4 -0.832667E-16 0.166533E-15 -0.325319E-16 0.00000 0.00000 5 -0.180411E-15 0.444089E-15 0.222045E-15 -0.222045E-15 0.444089E-15 SVD_TEST SVD computes the singular value decomposition of a real general matrix. Matrix order = 4 The matrix A: Col 1 2 3 4 Row 1 0.990000 0.200000E-02 0.600000E-02 0.200000E-02 2 0.200000E-02 0.990000 0.200000E-02 0.600000E-02 3 0.600000E-02 0.200000E-02 0.990000 0.200000E-02 4 0.200000E-02 0.600000E-02 0.200000E-02 0.990000 The singular values S 1 1.0000000 2 0.99200000 3 0.98400000 4 0.98400000 The U matrix: Col 1 2 3 4 Row 1 -0.500000 0.500000 0.706299 0.337928E-01 2 -0.500000 -0.500000 -0.337928E-01 0.706299 3 -0.500000 0.500000 -0.706299 -0.337928E-01 4 -0.500000 -0.500000 0.337928E-01 -0.706299 The V matrix: Col 1 2 3 4 Row 1 -0.500000 0.500000 0.706299 0.337928E-01 2 -0.500000 -0.500000 -0.337928E-01 0.706299 3 -0.500000 0.500000 -0.706299 -0.337928E-01 4 -0.500000 -0.500000 0.337928E-01 -0.706299 The product U * S * Transpose(V): Col 1 2 3 4 Row 1 0.990000 0.200000E-02 0.600000E-02 0.200000E-02 2 0.200000E-02 0.990000 0.200000E-02 0.600000E-02 3 0.600000E-02 0.200000E-02 0.990000 0.200000E-02 4 0.200000E-02 0.600000E-02 0.200000E-02 0.990000 TQL1_TEST TQL1 computes the eigenvalues of a real symmetric tridiagonal matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 TRIDIB_TEST TRIDIB computes some eigenvalues of a real symmetric tridiagonal matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 TSTURM_TEST TSTURM computes some eigenvalues and eigenvectors of a real symmetric tridiagonal matrix. Matrix order = 5 The matrix A: Col 1 2 3 4 5 Row 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 The eigenvalues Lambda: 1 0.26794919 2 1.0000000 3 2.0000000 4 3.0000000 5 3.7320508 The eigenvector matrix Z: Col 1 2 3 4 5 Row 1 0.288675 0.500000 -0.577350 0.500000 -0.288675 2 0.500000 0.500000 -0.102558E-14 -0.500000 0.500000 3 0.577350 -0.133227E-14 0.577350 -0.271424E-14 -0.577350 4 0.500000 -0.500000 -0.102558E-14 0.500000 0.500000 5 0.288675 -0.500000 -0.577350 -0.500000 -0.288675 The residual (A-Lambda*I)*X: Col 1 2 3 4 5 Row 1 0.130451E-14 0.555112E-16 0.155431E-14 0.199840E-14 0.444089E-15 2 0.471845E-15 0.832667E-15 0.160707E-14 0.532907E-14 0.666134E-15 3 -0.180411E-14 0.888178E-15 0.177636E-14 0.170342E-14 0.133227E-14 4 0.133227E-14 0.888178E-15 0.160707E-14 0.177636E-14 0.666134E-15 5 0.136002E-14 0.888178E-15 0.155431E-14 0.00000 0.444089E-15 EISPACK_TEST Normal end of execution. 22 January 2020 8:44:21.340 AM