rs.d

fortran版本的eslack,好不容易找到的。里面有程序的索引。主要包括了常用了矩阵计算代码。

1                    EXPLANATION OF COLUMN ENTRIES FOR THE SUMMARY STATISTICS

 -----------------------------------------------------------------------------------------------
 ORDER TQL2   TQLRAT IMTQL2 IMTQL1    LB      UB    M  IMTQLV   TSTURM   BISECT  M1 NO  TRIDIB  
 -----------------------------------------------------------------------------------------------

 UNDER 'ORDER' IS THE ORDER OF EACH TEST MATRIX.

 UNDER 'TQL2   TQLRAT' ARE THREE NUMBERS.  THE FIRST NUMBER, AN INTEGER, IS THE ABSOLUTE SUM OF
 THE ERROR FLAGS RETURNED SEPARATELY FROM  TQL2  AND  TQLRAT. THE SECOND NUMBER IS THE MEASURE
 OF PERFORMANCE BASED UPON THE RESIDUAL COMPUTED FOR THE  TQL2  PATH.  THE THIRD NUMBER
 MEASURES THE AGREEMENT OF THE EIGENVALUES FROM THE  TQL2  AND  TQLRAT  PATHS.

 UNDER 'IMTQL2 IMTQL1' ARE THREE NUMBERS WITH MEANING LIKE THOSE UNDER  'TQL2   TQLRAT'.       

 UNDER 'LB' AND 'UB' ARE THE INPUT VARIABLES SPECIFYING THE INTERVAL TO  BISECT  AND  TSTURM.  

 UNDER 'M' IS THE NUMBER OF EIGENVALUES DETERMINED BY  BISECT AND  TSTURM  THAT LIE IN THE
 INTERVAL (LB,UB).

 UNDER EACH OF 'IMTQLV', 'TSTURM', 'BISECT', AND 'TRIDIB' ARE TWO NUMBERS.  THE FIRST NUMBER,  
 AN INTEGER, IS THE ABSOLUTE SUM OF THE ERROR FLAGS RETURNED FROM THE RESPECTIVE PATH.         
 THE SECOND NUMBER IS THE MEASURE OF PERFORMANCE BASED UPON THE RESIDUAL COMPUTED FOR THE PATH.

 UNDER 'M1' AND 'NO' ARE THE VARIABLES SPECIFYING THE LOWER BOUNDARY INDEX AND THE NUMBER      
 OF EIGENVALUES TO  TRIDIB. 

 -1.0  AS THE MEASURE OF PERFORMANCE IS PRINTED IF AN ERROR IN THE CORRESPONDING PATH HAS
 PREVENTED THE COMPUTATION OF THE EIGENVECTORS.

 THE  TQL2    PATH USES THE EISPACK CODES   TRED2-TQL2  ,     
 AS CALLED FROM DRIVER SUBROUTINE  RS.
 THE  TQLRAT  PATH USES THE EISPACK CODES   TRED1-TQLRAT,     
 AS CALLED FROM DRIVER SUBROUTINE  RS.
 THE  IMTQL2  PATH USES THE EISPACK CODES   TRED2-IMTQL2.     
 THE  IMTQL1  PATH USES THE EISPACK CODES   TRED1-IMTQL1.     
 THE  IMTQLV  PATH USES THE EISPACK CODES   TRED1-IMTQLV-TINVIT-TRBAK ,
 AS CALLED FROM DRIVER SUBROUTINE RSM.
 THE  TSTURM  PATH USES THE EISPACK CODES   TRED1-TSTURM-TRBAK1.
 THE  BISECT  PATH USES THE EISPACK CODES   TRED1-BISECT-TINVIT-TRBAK1.
 THE  TRIDIB  PATH USES THE EISPACK CODES   TRED1-TRIDIB-TINVIT-TRBAK1.

 D.P. VERSION 04/15/83
1       TABULATION OF THE ERROR FLAG  ERROR  AND THE MEASURE OF PERFORMANCE  Y  FOR 
     THE  EISPACK  CODES.  THIS RUN DISPLAYS THESE STATISTICS FOR REAL SYMMETRIC MATRICES.    
0ORDER TQL2   TQLRAT IMTQL2 IMTQL1    LB      UB    M  IMTQLV   TSTURM   BISECT  M1 NO  TRIDIB 
   3  0  .014  .020  0  .026  .000 -1.e+30  1.e+30  3  0  .017  0  .030  0  .030  1  3  0  .030
   8  0  .025  .019  0  .025  .000 -1.e+30  1.e+30  8  0  .073  0  .032  0  .032  1  8  0  .032
  15  0  .039  .048  0  .025  .000 -1.e+30  1.e+30 15  0  .095  0  .023  0  .023  1 15  0  .023
  20  0  .038  .027  0  .028  .000 -1.e+30  1.e+30 20  0  .336  0  .133  0  .134  1 20  0  .134
   4  0  .018  .029  0  .018  .000 -1.e+30  1.e+30  4  0  .041  0  .041  0  .041  1  4  0  .041
  12  0  .021  .020  0  .018  .000 -1.e+30  1.e+30 12  0  .086  0  .032  0  .032  1 12  0  .032
   5  0  .024  .026  0  .016  .000 -1.e+30  1.e+30  5  0  .028  0  .015  0  .015  1  5  0  .015
  10  0  .000  .000  0  .000  .000 -1.e+30  1.e+30 10  0  .000  0  .000  0  .000  1 10  0  .000
   6  0  .010  .011  0  .006  .000 -1.e+30  1.e+30  6  0  .031  0  .031  0  .031  1  6  0  .031
   5  0  .000  .000  0  .000  .000 -1.e+30  1.e+30  5  0  .000  0  .000  0  .000  1  5  0  .000
   3  0  .000  .000  0  .000  .000 -1.e+30  1.e+30  3  0  .000  0  .000  0  .000  1  3  0  .000
  20  0  .000  .000  0  .000  .000 -1.e+30  1.e+30 20  0  .000  0  .000  0  .000  1 20  0  .000
   3  0  .020  .044  0  .012  .000 -1.e+30  1.e+30  3  0  .063  0  .063  0  .063  1  3  0  .063
   5  0  .055  .034  0  .053  .000 -1.e+30  1.e+30  5  0  .052  0  .052  0  .052  1  5  0  .052
   8  0  .030  .011  0  .023  .000 -1.e+30  1.e+30  8  0  .074  0  .033  0  .033  1  8  0  .033
   3  0  .049  .058  0  .038  .000 -1.e+30  1.e+30  3  0  .106  0  .106  0  .106  1  3  0  .106
   4  0  .010  .015  0  .035  .000 -1.e+30  1.e+30  4  0  .043  0  .043  0  .043  1  4  0  .043
   6  0  .030  .011  0  .035  .000 -1.e+30  1.e+30  6  0  .053  0  .026  0  .026  1  6  0  .026
   8  0  .016  .001  0  .016  .000 -1.e+30  1.e+30  8  0  .007  0  .022  0  .022  1  8  0  .022
   4  0  .041  .017  0  .059  .000 -1.e+30  1.e+30  4  0  .112  0  .046  0  .046  1  4  0  .046
   5  0  .038  .029  0  .046  .000 -1.e+30  1.e+30  5  0  .045  0  .070  0  .070  1  5  0  .070
0END OF DATA FOR SUBROUTINE  RMATIN(RGREADI). 
1