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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > 生物技术 > 查看源码
MyUtils.cpp
资源名称：ClusterTreeViewCode.zip　[点击查看]
上传用户：szb0815
上传日期：2007-06-13
资源大小：338k
文件大小：9k
源码类别：

生物技术

开发平台：

C++ Builder

MyUtils.cpp：源码内容
  1. //---------------------------------------------------------------------------
  2. #include <vcl.h>
  3. #pragma hdrstop
  5. #include "MyUtils.h"
  7. //---------------------------------------------------------------------------
  8. #pragma package(smart_init)
  11. float** __fastcall newmatrix(int n, int m)
  12. {
  13.     float **TempPointer = new float*[n];
  15.     for (int i=0;i<n;i++)
  16.     {
  17.         TempPointer[i] = new float[m];
  18.         for (int j=0; j<m; j++)
  19.         {
  20.             TempPointer[i][j] = 0;
  21.         }
  22.     }
  24.     return TempPointer;
  25. }
  27. void __fastcall deletematrix(float** Matrix, int n)
  28. {
  29.     for (int i=0; i<n; i++)
  30.     {
  31.         delete Matrix[i];
  32.     }
  34.     delete Matrix;
  35. }
  37. int** __fastcall newimatrix(int n, int m)
  38. {
  39.     int **TempPointer = new int*[n];
  41.     for (int i=0;i<n;i++)
  42.     {
  43.         TempPointer[i] = new int[m];
  44.         for (int j=0; j<m; j++)
  45.         {
  46.             TempPointer[i][j] = 0;
  47.         }
  48.     }
  50.     return TempPointer;
  51. }
  53. void __fastcall deleteimatrix(int** Matrix, int n)
  54. {
  55.     for (int i=0; i<n; i++)
  56.     {
  57.         delete Matrix[i];
  58.     }
  60.     delete Matrix;
  61. }
  63. float * __fastcall newvector(int n)
  64. {
  65.     float *TempVector = new float[n];
  66.     for (int i=0; i<n; i++)
  67.     {
  68.         TempVector[i] = 0;
  69.     }
  71.     return TempVector;
  72. }
  74. int * __fastcall newivector(int n)
  75. {
  76.     int *TempVector = new int[n];
  77.     for (int i=0; i<n; i++)
  78.     {
  79.         TempVector[i] = 0;
  80.     }
  82.     return TempVector;
  83. }
  85. #define NRANSI
  86. #define SWAP(a,b) temp=(a);(a)=(b);(b)=temp; #define M 7 #define NSTACK 50
  88. void sort2fi(unsigned long n, float arr[], int brr[])
  89. { unsigned long i,ir=n,j,k,l=1; int *istack,jstack=0; float a,temp;         int b; istack=ivector(1,NSTACK); for (;;) { if (ir-l < M) { for (j=l+1;j<=ir;j++) { a=arr[j]; b=brr[j]; for (i=j-1;i>=l;i--) { if (arr[i] <= a) break; arr[i+1]=arr[i]; brr[i+1]=brr[i]; } arr[i+1]=a; brr[i+1]=b; } if (!jstack) { free_ivector(istack,1,NSTACK); return; } ir=istack[jstack]; l=istack[jstack-1]; jstack -= 2; } else { k=(l+ir) >> 1; SWAP(arr[k],arr[l+1]) SWAP(brr[k],brr[l+1]) if (arr[l] > arr[ir]) { SWAP(arr[l],arr[ir]) SWAP(brr[l],brr[ir]) } if (arr[l+1] > arr[ir]) { SWAP(arr[l+1],arr[ir]) SWAP(brr[l+1],brr[ir]) } if (arr[l] > arr[l+1]) { SWAP(arr[l],arr[l+1]) SWAP(brr[l],brr[l+1]) } i=l+1; j=ir; a=arr[l+1]; b=brr[l+1]; for (;;) { do i++; while (arr[i] < a); do j--; while (arr[j] > a); if (j < i) break; SWAP(arr[i],arr[j]) SWAP(brr[i],brr[j]) } arr[l+1]=arr[j]; arr[j]=a; brr[l+1]=brr[j]; brr[j]=b; jstack += 2; if (jstack > NSTACK) nrerror("NSTACK too small in sort2."); if (ir-i+1 >= j-l) { istack[jstack]=ir; istack[jstack-1]=i; ir=j-1; } else { istack[jstack]=j-1; istack[jstack-1]=l; l=i; } } } } #undef M #undef NSTACK #undef SWAP #undef NRANSI
  91. #define NRANSI
  92. #define SWAP(a,b) temp=(a);(a)=(b);(b)=temp;
  93. #define M 7
  94. #define NSTACK 50
  96. void sort2i(unsigned long n, int arr[], int brr[])
  97. {
  98. unsigned long i,ir=n,j,k,l=1;
  99. int *istack,jstack=0;
  100. int a,b,temp;
  102. istack=ivector(1,NSTACK);
  103. for (;;) {
  104. if (ir-l < M) {
  105. for (j=l+1;j<=ir;j++) {
  106. a=arr[j];
  107. b=brr[j];
  108. for (i=j-1;i>=l;i--) {
  109. if (arr[i] <= a) break;
  110. arr[i+1]=arr[i];
  111. brr[i+1]=brr[i];
  112. }
  113. arr[i+1]=a;
  114. brr[i+1]=b;
  115. }
  116. if (!jstack) {
  117. free_ivector(istack,1,NSTACK);
  118. return;
  119. }
  120. ir=istack[jstack];
  121. l=istack[jstack-1];
  122. jstack -= 2;
  123. } else {
  124. k=(l+ir) >> 1;
  125. SWAP(arr[k],arr[l+1])
  126. SWAP(brr[k],brr[l+1])
  127. if (arr[l] > arr[ir]) {
  128. SWAP(arr[l],arr[ir])
  129. SWAP(brr[l],brr[ir])
  130. }
  131. if (arr[l+1] > arr[ir]) {
  132. SWAP(arr[l+1],arr[ir])
  133. SWAP(brr[l+1],brr[ir])
  134. }
  135. if (arr[l] > arr[l+1]) {
  136. SWAP(arr[l],arr[l+1])
  137. SWAP(brr[l],brr[l+1])
  138. }
  139. i=l+1;
  140. j=ir;
  141. a=arr[l+1];
  142. b=brr[l+1];
  143. for (;;) {
  144. do i++; while (arr[i] < a);
  145. do j--; while (arr[j] > a);
  146. if (j < i) break;
  147. SWAP(arr[i],arr[j])
  148. SWAP(brr[i],brr[j])
  149. }
  150. arr[l+1]=arr[j];
  151. arr[j]=a;
  152. brr[l+1]=brr[j];
  153. brr[j]=b;
  154. jstack += 2;
  155. if (jstack > NSTACK) nrerror("NSTACK too small in sort2.");
  156. if (ir-i+1 >= j-l) {
  157. istack[jstack]=ir;
  158. istack[jstack-1]=i;
  159. ir=j-1;
  160. } else {
  161. istack[jstack]=j-1;
  162. istack[jstack-1]=l;
  163. l=i;
  164. }
  165. }
  166. }
  167. }
  168. #undef M
  169. #undef NSTACK
  170. #undef SWAP
  171. #undef NRANSI
  174. #define NRANSI
  175. #define SWAP(a,b) temp=(a);(a)=(b);(b)=temp; #define M 7 #define NSTACK 50
  177. void sort3fff(unsigned long n, float arr[], float brr[], float crr[])
  178. { unsigned long i,ir=n,j,k,l=1; int *istack,jstack=0; float a,temp;         float b,c; istack=ivector(1,NSTACK); for (;;) { if (ir-l < M) { for (j=l+1;j<=ir;j++) { a=arr[j]; b=brr[j];                                 c=crr[j]; for (i=j-1;i>=l;i--) { if (arr[i] <= a) break; arr[i+1]=arr[i]; brr[i+1]=brr[i];                                         crr[i+1]=crr[i]; } arr[i+1]=a; brr[i+1]=b;                                 crr[i+1]=c; } if (!jstack) { free_ivector(istack,1,NSTACK); return; } ir=istack[jstack]; l=istack[jstack-1]; jstack -= 2; } else { k=(l+ir) >> 1; SWAP(arr[k],arr[l+1]) SWAP(brr[k],brr[l+1])                         SWAP(crr[k],crr[l+1]) if (arr[l] > arr[ir]) { SWAP(arr[l],arr[ir]) SWAP(brr[l],brr[ir])                                 SWAP(crr[l],crr[ir]) } if (arr[l+1] > arr[ir]) { SWAP(arr[l+1],arr[ir]) SWAP(brr[l+1],brr[ir]) SWAP(crr[l+1],crr[ir]) } if (arr[l] > arr[l+1]) { SWAP(arr[l],arr[l+1]) SWAP(brr[l],brr[l+1]) SWAP(crr[l],crr[l+1]) } i=l+1; j=ir; a=arr[l+1]; b=brr[l+1];                         c=crr[l+1]; for (;;) { do i++; while (arr[i] < a); do j--; while (arr[j] > a); if (j < i) break; SWAP(arr[i],arr[j]) SWAP(brr[i],brr[j]) SWAP(crr[i],crr[j]) } arr[l+1]=arr[j]; arr[j]=a; brr[l+1]=brr[j]; brr[j]=b; crr[l+1]=crr[j]; crr[j]=c; jstack += 2; if (jstack > NSTACK) nrerror("NSTACK too small in sort2."); if (ir-i+1 >= j-l) { istack[jstack]=ir; istack[jstack-1]=i; ir=j-1; } else { istack[jstack]=j-1; istack[jstack-1]=l; l=i; } } } } #undef M #undef NSTACK #undef SWAP #undef NRANSI
  180. #define NRANSI
  181. #define SWAP(a,b) temp=(a);(a)=(b);(b)=temp; #define M 7 #define NSTACK 50
  183. void sort2fi_0(unsigned long n, float arr[], int brr[])
  184. { unsigned long i,ir=n-1,j,k,l=0; int *istack,jstack=0; float a,temp;         int b; istack=ivector(1,NSTACK); for (;;) { if (ir-l < M) { for (j=l+1;j<ir;j++) { a=arr[j]; b=brr[j]; for (i=j-1;i>=l;i--) { if (arr[i] <= a) break; arr[i+1]=arr[i]; brr[i+1]=brr[i]; } arr[i+1]=a; brr[i+1]=b; } if (!jstack) { free_ivector(istack,1,NSTACK); return; } ir=istack[jstack]; l=istack[jstack-1]; jstack -= 2; } else { k=(l+ir) >> 1; SWAP(arr[k],arr[l+1]) SWAP(brr[k],brr[l+1]) if (arr[l] > arr[ir]) { SWAP(arr[l],arr[ir]) SWAP(brr[l],brr[ir]) } if (arr[l+1] > arr[ir]) { SWAP(arr[l+1],arr[ir]) SWAP(brr[l+1],brr[ir]) } if (arr[l] > arr[l+1]) { SWAP(arr[l],arr[l+1]) SWAP(brr[l],brr[l+1]) } i=l+1; j=ir; a=arr[l+1]; b=brr[l+1]; for (;;) { do i++; while (arr[i] < a); do j--; while (arr[j] > a); if (j < i) break; SWAP(arr[i],arr[j]) SWAP(brr[i],brr[j]) } arr[l+1]=arr[j]; arr[j]=a; brr[l+1]=brr[j]; brr[j]=b; jstack += 2; if (jstack > NSTACK) nrerror("NSTACK too small in sort2."); if (ir-i+1 >= j-l) { istack[jstack]=ir; istack[jstack-1]=i; ir=j-1; } else { istack[jstack]=j-1; istack[jstack-1]=l; l=i; } } } } #undef M #undef NSTACK #undef SWAP #undef NRANSI
  186. #define NRANSI
  187. #define SWAP(a,b) temp=(a);(a)=(b);(b)=temp; #define M 7 #define NSTACK 50
  189. void sort2_0(unsigned long n, float arr[], float brr[])
  190. { unsigned long i,ir=n-1,j,k,l=0; int *istack,jstack=0; float a,temp;         float b; istack=new int[NSTACK]; for (;;) { if (ir-l < M) { for (j=l+1;j<ir;j++) { a=arr[j]; b=brr[j]; for (i=j-1;i>=l;i--) { if (arr[i] <= a) break; arr[i+1]=arr[i]; brr[i+1]=brr[i]; } arr[i+1]=a; brr[i+1]=b; } if (!jstack) { delete istack; return; } ir=istack[jstack]; l=istack[jstack-1]; jstack -= 2; } else { k=(l+ir) >> 1; SWAP(arr[k],arr[l+1]) SWAP(brr[k],brr[l+1]) if (arr[l] > arr[ir]) { SWAP(arr[l],arr[ir]) SWAP(brr[l],brr[ir]) } if (arr[l+1] > arr[ir]) { SWAP(arr[l+1],arr[ir]) SWAP(brr[l+1],brr[ir]) } if (arr[l] > arr[l+1]) { SWAP(arr[l],arr[l+1]) SWAP(brr[l],brr[l+1]) } i=l+1; j=ir; a=arr[l+1]; b=brr[l+1]; for (;;) { do i++; while (arr[i] < a); do j--; while (arr[j] > a); if (j < i) break; SWAP(arr[i],arr[j]) SWAP(brr[i],brr[j]) } arr[l+1]=arr[j]; arr[j]=a; brr[l+1]=brr[j]; brr[j]=b; jstack += 2; if (jstack > NSTACK) nrerror("NSTACK too small in sort2."); if (ir-i+1 >= j-l) { istack[jstack]=ir; istack[jstack-1]=i; ir=j-1; } else { istack[jstack]=j-1; istack[jstack-1]=l; l=i; } } } } #undef M #undef NSTACK #undef SWAP #undef NRANSI