NAG Library Function Document

nag_gen_real_mat_print_comp (x04cbc)

1  Purpose

2  Specification

3  Description

4  References

5  Arguments

1:     order – Nag_OrderTypeInput

On entry: the order argument specifies the two-dimensional storage scheme being used, i.e., row-major ordering or column-major ordering. C language defined storage is specified by $\mathbf{order}=\mathrm{Nag\_RowMajor}$. See Section 3.2.1.3 in the Essential Introduction for a more detailed explanation of the use of this argument.

Constraint: $\mathbf{order}=\mathrm{Nag\_RowMajor}$ or $\mathrm{Nag\_ColMajor}$.

2:     matrix – Nag_MatrixTypeInput

On entry: indicates the part of the matrix to be printed.

$\mathbf{matrix}=\mathrm{Nag\_GeneralMatrix}$

The whole of the rectangular matrix.

$\mathbf{matrix}=\mathrm{Nag\_LowerMatrix}$

The lower triangle of the matrix, or the lower trapezium if the matrix has more rows than columns.

$\mathbf{matrix}=\mathrm{Nag\_UpperMatrix}$

The upper triangle of the matrix, or the upper trapezium if the matrix has more columns than rows.

Constraint: $\mathbf{matrix}=\mathrm{Nag\_GeneralMatrix}$, $\mathrm{Nag\_LowerMatrix}$ or $\mathrm{Nag\_UpperMatrix}$.

3:     diag – Nag_DiagTypeInput

On entry: indicates whether the diagonal elements of the matrix are to be printed.

$\mathbf{diag}=\mathrm{Nag\_NonRefDiag}$

The diagonal elements of the matrix are not referenced and not printed.

$\mathbf{diag}=\mathrm{Nag\_UnitDiag}$

The diagonal elements of the matrix are not referenced, but are assumed all to be unity, and are printed as such.

$\mathbf{diag}=\mathrm{Nag\_NonUnitDiag}$

The diagonal elements of the matrix are referenced and printed.

If $\mathbf{matrix}=\mathrm{Nag\_GeneralMatrix}$, then diag must be set to $\mathrm{Nag\_NonUnitDiag}$.

Constraints:

• if $\mathbf{matrix}\ne \mathrm{Nag\_GeneralMatrix}$, $\mathbf{diag}=\mathrm{Nag\_NonRefDiag}$, $\mathrm{Nag\_UnitDiag}$ or $\mathrm{Nag\_NonUnitDiag}$;
• if $\mathbf{matrix}=\mathrm{Nag\_GeneralMatrix}$, $\mathbf{diag}=\mathrm{Nag\_NonUnitDiag}$.

4:     m – IntegerInput

5:     n – IntegerInput

On entry: the number of rows and columns of the matrix, respectively, to be printed.

If either m or n is less than $1$, nag_gen_real_mat_print_comp (x04cbc) will exit immediately after printing title; no row or column labels are printed.

6:     a[$\mathit{dim}$] – const doubleInput

Note: the dimension, dim, of the array a must be at least

• $\max \left(1,\mathbf{pda}\times \mathbf{n}\right)$ when $\mathbf{order}=\mathrm{Nag\_ColMajor}$;
• $\max \left(1,\mathbf{m}\times \mathbf{pda}\right)$ when $\mathbf{order}=\mathrm{Nag\_RowMajor}$.

The $\left(i,j\right)$th element of the matrix $A$ is stored in

• $\mathbf{a}\left[\left(j-1\right)\times \mathbf{pda}+i-1\right]$ when $\mathbf{order}=\mathrm{Nag\_ColMajor}$;
• $\mathbf{a}\left[\left(i-1\right)\times \mathbf{pda}+j-1\right]$ when $\mathbf{order}=\mathrm{Nag\_RowMajor}$.

On entry: the matrix to be printed. Only the elements that will be referred to, as specified by arguments matrix and diag, need be set.

7:     pda – IntegerInput

On entry: the stride separating row or column elements (depending on the value of order) in the array a.

Constraints:

• if $\mathbf{order}=\mathrm{Nag\_ColMajor}$, $\mathbf{pda}\ge \max \left(1,\mathbf{m}\right)$;
• if $\mathbf{order}=\mathrm{Nag\_RowMajor}$, $\mathbf{pda}\ge \max \left(1,\mathbf{n}\right)$.

8:     form – const char *Input

On entry: a valid C format code. This should be of the form $\%\left[\mathit{flag}\right]ww.pp\left[\mathit{format}\text{ }\mathit{indicator}\right]$, where $\mathit{ww}.pp$ indicates that up to two digits may be used to specify the field width and precision respectively. Only % and $\mathit{format}\text{ }\mathit{indicator}$ must be present. $\mathit{flag}$ can be one of $-$, $+$, $<\text{space}>$ or # and $\mathit{format}\text{ }\mathit{indicator}$ can be e, E, f, g or G. Thus, possible formats include %f, %+23.15G, %.6e. form is used to print elements of the matrix $A$.

In addition, nag_gen_real_mat_print_comp (x04cbc) chooses its own format code when form is NULL or $\mathbf{form}=\text{'*'}$.

$\mathbf{form}=\mathbf{NULL}$

nag_gen_real_mat_print_comp (x04cbc) will choose a format code such that numbers will be printed with either a %8.4f, a %11.4f or a %13.4e format. The %8.4f code is chosen if the sizes of all the matrix elements to be printed lie between $0.001$ and $1.0$. The %11.4f code is chosen if the sizes of all the matrix elements to be printed lie between $0.001$ and $9999.9999$. Otherwise the %13.4e code is chosen.

$\mathbf{form}=\text{'*'}$

nag_gen_real_mat_print_comp (x04cbc) will choose a format code such that numbers will be printed to as many significant digits as are necessary to distinguish between neighbouring machine numbers. Thus any two numbers that are stored with different internal representations should look different on output.

Constraint: form must be of the form $\%\left[\mathit{flag}\right]ww.pp\left[\mathit{format}\text{ }\mathit{indicator}\right]$.

9:     title – const char *Input

On entry: a title to be printed above the matrix, or name of the matrix.

If $\mathbf{title}=\mathbf{NULL}$, no title (and no blank line) will be printed.

If title contains more than ncols characters, the contents of title will be wrapped onto more than one line, with the break after ncols characters.

Any trailing blank characters in title are ignored.

10:   labrow – Nag_LabelTypeInput

On entry: indicates the type of labelling to be applied to the rows of the matrix.

$\mathbf{labrow}=\mathrm{Nag\_NoLabels}$

Prints no row labels.

$\mathbf{labrow}=\mathrm{Nag\_IntegerLabels}$

Prints integer row labels.

$\mathbf{labrow}=\mathrm{Nag\_CharacterLabels}$

Prints character labels, which must be supplied in array rlabs.

Constraint: $\mathbf{labrow}=\mathrm{Nag\_NoLabels}$, $\mathrm{Nag\_IntegerLabels}$ or $\mathrm{Nag\_CharacterLabels}$.

11:   rlabs[$\mathit{dim}$] – const char *Input

Note: the dimension, dim, of the array rlabs must be at least

• $\mathbf{m}$ when $\mathbf{labrow}=\mathrm{Nag\_CharacterLabels}$;
• otherwise rlabs may be NULL.

On entry: if $\mathbf{labrow}=\mathrm{Nag\_CharacterLabels}$, rlabs must contain labels for the rows of the matrix; otherwise rlabs is not referenced and may be NULL.

Labels are right-justified when output, in a field which is as wide as necessary to hold the longest row label. Note that this field width is subtracted from the number of usable columns, ncols.

12:   labcol – Nag_LabelTypeInput

On entry: indicates the type of labelling to be applied to the columns of the matrix.

$\mathbf{labcol}=\mathrm{Nag\_NoLabels}$

Prints no column labels.

$\mathbf{labcol}=\mathrm{Nag\_IntegerLabels}$

Prints integer column labels.

$\mathbf{labcol}=\mathrm{Nag\_CharacterLabels}$

Prints character labels, which must be supplied in array clabs.

Constraint: $\mathbf{labcol}=\mathrm{Nag\_NoLabels}$, $\mathrm{Nag\_IntegerLabels}$ or $\mathrm{Nag\_CharacterLabels}$.

13:   clabs[$\mathit{dim}$] – const char *Input

Note: the dimension, dim, of the array clabs must be at least

• $\mathbf{n}$ when $\mathbf{labcol}=\mathrm{Nag\_CharacterLabels}$;
• otherwise clabs may be NULL.

On entry: if $\mathbf{labcol}=\mathrm{Nag\_CharacterLabels}$, clabs must contain labels for the columns of the matrix; otherwise clabs is not referenced and may be NULL.

Labels are right-justified when output. Any label that is too long for the column width, which is determined by form, is truncated.

14:   ncols – IntegerInput

On entry: the maximum output record length. If the number of columns of the matrix is too large to be accommodated in ncols characters, the matrix will be printed in parts, containing the largest possible number of matrix columns, and each part separated by a blank line.

ncols must be large enough to hold at least one column of the matrix using the format specifier in form. If a value less than or equal to 0 or greater than $132$ is supplied for ncols, then the value $80$ is used instead.

15:   indent – IntegerInput

On entry: the number of columns by which the matrix (and any title and labels) should be indented. The effective value of ncols is reduced by indent columns. If a value less than $0$ or greater than ncols is supplied for indent, the value $0$ is used instead.

16:   outfile – const char *Input

On entry: the name of a file to which output will be directed. If outfile is NULL the output will be directed to standard output.

17:   fail – NagError *Input/Output

The NAG error argument (see Section 3.6 in the Essential Introduction).

6  Error Indicators and Warnings

NE_ALLOC_FAIL

Memory allocation failed.

NE_BAD_PARAM

On entry, argument $⟨\mathit{\text{value}}⟩$ had an illegal value.

NE_COL_WIDTH

$⟨\mathit{\text{value}}⟩$ is not wide enough to hold at least one matrix column. $\mathbf{ncols}=⟨\mathit{\text{value}}⟩$ and $\mathbf{indent}=⟨\mathit{\text{value}}⟩$.

NE_INTERNAL_ERROR

An internal error has occurred in this function. Check the function call and any array sizes. If the call is correct then please contact NAG for assistance.

NE_INVALID_FORMAT

The string  has not been recognized as a valid format.

NE_NOT_APPEND_FILE

Cannot open file  for appending.

NE_NOT_CLOSE_FILE

Cannot close file .

NE_NOT_WRITE_FILE

Cannot open file  for writing.

7  Accuracy

8  Parallelism and Performance

9  Further Comments

#include <nag.h>
#include <nagx04.h>
#include <nag_stdlib.h>    
double *a = 0;
Integer n = 4;
if ( !(a = NAG_ALLOC(n, double)) )
  {
    Vprintf("Allocation failure\n");
    return -1;
  }
/* Read A from data file */
for (i = 0; i < n; ++i)
  Vscanf("%lf", &a[i]);
/* Print vector A as a row vector */
x04cbc(Nag_RowMajor, Nag_GeneralMatrix, Nag_NonUnitDiag,
1, n, a, n, 0, 0,  Nag_NoLabels, 0, Nag_IntegerLabels, 0,
0, 0, 0, NAGERR_DEFAULT);
/* Print vector A as a column vector */
x04cbc(Nag_RowMajor, Nag_GeneralMatrix, Nag_NonUnitDiag,
n, 1, a, 1, 0, 0,  Nag_IntegerLabels, 0, Nag_NoLabels, 0,
0, 0, 0, NAGERR_DEFAULT);

10  Example