lao2

求矩阵奇异分解svd算法

 ... INPUT MATRIX FROM ASCII FILE (not transposed)
 ... SINGULAR VECTORS STORED IN CORE
 ... LANCZOS VECTORS STORED IN CORE
 ... 
 ... SOLVE THE [A^TA]   EIGENPROBLEM
 ... NO. OF EQUATIONS          =   82
 ... MAX. NO. OF LANCZOS STEPS =   44
 ... MAX. NO. OF EIGENPAIRS    =   10
 ... LEFT  END OF THE INTERVAL =  0.00000000000000E+00
 ... RIGHT END OF THE INTERVAL =  1.00000000000000E-30
 ... WANT S-VECTORS?   [T/F]   =   T
 ... KAPPA                     =  1.00000000000000E-06
 Bellcore ADI Linguistics Data                                           belladit

           'belladit'
 ... NO. OF DOCUMENTS (ROWS)   =      374
 ... NO. OF TERMS     (COLS)   =       82
 ... ORDER OF MATRIX A         =       82
 ... 

 ...... ALLOCATED MEMORY (BYTES)=   1.04E+05
 ...... LANSO EXECUTION TIME= 2.592E-01
 ...... 
 ...... NUMBER OF LANCZOS STEPS =  44       NEIG =  10
 ...... 
 ......         COMPUTED RITZ VALUES  (ERROR BNDS)
 ...... 
 ......   1     3.55049373003382E+00  (   4.07E-07)
 ......   2     4.42887450180877E+00  (   1.72E-03)
 ......   3     5.12894426260373E+00  (   4.07E-03)
 ......   4     5.83039380671871E+00  (   5.46E-02)
 ......   5     6.99356947668029E+00  (   2.52E-01)
 ......   6     7.57252924696085E+00  (   1.40E+00)
 ......   7     8.34706138657497E+00  (   9.13E-01)
 ......   8     9.71733103701597E+00  (   4.65E-01)
 ......   9     1.09437372843391E+01  (   7.97E-01)
 ......  10     1.28429426794953E+01  (   1.39E+00)
 ......  11     1.43882899529109E+01  (   1.76E+00)
 ......  12     1.59406033877393E+01  (   2.28E+00)
 ......  13     1.76787598938949E+01  (   1.75E+00)
 ......  14     1.91448843228265E+01  (   4.06E+00)
 ......  15     2.08869330825166E+01  (   2.20E+00)
 ......  16     2.28846237250649E+01  (   2.98E+00)
 ......  17     2.48577102798676E+01  (   2.02E+00)
 ......  18     2.71953122304238E+01  (   2.50E+00)
 ......  19     2.87417530180587E+01  (   2.25E+00)
 ......  20     3.13642964784604E+01  (   1.94E+00)
 ......  21     3.29352918380765E+01  (   2.23E+00)
 ......  22     3.54845617295022E+01  (   2.49E+00)
 ......  23     3.66178781948609E+01  (   1.95E+00)
 ......  24     3.86858273595835E+01  (   2.14E+00)
 ......  25     4.01917894384419E+01  (   4.85E-01)
 ......  26     4.26648674046532E+01  (   4.09E-01)
 ......  27     4.42017220045778E+01  (   7.33E-02)
 ......  28     4.63925620005627E+01  (   1.30E-02)
 ......  29     4.81151289218662E+01  (   6.76E-04)
 ......  30     4.94033913980903E+01  (   1.44E-05)
 ......  31     5.27222276195975E+01  (   7.15E-07)
 ......  32     5.39926489888016E+01  (   2.90E-07)
 ......  33     5.47757915513079E+01  (   9.88E-08)
 ......  34     5.66699300278199E+01  (   2.70E-11)
 ......  35     6.20425301672993E+01  (   1.19E-15)
 ......  36     6.65518872251679E+01  (   6.72E-21)
 ......  37     6.98375125671880E+01  (   7.29E-24)
 ......  38     7.50088142326432E+01  (   7.11E-29)
 ......  39     7.93411862500450E+01  (   1.34E-30)
 ......  40     9.24374549973012E+01  (   1.31E-39)
 ......  41     1.11595482714292E+02  (   3.32E-54)
 ......  42     1.67949416448475E+02  (   1.80E-73)
 ......  43     2.02644907957103E+02  (   1.91E-81)
 ......  44     4.32089393862306E+02  (   0.00E+00)

 ...... LASVD EXECUTION TIME= 2.919E-01
 ...... 
 ......                      NSIG = 10
 ...... 
 ......         COMPUTED S-VALUES     (RES. NORMS)
 ...... 
 ......   1     7.87670807935011E+00  (   9.30E-09)
 ......   2     8.15793400470780E+00  (   2.16E-11)
 ......   3     8.35688414226187E+00  (   8.79E-13)
 ......   4     8.66076291285261E+00  (   3.13E-11)
 ......   5     8.90736696504896E+00  (   2.21E-11)
 ......   6     9.61443992114473E+00  (   7.71E-12)
 ......   7     1.05638763110088E+01  (   3.77E-12)
 ......   8     1.29595299470496E+01  (   1.39E-12)
 ......   9     1.42353401068293E+01  (   6.48E-13)
 ......  10     2.07867600616909E+01  (   9.75E-15)