NAG Library Routine Document

f06thf (zmload)

 Contents

    1  Purpose
    7  Accuracy
    10  Example

1
Purpose

f06thf forms the complex m by n rectangular or trapezoidal matrix A given by
aij = diag if i=j const if ij .  

2
Specification

Fortran Interface
Subroutine f06thf ( matrix, m, n, con, diag, a, lda)
Integer, Intent (In):: m, n, lda
Complex (Kind=nag_wp), Intent (In):: con, diag
Complex (Kind=nag_wp), Intent (Inout):: a(lda,*)
Character (1), Intent (In):: matrix
C Header Interface
#include nagmk26.h
void  f06thf_ (const char *matrix, const Integer *m, const Integer *n, const Complex *con, const Complex *diag, Complex a[], const Integer *lda, const Charlen length_matrix)

3
Description

None.

4
References

None.

5
Arguments

1:     matrix – Character(1)Input
On entry: the matrix type.
matrix='G'
General matrix.
matrix='U'
Upper trapezoidal matrix (upper triangular if m=n).
matrix='L'
Lower trapezoidal matrix (lower triangular if m=n).
Constraint: matrix='G', 'U' or 'L'.
2:     m – IntegerInput
On entry: m, the number of rows of the matrix A.
Constraint: m0.
3:     n – IntegerInput
On entry: n, the number of columns of the matrix A.
Constraint: n0.
4:     con – Complex (Kind=nag_wp)Input
On entry: the value to be assigned to the off-diagonal elements of A.
5:     diag – Complex (Kind=nag_wp)Input
On entry: the value to be assigned to the diagonal elements of A.
6:     alda* – Complex (Kind=nag_wp) arrayOutput
Note: the second dimension of the array a must be at least n.
On exit: the m by n general or trapezoidal matrix A.
  • If matrix='U', A is upper trapezoidal and the elements of the array below the diagonal are not referenced.
  • If matrix='L', A is lower trapezoidal and the elements of the array above the diagonal are not referenced.
7:     lda – IntegerInput
On entry: the first dimension of the array a as declared in the (sub)program from which f06thf is called.
Constraint: lda max1,m .

6
Error Indicators and Warnings

None.

7
Accuracy

Not applicable.

8
Parallelism and Performance

f06thf is not threaded in any implementation.

9
Further Comments

None.

10
Example

None.
© The Numerical Algorithms Group Ltd, Oxford, UK. 2017