·  Home
+   man pages
 -> Linux -> FreeBSD -> OpenBSD -> NetBSD -> Tru64 Unix -> HP-UX 11i -> IRIX
·  Linux HOWTOs
·  FreeBSD Tips
·  *niX Forums

man pages->IRIX man pages -> complib/CPOCO (3)
 Title
 Content
 Arch
 Section All Sections 1 - General Commands 2 - System Calls 3 - Subroutines 4 - Special Files 5 - File Formats 6 - Games 7 - Macros and Conventions 8 - Maintenance Commands 9 - Kernel Interface n - New Commands

### Contents

```
CPOCO(3F)							     CPOCO(3F)

```

### NAME[Toc][Back]

```     CPOCO   - CPOCO factors a complex Hermitian positive definite matrix and
estimates the condition of	the matrix.

If	 RCOND	is not needed, CPOFA is	slightly faster.  To solve  A*X	= B ,
follow CPOCO by CPOSL.  To	compute	 INVERSE(A)*C ,	follow CPOCO by	CPOSL.
To	compute	 DETERMINANT(A)	, follow CPOCO by CPODI.  To compute
INVERSE(A)	, follow CPOCO by CPODI.

```

### SYNOPSYS[Toc][Back]

```      SUBROUTINE CPOCO(A,LDA,N,RCOND,Z,INFO)
```

### DESCRIPTION[Toc][Back]

```     On	Entry

A COMPLEX(LDA, N)
the Hermitian matrix to	be factored.  Only the
diagonal and upper triangle are	used.

LDA INTEGER
the leading dimension of the array  A .

N INTEGER
the order of the matrix	 A .  On Return

A an upper	triangular matrix  R  so that  A =
CTRANS(R)*R where  CTRANS(R)  is the conjugate
transpose.  The	strict lower triangle is unaltered.
If  INFO .NE. 0	, the factorization is not complete.

RCOND REAL
an estimate of the reciprocal condition	of  A .
For the	system	A*X = B	, relative perturbations
in  A  and  B  of size	EPSILON	 may cause
relative perturbations in  X  of size  EPSILON/RCOND .
If  RCOND  is so small that the	logical	expression
1.0 + RCOND .EQ. 1.0
is true, then  A  may be singular to working
precision.  In particular,  RCOND  is zero  if
exact singularity is detected or the estimate
underflows.  If	INFO .NE. 0 , RCOND is unchanged.

Z COMPLEX(N)
a work vector whose contents are usually unimportant.
If  A  is close	to a singular matrix, then  Z  is
an approximate null vector in the sense	that
NORM(A*Z) = RCOND*NORM(A)*NORM(Z) .
If  INFO .NE. 0	, Z  is	unchanged.

INFO INTEGER
= 0  for normal	return.

Page 1

CPOCO(3F)							     CPOCO(3F)

= K  signals an	error condition.  The leading minor
of order  K  is	not positive definite.	LINPACK.  This version dated
08/14/78 .	 Cleve Moler, University of New	Mexico,	Argonne	National Lab.
Subroutines and Functions LINPACK CPOFA BLAS CAXPY,CDOTC,CSSCAL,SCASUM
Fortran ABS,AIMAG,AMAX1,CMPLX,CONJG,REAL

PPPPaaaaggggeeee 2222```
[ Back ]
Similar pages