dcoef.f
SUBROUTINE DCOEF (YH, YP, NCOMP, NROWB, NFC, NIC, B, BETA, COEF,
+ INHOMO, RE, AE, BY, CVEC, WORK, IWORK, IFLAG, NFCC)
C***BEGIN PROLOGUE DCOEF
....
....
Warning: this routine is not intended to be user-callable.
....
....
C***SUBSIDIARY
C***PURPOSE Subsidiary to DBVSUP
C***LIBRARY SLATEC
C***TYPE DOUBLE PRECISION (SCOEF-S, DCOEF-D)
C***AUTHOR Watts, H. A., (SNLA)
C***DESCRIPTION
C
C **********************************************************************
C INPUT to DCOEF
C **********************************************************************
C
C YH = matrix of homogeneous solutions.
C YP = vector containing particular solution.
C NCOMP = number of components per solution vector.
C NROWB = first dimension of B in calling program.
C NFC = number of base solution vectors.
C NFCC = 2*NFC for the special treatment of COMPLEX*16 valued
C equations. Otherwise, NFCC=NFC.
C NIC = number of specified initial conditions.
C B = boundary condition matrix at X = XFINAL.
C BETA = vector of nonhomogeneous boundary conditions at X = XFINAL.
C 1 - nonzero particular solution
C INHOMO = 2 - zero particular solution
C 3 - eigenvalue problem
C RE = relative error tolerance.
C AE = absolute error tolerance.
C BY = storage space for the matrix B*YH
C CVEC = storage space for the vector BETA-B*YP
C WORK = double precision array of internal storage. Dimension must
C be GE
C NFCC*(NFCC+4)
C IWORK = integer array of internal storage. Dimension must be GE
C 3+NFCC
C
C **********************************************************************
C OUTPUT from DCOEF
C **********************************************************************
C
C COEF = array containing superposition constants.
C IFLAG = indicator of success from DSUDS in solving the
C boundary equations.
C = 0 boundary equations are solved.
C = 1 boundary equations appear to have many solutions.
C = 2 boundary equations appear to be inconsistent.
C = 3 for this value of an eigenparameter, the boundary
C equations have only the zero solution.
C
C **********************************************************************
C
C Subroutine DCOEF solves for the superposition constants from the
C linear equations defined by the boundary conditions at X = XFINAL.
C
C B*YP + B*YH*COEF = BETA
C
C **********************************************************************
C
C***SEE ALSO DBVSUP
C***ROUTINES CALLED DDOT, DSUDS, XGETF, XSETF
C***COMMON BLOCKS DML5MC
C***REVISION HISTORY (YYMMDD)
C 750601 DATE WRITTEN
C 890531 Changed all specific intrinsics to generic. (WRB)
C 890831 Modified array declarations. (WRB)
C 890921 Realigned order of variables in certain COMMON blocks.
C (WRB)
C 890921 REVISION DATE from Version 3.2
C 891214 Prologue converted to Version 4.0 format. (BAB)
C 900328 Added TYPE section. (WRB)
C 910722 Updated AUTHOR section. (ALS)
C***END PROLOGUE DCOEF