37 #ifndef TEMPLATE_LAPACK_POCON_HEADER 38 #define TEMPLATE_LAPACK_POCON_HEADER 43 lda,
const Treal *anorm, Treal *rcond, Treal *work,
integer *
105 integer a_dim1, a_offset, i__1;
120 a_offset = 1 + a_dim1 * 1;
134 }
else if (*anorm < 0.) {
149 }
else if (*anorm == 0.) {
158 *(
unsigned char *)normin =
'N';
167 lda, &work[1], &scalel, &work[(*n << 1) + 1], info);
168 *(
unsigned char *)normin =
'Y';
173 a_offset], lda, &work[1], &scaleu, &work[(*n << 1) + 1],
180 a_offset], lda, &work[1], &scalel, &work[(*n << 1) + 1],
182 *(
unsigned char *)normin =
'Y';
187 lda, &work[1], &scaleu, &work[(*n << 1) + 1], info);
192 scale = scalel * scaleu;
195 if (scale < (d__1 = work[ix],
absMACRO(d__1)) * smlnum || scale == 0.)
207 *rcond = 1. / ainvnm / *anorm;
#define absMACRO(x)
Definition: template_blas_common.h:47
integer template_blas_idamax(const integer *n, const Treal *dx, const integer *incx)
Definition: template_blas_idamax.h:42
int template_lapack_pocon(const char *uplo, const integer *n, const Treal *a, const integer *lda, const Treal *anorm, Treal *rcond, Treal *work, integer *iwork, integer *info)
Definition: template_lapack_pocon.h:42
int integer
Definition: template_blas_common.h:40
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
int template_lapack_lacon(const integer *n, Treal *v, Treal *x, integer *isgn, Treal *est, integer *kase)
Definition: template_lapack_lacon.h:42
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
int template_lapack_latrs(const char *uplo, const char *trans, const char *diag, const char *normin, const integer *n, const Treal *a, const integer *lda, Treal *x, Treal *scale, Treal *cnorm, integer *info)
Definition: template_lapack_latrs.h:42
Treal template_lapack_lamch(const char *cmach, Treal dummyReal)
Definition: template_lapack_lamch.h:202
int template_lapack_rscl(const integer *n, const Treal *sa, Treal *sx, const integer *incx)
Definition: template_lapack_rscl.h:42
bool logical
Definition: template_blas_common.h:41
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46