rsg.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 SYSTEM.

 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 OR VALUES.

 THE  TQL2    PATH USES THE EISPACK CODES   REDUC-TRED2-TQL2  -REBAK , 
 AS CALLED FROM DRIVER SUBROUTINE  RSG.
 THE  TQLRAT  PATH USES THE EISPACK CODES   REDUC-TRED1-TQLRAT, 
 AS CALLED FROM DRIVER SUBROUTINE  RSG.
 THE  IMTQL2  PATH USES THE EISPACK CODES   REDUC-TRED2-IMTQL2-REBAK . 
 THE  IMTQL1  PATH USES THE EISPACK CODES   REDUC-TRED1-IMTQL1. 
 THE  IMTQLV  PATH USES THE EISPACK CODES   REDUC-TRED1-IMTQLV-TINVIT-TRBAK1-REBAK.  
 THE  TSTURM  PATH USES THE EISPACK CODES   REDUC-TRED1-TSTURM-TRBAK1-REBAK . 
 THE  BISECT  PATH USES THE EISPACK CODES   REDUC-TRED1-BISECT-TINVIT-TRBAK1-REBAK.  
 THE  TRIDIB  PATH USES THE EISPACK CODES   REDUC-TRED1-TRIDIB-TINVIT-TRBAK1-REBAK.  
 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 THE REAL SYMMETRIC
     GENERALIZED EIGENPROBLEM  A*X = (LAMBDA)*B*X .
0ORDER TQL2   TQLRAT IMTQL2 IMTQL1    LB      UB    M  IMTQLV   TSTURM   BISECT  M1 NO  TRIDIB 
   3  0  .069  .042  0  .069  .000 -1.e+30  1.e+30  3  0  .052  0  .052  0  .052  1  3  0  .052
   8  0  .126  .003  0  .161  .000 -1.e+30  1.e+30  8  0  .161  0  .228  0  .228  1  8  0  .228
  15  0  .050  .024  0  .040  .000 -1.e+30  1.e+30 15  0  .125  0  .104  0  .104  1 15  0  .104
  20  0  .031  .030  0  .036  .000 -1.e+30  1.e+30 20  0  .132  0  .045  0  .045  1 20  0  .045
   4  0  .063  .014  0  .043  .000 -1.e+30  1.e+30  4  0  .043  0  .088  0  .088  1  4  0  .088
  12  0  .035  .010  0  .036  .000 -1.e+30  1.e+30 12  0  .108  0  .108  0  .108  1 12  0  .108
   5  0  .019  .015  0  .025  .000 -1.e+30  1.e+30  5  0  .052  0  .033  0  .033  1  5  0  .033
  10  0  .068  .023  0  .062  .000 -1.e+30  1.e+30 10  0  .165  0  .110  0  .110  1 10  0  .110
   6  0  .053  .014  0  .050  .000 -1.e+30  1.e+30  6  0  .055  0  .070  0  .070  1  6  0  .070
   5  0  .065  .014  0  .036  .000 -1.e+30  1.e+30  5  0  .057  0  .111  0  .111  1  5  0  .111
   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  .050  .006  0  .040  .000 -1.e+30  1.e+30 20  0  .160  0  .117  0  .117  1 20  0  .117
   3  0  .149  .056  0  .092  .000 -1.e+30  1.e+30  3  0  .113  0  .141  0  .141  1  3  0  .141
   5  0  .024  .034  0  .029  .000 -1.e+30  1.e+30  5  0  .026  0  .026  0  .026  1  5  0  .026
   8  0  .034  .003  0  .049  .000 -1.e+30  1.e+30  8  0  .049  0  .073  0  .073  1  8  0  .073
   3  0  .051  .084  0  .049  .000 -1.e+30  1.e+30  3  0  .077  0  .043  0  .043  1  3  0  .043
   4  0  .020  .001  0  .018  .000 -1.e+30  1.e+30  4  0  .016  0  .016  0  .016  1  4  0  .016
   6  0  .149  .065  0  .066  .000 -1.e+30  1.e+30  6  0  .111  0  .055  0  .055  1  6  0  .055
   8  0  .015  .001  0  .015  .000 -1.e+30  1.e+30  8  0  .011  0  .017  0  .017  1  8  0  .017
   4  0  .051  .018  0  .060  .000 -1.e+30  1.e+30  4  0  .076  0  .129  0  .129  1  4  0  .129
   5  0  .031  .020  0  .041  .000 -1.e+30  1.e+30  5  0  .042  0  .042  0  .042  1  5  0  .042
0END OF DATA FOR SUBROUTINE  RMATIN(RSGREADI). 
1