Página inicial Explorar Ajuda
Registrar Entrar
lij
/
PS_Method_CCS_test
1
0
Fork 0
Arquivos Issues 0 Pull Requests 0 Wiki
Tree: 6f9ff47688
Branches Tags
PS_Method_CC...
/
libarmadillo
/
include
/
armadillo_bits
/
def_arpack.hpp

def_arpack.hpp 10 KB
Histórico Raw

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107 
  1. // Copyright 2008-2016 Conrad Sanderson (http://conradsanderson.id.au)
    2. 

  2. // Copyright 2008-2016 National ICT Australia (NICTA)
    3. 

  3. // 
    4. 

  4. // Licensed under the Apache License, Version 2.0 (the "License");
    5. 

  5. // you may not use this file except in compliance with the License.
    6. 

  6. // You may obtain a copy of the License at
    7. 

  7. // http://www.apache.org/licenses/LICENSE-2.0
    8. 

  8. // 
    9. 

  9. // Unless required by applicable law or agreed to in writing, software
    10. 

  10. // distributed under the License is distributed on an "AS IS" BASIS,
    11. 

  11. // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
    12. 

  12. // See the License for the specific language governing permissions and
    13. 

  13. // limitations under the License.
    14. 

  14. // ------------------------------------------------------------------------
    15. 

  15. 

  16. 

  17. #ifdef ARMA_USE_ARPACK
    18. 

  18. 

  19. // I'm not sure this is necessary.
    20. 

  20. #if !defined(ARMA_BLAS_CAPITALS)
    21. 

  21.   
    22. 

  22.   #define arma_snaupd snaupd
    23. 

  23.   #define arma_dnaupd dnaupd
    24. 

  24.   #define arma_cnaupd cnaupd
    25. 

  25.   #define arma_znaupd znaupd
    26. 

  26.   
    27. 

  27.   #define arma_sneupd sneupd
    28. 

  28.   #define arma_dneupd dneupd
    29. 

  29.   #define arma_cneupd cneupd
    30. 

  30.   #define arma_zneupd zneupd
    31. 

  31.   
    32. 

  32.   #define arma_ssaupd ssaupd
    33. 

  33.   #define arma_dsaupd dsaupd
    34. 

  34.   
    35. 

  35.   #define arma_sseupd sseupd
    36. 

  36.   #define arma_dseupd dseupd
    37. 

  37. 

  38. #else
    39. 

  39.   
    40. 

  40.   #define arma_snaupd SNAUPD
    41. 

  41.   #define arma_dnaupd DNAUPD
    42. 

  42.   #define arma_cnaupd CNAUPD
    43. 

  43.   #define arma_znaupd ZNAUPD
    44. 

  44.   
    45. 

  45.   #define arma_sneupd SNEUPD
    46. 

  46.   #define arma_dneupd DNEUPD
    47. 

  47.   #define arma_cneupd CNEUPD
    48. 

  48.   #define arma_zneupd ZNEUPD
    49. 

  49.   
    50. 

  50.   #define arma_ssaupd SSAUPD
    51. 

  51.   #define arma_dsaupd DSAUPD
    52. 

  52.   
    53. 

  53.   #define arma_sseupd SSEUPD
    54. 

  54.   #define arma_dseupd DSEUPD
    55. 

  55.   
    56. 

  56. #endif
    57. 

  57. 

  58. extern "C"
    59. 

  59. {
    60. 

  60. #if defined(ARMA_USE_FORTRAN_HIDDEN_ARGS)
    61. 

  61.   
    62. 

  62.   // eigendecomposition of non-symmetric positive semi-definite matrices
    63. 

  63.   void arma_fortran(arma_snaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,  float* resid, blas_int* ncv,  float* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,  float* workd,  float* workl, blas_int* lworkl, blas_int* info, blas_len bmat_len, blas_len which_len);
    64. 

  64.   void arma_fortran(arma_dnaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol, double* resid, blas_int* ncv, double* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr, double* workd, double* workl, blas_int* lworkl, blas_int* info, blas_len bmat_len, blas_len which_len);
    65. 

  65.   void arma_fortran(arma_cnaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,   void* resid, blas_int* ncv,   void* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,   void* workd,   void* workl, blas_int* lworkl,  float* rwork, blas_int* info, blas_len bmat_len, blas_len which_len);
    66. 

  66.   void arma_fortran(arma_znaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol,   void* resid, blas_int* ncv,   void* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,   void* workd,   void* workl, blas_int* lworkl, double* rwork, blas_int* info, blas_len bmat_len, blas_len which_len);
    67. 

  67.   
    68. 

  68.   // recovery of eigenvectors after naupd(); uses blas_int for LOGICAL types
    69. 

  69.   void arma_fortran(arma_sneupd)(blas_int* rvec, char* howmny, blas_int* select,  float* dr,  float* di,  float* z, blas_int* ldz,  float* sigmar,  float* sigmai,  float* workev, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,  float* resid, blas_int* ncv,  float* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,  float* workd,  float* workl, blas_int* lworkl, blas_int* info, blas_len howmny_len, blas_len bmat_len, blas_len which_len);
    70. 

  70.   void arma_fortran(arma_dneupd)(blas_int* rvec, char* howmny, blas_int* select, double* dr, double* di, double* z, blas_int* ldz, double* sigmar, double* sigmai, double* workev, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol, double* resid, blas_int* ncv, double* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr, double* workd, double* workl, blas_int* lworkl, blas_int* info, blas_len howmny_len, blas_len bmat_len, blas_len which_len);
    71. 

  71.   void arma_fortran(arma_cneupd)(blas_int* rvec, char* howmny, blas_int* select, void* d,   void* z, blas_int* ldz,   void* sigma,   void* workev, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,   void* resid, blas_int* ncv,   void* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,   void* workd, void* workl, blas_int* lworkl,  float* rwork, blas_int* info, blas_len howmny_len, blas_len bmat_len, blas_len which_len);
    72. 

  72.   void arma_fortran(arma_zneupd)(blas_int* rvec, char* howmny, blas_int* select, void* d,   void* z, blas_int* ldz,   void* sigma,   void* workev, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol,   void* resid, blas_int* ncv,   void* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,   void* workd, void* workl, blas_int* lworkl, double* rwork, blas_int* info, blas_len howmny_len, blas_len bmat_len, blas_len which_len);
    73. 

  73.   
    74. 

  74.   // eigendecomposition of symmetric positive semi-definite matrices
    75. 

  75.   void arma_fortran(arma_ssaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,  float* resid, blas_int* ncv,  float* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,  float* workd,  float* workl, blas_int* lworkl, blas_int* info, blas_len bmat_len, blas_len which_len);
    76. 

  76.   void arma_fortran(arma_dsaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol, double* resid, blas_int* ncv, double* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr, double* workd, double* workl, blas_int* lworkl, blas_int* info, blas_len bmat_len, blas_len which_len);
    77. 

  77.   
    78. 

  78.   // recovery of eigenvectors after saupd(); uses blas_int for LOGICAL types
    79. 

  79.   void arma_fortran(arma_sseupd)(blas_int* rvec, char* howmny, blas_int* select,  float* d,  float* z, blas_int* ldz,  float* sigma, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,  float* resid, blas_int* ncv,  float* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,  float* workd,  float* workl, blas_int* lworkl, blas_int* info, blas_len howmny_len, blas_len bmat_len, blas_len which_len);
    80. 

  80.   void arma_fortran(arma_dseupd)(blas_int* rvec, char* howmny, blas_int* select, double* d, double* z, blas_int* ldz, double* sigma, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol, double* resid, blas_int* ncv, double* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr, double* workd, double* workl, blas_int* lworkl, blas_int* info, blas_len howmny_len, blas_len bmat_len, blas_len which_len);
    81. 

  81.   
    82. 

  82. #else
    83. 

  83.   
    84. 

  84.   // eigendecomposition of non-symmetric positive semi-definite matrices
    85. 

  85.   void arma_fortran(arma_snaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,  float* resid, blas_int* ncv,  float* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,  float* workd,  float* workl, blas_int* lworkl, blas_int* info);
    86. 

  86.   void arma_fortran(arma_dnaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol, double* resid, blas_int* ncv, double* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr, double* workd, double* workl, blas_int* lworkl, blas_int* info);
    87. 

  87.   void arma_fortran(arma_cnaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,   void* resid, blas_int* ncv,   void* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,   void* workd,   void* workl, blas_int* lworkl,  float* rwork, blas_int* info);
    88. 

  88.   void arma_fortran(arma_znaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol,   void* resid, blas_int* ncv,   void* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,   void* workd,   void* workl, blas_int* lworkl, double* rwork, blas_int* info);
    89. 

  89.   
    90. 

  90.   // recovery of eigenvectors after naupd(); uses blas_int for LOGICAL types
    91. 

  91.   void arma_fortran(arma_sneupd)(blas_int* rvec, char* howmny, blas_int* select,  float* dr,  float* di,  float* z, blas_int* ldz,  float* sigmar,  float* sigmai,  float* workev, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,  float* resid, blas_int* ncv,  float* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,  float* workd,  float* workl, blas_int* lworkl, blas_int* info);
    92. 

  92.   void arma_fortran(arma_dneupd)(blas_int* rvec, char* howmny, blas_int* select, double* dr, double* di, double* z, blas_int* ldz, double* sigmar, double* sigmai, double* workev, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol, double* resid, blas_int* ncv, double* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr, double* workd, double* workl, blas_int* lworkl, blas_int* info);
    93. 

  93.   void arma_fortran(arma_cneupd)(blas_int* rvec, char* howmny, blas_int* select, void* d,   void* z, blas_int* ldz,   void* sigma,   void* workev, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,   void* resid, blas_int* ncv,   void* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,   void* workd, void* workl, blas_int* lworkl,  float* rwork, blas_int* info);
    94. 

  94.   void arma_fortran(arma_zneupd)(blas_int* rvec, char* howmny, blas_int* select, void* d,   void* z, blas_int* ldz,   void* sigma,   void* workev, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol,   void* resid, blas_int* ncv,   void* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,   void* workd, void* workl, blas_int* lworkl, double* rwork, blas_int* info);
    95. 

  95.   
    96. 

  96.   // eigendecomposition of symmetric positive semi-definite matrices
    97. 

  97.   void arma_fortran(arma_ssaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,  float* resid, blas_int* ncv,  float* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,  float* workd,  float* workl, blas_int* lworkl, blas_int* info);
    98. 

  98.   void arma_fortran(arma_dsaupd)(blas_int* ido, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol, double* resid, blas_int* ncv, double* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr, double* workd, double* workl, blas_int* lworkl, blas_int* info);
    99. 

  99.   
    100. 

  100.   // recovery of eigenvectors after saupd(); uses blas_int for LOGICAL types
    101. 

  101.   void arma_fortran(arma_sseupd)(blas_int* rvec, char* howmny, blas_int* select,  float* d,  float* z, blas_int* ldz,  float* sigma, char* bmat, blas_int* n, char* which, blas_int* nev,  float* tol,  float* resid, blas_int* ncv,  float* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr,  float* workd,  float* workl, blas_int* lworkl, blas_int* info);
    102. 

  102.   void arma_fortran(arma_dseupd)(blas_int* rvec, char* howmny, blas_int* select, double* d, double* z, blas_int* ldz, double* sigma, char* bmat, blas_int* n, char* which, blas_int* nev, double* tol, double* resid, blas_int* ncv, double* v, blas_int* ldv, blas_int* iparam, blas_int* ipntr, double* workd, double* workl, blas_int* lworkl, blas_int* info);
    103. 

  103.   
    104. 

  104. #endif
    105. 

  105. }
    106. 

  106. 

  107. #endif
    108.