NAME

SYNOPSIS

  SUBROUTINE DCHUD( A, LDA, N, X, Z, LDZ, NZ, Y, RHO, COS, SIN)
  INTEGER LDA, N, LDZ, NZ
  DOUBLE PRECISION A(LDA,*), X(*), Z(LDZ,*), Y(*), RHO(*), COS(*), SIN(*)
 
  SUBROUTINE DCHUD_64( A, LDA, N, X, Z, LDZ, NZ, Y, RHO, COS, SIN)
  INTEGER*8 LDA, N, LDZ, NZ
  DOUBLE PRECISION A(LDA,*), X(*), Z(LDZ,*), Y(*), RHO(*), COS(*), SIN(*)
 

C INTERFACE

void dchud(double *a, int lda, int n, double *x, double *z, int ldz, int nz, double *y, double *rho, double *cos, double *sin);

void dchud_64(double *a, long lda, long n, double *x, double *z, long ldz, long nz, double *y, double *rho, double *cos, double *sin);

ARGUMENTS

* A (input/output)

On entry, the upper triangular matrix A. On exit, A has been updated. The strict lower triangle of A is not referenced.

* LDA (input)

Leading dimension of the array A as specified in a dimension or type statement. LDA >= max(1,N).

* N (input)

Order of the matrix A. N >= 0.

* X (input)

Row to be added to A.

* Z (input/output)

Vectors to be updated with A.

* LDZ (input)

Leading dimension on the array Z as specified in a dimension or type statement. LDZ >= max(1,N).

* NZ (input)

Number of vectors to be updated with A. NZ >= 0. If NZ = 0 then Z, Y, and RHO are not used.

* Y (input)

Scalars for updating the vectors in Z.

* RHO (input/output)

On entry, the norms of the residual vectors that are to be updated. On exit, RHO has been updated. If RHO(i) is negative on entry then it is not changed.

* COS (input)

Cosines of the transforming rotations.

* SIN (input)

Sines of the transforming rotations.