SVMF32.cpp 6.7 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259
  1. #include "SVMF32.h"
  2. #include "Error.h"
  3. void SVMF32::test_svm_linear_predict_f32()
  4. {
  5. const float32_t *inp = samples.ptr();
  6. int32_t *refp = ref.ptr();
  7. int32_t *outp = output.ptr();
  8. int32_t *result;
  9. result=outp;
  10. for(int i =0; i < this->nbTestSamples; i++)
  11. {
  12. arm_svm_linear_predict_f32(&this->linear,inp,result);
  13. result++;
  14. inp += this->vecDim;
  15. }
  16. ASSERT_EQ(ref,output);
  17. }
  18. void SVMF32::test_svm_polynomial_predict_f32()
  19. {
  20. const float32_t *inp = samples.ptr();
  21. int32_t *refp = ref.ptr();
  22. int32_t *outp = output.ptr();
  23. int32_t *result;
  24. result=outp;
  25. for(int i =0; i < this->nbTestSamples; i++)
  26. {
  27. arm_svm_polynomial_predict_f32(&this->poly,inp,result);
  28. result++;
  29. inp += this->vecDim;
  30. }
  31. ASSERT_EQ(ref,output);
  32. }
  33. void SVMF32::test_svm_rbf_predict_f32()
  34. {
  35. const float32_t *inp = samples.ptr();
  36. int32_t *refp = ref.ptr();
  37. int32_t *outp = output.ptr();
  38. int32_t *result;
  39. result=outp;
  40. for(int i =0; i < this->nbTestSamples; i++)
  41. {
  42. arm_svm_rbf_predict_f32(&this->rbf,inp,result);
  43. result++;
  44. inp += this->vecDim;
  45. }
  46. ASSERT_EQ(ref,output);
  47. }
  48. void SVMF32::test_svm_sigmoid_predict_f32()
  49. {
  50. const float32_t *inp = samples.ptr();
  51. int32_t *refp = ref.ptr();
  52. int32_t *outp = output.ptr();
  53. int32_t *result;
  54. result=outp;
  55. for(int i =0; i < this->nbTestSamples; i++)
  56. {
  57. arm_svm_sigmoid_predict_f32(&this->sigmoid,inp,result);
  58. result++;
  59. inp += this->vecDim;
  60. }
  61. ASSERT_EQ(ref,output);
  62. }
  63. void SVMF32::setUp(Testing::testID_t id,std::vector<Testing::param_t>& testparams,Client::PatternMgr *mgr)
  64. {
  65. int kind;
  66. Testing::nbSamples_t nb=MAX_NB_SAMPLES;
  67. switch(id)
  68. {
  69. case SVMF32::TEST_SVM_LINEAR_PREDICT_F32_1:
  70. {
  71. samples.reload(SVMF32::SAMPLES1_F32_ID,mgr,nb);
  72. params.reload(SVMF32::PARAMS1_F32_ID,mgr,nb);
  73. dims.reload(SVMF32::DIMS1_S16_ID,mgr,nb);
  74. ref.reload(SVMF32::REF1_S32_ID,mgr,nb);
  75. }
  76. break;
  77. case SVMF32::TEST_SVM_POLYNOMIAL_PREDICT_F32_2:
  78. {
  79. samples.reload(SVMF32::SAMPLES2_F32_ID,mgr,nb);
  80. params.reload(SVMF32::PARAMS2_F32_ID,mgr,nb);
  81. dims.reload(SVMF32::DIMS2_S16_ID,mgr,nb);
  82. ref.reload(SVMF32::REF2_S32_ID,mgr,nb);
  83. }
  84. break;
  85. case SVMF32::TEST_SVM_RBF_PREDICT_F32_3:
  86. {
  87. samples.reload(SVMF32::SAMPLES3_F32_ID,mgr,nb);
  88. params.reload(SVMF32::PARAMS3_F32_ID,mgr,nb);
  89. dims.reload(SVMF32::DIMS3_S16_ID,mgr,nb);
  90. ref.reload(SVMF32::REF3_S32_ID,mgr,nb);
  91. }
  92. break;
  93. case SVMF32::TEST_SVM_SIGMOID_PREDICT_F32_4:
  94. {
  95. samples.reload(SVMF32::SAMPLES4_F32_ID,mgr,nb);
  96. params.reload(SVMF32::PARAMS4_F32_ID,mgr,nb);
  97. dims.reload(SVMF32::DIMS4_S16_ID,mgr,nb);
  98. ref.reload(SVMF32::REF4_S32_ID,mgr,nb);
  99. }
  100. break;
  101. case SVMF32::TEST_SVM_RBF_PREDICT_F32_5:
  102. {
  103. samples.reload(SVMF32::SAMPLES5_F32_ID,mgr,nb);
  104. params.reload(SVMF32::PARAMS5_F32_ID,mgr,nb);
  105. dims.reload(SVMF32::DIMS5_S16_ID,mgr,nb);
  106. ref.reload(SVMF32::REF5_S32_ID,mgr,nb);
  107. }
  108. break;
  109. }
  110. const int16_t *dimsp = dims.ptr();
  111. const float32_t *paramsp = params.ptr();
  112. kind = dimsp[0];
  113. this->classes[0] = dimsp[1];
  114. this->classes[1] = dimsp[2];
  115. this->nbTestSamples=dimsp[3];
  116. this->vecDim = dimsp[4];
  117. this->nbSupportVectors = dimsp[5];
  118. this->intercept=paramsp[this->vecDim*this->nbSupportVectors + this->nbSupportVectors];
  119. this->supportVectors=paramsp;
  120. this->dualCoefs=paramsp + (this->vecDim*this->nbSupportVectors);
  121. switch(kind)
  122. {
  123. case SVMF32::POLY:
  124. this->degree = dimsp[6];
  125. this->coef0 =paramsp[this->vecDim*this->nbSupportVectors + this->nbSupportVectors + 1] ;
  126. this->gamma=paramsp[this->vecDim*this->nbSupportVectors + this->nbSupportVectors + 2];
  127. break;
  128. case SVMF32::RBF:
  129. this->gamma=paramsp[this->vecDim*this->nbSupportVectors + this->nbSupportVectors + 1];
  130. break;
  131. case SVMF32::SIGMOID:
  132. this->coef0 =paramsp[this->vecDim*this->nbSupportVectors + this->nbSupportVectors + 1] ;
  133. this->gamma=paramsp[this->vecDim*this->nbSupportVectors + this->nbSupportVectors + 2];
  134. break;
  135. }
  136. switch(id)
  137. {
  138. case SVMF32::TEST_SVM_LINEAR_PREDICT_F32_1:
  139. {
  140. arm_svm_linear_init_f32(&linear,
  141. this->nbSupportVectors,
  142. this->vecDim,
  143. this->intercept,
  144. this->dualCoefs,
  145. this->supportVectors,
  146. this->classes);
  147. }
  148. break;
  149. case SVMF32::TEST_SVM_POLYNOMIAL_PREDICT_F32_2:
  150. {
  151. arm_svm_polynomial_init_f32(&poly,
  152. this->nbSupportVectors,
  153. this->vecDim,
  154. this->intercept,
  155. this->dualCoefs,
  156. this->supportVectors,
  157. this->classes,
  158. this->degree,
  159. this->coef0,
  160. this->gamma
  161. );
  162. }
  163. break;
  164. case SVMF32::TEST_SVM_RBF_PREDICT_F32_3:
  165. case SVMF32::TEST_SVM_RBF_PREDICT_F32_5:
  166. {
  167. arm_svm_rbf_init_f32(&rbf,
  168. this->nbSupportVectors,
  169. this->vecDim,
  170. this->intercept,
  171. this->dualCoefs,
  172. this->supportVectors,
  173. this->classes,
  174. this->gamma
  175. );
  176. }
  177. break;
  178. case SVMF32::TEST_SVM_SIGMOID_PREDICT_F32_4:
  179. {
  180. arm_svm_sigmoid_init_f32(&sigmoid,
  181. this->nbSupportVectors,
  182. this->vecDim,
  183. this->intercept,
  184. this->dualCoefs,
  185. this->supportVectors,
  186. this->classes,
  187. this->coef0,
  188. this->gamma
  189. );
  190. }
  191. break;
  192. }
  193. output.create(ref.nbSamples(),SVMF32::OUT_S32_ID,mgr);
  194. }
  195. void SVMF32::tearDown(Testing::testID_t id,Client::PatternMgr *mgr)
  196. {
  197. output.dump(mgr);
  198. }