MA02CD

Pertransposing the central band of a square matrix

[Specification] [Arguments] [Method] [References] [Comments] [Example]

Purpose

  To compute the pertranspose of a central band of a square matrix.

Specification
      SUBROUTINE MA02CD( N, KL, KU, A, LDA )
C     .. Scalar Arguments ..
      INTEGER          KL, KU, LDA, N
C     .. Array Arguments ..
      DOUBLE PRECISION A(LDA,*)

Arguments

Input/Output Parameters

  N       (input) INTEGER
          The order of the square matrix A.  N >= 0.

  KL      (input) INTEGER
          The number of subdiagonals of A to be pertransposed.
          0 <= KL <= N-1.

  KU      (input) INTEGER
          The number of superdiagonals of A to be pertransposed.
          0 <= KU <= N-1.

  A       (input/output) DOUBLE PRECISION array, dimension (LDA,N)
          On entry, the leading N-by-N part of this array must 
          contain a square matrix whose central band formed from 
          the KL subdiagonals, the main diagonal and the KU   
          superdiagonals will be pertransposed.
          On exit, the leading N-by-N part of this array contains 
          the matrix A with its central band (the KL subdiagonals,
          the main diagonal and the KU superdiagonals) pertransposed 
          (that is the elements of each antidiagonal appear in 
          reversed order). This is equivalent to forming P*B'*P,
          where B is the matrix formed from the central band of A
          and P is a permutation matrix with ones down the secondary
          diagonal.

  LDA     INTEGER
          The leading dimension of the array A.  LDA >= max(1,N).

Further Comments
  None
Example

Program Text

  None
Program Data
  None
Program Results
  None

Return to Supporting Routines index