| /* blas/gsl_blas.h |
| * |
| * Copyright (C) 1996, 1997, 1998, 1999, 2000 Gerard Jungman |
| * |
| * This program is free software; you can redistribute it and/or modify |
| * it under the terms of the GNU General Public License as published by |
| * the Free Software Foundation; either version 2 of the License, or (at |
| * your option) any later version. |
| * |
| * This program is distributed in the hope that it will be useful, but |
| * WITHOUT ANY WARRANTY; without even the implied warranty of |
| * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
| * General Public License for more details. |
| * |
| * You should have received a copy of the GNU General Public License |
| * along with this program; if not, write to the Free Software |
| * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. |
| */ |
| |
| /* |
| * Author: G. Jungman |
| */ |
| #ifndef __GSL_BLAS_H__ |
| #define __GSL_BLAS_H__ |
| |
| #include <gsl/gsl_vector.h> |
| #include <gsl/gsl_matrix.h> |
| |
| #include <gsl/gsl_blas_types.h> |
| |
| |
| #undef __BEGIN_DECLS |
| #undef __END_DECLS |
| #ifdef __cplusplus |
| # define __BEGIN_DECLS extern "C" { |
| # define __END_DECLS } |
| #else |
| # define __BEGIN_DECLS /* empty */ |
| # define __END_DECLS /* empty */ |
| #endif |
| |
| __BEGIN_DECLS |
| |
| |
| /* ======================================================================== |
| * Level 1 |
| * ======================================================================== |
| */ |
| |
| int gsl_blas_sdsdot (float alpha, |
| const gsl_vector_float * X, |
| const gsl_vector_float * Y, |
| float * result |
| ); |
| |
| int gsl_blas_dsdot (const gsl_vector_float * X, |
| const gsl_vector_float * Y, |
| double * result |
| ); |
| |
| int gsl_blas_sdot (const gsl_vector_float * X, |
| const gsl_vector_float * Y, |
| float * result |
| ); |
| |
| int gsl_blas_ddot (const gsl_vector * X, |
| const gsl_vector * Y, |
| double * result |
| ); |
| |
| |
| int gsl_blas_cdotu (const gsl_vector_complex_float * X, |
| const gsl_vector_complex_float * Y, |
| gsl_complex_float * dotu); |
| |
| int gsl_blas_cdotc (const gsl_vector_complex_float * X, |
| const gsl_vector_complex_float * Y, |
| gsl_complex_float * dotc); |
| |
| int gsl_blas_zdotu (const gsl_vector_complex * X, |
| const gsl_vector_complex * Y, |
| gsl_complex * dotu); |
| |
| int gsl_blas_zdotc (const gsl_vector_complex * X, |
| const gsl_vector_complex * Y, |
| gsl_complex * dotc); |
| |
| |
| float gsl_blas_snrm2 (const gsl_vector_float * X); |
| float gsl_blas_sasum (const gsl_vector_float * X); |
| double gsl_blas_dnrm2 (const gsl_vector * X); |
| double gsl_blas_dasum (const gsl_vector * X); |
| float gsl_blas_scnrm2 (const gsl_vector_complex_float * X); |
| float gsl_blas_scasum (const gsl_vector_complex_float * X); |
| double gsl_blas_dznrm2 (const gsl_vector_complex * X); |
| double gsl_blas_dzasum (const gsl_vector_complex * X); |
| |
| |
| CBLAS_INDEX_t gsl_blas_isamax (const gsl_vector_float * X); |
| CBLAS_INDEX_t gsl_blas_idamax (const gsl_vector * X); |
| CBLAS_INDEX_t gsl_blas_icamax (const gsl_vector_complex_float * X); |
| CBLAS_INDEX_t gsl_blas_izamax (const gsl_vector_complex * X); |
| |
| |
| int gsl_blas_sswap (gsl_vector_float * X, |
| gsl_vector_float * Y); |
| |
| int gsl_blas_scopy (const gsl_vector_float * X, |
| gsl_vector_float * Y); |
| |
| int gsl_blas_saxpy (float alpha, |
| const gsl_vector_float * X, |
| gsl_vector_float * Y); |
| |
| int gsl_blas_dswap (gsl_vector * X, |
| gsl_vector * Y); |
| |
| int gsl_blas_dcopy (const gsl_vector * X, |
| gsl_vector * Y); |
| |
| int gsl_blas_daxpy (double alpha, |
| const gsl_vector * X, |
| gsl_vector * Y); |
| |
| int gsl_blas_cswap (gsl_vector_complex_float * X, |
| gsl_vector_complex_float * Y); |
| |
| int gsl_blas_ccopy (const gsl_vector_complex_float * X, |
| gsl_vector_complex_float * Y); |
| |
| int gsl_blas_caxpy (const gsl_complex_float alpha, |
| const gsl_vector_complex_float * X, |
| gsl_vector_complex_float * Y); |
| |
| int gsl_blas_zswap (gsl_vector_complex * X, |
| gsl_vector_complex * Y); |
| |
| int gsl_blas_zcopy (const gsl_vector_complex * X, |
| gsl_vector_complex * Y); |
| |
| int gsl_blas_zaxpy (const gsl_complex alpha, |
| const gsl_vector_complex * X, |
| gsl_vector_complex * Y); |
| |
| |
| int gsl_blas_srotg (float a[], float b[], float c[], float s[]); |
| |
| int gsl_blas_srotmg (float d1[], float d2[], float b1[], float b2, float P[]); |
| |
| int gsl_blas_srot (gsl_vector_float * X, |
| gsl_vector_float * Y, |
| float c, float s); |
| |
| int gsl_blas_srotm (gsl_vector_float * X, |
| gsl_vector_float * Y, |
| const float P[]); |
| |
| int gsl_blas_drotg (double a[], double b[], double c[], double s[]); |
| |
| int gsl_blas_drotmg (double d1[], double d2[], double b1[], |
| double b2, double P[]); |
| |
| int gsl_blas_drot (gsl_vector * X, |
| gsl_vector * Y, |
| const double c, const double s); |
| |
| int gsl_blas_drotm (gsl_vector * X, |
| gsl_vector * Y, |
| const double P[]); |
| |
| |
| void gsl_blas_sscal (float alpha, gsl_vector_float * X); |
| void gsl_blas_dscal (double alpha, gsl_vector * X); |
| void gsl_blas_cscal (const gsl_complex_float alpha, gsl_vector_complex_float * X); |
| void gsl_blas_zscal (const gsl_complex alpha, gsl_vector_complex * X); |
| void gsl_blas_csscal (float alpha, gsl_vector_complex_float * X); |
| void gsl_blas_zdscal (double alpha, gsl_vector_complex * X); |
| |
| |
| /* =========================================================================== |
| * Level 2 |
| * =========================================================================== |
| */ |
| |
| /* |
| * Routines with standard 4 prefixes (S, D, C, Z) |
| */ |
| int gsl_blas_sgemv (CBLAS_TRANSPOSE_t TransA, |
| float alpha, |
| const gsl_matrix_float * A, |
| const gsl_vector_float * X, |
| float beta, |
| gsl_vector_float * Y); |
| |
| int gsl_blas_strmv (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, |
| const gsl_matrix_float * A, |
| gsl_vector_float * X); |
| |
| int gsl_blas_strsv (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, |
| const gsl_matrix_float * A, |
| gsl_vector_float * X); |
| |
| int gsl_blas_dgemv (CBLAS_TRANSPOSE_t TransA, |
| double alpha, |
| const gsl_matrix * A, |
| const gsl_vector * X, |
| double beta, |
| gsl_vector * Y); |
| |
| int gsl_blas_dtrmv (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, |
| const gsl_matrix * A, |
| gsl_vector * X); |
| |
| int gsl_blas_dtrsv (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, |
| const gsl_matrix * A, |
| gsl_vector * X); |
| |
| int gsl_blas_cgemv (CBLAS_TRANSPOSE_t TransA, |
| const gsl_complex_float alpha, |
| const gsl_matrix_complex_float * A, |
| const gsl_vector_complex_float * X, |
| const gsl_complex_float beta, |
| gsl_vector_complex_float * Y); |
| |
| int gsl_blas_ctrmv (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, |
| const gsl_matrix_complex_float * A, |
| gsl_vector_complex_float * X); |
| |
| int gsl_blas_ctrsv (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, |
| const gsl_matrix_complex_float * A, |
| gsl_vector_complex_float * X); |
| |
| int gsl_blas_zgemv (CBLAS_TRANSPOSE_t TransA, |
| const gsl_complex alpha, |
| const gsl_matrix_complex * A, |
| const gsl_vector_complex * X, |
| const gsl_complex beta, |
| gsl_vector_complex * Y); |
| |
| int gsl_blas_ztrmv (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, |
| const gsl_matrix_complex * A, |
| gsl_vector_complex * X); |
| |
| int gsl_blas_ztrsv (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t TransA, CBLAS_DIAG_t Diag, |
| const gsl_matrix_complex * A, |
| gsl_vector_complex *X); |
| |
| /* |
| * Routines with S and D prefixes only |
| */ |
| int gsl_blas_ssymv (CBLAS_UPLO_t Uplo, |
| float alpha, |
| const gsl_matrix_float * A, |
| const gsl_vector_float * X, |
| float beta, |
| gsl_vector_float * Y); |
| |
| int gsl_blas_sger (float alpha, |
| const gsl_vector_float * X, |
| const gsl_vector_float * Y, |
| gsl_matrix_float * A); |
| |
| int gsl_blas_ssyr (CBLAS_UPLO_t Uplo, |
| float alpha, |
| const gsl_vector_float * X, |
| gsl_matrix_float * A); |
| |
| int gsl_blas_ssyr2 (CBLAS_UPLO_t Uplo, |
| float alpha, |
| const gsl_vector_float * X, |
| const gsl_vector_float * Y, |
| gsl_matrix_float * A); |
| |
| int gsl_blas_dsymv (CBLAS_UPLO_t Uplo, |
| double alpha, |
| const gsl_matrix * A, |
| const gsl_vector * X, |
| double beta, |
| gsl_vector * Y); |
| int gsl_blas_dger (double alpha, |
| const gsl_vector * X, |
| const gsl_vector * Y, |
| gsl_matrix * A); |
| |
| int gsl_blas_dsyr (CBLAS_UPLO_t Uplo, |
| double alpha, |
| const gsl_vector * X, |
| gsl_matrix * A); |
| |
| int gsl_blas_dsyr2 (CBLAS_UPLO_t Uplo, |
| double alpha, |
| const gsl_vector * X, |
| const gsl_vector * Y, |
| gsl_matrix * A); |
| |
| /* |
| * Routines with C and Z prefixes only |
| */ |
| |
| int gsl_blas_chemv (CBLAS_UPLO_t Uplo, |
| const gsl_complex_float alpha, |
| const gsl_matrix_complex_float * A, |
| const gsl_vector_complex_float * X, |
| const gsl_complex_float beta, |
| gsl_vector_complex_float * Y); |
| |
| int gsl_blas_cgeru (const gsl_complex_float alpha, |
| const gsl_vector_complex_float * X, |
| const gsl_vector_complex_float * Y, |
| gsl_matrix_complex_float * A); |
| |
| int gsl_blas_cgerc (const gsl_complex_float alpha, |
| const gsl_vector_complex_float * X, |
| const gsl_vector_complex_float * Y, |
| gsl_matrix_complex_float * A); |
| |
| int gsl_blas_cher (CBLAS_UPLO_t Uplo, |
| float alpha, |
| const gsl_vector_complex_float * X, |
| gsl_matrix_complex_float * A); |
| |
| int gsl_blas_cher2 (CBLAS_UPLO_t Uplo, |
| const gsl_complex_float alpha, |
| const gsl_vector_complex_float * X, |
| const gsl_vector_complex_float * Y, |
| gsl_matrix_complex_float * A); |
| |
| int gsl_blas_zhemv (CBLAS_UPLO_t Uplo, |
| const gsl_complex alpha, |
| const gsl_matrix_complex * A, |
| const gsl_vector_complex * X, |
| const gsl_complex beta, |
| gsl_vector_complex * Y); |
| |
| int gsl_blas_zgeru (const gsl_complex alpha, |
| const gsl_vector_complex * X, |
| const gsl_vector_complex * Y, |
| gsl_matrix_complex * A); |
| |
| int gsl_blas_zgerc (const gsl_complex alpha, |
| const gsl_vector_complex * X, |
| const gsl_vector_complex * Y, |
| gsl_matrix_complex * A); |
| |
| int gsl_blas_zher (CBLAS_UPLO_t Uplo, |
| double alpha, |
| const gsl_vector_complex * X, |
| gsl_matrix_complex * A); |
| |
| int gsl_blas_zher2 (CBLAS_UPLO_t Uplo, |
| const gsl_complex alpha, |
| const gsl_vector_complex * X, |
| const gsl_vector_complex * Y, |
| gsl_matrix_complex * A); |
| |
| /* |
| * =========================================================================== |
| * Prototypes for level 3 BLAS |
| * =========================================================================== |
| */ |
| |
| /* |
| * Routines with standard 4 prefixes (S, D, C, Z) |
| */ |
| int gsl_blas_sgemm (CBLAS_TRANSPOSE_t TransA, |
| CBLAS_TRANSPOSE_t TransB, |
| float alpha, |
| const gsl_matrix_float * A, |
| const gsl_matrix_float * B, |
| float beta, |
| gsl_matrix_float * C); |
| |
| int gsl_blas_ssymm (CBLAS_SIDE_t Side, CBLAS_UPLO_t Uplo, |
| float alpha, |
| const gsl_matrix_float * A, |
| const gsl_matrix_float * B, |
| float beta, |
| gsl_matrix_float * C); |
| |
| int gsl_blas_ssyrk (CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t Trans, |
| float alpha, |
| const gsl_matrix_float * A, |
| float beta, |
| gsl_matrix_float * C); |
| |
| int gsl_blas_ssyr2k (CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t Trans, |
| float alpha, |
| const gsl_matrix_float * A, |
| const gsl_matrix_float * B, |
| float beta, |
| gsl_matrix_float * C); |
| |
| int gsl_blas_strmm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, |
| CBLAS_DIAG_t Diag, |
| float alpha, |
| const gsl_matrix_float * A, |
| gsl_matrix_float * B); |
| |
| int gsl_blas_strsm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, |
| CBLAS_DIAG_t Diag, |
| float alpha, |
| const gsl_matrix_float * A, |
| gsl_matrix_float * B); |
| |
| int gsl_blas_dgemm (CBLAS_TRANSPOSE_t TransA, |
| CBLAS_TRANSPOSE_t TransB, |
| double alpha, |
| const gsl_matrix * A, |
| const gsl_matrix * B, |
| double beta, |
| gsl_matrix * C); |
| |
| int gsl_blas_dsymm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, |
| double alpha, |
| const gsl_matrix * A, |
| const gsl_matrix * B, |
| double beta, |
| gsl_matrix * C); |
| |
| int gsl_blas_dsyrk (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t Trans, |
| double alpha, |
| const gsl_matrix * A, |
| double beta, |
| gsl_matrix * C); |
| |
| int gsl_blas_dsyr2k (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t Trans, |
| double alpha, |
| const gsl_matrix * A, |
| const gsl_matrix * B, |
| double beta, |
| gsl_matrix * C); |
| |
| int gsl_blas_dtrmm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, |
| CBLAS_DIAG_t Diag, |
| double alpha, |
| const gsl_matrix * A, |
| gsl_matrix * B); |
| |
| int gsl_blas_dtrsm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, |
| CBLAS_DIAG_t Diag, |
| double alpha, |
| const gsl_matrix * A, |
| gsl_matrix * B); |
| |
| int gsl_blas_cgemm (CBLAS_TRANSPOSE_t TransA, |
| CBLAS_TRANSPOSE_t TransB, |
| const gsl_complex_float alpha, |
| const gsl_matrix_complex_float * A, |
| const gsl_matrix_complex_float * B, |
| const gsl_complex_float beta, |
| gsl_matrix_complex_float * C); |
| |
| int gsl_blas_csymm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, |
| const gsl_complex_float alpha, |
| const gsl_matrix_complex_float * A, |
| const gsl_matrix_complex_float * B, |
| const gsl_complex_float beta, |
| gsl_matrix_complex_float * C); |
| |
| int gsl_blas_csyrk (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t Trans, |
| const gsl_complex_float alpha, |
| const gsl_matrix_complex_float * A, |
| const gsl_complex_float beta, |
| gsl_matrix_complex_float * C); |
| |
| int gsl_blas_csyr2k (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t Trans, |
| const gsl_complex_float alpha, |
| const gsl_matrix_complex_float * A, |
| const gsl_matrix_complex_float * B, |
| const gsl_complex_float beta, |
| gsl_matrix_complex_float * C); |
| |
| int gsl_blas_ctrmm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, |
| CBLAS_DIAG_t Diag, |
| const gsl_complex_float alpha, |
| const gsl_matrix_complex_float * A, |
| gsl_matrix_complex_float * B); |
| |
| int gsl_blas_ctrsm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, |
| CBLAS_DIAG_t Diag, |
| const gsl_complex_float alpha, |
| const gsl_matrix_complex_float * A, |
| gsl_matrix_complex_float * B); |
| |
| int gsl_blas_zgemm (CBLAS_TRANSPOSE_t TransA, |
| CBLAS_TRANSPOSE_t TransB, |
| const gsl_complex alpha, |
| const gsl_matrix_complex * A, |
| const gsl_matrix_complex * B, |
| const gsl_complex beta, |
| gsl_matrix_complex * C); |
| |
| int gsl_blas_zsymm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, |
| const gsl_complex alpha, |
| const gsl_matrix_complex * A, |
| const gsl_matrix_complex * B, |
| const gsl_complex beta, |
| gsl_matrix_complex * C); |
| |
| int gsl_blas_zsyrk (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t Trans, |
| const gsl_complex alpha, |
| const gsl_matrix_complex * A, |
| const gsl_complex beta, |
| gsl_matrix_complex * C); |
| |
| int gsl_blas_zsyr2k (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t Trans, |
| const gsl_complex alpha, |
| const gsl_matrix_complex * A, |
| const gsl_matrix_complex * B, |
| const gsl_complex beta, |
| gsl_matrix_complex *C); |
| |
| int gsl_blas_ztrmm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, |
| CBLAS_DIAG_t Diag, |
| const gsl_complex alpha, |
| const gsl_matrix_complex * A, |
| gsl_matrix_complex * B); |
| |
| int gsl_blas_ztrsm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, CBLAS_TRANSPOSE_t TransA, |
| CBLAS_DIAG_t Diag, |
| const gsl_complex alpha, |
| const gsl_matrix_complex * A, |
| gsl_matrix_complex * B); |
| |
| /* |
| * Routines with prefixes C and Z only |
| */ |
| int gsl_blas_chemm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, |
| const gsl_complex_float alpha, |
| const gsl_matrix_complex_float * A, |
| const gsl_matrix_complex_float * B, |
| const gsl_complex_float beta, |
| gsl_matrix_complex_float * C); |
| |
| int gsl_blas_cherk (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t Trans, |
| float alpha, |
| const gsl_matrix_complex_float * A, |
| float beta, |
| gsl_matrix_complex_float * C); |
| |
| int gsl_blas_cher2k (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t Trans, |
| const gsl_complex_float alpha, |
| const gsl_matrix_complex_float * A, |
| const gsl_matrix_complex_float * B, |
| float beta, |
| gsl_matrix_complex_float * C); |
| |
| int gsl_blas_zhemm (CBLAS_SIDE_t Side, |
| CBLAS_UPLO_t Uplo, |
| const gsl_complex alpha, |
| const gsl_matrix_complex * A, |
| const gsl_matrix_complex * B, |
| const gsl_complex beta, |
| gsl_matrix_complex * C); |
| |
| int gsl_blas_zherk (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t Trans, |
| double alpha, |
| const gsl_matrix_complex * A, |
| double beta, |
| gsl_matrix_complex * C); |
| |
| int gsl_blas_zher2k (CBLAS_UPLO_t Uplo, |
| CBLAS_TRANSPOSE_t Trans, |
| const gsl_complex alpha, |
| const gsl_matrix_complex * A, |
| const gsl_matrix_complex * B, |
| double beta, |
| gsl_matrix_complex * C); |
| |
| |
| __END_DECLS |
| |
| #endif /* __GSL_BLAS_H__ */ |