--- lapack.h.orig 2009-03-03 16:06:07.000000000 +0000 +++ lapack.h 2009-03-03 16:07:13.000000000 +0000 @@ -11,6 +11,7 @@ * */ /* Modified to include xGEQP3 by Kian Ming A. Chai (14 May 2008) */ +/* Modified to include ppcon/trcon by Kian Ming A. Chai (3 March 2009) */ #ifndef BOOST_NUMERIC_BINDINGS_LAPACK_LAPACK_H #define BOOST_NUMERIC_BINDINGS_LAPACK_LAPACK_H @@ -764,7 +765,41 @@ dcomplex_t *a, const int *lda, dcomplex_t *b, const int *ldb, double *s, const double *rcond, int *rank, dcomplex_t *work, const int *lwork, double *rwork, int *iwork, int *info); + + + /*********************************************************************/ + /* Get condition numbers */ + /*********************************************************************/ + + void LAPACK_STRCON( const char *norm, const char *uplo, const char *diag, + const int *n, const float *a, const int *lda, + float *rcond, float *work, int *iwork, int *info); + + void LAPACK_DTRCON( const char *norm, const char *uplo, const char *diag, + const int *n, const double *a, const int *lda, + double *rcond, double *work, int *iwork, int *info); + + void LAPACK_CTRCON( const char *norm, const char *uplo, const char *diag, + const int *n, const fcomplex_t *a, const int *lda, + float *rcond, fcomplex_t *work, float *rwork, int *info); + + void LAPACK_ZTRCON( const char *norm, const char *uplo, const char *diag, + const int *n, const dcomplex_t *a, const int *lda, + double *rcond, dcomplex_t *work, double *rwork, int *info); + + + + void LAPACK_SPPCON( const char *uplo, const int *n, const float *ap, const float *anorm, + float *rcond, float *work, int *iwork, int *info); + + void LAPACK_DPPCON( const char *uplo, const int *n, const double *ap, const double *anorm, + double *rcond, double *work, int *iwork, int *info); + + void LAPACK_CPPCON( const char *uplo, const int *n, const fcomplex_t *ap, const float *anorm, + float *rcond, fcomplex_t *work, float *rwork, int *info); + void LAPACK_ZPPCON( const char *uplo, const int *n, const dcomplex_t *ap, const double *anorm, + double *rcond, dcomplex_t *work, double *rwork, int *info); /********************************************************************/