20130111 081751.627 TOMS358_PRB: FORTRAN77 version Test TOMS358 library. TEST01 Call ACM TOMS Algorithm 358 for the singular value decomposition: A = U S V* of an M by N complex matrix. Matrix row order M = 5 Matrix column order N = 2 Matrix A: 0.4499 -0.1267 -0.2361 0.0775 -0.8432 -0.3443 0.0186 -0.6332 0.5896 0.2601 0.8928 0.0103 0.3911 0.3234 -0.5605 0.7638 -0.1395 -0.1561 0.3064 0.0263 Singular values: 1.75208 0.973251 U: -0.1266 -0.0614 -0.4090 0.2535 0.0002 0.7742 -0.2394 -0.1659 0.2383 -0.0877 0.5252 0.2924 0.3566 0.0346 0.0885 0.4135 0.4869 -0.2559 -0.0419 0.1650 -0.4433 0.2315 -0.1410 -0.6706 0.3780 0.2143 0.2607 0.1109 -0.0243 0.1032 -0.2728 -0.5167 -0.1239 0.3124 0.1092 -0.0869 0.7119 -0.0823 0.0956 -0.0585 -0.0394 0.1605 0.2388 -0.0240 -0.0292 -0.1127 0.0318 0.1334 0.9393 -0.0205 V: -0.6206 0.0000 -0.7841 0.0000 -0.4528 0.6402 0.3583 -0.5067 Matrix U S V* (should equal the original A): 0.4499 -0.1267 -0.2361 0.0775 -0.8432 -0.3443 0.0186 -0.6332 0.5896 0.2601 0.8929 0.0103 0.3911 0.3234 -0.5605 0.7638 -0.1395 -0.1561 0.3064 0.0263 TEST02 CSVDC is the LINPACK routine for the singular value decomposition: A = U S V* of an M by N complex matrix. Matrix row order M = 5 Matrix column order N = 2 Matrix A: 0.4499 -0.1267 -0.2361 0.0775 -0.8432 -0.3443 0.0186 -0.6332 0.5896 0.2601 0.8928 0.0103 0.3911 0.3234 -0.5605 0.7638 -0.1395 -0.1561 0.3064 0.0263 Singular values: 1.75208 0.00000 0.973251 0.00000 U: 0.1266 0.0614 -0.4090 0.2535 0.0002 0.7742 -0.2394 -0.1659 0.2383 -0.0877 -0.5252 -0.2924 0.3566 0.0346 0.0885 0.4135 0.4869 -0.2559 -0.0419 0.1650 0.4433 -0.2315 -0.1410 -0.6706 0.3780 0.2143 0.2607 0.1109 -0.0243 0.1032 0.2728 0.5167 -0.1239 0.3124 0.1092 -0.0869 0.7119 -0.0823 0.0956 -0.0585 0.0394 -0.1605 0.2388 -0.0240 -0.0292 -0.1127 0.0318 0.1334 0.9393 -0.0205 V: 0.6206 0.0000 -0.7841 0.0000 0.4528 -0.6402 0.3583 -0.5067 Matrix U S V* (should equal the original A): 0.4499 -0.1267 -0.2361 0.0775 -0.8432 -0.3443 0.0186 -0.6332 0.5896 0.2601 0.8928 0.0103 0.3911 0.3234 -0.5605 0.7638 -0.1395 -0.1561 0.3064 0.0263 TEST01 Call ACM TOMS Algorithm 358 for the singular value decomposition: A = U S V* of an M by N complex matrix. Matrix row order M = 6 Matrix column order N = 4 Matrix A: 0.4499 -0.1267 0.0186 -0.6332 0.4350 -0.2666 0.0314 -0.4336 -0.8432 -0.3443 0.8928 0.0103 -0.2009 0.2707 0.2491 0.5787 0.5896 0.2601 -0.5605 0.7638 -0.0975 0.9019 -0.5084 0.5621 0.3911 0.3234 0.3064 0.0263 -0.7702 -0.3143 -0.5066 0.6005 -0.1395 -0.1561 0.5008 -0.7799 -0.8892 0.2657 -0.1022 -0.4500 -0.2361 0.0775 0.3505 0.0166 -0.7799 -0.5512 -0.1046 0.3267 Singular values: 2.18690 2.02082 1.17274 0.726572 U: -0.0874 0.0785 0.2113 -0.2947 0.4747 -0.1990 -0.4352 0.0507 -0.4813 -0.3792 0.1438 -0.0394 -0.3294 -0.2841 0.0462 0.4345 -0.5446 -0.0702 -0.2664 -0.0385 -0.3957 -0.2176 -0.1966 0.0580 0.6394 0.1111 -0.1603 0.3537 0.0598 0.1787 -0.0076 -0.2420 -0.4508 0.0966 0.2205 0.2695 -0.0528 0.0133 -0.5044 0.2021 0.3239 0.3169 -0.5142 0.1930 0.1200 0.1441 -0.3858 -0.1041 -0.4521 0.2903 -0.1047 0.2017 0.1870 -0.3645 -0.0957 -0.5333 0.0578 0.3710 0.1363 0.1966 -0.2807 -0.0737 -0.4121 0.1000 -0.0498 0.1494 0.0217 0.2851 0.0587 -0.1685 0.7729 0.0147 V: 0.3631 0.0000 -0.0507 0.0000 0.9271 0.0000 -0.0783 0.0000 -0.5442 -0.4928 0.0469 0.2852 0.1888 0.1697 -0.3180 -0.4602 0.3399 -0.3167 0.6687 -0.4119 -0.1282 0.1010 -0.3750 -0.0063 -0.2699 -0.2013 0.5198 0.1639 0.1877 0.1192 0.6339 0.3720 Matrix U S V* (should equal the original A): 0.4499 -0.1267 0.0186 -0.6332 0.4350 -0.2666 0.0314 -0.4336 -0.8432 -0.3443 0.8928 0.0103 -0.2009 0.2707 0.2491 0.5787 0.5896 0.2601 -0.5605 0.7638 -0.0975 0.9019 -0.5084 0.5621 0.3911 0.3234 0.3064 0.0263 -0.7702 -0.3143 -0.5066 0.6005 -0.1395 -0.1561 0.5008 -0.7799 -0.8892 0.2657 -0.1022 -0.4500 -0.2361 0.0775 0.3505 0.0166 -0.7799 -0.5512 -0.1046 0.3267 TEST02 CSVDC is the LINPACK routine for the singular value decomposition: A = U S V* of an M by N complex matrix. Matrix row order M = 6 Matrix column order N = 4 Matrix A: 0.4499 -0.1267 0.0186 -0.6332 0.4350 -0.2666 0.0314 -0.4336 -0.8432 -0.3443 0.8928 0.0103 -0.2009 0.2707 0.2491 0.5787 0.5896 0.2601 -0.5605 0.7638 -0.0975 0.9019 -0.5084 0.5621 0.3911 0.3234 0.3064 0.0263 -0.7702 -0.3143 -0.5066 0.6005 -0.1395 -0.1561 0.5008 -0.7799 -0.8892 0.2657 -0.1022 -0.4500 -0.2361 0.0775 0.3505 0.0166 -0.7799 -0.5512 -0.1046 0.3267 Singular values: 2.18690 0.00000 2.02082 0.00000 1.17274 -0.00000 0.726572 -0.00000 U: 0.0874 -0.0785 0.2113 -0.2947 0.4747 -0.1990 0.4352 -0.0507 -0.4813 -0.3792 0.1438 -0.0394 0.3294 0.2841 0.0462 0.4345 -0.5446 -0.0702 0.2664 0.0385 -0.3957 -0.2176 -0.1966 0.0580 -0.6394 -0.1111 -0.1603 0.3537 0.0598 0.1787 0.0076 0.2420 -0.4508 0.0966 0.2205 0.2695 0.0528 -0.0133 -0.5044 0.2021 0.3239 0.3169 0.5142 -0.1930 0.1200 0.1441 -0.3858 -0.1041 0.4521 -0.2903 -0.1047 0.2017 0.1870 -0.3645 0.0957 0.5333 0.0578 0.3710 0.1363 0.1966 0.2807 0.0737 -0.4121 0.1000 -0.0498 0.1494 -0.0217 -0.2851 0.0587 -0.1685 0.7729 0.0147 V: -0.3631 0.0000 -0.0507 0.0000 0.9271 -0.0000 0.0783 -0.0000 0.5442 0.4928 0.0469 0.2852 0.1888 0.1697 0.3180 0.4602 -0.3399 0.3167 0.6687 -0.4119 -0.1282 0.1010 0.3750 0.0063 0.2699 0.2013 0.5198 0.1638 0.1877 0.1192 -0.6339 -0.3720 Matrix U S V* (should equal the original A): 0.4499 -0.1267 0.0186 -0.6332 0.4350 -0.2666 0.0314 -0.4336 -0.8432 -0.3443 0.8928 0.0103 -0.2009 0.2707 0.2491 0.5787 0.5896 0.2601 -0.5605 0.7638 -0.0975 0.9019 -0.5084 0.5621 0.3911 0.3234 0.3064 0.0263 -0.7702 -0.3143 -0.5066 0.6005 -0.1395 -0.1561 0.5008 -0.7799 -0.8892 0.2657 -0.1022 -0.4500 -0.2361 0.0775 0.3505 0.0166 -0.7799 -0.5512 -0.1046 0.3267 TEST01 Call ACM TOMS Algorithm 358 for the singular value decomposition: A = U S V* of an M by N complex matrix. Matrix row order M = 5 Matrix column order N = 5 Matrix A: 0.4499 -0.1267 -0.2361 0.0775 0.5008 -0.7799 -0.7702 -0.3143 -0.5084 0.5621 -0.8432 -0.3443 0.0186 -0.6332 0.3505 0.0166 -0.8892 0.2657 -0.5066 0.6005 0.5896 0.2601 0.8928 0.0103 0.4350 -0.2666 -0.7799 -0.5512 -0.1022 -0.4500 0.3911 0.3234 -0.5605 0.7638 -0.2009 0.2707 0.0314 -0.4336 -0.1046 0.3267 -0.1395 -0.1561 0.3064 0.0263 -0.0975 0.9019 0.2491 0.5787 0.4106 -0.8102 Singular values: 2.43285 1.94345 1.15708 0.705795 0.248448 U: 0.2242 -0.4836 -0.2151 -0.3438 0.0537 0.0352 -0.3236 0.2017 0.3038 0.5564 -0.0470 -0.5195 0.3621 0.0569 -0.7086 -0.0724 0.0222 -0.2600 0.0696 -0.1153 0.5033 -0.1529 -0.2687 0.2132 0.1156 0.5523 0.0944 -0.3937 0.2180 -0.2752 0.0390 0.1680 -0.3397 -0.3904 -0.1128 -0.3801 -0.4899 -0.2989 0.1194 -0.4491 -0.0944 0.3583 0.1730 0.5346 -0.1006 0.0531 -0.5220 -0.1440 0.4013 0.2825 V: 0.2731 0.0000 -0.4363 0.0000 0.5986 0.0000 -0.4818 0.0000 -0.3803 0.0000 0.3175 -0.0395 -0.1459 0.5963 -0.1079 0.3004 0.3116 0.2061 -0.1692 -0.5007 0.4497 -0.0720 0.2791 -0.1763 -0.3194 0.4403 -0.3277 0.4089 -0.0848 0.3256 -0.1286 0.6084 0.2933 -0.2009 0.2858 -0.0985 0.2001 0.5267 -0.2325 -0.1550 -0.3842 0.2852 -0.3889 0.2152 -0.1289 0.3687 0.1898 -0.0146 -0.2733 0.5568 Matrix U S V* (should equal the original A): 0.4499 -0.1267 -0.2361 0.0775 0.5008 -0.7799 -0.7702 -0.3143 -0.5084 0.5621 -0.8432 -0.3443 0.0186 -0.6332 0.3505 0.0166 -0.8892 0.2657 -0.5066 0.6005 0.5896 0.2601 0.8928 0.0103 0.4350 -0.2666 -0.7799 -0.5512 -0.1022 -0.4500 0.3911 0.3234 -0.5605 0.7638 -0.2009 0.2707 0.0314 -0.4336 -0.1046 0.3267 -0.1395 -0.1561 0.3064 0.0263 -0.0975 0.9019 0.2491 0.5787 0.4106 -0.8102 TEST02 CSVDC is the LINPACK routine for the singular value decomposition: A = U S V* of an M by N complex matrix. Matrix row order M = 5 Matrix column order N = 5 Matrix A: 0.4499 -0.1267 -0.2361 0.0775 0.5008 -0.7799 -0.7702 -0.3143 -0.5084 0.5621 -0.8432 -0.3443 0.0186 -0.6332 0.3505 0.0166 -0.8892 0.2657 -0.5066 0.6005 0.5896 0.2601 0.8928 0.0103 0.4350 -0.2666 -0.7799 -0.5512 -0.1022 -0.4500 0.3911 0.3234 -0.5605 0.7638 -0.2009 0.2707 0.0314 -0.4336 -0.1046 0.3267 -0.1395 -0.1561 0.3064 0.0263 -0.0975 0.9019 0.2491 0.5787 0.4106 -0.8102 Singular values: 2.43285 0.00000 1.94346 0.00000 1.15708 0.00000 0.705795 0.00000 0.248448 0.00000 U: 0.2242 -0.4836 -0.2151 -0.3438 0.0537 0.0352 -0.3236 0.2017 -0.3038 -0.5564 -0.0470 -0.5195 0.3621 0.0569 -0.7086 -0.0724 0.0222 -0.2600 -0.0696 0.1153 0.5033 -0.1529 -0.2687 0.2132 0.1156 0.5523 0.0944 -0.3937 -0.2180 0.2752 0.0390 0.1680 -0.3397 -0.3904 -0.1128 -0.3801 -0.4899 -0.2989 -0.1194 0.4491 -0.0944 0.3583 0.1730 0.5346 -0.1006 0.0531 -0.5220 -0.1440 -0.4013 -0.2825 V: 0.2731 0.0000 -0.4363 0.0000 0.5986 0.0000 -0.4818 0.0000 0.3803 0.0000 0.3175 -0.0395 -0.1459 0.5963 -0.1079 0.3004 0.3116 0.2061 0.1692 0.5007 0.4497 -0.0720 0.2791 -0.1763 -0.3194 0.4403 -0.3277 0.4089 0.0848 -0.3256 -0.1286 0.6084 0.2933 -0.2009 0.2858 -0.0985 0.2001 0.5267 0.2325 0.1550 -0.3842 0.2852 -0.3889 0.2152 -0.1289 0.3687 0.1898 -0.0146 0.2733 -0.5568 Matrix U S V* (should equal the original A): 0.4499 -0.1267 -0.2361 0.0775 0.5008 -0.7799 -0.7702 -0.3143 -0.5084 0.5621 -0.8432 -0.3443 0.0186 -0.6332 0.3505 0.0166 -0.8892 0.2657 -0.5066 0.6005 0.5896 0.2601 0.8928 0.0103 0.4350 -0.2666 -0.7799 -0.5512 -0.1022 -0.4500 0.3911 0.3234 -0.5605 0.7638 -0.2009 0.2707 0.0314 -0.4336 -0.1046 0.3267 -0.1395 -0.1561 0.3064 0.0263 -0.0975 0.9019 0.2491 0.5787 0.4106 -0.8102 TEST02 CSVDC is the LINPACK routine for the singular value decomposition: A = U S V* of an M by N complex matrix. Matrix row order M = 3 Matrix column order N = 5 Matrix A: 0.4499 -0.1267 0.3911 0.3234 0.0186 -0.6332 0.3064 0.0263 0.4350 -0.2666 -0.8432 -0.3443 -0.1395 -0.1561 0.8928 0.0103 0.5008 -0.7799 -0.2009 0.2707 0.5896 0.2601 -0.2361 0.0775 -0.5605 0.7638 0.3505 0.0166 -0.0975 0.9019 Singular values: 1.86953 0.00000 1.44684 -0.00000 0.779459 0.00000 U: -0.0689 0.0675 -0.5527 0.3111 0.7200 0.2645 -0.5653 -0.4990 0.4941 -0.1448 0.3823 0.1420 0.6343 -0.1404 0.1078 -0.5666 0.4646 -0.1717 V: 0.5062 0.0000 -0.5105 0.0000 0.1904 0.0000 0.6685 0.0000 0.0000 0.0000 -0.0048 0.0075 -0.1598 0.3616 0.2165 -0.1090 -0.1801 0.3015 -0.5443 0.6044 -0.5438 -0.4749 -0.1803 -0.1682 -0.2602 0.4171 0.3483 0.1124 -0.1025 0.1783 0.1641 -0.3894 0.1573 0.1536 0.6006 0.4664 -0.1752 0.2794 0.2629 -0.1300 -0.1380 -0.1573 -0.6797 -0.1097 0.0054 -0.2915 -0.4161 0.1184 0.4546 0.0586 Matrix U S V* (should equal the original A): 0.4499 -0.1267 0.3911 0.3234 0.0186 -0.6332 0.3064 0.0263 0.4350 -0.2666 -0.8432 -0.3443 -0.1395 -0.1561 0.8928 0.0103 0.5008 -0.7799 -0.2009 0.2707 0.5896 0.2601 -0.2361 0.0775 -0.5605 0.7638 0.3505 0.0166 -0.0975 0.9019 TOMS358_PRB: Normal end of execution. 20130111 081751.648