开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 图形图像处理 > 查看源码
utils.cc
资源名称：contourletHMT.rar　[点击查看]
上传用户：l56789
上传日期：2022-02-25
资源大小：2422k
文件大小：2k
源码类别：

图形图像处理

开发平台：

Matlab

utils.cc：源码内容
  1. //-----------------------------------------------------------------------------
  2. // utils.cc
  3. //-----------------------------------------------------------------------------
  5. #include <stdio.h>
  6. #include <string.h>
  7. #include <math.h>
  8. #include "utils.hh"
  11. /*---------------------------------------------------------------------------*/
  12. int ipow(int n, int j)
  13. {
  14.     register int i, res;
  15.   
  16.     res = 1;
  17.     for (i = 0; i < j; i++)
  18. res *= n;
  20.     return res;
  21. }
  24. /*-----------------------------------------------------------------------*/
  25. // For random generator (Numerical Recipe book)
  26. #define IA 16807
  27. #define IM 2147483647
  28. #define AM (1.0/IM)
  29. #define IQ 127773
  30. #define IR 2836
  31. #define NTAB 32
  32. #define NDIV (1+(IM-1)/NTAB)
  33. #define EPS 1.2e-7
  34. #define RNMX (1.0-EPS)
  36. double ran1(int& idum)
  37. {
  38.     int j;
  39.     int k;
  40.     static int iy=0;
  41.     static int iv[NTAB];
  42.     double temp;
  44.     if (idum <= 0 || !iy) {
  45. if (-(idum) < 1) idum=1;
  46. else idum = -(idum);
  47. for (j=NTAB+7;j>=0;j--) {
  48.     k=(idum)/IQ;
  49.     idum=IA*(idum-k*IQ)-IR*k;
  50.     if (idum < 0) idum += IM;
  51.     if (j < NTAB) iv[j] = idum;
  52. }
  53. iy=iv[0];
  54.     }
  55.     k=(idum)/IQ;
  56.     idum=IA*(idum-k*IQ)-IR*k;
  57.     if (idum < 0) idum += IM;
  58.     j=iy/NDIV;
  59.     iy=iv[j];
  60.     iv[j] = idum;
  61.     if ((temp=AM*iy) > RNMX) return RNMX;
  62.     else return temp;
  63. }
  66. /*-----------------------------------------------------------------------*/
  67. double rangas(int& idum)
  68. {
  69.     static int iset=0;
  70.     static double gset;
  71.     double fac,rsq,v1,v2;
  73.     if  (iset == 0) {
  74. do {
  75.     v1=2.0*ran1(idum)-1.0;
  76.     v2=2.0*ran1(idum)-1.0;
  77.     rsq=v1*v1+v2*v2;
  78. } while (rsq >= 1.0 || rsq == 0.0);
  79. fac=sqrt(-2.0*log(rsq)/rsq);
  80. gset=v1*fac;
  81. iset=1;
  82. return v2*fac;
  83.     } else {
  84. iset=0;
  85. return gset;
  86.     }
  87. }
  90. /*---------------------------------------------------------------------------*/
  91. void ranprobs(vector<double>& v, int& idum)
  92. {
  93.     register double sum = 0.0;
  94.     register int i, size = v.size();
  96.     for (i = 0; i < size; i++) {
  97. v[i] = ran1(idum);
  98. sum += v[i];
  99.     }
  101.     // Normalize so add to one
  102.     for (i = 0; i < size; i++)
  103. v[i] = v[i] / sum;
  104. }
  107. /*---------------------------------------------------------------------------*/
  108. int ranind(const vector<double>& vprob, int& idum)
  109. {
  110.     register int k = 0;
  111.     register double cum = 0.0;
  112.     register double rn = ran1(idum);
  114.     while (cum < rn)
  115. cum += vprob[k++];
  117.     return (k-1);
  118. }