X04CEF (PDF version)
X04 Chapter Contents
X04 Chapter Introduction
NAG Library Manual

NAG Library Routine Document

X04CEF

Note:  before using this routine, please read the Users' Note for your implementation to check the interpretation of bold italicised terms and other implementation-dependent details.

 Contents

    1  Purpose
    7  Accuracy

1  Purpose

X04CEF is an easy-to-use routine to print a real band matrix stored in a packed two-dimensional array.

2  Specification

SUBROUTINE X04CEF ( M, N, KL, KU, A, LDA, TITLE, IFAIL)
INTEGER  M, N, KL, KU, LDA, IFAIL
REAL (KIND=nag_wp)  A(LDA,*)
CHARACTER(*)  TITLE

3  Description

X04CEF prints a real band matrix stored in a packed two-dimensional array. It is an easy-to-use driver for X04CFF. The routine uses default values for the format in which numbers are printed, for labelling the rows and columns, and for output record length.
X04CEF will choose a format code such that numbers will be printed with an F8.4, an F11.4 or a 1PE13.4 format . The F8.4 code is chosen if the sizes of all the matrix elements to be printed lie between 0.001 and 1.0. The F11.4 code is chosen if the sizes of all the matrix elements to be printed lie between 0.001 and 9999.9999. Otherwise the 1PE13.4 code is chosen.
The matrix is printed with integer row and column labels, and with a maximum record length of 80.
The matrix is output to the unit defined by X04ABF.

4  References

None.

5  Parameters

1:     M – INTEGERInput
2:     N – INTEGERInput
On entry: the number of rows and columns of the band matrix, respectively, to be printed.
If either M or N is less than 1, X04CEF will exit immediately after printing TITLE; no row or column labels are printed.
3:     KL – INTEGERInput
On entry: the number of subdiagonals of the band matrix A.
Constraint: KL0.
4:     KU – INTEGERInput
On entry: the number of superdiagonals of the band matrix A.
Constraint: KU0.
5:     ALDA* – REAL (KIND=nag_wp) arrayInput
Note: the second dimension of the array A must be at least max1,minM+KU,N.
On entry: the band matrix to be printed.
The matrix is stored in rows 1 to kl+ku+1, more precisely, the element Aij must be stored in
Aku+1+i-jj  for ​max1,j-kuiminm,j+kl. 
6:     LDA – INTEGERInput
On entry: the first dimension of the array A as declared in the (sub)program from which X04CEF is called.
Constraint: LDAKL+KU+1.
7:     TITLE – CHARACTER(*)Input
On entry: a title to be printed above the matrix.
If TITLE=' ', no title (and no blank line) will be printed.
If TITLE contains more than 80 characters, the contents of TITLE will be wrapped onto more than one line, with the break after 80 characters.
Any trailing blank characters in TITLE are ignored.
8:     IFAIL – INTEGERInput/Output
On entry: IFAIL must be set to 0, -1​ or ​1. If you are unfamiliar with this parameter you should refer to Section 3.3 in the Essential Introduction for details.
For environments where it might be inappropriate to halt program execution when an error is detected, the value -1​ or ​1 is recommended. If the output of error messages is undesirable, then the value 1 is recommended. Otherwise, if you are not familiar with this parameter, the recommended value is 0. When the value -1​ or ​1 is used it is essential to test the value of IFAIL on exit.
On exit: IFAIL=0 unless the routine detects an error or a warning has been flagged (see Section 6).

6  Error Indicators and Warnings

If on entry IFAIL=0 or -1, explanatory error messages are output on the current error message unit (as defined by X04AAF).
Errors or warnings detected by the routine:
IFAIL=1
On entry,KL<0.
IFAIL=2
On entry,KU<0.
IFAIL=3
On entry,LDA<KL+KU+1.
IFAIL=-99
An unexpected error has been triggered by this routine. Please contact NAG.
See Section 3.8 in the Essential Introduction for further information.
IFAIL=-399
Your licence key may have expired or may not have been installed correctly.
See Section 3.7 in the Essential Introduction for further information.
IFAIL=-999
Dynamic memory allocation failed.
See Section 3.6 in the Essential Introduction for further information.

7  Accuracy

Not applicable.

8  Parallelism and Performance

Not applicable.

9  Further Comments

A call to X04CEF is equivalent to a call to X04CFF with the following argument values:

NCOLS = 80
INDENT = 0
LABROW = 'I'
LABCOL = 'I'
FORM = ' '

10  Example

The example program calls X04CEF to print a 5 by 5 band matrix with one subdiagonal and one superdiagonal.

10.1  Program Text

Program Text (x04cefe.f90)

10.2  Program Data

None.

10.3  Program Results

Program Results (x04cefe.r)


X04CEF (PDF version)
X04 Chapter Contents
X04 Chapter Introduction
NAG Library Manual

© The Numerical Algorithms Group Ltd, Oxford, UK. 2015