![]() |
LAPACK
3.4.1
LAPACK: Linear Algebra PACKage
|
00001 *> \brief \b DCHKRFP 00002 * 00003 * =========== DOCUMENTATION =========== 00004 * 00005 * Online html documentation available at 00006 * http://www.netlib.org/lapack/explore-html/ 00007 * 00008 * Definition: 00009 * =========== 00010 * 00011 * PROGRAM DCHKRFP 00012 * 00013 * 00014 *> \par Purpose: 00015 * ============= 00016 *> 00017 *> \verbatim 00018 *> 00019 *> DCHKRFP is the main test program for the DOUBLE PRECISION linear 00020 *> equation routines with RFP storage format 00021 *> 00022 *> \endverbatim 00023 * 00024 * Arguments: 00025 * ========== 00026 * 00027 *> \verbatim 00028 *> MAXIN INTEGER 00029 *> The number of different values that can be used for each of 00030 *> M, N, or NB 00031 *> 00032 *> MAXRHS INTEGER 00033 *> The maximum number of right hand sides 00034 *> 00035 *> NTYPES INTEGER 00036 *> 00037 *> NMAX INTEGER 00038 *> The maximum allowable value for N. 00039 *> 00040 *> NIN INTEGER 00041 *> The unit number for input 00042 *> 00043 *> NOUT INTEGER 00044 *> The unit number for output 00045 *> \endverbatim 00046 * 00047 * Authors: 00048 * ======== 00049 * 00050 *> \author Univ. of Tennessee 00051 *> \author Univ. of California Berkeley 00052 *> \author Univ. of Colorado Denver 00053 *> \author NAG Ltd. 00054 * 00055 *> \date April 2012 00056 * 00057 *> \ingroup double_lin 00058 * 00059 * ===================================================================== 00060 PROGRAM DCHKRFP 00061 * 00062 * -- LAPACK test routine (version 3.4.1) -- 00063 * -- LAPACK is a software package provided by Univ. of Tennessee, -- 00064 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- 00065 * April 2012 00066 * 00067 * ===================================================================== 00068 * 00069 * .. Parameters .. 00070 INTEGER MAXIN 00071 PARAMETER ( MAXIN = 12 ) 00072 INTEGER NMAX 00073 PARAMETER ( NMAX = 50 ) 00074 INTEGER MAXRHS 00075 PARAMETER ( MAXRHS = 16 ) 00076 INTEGER NTYPES 00077 PARAMETER ( NTYPES = 9 ) 00078 INTEGER NIN, NOUT 00079 PARAMETER ( NIN = 5, NOUT = 6 ) 00080 * .. 00081 * .. Local Scalars .. 00082 LOGICAL FATAL, TSTERR 00083 INTEGER VERS_MAJOR, VERS_MINOR, VERS_PATCH 00084 INTEGER I, NN, NNS, NNT 00085 DOUBLE PRECISION EPS, S1, S2, THRESH 00086 00087 * .. 00088 * .. Local Arrays .. 00089 INTEGER NVAL( MAXIN ), NSVAL( MAXIN ), NTVAL( NTYPES ) 00090 DOUBLE PRECISION WORKA( NMAX, NMAX ) 00091 DOUBLE PRECISION WORKASAV( NMAX, NMAX ) 00092 DOUBLE PRECISION WORKB( NMAX, MAXRHS ) 00093 DOUBLE PRECISION WORKXACT( NMAX, MAXRHS ) 00094 DOUBLE PRECISION WORKBSAV( NMAX, MAXRHS ) 00095 DOUBLE PRECISION WORKX( NMAX, MAXRHS ) 00096 DOUBLE PRECISION WORKAFAC( NMAX, NMAX ) 00097 DOUBLE PRECISION WORKAINV( NMAX, NMAX ) 00098 DOUBLE PRECISION WORKARF( (NMAX*(NMAX+1))/2 ) 00099 DOUBLE PRECISION WORKAP( (NMAX*(NMAX+1))/2 ) 00100 DOUBLE PRECISION WORKARFINV( (NMAX*(NMAX+1))/2 ) 00101 DOUBLE PRECISION D_WORK_DLATMS( 3 * NMAX ) 00102 DOUBLE PRECISION D_WORK_DPOT01( NMAX ) 00103 DOUBLE PRECISION D_TEMP_DPOT02( NMAX, MAXRHS ) 00104 DOUBLE PRECISION D_TEMP_DPOT03( NMAX, NMAX ) 00105 DOUBLE PRECISION D_WORK_DLANSY( NMAX ) 00106 DOUBLE PRECISION D_WORK_DPOT02( NMAX ) 00107 DOUBLE PRECISION D_WORK_DPOT03( NMAX ) 00108 * .. 00109 * .. External Functions .. 00110 DOUBLE PRECISION DLAMCH, DSECND 00111 EXTERNAL DLAMCH, DSECND 00112 * .. 00113 * .. External Subroutines .. 00114 EXTERNAL ILAVER, DDRVRFP, DDRVRF1, DDRVRF2, DDRVRF3, 00115 + DDRVRF4 00116 * .. 00117 * .. Executable Statements .. 00118 * 00119 S1 = DSECND( ) 00120 FATAL = .FALSE. 00121 * 00122 * Read a dummy line. 00123 * 00124 READ( NIN, FMT = * ) 00125 * 00126 * Report LAPACK version tag (e.g. LAPACK-3.2.0) 00127 * 00128 CALL ILAVER( VERS_MAJOR, VERS_MINOR, VERS_PATCH ) 00129 WRITE( NOUT, FMT = 9994 ) VERS_MAJOR, VERS_MINOR, VERS_PATCH 00130 * 00131 * Read the values of N 00132 * 00133 READ( NIN, FMT = * )NN 00134 IF( NN.LT.1 ) THEN 00135 WRITE( NOUT, FMT = 9996 )' NN ', NN, 1 00136 NN = 0 00137 FATAL = .TRUE. 00138 ELSE IF( NN.GT.MAXIN ) THEN 00139 WRITE( NOUT, FMT = 9995 )' NN ', NN, MAXIN 00140 NN = 0 00141 FATAL = .TRUE. 00142 END IF 00143 READ( NIN, FMT = * )( NVAL( I ), I = 1, NN ) 00144 DO 10 I = 1, NN 00145 IF( NVAL( I ).LT.0 ) THEN 00146 WRITE( NOUT, FMT = 9996 )' M ', NVAL( I ), 0 00147 FATAL = .TRUE. 00148 ELSE IF( NVAL( I ).GT.NMAX ) THEN 00149 WRITE( NOUT, FMT = 9995 )' M ', NVAL( I ), NMAX 00150 FATAL = .TRUE. 00151 END IF 00152 10 CONTINUE 00153 IF( NN.GT.0 ) 00154 $ WRITE( NOUT, FMT = 9993 )'N ', ( NVAL( I ), I = 1, NN ) 00155 * 00156 * Read the values of NRHS 00157 * 00158 READ( NIN, FMT = * )NNS 00159 IF( NNS.LT.1 ) THEN 00160 WRITE( NOUT, FMT = 9996 )' NNS', NNS, 1 00161 NNS = 0 00162 FATAL = .TRUE. 00163 ELSE IF( NNS.GT.MAXIN ) THEN 00164 WRITE( NOUT, FMT = 9995 )' NNS', NNS, MAXIN 00165 NNS = 0 00166 FATAL = .TRUE. 00167 END IF 00168 READ( NIN, FMT = * )( NSVAL( I ), I = 1, NNS ) 00169 DO 30 I = 1, NNS 00170 IF( NSVAL( I ).LT.0 ) THEN 00171 WRITE( NOUT, FMT = 9996 )'NRHS', NSVAL( I ), 0 00172 FATAL = .TRUE. 00173 ELSE IF( NSVAL( I ).GT.MAXRHS ) THEN 00174 WRITE( NOUT, FMT = 9995 )'NRHS', NSVAL( I ), MAXRHS 00175 FATAL = .TRUE. 00176 END IF 00177 30 CONTINUE 00178 IF( NNS.GT.0 ) 00179 $ WRITE( NOUT, FMT = 9993 )'NRHS', ( NSVAL( I ), I = 1, NNS ) 00180 * 00181 * Read the matrix types 00182 * 00183 READ( NIN, FMT = * )NNT 00184 IF( NNT.LT.1 ) THEN 00185 WRITE( NOUT, FMT = 9996 )' NMA', NNT, 1 00186 NNT = 0 00187 FATAL = .TRUE. 00188 ELSE IF( NNT.GT.NTYPES ) THEN 00189 WRITE( NOUT, FMT = 9995 )' NMA', NNT, NTYPES 00190 NNT = 0 00191 FATAL = .TRUE. 00192 END IF 00193 READ( NIN, FMT = * )( NTVAL( I ), I = 1, NNT ) 00194 DO 320 I = 1, NNT 00195 IF( NTVAL( I ).LT.0 ) THEN 00196 WRITE( NOUT, FMT = 9996 )'TYPE', NTVAL( I ), 0 00197 FATAL = .TRUE. 00198 ELSE IF( NTVAL( I ).GT.NTYPES ) THEN 00199 WRITE( NOUT, FMT = 9995 )'TYPE', NTVAL( I ), NTYPES 00200 FATAL = .TRUE. 00201 END IF 00202 320 CONTINUE 00203 IF( NNT.GT.0 ) 00204 $ WRITE( NOUT, FMT = 9993 )'TYPE', ( NTVAL( I ), I = 1, NNT ) 00205 * 00206 * Read the threshold value for the test ratios. 00207 * 00208 READ( NIN, FMT = * )THRESH 00209 WRITE( NOUT, FMT = 9992 )THRESH 00210 * 00211 * Read the flag that indicates whether to test the error exits. 00212 * 00213 READ( NIN, FMT = * )TSTERR 00214 * 00215 IF( FATAL ) THEN 00216 WRITE( NOUT, FMT = 9999 ) 00217 STOP 00218 END IF 00219 * 00220 IF( FATAL ) THEN 00221 WRITE( NOUT, FMT = 9999 ) 00222 STOP 00223 END IF 00224 * 00225 * Calculate and print the machine dependent constants. 00226 * 00227 EPS = DLAMCH( 'Underflow threshold' ) 00228 WRITE( NOUT, FMT = 9991 )'underflow', EPS 00229 EPS = DLAMCH( 'Overflow threshold' ) 00230 WRITE( NOUT, FMT = 9991 )'overflow ', EPS 00231 EPS = DLAMCH( 'Epsilon' ) 00232 WRITE( NOUT, FMT = 9991 )'precision', EPS 00233 WRITE( NOUT, FMT = * ) 00234 * 00235 * Test the error exit of: 00236 * 00237 IF( TSTERR ) 00238 $ CALL DERRRFP( NOUT ) 00239 * 00240 * Test the routines: dpftrf, dpftri, dpftrs (as in DDRVPO). 00241 * This also tests the routines: dtfsm, dtftri, dtfttr, dtrttf. 00242 * 00243 CALL DDRVRFP( NOUT, NN, NVAL, NNS, NSVAL, NNT, NTVAL, THRESH, 00244 $ WORKA, WORKASAV, WORKAFAC, WORKAINV, WORKB, 00245 $ WORKBSAV, WORKXACT, WORKX, WORKARF, WORKARFINV, 00246 $ D_WORK_DLATMS, D_WORK_DPOT01, D_TEMP_DPOT02, 00247 $ D_TEMP_DPOT03, D_WORK_DLANSY, D_WORK_DPOT02, 00248 $ D_WORK_DPOT03 ) 00249 * 00250 * Test the routine: dlansf 00251 * 00252 CALL DDRVRF1( NOUT, NN, NVAL, THRESH, WORKA, NMAX, WORKARF, 00253 + D_WORK_DLANSY ) 00254 * 00255 * Test the convertion routines: 00256 * dtfttp, dtpttf, dtfttr, dtrttf, dtrttp and dtpttr. 00257 * 00258 CALL DDRVRF2( NOUT, NN, NVAL, WORKA, NMAX, WORKARF, 00259 + WORKAP, WORKASAV ) 00260 * 00261 * Test the routine: dtfsm 00262 * 00263 CALL DDRVRF3( NOUT, NN, NVAL, THRESH, WORKA, NMAX, WORKARF, 00264 + WORKAINV, WORKAFAC, D_WORK_DLANSY, 00265 + D_WORK_DPOT03, D_WORK_DPOT01 ) 00266 * 00267 * 00268 * Test the routine: dsfrk 00269 * 00270 CALL DDRVRF4( NOUT, NN, NVAL, THRESH, WORKA, WORKAFAC, NMAX, 00271 + WORKARF, WORKAINV, NMAX, D_WORK_DLANSY) 00272 * 00273 CLOSE ( NIN ) 00274 S2 = DSECND( ) 00275 WRITE( NOUT, FMT = 9998 ) 00276 WRITE( NOUT, FMT = 9997 )S2 - S1 00277 * 00278 9999 FORMAT( / ' Execution not attempted due to input errors' ) 00279 9998 FORMAT( / ' End of tests' ) 00280 9997 FORMAT( ' Total time used = ', F12.2, ' seconds', / ) 00281 9996 FORMAT( ' !! Invalid input value: ', A4, '=', I6, '; must be >=', 00282 $ I6 ) 00283 9995 FORMAT( ' !! Invalid input value: ', A4, '=', I6, '; must be <=', 00284 $ I6 ) 00285 9994 FORMAT( / ' Tests of the DOUBLE PRECISION LAPACK RFP routines ', 00286 $ / ' LAPACK VERSION ', I1, '.', I1, '.', I1, 00287 $ / / ' The following parameter values will be used:' ) 00288 9993 FORMAT( 4X, A4, ': ', 10I6, / 11X, 10I6 ) 00289 9992 FORMAT( / ' Routines pass computational tests if test ratio is ', 00290 $ 'less than', F8.2, / ) 00291 9991 FORMAT( ' Relative machine ', A, ' is taken to be', D16.6 ) 00292 * 00293 * End of DCHKRFP 00294 * 00295 END