dsifa
dsifa - (obsolete) compute the UDU factorization of a symmetric matrix A. It is typical to
follow a call to SSIFA with a call to SSISL to solve Ax = b or to SSIDI
to compute the determinant, inverse, and inertia of A.
SUBROUTINE DSIFA( A, LDA, N, IPIVOT, INFO)
INTEGER LDA, N, INFO
INTEGER IPIVOT(*)
DOUBLE PRECISION A(LDA,*)
SUBROUTINE DSIFA_64( A, LDA, N, IPIVOT, INFO)
INTEGER*8 LDA, N, INFO
INTEGER*8 IPIVOT(*)
DOUBLE PRECISION A(LDA,*)
#include <sunperf.h>
void dsifa(double *a, int lda, int n, int *ipivot, int *info);
void dsifa_64(double *a, long lda, long n, long *ipivot, long *info);
-
* A (input/output)
-
On entry, the upper triangle of the matrix A.
On exit, a UDU factorization of the matrix A. 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.
-
* IPIVOT (output)
-
On exit, a vector of pivot indices.
-
* INFO (output)
-
On exit:
INFO = 0 Subroutine completed normally.
INFO > 0 Returns a value k if the kth pivot block is singular to
indicate that SSISL or SSIDI will divide by zero if called.