The NTRA argument can be equal to -1 if one does not want a snapshot file
(and this is the case with sample data {s,d,c,z}in3).
The routines {S,D,C,Z}PRCN3 will try to use their first argument as an output
unit number, so we avoid calling them when NTRA < 0.
Patch originally written by Camm Maguire.
tags/v0.2.16.rc1
| @@ -1365,8 +1365,9 @@ | |||||
| * | * | ||||
| 150 CONTINUE | 150 CONTINUE | ||||
| WRITE( NOUT, FMT = 9996 )SNAME | WRITE( NOUT, FMT = 9996 )SNAME | ||||
| CALL CPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| IF( TRACE ) | |||||
| $ CALL CPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| * | * | ||||
| 160 CONTINUE | 160 CONTINUE | ||||
| RETURN | RETURN | ||||
| @@ -1365,8 +1365,9 @@ | |||||
| * | * | ||||
| 150 CONTINUE | 150 CONTINUE | ||||
| WRITE( NOUT, FMT = 9996 )SNAME | WRITE( NOUT, FMT = 9996 )SNAME | ||||
| CALL CPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| IF( TRACE ) | |||||
| $ CALL CPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| * | * | ||||
| 160 CONTINUE | 160 CONTINUE | ||||
| RETURN | RETURN | ||||
| @@ -1335,8 +1335,9 @@ | |||||
| * | * | ||||
| 150 CONTINUE | 150 CONTINUE | ||||
| WRITE( NOUT, FMT = 9996 )SNAME | WRITE( NOUT, FMT = 9996 )SNAME | ||||
| CALL DPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| IF( TRACE ) | |||||
| $ CALL DPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| * | * | ||||
| 160 CONTINUE | 160 CONTINUE | ||||
| RETURN | RETURN | ||||
| @@ -1339,8 +1339,9 @@ | |||||
| * | * | ||||
| 150 CONTINUE | 150 CONTINUE | ||||
| WRITE( NOUT, FMT = 9996 )SNAME | WRITE( NOUT, FMT = 9996 )SNAME | ||||
| CALL SPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| IF( TRACE ) | |||||
| $ CALL SPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| * | * | ||||
| 160 CONTINUE | 160 CONTINUE | ||||
| RETURN | RETURN | ||||
| @@ -1366,8 +1366,9 @@ | |||||
| * | * | ||||
| 150 CONTINUE | 150 CONTINUE | ||||
| WRITE( NOUT, FMT = 9996 )SNAME | WRITE( NOUT, FMT = 9996 )SNAME | ||||
| CALL ZPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| IF( TRACE ) | |||||
| $ CALL ZPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| * | * | ||||
| 160 CONTINUE | 160 CONTINUE | ||||
| RETURN | RETURN | ||||
| @@ -1366,8 +1366,9 @@ | |||||
| * | * | ||||
| 150 CONTINUE | 150 CONTINUE | ||||
| WRITE( NOUT, FMT = 9996 )SNAME | WRITE( NOUT, FMT = 9996 )SNAME | ||||
| CALL ZPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| IF( TRACE ) | |||||
| $ CALL ZPRCN3( NTRA, NC, SNAME, IORDER, SIDE, UPLO, TRANSA, DIAG, | |||||
| $ M, N, ALPHA, LDA, LDB) | |||||
| * | * | ||||
| 160 CONTINUE | 160 CONTINUE | ||||
| RETURN | RETURN | ||||