dlpdp.f
SUBROUTINE DLPDP (A, MDA, M, N1, N2, PRGOPT, X, WNORM, MODE, WS,
+ IS)
C***BEGIN PROLOGUE DLPDP
....
....
Warning: this routine is not intended to be user-callable.
....
....
C***SUBSIDIARY
C***PURPOSE Subsidiary to DLSEI
C***LIBRARY SLATEC
C***TYPE DOUBLE PRECISION (LPDP-S, DLPDP-D)
C***AUTHOR Hanson, R. J., (SNLA)
C Haskell, K. H., (SNLA)
C***DESCRIPTION
C
C **** Double Precision version of LPDP ****
C DIMENSION A(MDA,N+1),PRGOPT(*),X(N),WS((M+2)*(N+7)),IS(M+N+1),
C where N=N1+N2. This is a slight overestimate for WS(*).
C
C Determine an N1-vector W, and
C an N2-vector Z
C which minimizes the Euclidean length of W
C subject to G*W+H*Z .GE. Y.
C This is the least projected distance problem, LPDP.
C The matrices G and H are of respective
C dimensions M by N1 and M by N2.
C
C Called by subprogram DLSI( ).
C
C The matrix
C (G H Y)
C
C occupies rows 1,...,M and cols 1,...,N1+N2+1 of A(*,*).
C
C The solution (W) is returned in X(*).
C (Z)
C
C The value of MODE indicates the status of
C the computation after returning to the user.
C
C MODE=1 The solution was successfully obtained.
C
C MODE=2 The inequalities are inconsistent.
C
C***SEE ALSO DLSEI
C***ROUTINES CALLED DCOPY, DDOT, DNRM2, DSCAL, DWNNLS
C***REVISION HISTORY (YYMMDD)
C 790701 DATE WRITTEN
C 890531 Changed all specific intrinsics to generic. (WRB)
C 891214 Prologue converted to Version 4.0 format. (BAB)
C 900328 Added TYPE section. (WRB)
C 910408 Updated the AUTHOR section. (WRB)
C***END PROLOGUE DLPDP