PageRenderTime 58ms CodeModel.GetById 31ms RepoModel.GetById 1ms app.codeStats 0ms

/Core/Lapack.cpp

https://github.com/leppert/fire
C++ | 171 lines | 130 code | 40 blank | 1 comment | 0 complexity | 08dd335457cd4e187e0800366298489f MD5 | raw file
    

  1. #include "Lapack.hh"
    2. 

  2. 

  3. #define F77NAME(x) x##_
    4. 

  4. 

  5. // Type dependend Lapack definitions:
    6. 

  6. 

  7. namespace Lapack{
    8. 

  8. 

  9.   extern "C" {
    10. 

  10. 

  11.     void F77NAME(dgelsd)(long* m, long* n, long* nrhs, double* a, long* lda, double* b, long* ldb, 
    12. 

  12.                          double* s, double* rcond, long* rank, double* work, long* lwork, long* iwork,
    13. 

  13.                          long* info);
    14. 

  14. 

  15. 			 
    16. 

  16.   
    17. 

  17.     void F77NAME(dggevx)(char* balanc, char* jobvl, char* jobvr, char* sense,
    18. 

  18.                          long* n, double* a, long* lda, double* b, long* ldb,
    19. 

  19.                          double* alphar, double* alphai, double* beta, double* vl, long* ldvl, double* vr, long* ldvr,
    20. 

  20.                          long* ilo, long* ihi, double* lscale, double* rscale, double* abnrm, double* bbnrm, 
    21. 

  21.                          double* rconde, double* rcondv, double* work, long* lwork, long* iwork, bool* bwork, long* info);
    22. 

  22.   
    23. 

  23.   
    24. 

  24.     void F77NAME(sggev)(char *jobvl,char *jobvr,int *n,double *a,int *lda,double *b, int *ld,
    25. 

  25.                         double *alphar,double *alphai,double *beta,double *vl,
    26. 

  26.                         int *ldvl,double *vr,int *ldvr,double *work,int *lwork,
    27. 

  27.                         int *info) ;
    28. 

  28.     void F77NAME(dggev)(char *jobvl,char *jobvr,int *n,double *a,int *lda,double *b, int *ld,
    29. 

  29.                         double *alphar,double *alphai,double *beta,double *vl,
    30. 

  30.                         int *ldvl,double *vr,int *ldvr,double *work,int *lwork,
    31. 

  31.                         int *info) ;
    32. 

  32.   
    33. 

  33.   
    34. 

  34.     void F77NAME(dsygvd)(long* itype, char* jobz, char* uplo, long* n, double* a, long* lda, double* b, long* ldb,
    35. 

  35.                          double* w, double* work, long* lwork, long* iwork, long* liwork, long* info);
    36. 

  36. 

  37.     void F77NAME(dgels)(char* trans, long* m, long* n, long* nrhs, double* A, long* lda, 
    38. 

  38.                         double* B, long* ldb, double* work, long* lwork, long* info);
    39. 

  39. 

  40.     void F77NAME(sgels)(char* trans, long* m, long* n, long* nrhs, float* A, long* lda, 
    41. 

  41.                         float* B, long* ldb, float* work, long* lwork, long* info);
    42. 

  42. 

  43.     void F77NAME(dgelss)(long* m, long* n, long* nrhs, double* A, long* lda, double* B, long* ldb, 
    44. 

  44.                          double* s, double* rcond, long* rank, double* work, long* lwork, long* info);
    45. 

  45. 

  46.     void F77NAME(sgelss)(long* m, long* n, long* nrhs, float* A, long* lda, float* B, long* ldb, 
    47. 

  47.                          float* s, float* rcond, long* rank, float* work, long* lwork, long* info);
    48. 

  48. 

  49.     void F77NAME(dgebrd)(long* m, long *n, double* a, long* lda, double* d, double* e, 
    50. 

  50.                          double* tauq,double* taup, double* work, long* lwork,long* info);
    51. 

  51. 

  52.     void F77NAME(dbdsdc) (char* uplo, char* compq, long* n, double* d, double* e, double* u, 
    53. 

  53.                           long* ldu, double* vt, long* ldvt, double* q, long* iq, double* work, 
    54. 

  54.                           long* iwork, long* info );
    55. 

  55. 

  56.     void F77NAME(dgesdd) (char* jobz, long* m, long* n, double* a, long* lda, double* s, 
    57. 

  57.                           double* u, long *ldu, double* vt, long* ldvt, double* work, 
    58. 

  58.                           long* lwork, long* iwork, long* info );
    59. 

  59. 

  60.   }
    61. 

  61. 

  62.   void dgelsd(long* m, long* n, long* nrhs, double* a, long* lda, double* b, long* ldb, 
    63. 

  63.               double* s, double* rcond, long* rank, double* work, long* lwork, long* iwork, long* info){
    64. 

  64. #ifdef HAVE_LAPACK
    65. 

  65.     F77NAME(dgelsd)( m, n, nrhs, a, lda, b, ldb, s, rcond, rank, work,lwork, iwork, info);
    66. 

  66. #endif
    67. 

  67.   }
    68. 

  68.     
    69. 

  69. 

  70. 

  71.   void dggevx(char* balanc, char* jobvl, char* jobvr, char* sense,
    72. 

  72.               long* n, double* a, long* lda, double* b, long* ldb,
    73. 

  73.               double* alphar, double* alphai, double* beta, double* vl, long* ldvl, double* vr, long* ldvr,
    74. 

  74.               long* ilo, long* ihi, double* lscale, double* rscale, double* abnrm, double* bbnrm, 
    75. 

  75.               double* rconde, double* rcondv, double* work, long* lwork, long* iwork, bool* bwork, long* info)
    76. 

  76.   {
    77. 

  77. #ifdef HAVE_LAPACK
    78. 

  78.     F77NAME(dggevx)(balanc, jobvl, jobvr, sense, n, a, lda, b, ldb,
    79. 

  79.                     alphar, alphai, beta, vl, ldvl, vr, ldvr,
    80. 

  80.                     ilo, ihi, lscale, rscale, abnrm, bbnrm, 
    81. 

  81.                     rconde, rcondv, work, lwork, iwork, bwork, info);
    82. 

  82. #endif
    83. 

  83.   }
    84. 

  84. 

  85. 

  86.   void sggev(char *jobvl,char *jobvr,int *n,double *a,int *lda,double *b, int *ld,
    87. 

  87.              double *alphar,double *alphai,double *beta,double *vl,
    88. 

  88.              int *ldvl,double *vr,int *ldvr,double *work,int *lwork,
    89. 

  89.              int *info) {
    90. 

  90. #ifdef HAVE_LAPACK
    91. 

  91.     F77NAME(sggev)(jobvl,jobvr,n,a,lda,b,ld,alphar,alphai,beta,vl,ldvl,vr,ldvr,work,lwork,info);
    92. 

  92. #endif
    93. 

  93.   }
    94. 

  94.   void dggev(char *jobvl,char *jobvr,int *n,double *a,int *lda,double *b, int *ld,
    95. 

  95.              double *alphar,double *alphai,double *beta,double *vl,
    96. 

  96.              int *ldvl,double *vr,int *ldvr,double *work,int *lwork,
    97. 

  97.              int *info) {
    98. 

  98. #ifdef HAVE_LAPACK
    99. 

  99.     F77NAME(dggev)(jobvl,jobvr,n,a,lda,b,ld,alphar,alphai,beta,vl,ldvl,vr,ldvr,work,lwork,info);
    100. 

  100. #endif
    101. 

  101.   }
    102. 

  102. 

  103. 

  104.   void dsygvd(long* itype, char* jobz, char* uplo, long* n, double* a, long* lda, double* b, long* ldb,
    105. 

  105.               double* w, double* work, long* lwork, long* iwork, long* liwork, long* info)
    106. 

  106.   {
    107. 

  107. #ifdef HAVE_LAPACK
    108. 

  108.     F77NAME(dsygvd)(itype, jobz, uplo, n, a, lda, b, ldb,
    109. 

  109.                     w, work, lwork, iwork, liwork, info);
    110. 

  110. #endif
    111. 

  111.   }
    112. 

  112. 

  113. 

  114.   void gels(char* trans, long* m, long* n, long* nrhs, double* A, long* lda, 
    115. 

  115.             double* B, long* ldb, double* work, long* lwork, long* info)
    116. 

  116.   {
    117. 

  117. #ifdef HAVE_LAPACK
    118. 

  118.     F77NAME(dgels)(trans, m, n, nrhs, A, lda, B, ldb, work, lwork, info);
    119. 

  119. #endif
    120. 

  120.   }
    121. 

  121. 

  122.   void gels(char* trans, long* m, long* n, long* nrhs, float* A, long* lda, 
    123. 

  123.             float* B, long* ldb, float* work, long* lwork, long* info)
    124. 

  124.   {
    125. 

  125. #ifdef HAVE_LAPACK
    126. 

  126.     F77NAME(sgels)(trans, m, n, nrhs, A, lda, B, ldb, work, lwork, info);
    127. 

  127. #endif
    128. 

  128.   }
    129. 

  129. 

  130.   void gelss(long* m, long* n, long* nrhs, double* A, long* lda, double* B, long* ldb, 
    131. 

  131.              double* s, double* rcond, long* rank, double* work, long* lwork, long* info)
    132. 

  132.   {
    133. 

  133. #ifdef HAVE_LAPACK
    134. 

  134.     F77NAME(dgelss)(m, n, nrhs, A, lda, B, ldb, s, rcond, rank, work, lwork, info);
    135. 

  135. #endif
    136. 

  136.   }
    137. 

  137. 

  138.   void gelss(long* m, long* n, long* nrhs, float* A, long* lda, float* B, long* ldb, 
    139. 

  139.              float* s, float* rcond, long* rank, float* work, long* lwork, long* info)
    140. 

  140.   {
    141. 

  141. #ifdef HAVE_LAPACK
    142. 

  142.     F77NAME(sgelss)(m, n, nrhs, A, lda, B, ldb, s, rcond, rank, work, lwork, info);
    143. 

  143. #endif
    144. 

  144.   }  
    145. 

  145. 

  146. 

  147.   void dgebrd(long* m, long *n, double* a, long* lda, double* d, double* e, double* tauq, 
    148. 

  148.               double* taup, double* work, long* lwork,long* info){
    149. 

  149. 

  150. #ifdef HAVE_LAPACK
    151. 

  151.     F77NAME(dgebrd)( m, n, a,  lda, d, e, tauq, taup, work, lwork, info);
    152. 

  152. #endif
    153. 

  153.   }
    154. 

  154. 

  155.   void dbdsdc ( char* uplo, char* compq, long* n, double* d, double* e, double* u, long* ldu, 
    156. 

  156.                 double* vt, long* ldvt, double* q, long* iq, double* work, long* iwork, long* info ){
    157. 

  157. #ifdef HAVE_LAPACK
    158. 

  158.     F77NAME(dbdsdc) (uplo, compq, n, d, e, u, ldu, vt, ldvt, q, iq, work, iwork, info );
    159. 

  159. #endif
    160. 

  160.   }
    161. 

  161. 

  162.   void dgesdd( char* jobz, long* m, long* n, double* a, long* lda, double* s, double* u, 
    163. 

  163.                long *ldu, double* vt, long* ldvt, double* work, long* lwork, long* iwork, 
    164. 

  164.                long* info ){
    165. 

  165. 

  166. #ifdef HAVE_LAPACK
    167. 

  167.     F77NAME(dgesdd)(jobz,  m,  n,  a,  lda,  s,  u, ldu,  vt,  ldvt,  work,  lwork,  iwork, info );
    168. 

  168. #endif
    169. 

  169.   }
    170. 

  170. 

  171. } //namespace
    172. 

  172.