Trang chủ Khám phá Trợ giúp
Đăng nhập
RONCC
/
Geografic-Information-System
Xem 1
Star 0
Fork 0
Các tập tin Các vấn đề 0 Yêu cầu khéo về 0 Wiki
Tree: c458350789
Branches Tags
Geografic-In...
/
include
/
grass
/
defs
/
gmath.h

gmath.h 8.5 KB
Lịch sử Raw

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183 
  1. #ifndef GRASS_GMATHDEFS_H
    2. 

  2. #define GRASS_GMATHDEFS_H
    3. 

  3. 

  4. /* dalloc.c */
    5. 

  5. double *G_alloc_vector(size_t);
    6. 

  6. double **G_alloc_matrix(int, int);
    7. 

  7. float  *G_alloc_fvector(size_t);
    8. 

  8. float  **G_alloc_fmatrix(int, int);
    9. 

  9. void G_free_vector(double *);
    10. 

  10. void G_free_matrix(double **);
    11. 

  11. void G_free_fvector(float *);
    12. 

  12. void G_free_fmatrix(float **);
    13. 

  13. 

  14. /* ialloc.c */
    15. 

  15. int *G_alloc_ivector(size_t);
    16. 

  16. int **G_alloc_imatrix(int, int);
    17. 

  17. void G_free_ivector(int *);
    18. 

  18. void G_free_imatrix(int **);
    19. 

  19. 

  20. /* fft.c */
    21. 

  21. extern int fft(int, double *[2], int, int, int);
    22. 

  22. extern int fft2(int, double (*)[2], int, int, int);
    23. 

  23. 

  24. /* gauss.c */
    25. 

  25. extern double G_math_rand_gauss(double);
    26. 

  26. 

  27. /* max_pow2.c */
    28. 

  28. extern long G_math_max_pow2 (long n);
    29. 

  29. extern long G_math_min_pow2 (long n);
    30. 

  30. 

  31. /* rand1.c */
    32. 

  32. extern void G_math_srand(int);
    33. 

  33. extern int G_math_srand_auto(void);
    34. 

  34. extern float G_math_rand(void);
    35. 

  35. 

  36. /* del2g.c */
    37. 

  37. extern int del2g(double *[2], int, double);
    38. 

  38. 

  39. /* getg.c */
    40. 

  40. extern int getg(double, double *[2], int);
    41. 

  41. 

  42. /* eigen_tools.c */
    43. 

  43. extern int G_math_egvorder(double *, double **, long);
    44. 

  44. 

  45. /* mult.c */
    46. 

  46. extern int G_math_complex_mult (double *v1[2], int size1, double *v2[2], int size2, double *v3[2], int size3);
    47. 

  47. 

  48. /* lu.c*/
    49. 

  49. extern int G_ludcmp(double **, int, int *, double *);
    50. 

  50. extern void G_lubksb(double **a, int n, int *indx, double b[]);
    51. 

  51. 

  52. /* findzc.c */
    53. 

  53. extern int G_math_findzc(double conv[], int size, double zc[], double thresh, int num_orients);
    54. 

  54. 

  55. 

  56. /* *************************************************************** */
    57. 

  57. /* ***** WRAPPER FOR CCMATH FUNCTIONS USED IN GRASS ************** */
    58. 

  58. /* *************************************************************** */
    59. 

  59. extern int G_math_solv(double **,double *,int);
    60. 

  60. extern int G_math_solvps(double **,double *,int);
    61. 

  61. extern void G_math_solvtd(double *,double *,double *,double *,int);
    62. 

  62. extern int G_math_solvru(double **,double *,int);
    63. 

  63. extern int G_math_minv(double **,int);
    64. 

  64. extern int G_math_psinv(double **,int);
    65. 

  65. extern int G_math_ruinv(double **,int);
    66. 

  66. extern void G_math_eigval(double **,double *,int);
    67. 

  67. extern void G_math_eigen(double **,double *,int);
    68. 

  68. extern double G_math_evmax(double **,double *,int);
    69. 

  69. extern int G_math_svdval(double *,double **,int,int);
    70. 

  70. extern int G_math_sv2val(double *,double **,int,int);
    71. 

  71. extern int G_math_svduv(double *,double **,double **, int,double **,int);
    72. 

  72. extern int G_math_sv2uv(double *,double **,double **,int,double **,int);
    73. 

  73. extern int G_math_svdu1v(double *,double **,int,double **,int);
    74. 

  74. 

  75. /* *************************************************************** */
    76. 

  76. /* *************** LINEARE EQUATION SYSTEM PART ****************** */
    77. 

  77. /* *************************************************************** */
    78. 

  78. 

  79. /* Sparse matrix and sparse vector functions
    80. 

  80.  * */
    81. 

  81. extern G_math_spvector *G_math_alloc_spvector(int );
    82. 

  82. extern G_math_spvector **G_math_alloc_spmatrix(int );
    83. 

  83. extern void G_math_free_spmatrix(G_math_spvector ** , int );
    84. 

  84. extern void G_math_free_spvector(G_math_spvector * );
    85. 

  85. extern int G_math_add_spvector(G_math_spvector **, G_math_spvector * , int );
    86. 

  86. extern G_math_spvector **G_math_A_to_Asp(double **, int, double);
    87. 

  87. extern double **G_math_Asp_to_A(G_math_spvector **, int);
    88. 

  88. extern double **G_math_Asp_to_sband_matrix(G_math_spvector **, int, int);
    89. 

  89. extern G_math_spvector **G_math_sband_matrix_to_Asp(double **, int, int, double);
    90. 

  90. extern void G_math_print_spmatrix(G_math_spvector **, int);
    91. 

  91. extern void G_math_Ax_sparse(G_math_spvector **, double *, double *, int );
    92. 

  92. 

  93. /*Symmetric band matrix handling */
    94. 

  94. extern double **G_math_matrix_to_sband_matrix(double **, int, int);
    95. 

  95. extern double **G_math_sband_matrix_to_matrix(double **, int, int);
    96. 

  96. extern void G_math_Ax_sband(double ** A, double *x, double *y, int rows, int bandwidth);
    97. 

  97. 

  98. /*linear equation solver, most of them are multithreaded with OpenMP*/
    99. 

  99. extern int G_math_solver_gauss(double **, double *, double *, int );
    100. 

  100. extern int G_math_solver_lu(double **, double *, double *, int );
    101. 

  101. extern int G_math_solver_cholesky(double **, double *, double *, int , int );
    102. 

  102. extern void G_math_solver_cholesky_sband(double **, double *, double *, int, int);
    103. 

  103. extern int G_math_solver_jacobi(double **, double *, double *, int , int , double , double );
    104. 

  104. extern int G_math_solver_gs(double **, double *, double *, int , int , double , double );
    105. 

  105. 

  106. extern int G_math_solver_pcg(double **, double *, double *, int , int , double , int );
    107. 

  107. extern int G_math_solver_cg(double **, double *, double *, int , int , double );
    108. 

  108. extern int G_math_solver_cg_sband(double **, double *, double *, int, int, int, double);
    109. 

  109. extern int G_math_solver_bicgstab(double **, double *, double *, int , int , double );
    110. 

  110. extern int G_math_solver_sparse_jacobi(G_math_spvector **, double *, double *, int , int , double , double );
    111. 

  111. extern int G_math_solver_sparse_gs(G_math_spvector **, double *, double *, int , int , double , double );
    112. 

  112. extern int G_math_solver_sparse_pcg(G_math_spvector **, double *, double *, int , int , double , int );
    113. 

  113. extern int G_math_solver_sparse_cg(G_math_spvector **, double *, double *, int , int , double );
    114. 

  114. extern int G_math_solver_sparse_bicgstab(G_math_spvector **, double *, double *, int , int , double );
    115. 

  115. 

  116. /* solver algoithms and helper functions*/
    117. 

  117. extern void G_math_gauss_elimination(double **, double *, int );
    118. 

  118. extern void G_math_lu_decomposition(double **, double *, int );
    119. 

  119. extern int G_math_cholesky_decomposition(double **, int , int );
    120. 

  120. extern void G_math_cholesky_sband_decomposition(double **, double **, int, int);
    121. 

  121. extern void G_math_backward_substitution(double **, double *, double *, int );
    122. 

  122. extern void G_math_forward_substitution(double **, double *, double *, int );
    123. 

  123. extern void G_math_cholesky_sband_substitution(double **, double *, double *, int, int);
    124. 

  124. 

  125. /*BLAS like level 1,2 and 3 functions*/
    126. 

  126. 

  127. /*level 1 vector - vector grass implementation with OpenMP thread support*/
    128. 

  128. extern void G_math_d_x_dot_y(double *, double *, double *, int );
    129. 

  129. extern void G_math_d_asum_norm(double *, double *, int );
    130. 

  130. extern void G_math_d_euclid_norm(double *, double *, int );
    131. 

  131. extern void G_math_d_max_norm(double *, double *, int );
    132. 

  132. extern void G_math_d_ax_by(double *, double *, double *, double , double , int );
    133. 

  133. extern void G_math_d_copy(double *, double *, int );
    134. 

  134. 

  135. extern void G_math_f_x_dot_y(float *, float *, float *, int );
    136. 

  136. extern void G_math_f_asum_norm(float *, float *, int );
    137. 

  137. extern void G_math_f_euclid_norm(float *, float *, int );
    138. 

  138. extern void G_math_f_max_norm(float *, float *, int );
    139. 

  139. extern void G_math_f_ax_by(float *, float *, float *, float , float , int );
    140. 

  140. extern void G_math_f_copy(float *, float *, int );
    141. 

  141. 

  142. extern void G_math_i_x_dot_y(int *, int *,  double *, int );
    143. 

  143. extern void G_math_i_asum_norm(int *,  double *, int );
    144. 

  144. extern void G_math_i_euclid_norm(int *,  double *,int );
    145. 

  145. extern void G_math_i_max_norm(int *,  int *, int );
    146. 

  146. extern void G_math_i_ax_by(int *, int *, int *, int , int , int );
    147. 

  147. extern void G_math_i_copy(int *, int *, int );
    148. 

  148. 

  149. /*ATLAS blas level 1 wrapper*/
    150. 

  150. extern double G_math_ddot(double *, double *, int );
    151. 

  151. extern float G_math_sdot(float *, float *, int );
    152. 

  152. extern float G_math_sdsdot(float *, float *, float , int );
    153. 

  153. extern double G_math_dnrm2(double *, int );
    154. 

  154. extern double G_math_dasum(double *, int );
    155. 

  155. extern double G_math_idamax(double *, int );
    156. 

  156. extern float  G_math_snrm2(float *, int );
    157. 

  157. extern float  G_math_sasum(float *, int );
    158. 

  158. extern float  G_math_isamax(float *, int );
    159. 

  159. extern void G_math_dscal(double *, double , int );
    160. 

  160. extern void G_math_sscal(float *, float , int );
    161. 

  161. extern void G_math_dcopy(double *, double *, int );
    162. 

  162. extern void G_math_scopy(float *, float *, int );
    163. 

  163. extern void G_math_daxpy(double *, double *, double , int );
    164. 

  164. extern void G_math_saxpy(float *, float *, float , int );
    165. 

  165. 

  166. /*level 2 matrix - vector grass implementation with OpenMP thread support*/
    167. 

  167. extern void G_math_d_Ax(double **, double *, double *, int , int );
    168. 

  168. extern void G_math_f_Ax(float **, float *, float *, int , int );
    169. 

  169. extern void G_math_d_x_dyad_y(double *, double *, double **, int, int );
    170. 

  170. extern void G_math_f_x_dyad_y(float *, float *, float **, int, int );
    171. 

  171. extern void G_math_d_aAx_by(double **, double *, double *, double , double , double *, int , int );
    172. 

  172. extern void G_math_f_aAx_by(float **, float *, float *, float , float , float *, int , int );
    173. 

  173. extern int G_math_d_A_T(double **A, int rows);
    174. 

  174. extern int G_math_f_A_T(float **A, int rows);
    175. 

  175. 

  176. /*level 3 matrix - matrix grass implementation with OpenMP thread support*/
    177. 

  177. extern void G_math_d_aA_B(double **, double **, double , double **, int , int );
    178. 

  178. extern void G_math_f_aA_B(float **, float **, float , float **, int , int );
    179. 

  179. extern void G_math_d_AB(double **, double **, double **, int , int , int );
    180. 

  180. extern void G_math_f_AB(float **,  float **,  float **,  int , int , int );
    181. 

  181. 

  182. #endif /* GRASS_GMATHDEFS_H */
    183. 

  183.