du12us.f
SUBROUTINE DU12US (A, MDA, M, N, B, MDB, NB, MODE, KRANK, RNORM,
+ H, W, IR, IC)
C***BEGIN PROLOGUE DU12US
....
....
Warning: this routine is not intended to be user-callable.
....
....
C***SUBSIDIARY
C***PURPOSE Subsidiary to DULSIA
C***LIBRARY SLATEC
C***TYPE DOUBLE PRECISION (U12US-S, DU12US-D)
C***AUTHOR (UNKNOWN)
C***DESCRIPTION
C
C Given the Householder LQ factorization of A, this
C subroutine solves the system AX=B. If the system
C is of reduced rank, this routine returns a solution
C according to the selected mode.
C
C Note - If MODE.NE.2, W is never accessed.
C
C***SEE ALSO DULSIA
C***ROUTINES CALLED DAXPY, DDOT, DNRM2, DSWAP
C***REVISION HISTORY (YYMMDD)
C 810801 DATE WRITTEN
C 890531 Changed all specific intrinsics to generic. (WRB)
C 890831 Modified array declarations. (WRB)
C 891214 Prologue converted to Version 4.0 format. (BAB)
C 900328 Added TYPE section. (WRB)
C***END PROLOGUE DU12US