| Brian Silverman | 72890c2 | 2015-09-19 14:37:37 -0400 | [diff] [blame^] | 1 | //===================================================== |
| 2 | // File : blas_interface.hh |
| 3 | // Author : L. Plagne <laurent.plagne@edf.fr)> |
| 4 | // Copyright (C) EDF R&D, lun sep 30 14:23:28 CEST 2002 |
| 5 | //===================================================== |
| 6 | // |
| 7 | // This program is free software; you can redistribute it and/or |
| 8 | // modify it under the terms of the GNU General Public License |
| 9 | // as published by the Free Software Foundation; either version 2 |
| 10 | // of the License, or (at your option) any later version. |
| 11 | // |
| 12 | // This program is distributed in the hope that it will be useful, |
| 13 | // but WITHOUT ANY WARRANTY; without even the implied warranty of |
| 14 | // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
| 15 | // GNU General Public License for more details. |
| 16 | // You should have received a copy of the GNU General Public License |
| 17 | // along with this program; if not, write to the Free Software |
| 18 | // Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. |
| 19 | // |
| 20 | #ifndef blas_PRODUIT_MATRICE_VECTEUR_HH |
| 21 | #define blas_PRODUIT_MATRICE_VECTEUR_HH |
| 22 | |
| 23 | #include <c_interface_base.h> |
| 24 | #include <complex> |
| 25 | extern "C" |
| 26 | { |
| 27 | #include "blas.h" |
| 28 | |
| 29 | // Cholesky Factorization |
| 30 | // void spotrf_(const char* uplo, const int* n, float *a, const int* ld, int* info); |
| 31 | // void dpotrf_(const char* uplo, const int* n, double *a, const int* ld, int* info); |
| 32 | void ssytrd_(char *uplo, const int *n, float *a, const int *lda, float *d, float *e, float *tau, float *work, int *lwork, int *info ); |
| 33 | void dsytrd_(char *uplo, const int *n, double *a, const int *lda, double *d, double *e, double *tau, double *work, int *lwork, int *info ); |
| 34 | void sgehrd_( const int *n, int *ilo, int *ihi, float *a, const int *lda, float *tau, float *work, int *lwork, int *info ); |
| 35 | void dgehrd_( const int *n, int *ilo, int *ihi, double *a, const int *lda, double *tau, double *work, int *lwork, int *info ); |
| 36 | |
| 37 | // LU row pivoting |
| 38 | // void dgetrf_( int *m, int *n, double *a, int *lda, int *ipiv, int *info ); |
| 39 | // void sgetrf_(const int* m, const int* n, float *a, const int* ld, int* ipivot, int* info); |
| 40 | // LU full pivoting |
| 41 | void sgetc2_(const int* n, float *a, const int *lda, int *ipiv, int *jpiv, int*info ); |
| 42 | void dgetc2_(const int* n, double *a, const int *lda, int *ipiv, int *jpiv, int*info ); |
| 43 | #ifdef HAS_LAPACK |
| 44 | #endif |
| 45 | } |
| 46 | |
| 47 | #define MAKE_STRING2(S) #S |
| 48 | #define MAKE_STRING(S) MAKE_STRING2(S) |
| 49 | |
| 50 | #define CAT2(A,B) A##B |
| 51 | #define CAT(A,B) CAT2(A,B) |
| 52 | |
| 53 | |
| 54 | template<class real> class blas_interface; |
| 55 | |
| 56 | |
| 57 | static char notrans = 'N'; |
| 58 | static char trans = 'T'; |
| 59 | static char nonunit = 'N'; |
| 60 | static char lower = 'L'; |
| 61 | static char right = 'R'; |
| 62 | static char left = 'L'; |
| 63 | static int intone = 1; |
| 64 | |
| 65 | |
| 66 | |
| 67 | #define SCALAR float |
| 68 | #define SCALAR_PREFIX s |
| 69 | #include "blas_interface_impl.hh" |
| 70 | #undef SCALAR |
| 71 | #undef SCALAR_PREFIX |
| 72 | |
| 73 | |
| 74 | #define SCALAR double |
| 75 | #define SCALAR_PREFIX d |
| 76 | #include "blas_interface_impl.hh" |
| 77 | #undef SCALAR |
| 78 | #undef SCALAR_PREFIX |
| 79 | |
| 80 | #endif |
| 81 | |
| 82 | |
| 83 | |