apps.c 60 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985986987988989990991992993994995996997998999100010011002100310041005100610071008100910101011101210131014101510161017101810191020102110221023102410251026102710281029103010311032103310341035103610371038103910401041104210431044104510461047104810491050105110521053105410551056105710581059106010611062106310641065106610671068106910701071107210731074107510761077107810791080108110821083108410851086108710881089109010911092109310941095109610971098109911001101110211031104110511061107110811091110111111121113111411151116111711181119112011211122112311241125112611271128112911301131113211331134113511361137113811391140114111421143114411451146114711481149115011511152115311541155115611571158115911601161116211631164116511661167116811691170117111721173117411751176117711781179
  1. #include "apps.h"
  2. #include "scc_signals.h"
  3. #include "libfunctions.h"
  4. #include "my_rtrm.h"
  5. #define SWAP(a,b) {float tmp; tmp=a; a=b; b=tmp;}
  6. #define FFT_MAX 136192
  7. #define PAGE_SIZE 4096
  8. /*
  9. static float input_vector[25][10] = {
  10. {-0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119},
  11. {-0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467},
  12. {-0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565},
  13. {-0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533},
  14. {-0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941},
  15. {-1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784},
  16. {0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575},
  17. {0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371},
  18. {1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241},
  19. {1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752},
  20. {-0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115},
  21. {-0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876},
  22. {-0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559},
  23. {-0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523},
  24. {1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364},
  25. {-0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320},
  26. {3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579},
  27. {-0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210},
  28. {1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734},
  29. {-0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990},
  30. {-0.632187, -0.629597, -0.560808, -0.439143, -0.550906, -3.344673, -1.389866, 0.268649, 0.664319, 1.150327},
  31. {-0.205756, -0.210472, -0.151426, -0.121347, 0.196067, -3.136218, -2.621049, -0.026517, 0.358534, 0.714117},
  32. {-0.418011, -0.424854, -0.461205, -0.428858, -0.801747, 1.933860, -0.129047, -0.674498, -0.880092, -0.752953},
  33. {-0.625461, -0.633598, -0.651167, -0.621632, -0.312866, -3.908468, -2.380095, -0.118114, 0.233478, 0.722539},
  34. {-0.525633, -0.521436, -0.552314, -0.527505, -0.233392, -3.763046, -2.487090, -0.133160, 0.156544, 0.642195},
  35. };
  36. */
  37. /*
  38. static float input_vector[2][D_sv] = {
  39. {-0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  40. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  41. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  42. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  43. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  44. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  45. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  46. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  47. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  48. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752},
  49. {-0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  50. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  51. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  52. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  53. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  54. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  55. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  56. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  57. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  58. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990},
  59. };
  60. */
  61. /*
  62. static float input_vector[1][D_sv] = {
  63. {-0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  64. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  65. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  66. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  67. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  68. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  69. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  70. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  71. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  72. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  73. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  74. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  75. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  76. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  77. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  78. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  79. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  80. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  81. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  82. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  83. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  84. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  85. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  86. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  87. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  88. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  89. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  90. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  91. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  92. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  93. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  94. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  95. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  96. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  97. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  98. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  99. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  100. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  101. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  102. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  103. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  104. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  105. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  106. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  107. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  108. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  109. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  110. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  111. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  112. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  113. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  114. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  115. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  116. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  117. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  118. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  119. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  120. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  121. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  122. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  123. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  124. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  125. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  126. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  127. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  128. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  129. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  130. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  131. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  132. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  133. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  134. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  135. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  136. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  137. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  138. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  139. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  140. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  141. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  142. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  143. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  144. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  145. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  146. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  147. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  148. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  149. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  150. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  151. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  152. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  153. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  154. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  155. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  156. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  157. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  158. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  159. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  160. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  161. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  162. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990
  163. },
  164. };
  165. */
  166. /*
  167. static float input_vector[1][2100] = {
  168. {-0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  169. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  170. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  171. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  172. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  173. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  174. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  175. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  176. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  177. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  178. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  179. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  180. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  181. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  182. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  183. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  184. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  185. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  186. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  187. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  188. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  189. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  190. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  191. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  192. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  193. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  194. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  195. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  196. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  197. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  198. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  199. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  200. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  201. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  202. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  203. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  204. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  205. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  206. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  207. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  208. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  209. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  210. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  211. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  212. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  213. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  214. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  215. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  216. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  217. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  218. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  219. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  220. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  221. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  222. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  223. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  224. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  225. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  226. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  227. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  228. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  229. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  230. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  231. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  232. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  233. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  234. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  235. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  236. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  237. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  238. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  239. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  240. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  241. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  242. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  243. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  244. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  245. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  246. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  247. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  248. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  249. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  250. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  251. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  252. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  253. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  254. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  255. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  256. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  257. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  258. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  259. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  260. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  261. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  262. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  263. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  264. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  265. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  266. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  267. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  268. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  269. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  270. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  271. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  272. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  273. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  274. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  275. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  276. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  277. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  278. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  279. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  280. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  281. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  282. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  283. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  284. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  285. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  286. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  287. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  288. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  289. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  290. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  291. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  292. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  293. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  294. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  295. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  296. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  297. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  298. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  299. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  300. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  301. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  302. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  303. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  304. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  305. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  306. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  307. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  308. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  309. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  310. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  311. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  312. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  313. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  314. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  315. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  316. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  317. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  318. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  319. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  320. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  321. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  322. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  323. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  324. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  325. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  326. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  327. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  328. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  329. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  330. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  331. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  332. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  333. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  334. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  335. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  336. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  337. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  338. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  339. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  340. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  341. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  342. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  343. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  344. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  345. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  346. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  347. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  348. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  349. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  350. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  351. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  352. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  353. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  354. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  355. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  356. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  357. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752,
  358. -0.664593, -0.687087, -0.075215, 0.583554, 0.927597, 3.671647, 3.835722, 0.694133, -2.007198, -2.345115,
  359. -0.343777, -0.344074, -0.221304, -0.153775, -0.494637, 1.743759, -0.590775, -0.413656, -0.389204, -0.221876,
  360. -0.159571, -0.163080, -0.513428, -0.775665, -0.891170, -3.118664, -1.247074, 0.913688, 1.009656, 1.170559,
  361. -0.761538, -0.755972, -0.773967, -0.674253, -0.978258, 1.514730, -0.145977, -0.620700, -0.857082, -0.765523,
  362. 1.344641, 1.329151, 1.633256, 2.020196, 1.777348, -1.744907, -5.928727, -4.032057, -0.585966, 1.072364,
  363. -0.468882, -0.486856, 0.935241, 2.032263, 2.219547, -0.225632, -3.438104, -4.412505, -0.550391, 2.234320,
  364. 3.177341, 3.204259, 2.749250, 1.741244, 1.197538, -2.985864, -6.084715, -1.447878, 0.904210, 1.495579,
  365. -0.998953, -1.005918, -0.988911, -0.863153, -1.003750, 1.256336, -0.471785, -0.860056, -0.872804, -0.629210,
  366. 1.957319, 1.966453, 1.840960, 1.405216, 1.313205, -0.956540, -3.535391, -2.280320, -1.179478, -0.458734,
  367. -0.326054, -0.331524, -0.134208, -0.218622, -0.158037, 3.128528, 4.020623, -1.129257, -1.524952, -1.377990,
  368. -0.390695, -0.381094, -0.341754, -0.272043, -0.588159, 2.297114, -0.116822, -0.507031, -0.563574, -0.534119,
  369. -0.250457, -0.240673, -0.210423, -0.217804, -0.453731, 1.784484, -0.142005, -0.230148, -0.339593, -0.298467,
  370. -0.177671, -0.183887, -0.191644, -0.205313, -0.441613, 1.599526, 0.023297, -0.322458, -0.539978, -0.457565,
  371. -0.149024, -0.140364, -0.142723, -0.066651, -0.344126, 1.830863, 0.154307, -0.171898, -0.343950, -0.379533,
  372. -0.155498, -0.154145, -0.181875, -0.189322, -0.122244, -3.037579, -1.476675, -0.089878, 0.004249, 0.169941,
  373. -1.786535, -1.809749, -1.445913, -1.050310, -0.095173, 2.964720, 4.706277, -0.649227, -1.362017, -1.447784,
  374. 0.034026, 0.038159, -0.018497, -0.067303, 0.059710, -2.908602, -1.484229, 0.108804, 0.292757, 0.566575,
  375. 0.527511, 0.518252, 0.762342, 1.148066, 0.772846, -2.113671, -2.704303, -0.664257, 0.460104, 1.295371,
  376. 1.370790, 1.383803, 1.340450, 0.634441, -0.596868, -2.427651, -2.352223, 0.454767, 1.405430, 1.838241,
  377. 1.455314, 1.464866, 2.487861, 2.464788, 1.963237, 0.324097, -1.869693, -2.554004, -2.022594, 1.114752
  378. },
  379. };
  380. */
  381. static float **svm_vectors, *svm_coef;
  382. static int *vector, **matrix;
  383. static float input_vector[D_sv];
  384. //static float matr_speedup[NUM_OF_MATRICES][MAX_WORKERS_COUNT];
  385. //static int matr_times[NUM_OF_MATRICES][MAX_WORKERS_COUNT];
  386. static float Exec_Speedup[MAX_WORKERS_COUNT];
  387. static int Exec_Latencies[MAX_WORKERS_COUNT];
  388. //static float **vectors, *coef;
  389. //2*(N+rootN*pad_length)*sizeof(float)+PAGE_SIZE);
  390. static int P = 1; /* DEFAULT_P = 1 */
  391. static int M = 16; /* DEFAULT_M = 10 */
  392. static int N = 65536; /* N = 2^M */
  393. static int rootN = 256; /* rootN = sqrt(N) */
  394. static int num_cache_lines = 65536;
  395. #define PADLENGTH 2
  396. static float *x_local; /* x is the original time-domain data */
  397. static float *trans; /* trans is used as scratch space */
  398. static float *umain; /* umain is roots of unity for 1D FFTs */
  399. static float *umain2; /* umain2 is entire roots of unity matrix*/
  400. static float *upriv;
  401. void execute_workload_svm (int lower_bound, int upper_bound);
  402. void execute_workload_matrix (int lower_bound, int upper_bound);
  403. void matrix_transpose(int n1, float *src, float *dest, int node_id, int myFirst, int myLast, int pad_length);
  404. void FFT1D(int direction, int M, int N, float *x, float *scratch, float *upriv, float *umain2, int node_id, int myFirst, int myLast, int pad_length, int P);
  405. void copyColumn(int n1, float *src, float *dest);
  406. void single_FFT1D(int direction, int M, int N, float *u, float *x);
  407. void twiddle_Col(int direction, int n1, int N, int j, float *u, float *x, int pad_length);
  408. void reverse(int N, int M, float *x);
  409. int reverse_bit(int M, int k);
  410. void execute_workload_svm (int lower_bound, int upper_bound) {
  411. int i = 0, j = 0;
  412. float diff = 0, norma = 0, local_sum[N_sv];
  413. /* int vector_id = 0; Removed 16.02. Only one test vector */
  414. if (base_offset == -1) {
  415. base_offset = cur_agent.my_agent * N_sv;
  416. //fprintf(log_file, "My agent is %d. Calculated base_offset is %d\n",cur_agent.my_agent,base_offset);
  417. }
  418. for (i = lower_bound; i <= upper_bound; i++) {
  419. local_sum[i] = 0;
  420. scc_signals_check();
  421. for (j = 0; j < D_sv; j++){
  422. diff = svm_vectors[i][j] - input_vector[j];
  423. norma += diff*diff;
  424. }
  425. local_sum[i] += (float) (exp((double) (-gamma*norma))*svm_coef[i]);
  426. norma = 0;
  427. }
  428. for (i=lower_bound; i<=upper_bound; i++)
  429. manager_result_out[base_offset+i] = (int) local_sum[i];
  430. }
  431. void execute_workload_matrix (int lower_bound, int upper_bound) {
  432. int i, j, local_sum[MAX_ARRAY_SIZE];
  433. if (base_offset == -1) {
  434. //matrix_out = (int*) shmat (cur_agent.segment_id, NULL, 0);
  435. base_offset = cur_agent.my_agent * MAX_ARRAY_SIZE;
  436. }
  437. for (i=lower_bound; i<=upper_bound; i++) {
  438. local_sum[i] = 0;
  439. scc_signals_check();
  440. //signals_enable();
  441. for (j=0; j<cur_agent.array_size; j++)
  442. local_sum[i] += matrix[i][j] * vector[j];
  443. //signals_disable();
  444. }
  445. for (i=lower_bound; i<=upper_bound; i++)
  446. manager_result_out[base_offset+i] = local_sum[i];
  447. }
  448. void execute_workload_fft (int lower_bound, int upper_bound) {
  449. int work_id = 0, pad_length = PADLENGTH;
  450. if ((lower_bound == 0) && (upper_bound == FFT_MAX)) {
  451. P = 1;
  452. } else {
  453. P = 2;
  454. }
  455. /* FIXME works only because fft is restricted to two workers */
  456. if (lower_bound > 0) {
  457. work_id = 1;
  458. }
  459. FFT1D(1, M, N, x_local, trans, upriv, umain2, work_id, lower_bound, upper_bound, pad_length, P); //HACK node_id - 1 important!!
  460. }
  461. void execute_workload (int lower_bound, int upper_bound) {
  462. if (executed_app == MATRIX_MUL) {
  463. execute_workload_matrix (lower_bound, upper_bound);
  464. } else if (executed_app == SVM) {
  465. execute_workload_svm (lower_bound, upper_bound);
  466. } else if (executed_app == FFT) {
  467. execute_workload_fft (lower_bound, upper_bound);
  468. }
  469. }
  470. void init_speedup_structs (void) {
  471. if (executed_app == MATRIX_MUL) {
  472. if (MATRIX_ARRAY_SIZE == 1024) {
  473. #ifdef PLAT_SCC
  474. Exec_Speedup[0] = 1.0;
  475. Exec_Speedup[1] = 1.188;
  476. Exec_Speedup[2] = 2.264;
  477. Exec_Speedup[3] = 3.0;
  478. Exec_Speedup[4] = 3.429;
  479. Exec_Speedup[5] = 4.0;
  480. Exec_Speedup[6] = 8.0;
  481. Exec_Speedup[7] = 0.0;
  482. Exec_Latencies[0] = 120;//29352;
  483. Exec_Latencies[1] = 101;//15112;
  484. Exec_Latencies[2] = 53;//11194;
  485. Exec_Latencies[3] = 40;//10313;
  486. Exec_Latencies[4] = 35;//8645;
  487. Exec_Latencies[5] = 30;//7871;
  488. Exec_Latencies[6] = 15;//6715;
  489. #else
  490. Exec_Speedup[0] = 1.0;
  491. Exec_Speedup[1] = 1.065;
  492. Exec_Speedup[2] = 1.270;
  493. Exec_Speedup[3] = 0.0;
  494. Exec_Speedup[4] = 0.0;
  495. Exec_Speedup[5] = 0.0;
  496. Exec_Speedup[6] = 0.0;
  497. Exec_Speedup[7] = 0.0;
  498. Exec_Latencies[0] = 100000000;//29352;
  499. Exec_Latencies[1] = 31;//15112;
  500. Exec_Latencies[2] = 29;//11194;
  501. Exec_Latencies[3] = 24;//10313;
  502. Exec_Latencies[4] = 0;//8645;
  503. Exec_Latencies[5] = 0;//7871;
  504. Exec_Latencies[6] = 0;//6715;
  505. Exec_Latencies[7] = 0;//7014;
  506. #endif
  507. } else if (MATRIX_ARRAY_SIZE == 2048) {
  508. #ifdef PLAT_SCC
  509. Exec_Speedup[0] = 1.0;
  510. Exec_Speedup[1] = 1.091;
  511. Exec_Speedup[2] = 1.2;
  512. Exec_Speedup[3] = 1.491;
  513. Exec_Speedup[4] = 1.791;
  514. Exec_Speedup[5] = 2.824;
  515. Exec_Speedup[6] = 3.0;
  516. Exec_Latencies[0] = 240;//112276;
  517. Exec_Latencies[1] = 220;//58880;
  518. Exec_Latencies[2] = 200;//40305;
  519. Exec_Latencies[3] = 161;//31705;
  520. Exec_Latencies[4] = 134;//28309;
  521. Exec_Latencies[5] = 85;//24512;
  522. Exec_Latencies[6] = 80;//22239;
  523. //matr_times[1][7] = 23;//20332;
  524. #else
  525. Exec_Speedup[0] = 1.0;
  526. Exec_Speedup[1] = 1.331;
  527. Exec_Speedup[2] = 2.009;
  528. Exec_Speedup[3] = 2.315;
  529. Exec_Speedup[4] = 2.572;
  530. Exec_Speedup[5] = 0.0;
  531. Exec_Speedup[6] = 0.0;
  532. Exec_Speedup[7] = 0.0;//5.522;
  533. Exec_Latencies[0] = 100000000;//112276;
  534. Exec_Latencies[1] = 116;//58880;
  535. Exec_Latencies[2] = 87;//40305;
  536. Exec_Latencies[3] = 58;//31705;
  537. Exec_Latencies[4] = 50;//28309;
  538. Exec_Latencies[5] = 45;//24512;
  539. Exec_Latencies[6] = 0;//22239;
  540. Exec_Latencies[7] = 0;//20332;
  541. #endif
  542. } else if (MATRIX_ARRAY_SIZE == 4096) {
  543. #ifdef PLAT_SCC
  544. Exec_Speedup[0] = 1.0;
  545. Exec_Speedup[1] = 2.001;
  546. Exec_Speedup[2] = 2.976;
  547. Exec_Speedup[3] = 4.032;
  548. Exec_Speedup[4] = 5.034;
  549. Exec_Speedup[5] = 6.25;
  550. Exec_Speedup[6] = 6.678;
  551. Exec_Speedup[7] = 6.819;
  552. Exec_Latencies[0] = 750;//384005;
  553. Exec_Latencies[1] = 374;//231583;
  554. Exec_Latencies[2] = 252;//157966;
  555. Exec_Latencies[3] = 186;//121222;
  556. Exec_Latencies[4] = 149;//101208;
  557. Exec_Latencies[5] = 120;//87852;
  558. Exec_Latencies[6] = 110;//78093;
  559. #else
  560. Exec_Speedup[0] = 1.0;
  561. Exec_Speedup[1] = 1.517;
  562. Exec_Speedup[2] = 1.958;
  563. Exec_Speedup[3] = 2.112;
  564. Exec_Speedup[4] = 2.878;
  565. Exec_Speedup[5] = 3.338;
  566. Exec_Speedup[6] = 4.241;
  567. Exec_Speedup[7] = 0.0;//5.073;
  568. Exec_Latencies[0] = 100000000;//384005;
  569. Exec_Latencies[1] = 431;//231583;
  570. Exec_Latencies[2] = 284;//157966;
  571. Exec_Latencies[3] = 220;//121222;
  572. Exec_Latencies[4] = 204;//101208;
  573. Exec_Latencies[5] = 150;//87852;
  574. Exec_Latencies[6] = 129;//78093;
  575. Exec_Latencies[7] = 102;//75690;
  576. #endif
  577. } else {
  578. printf("Unknown array size\n");
  579. exit(0);
  580. }
  581. } else if (executed_app == SVM) {
  582. Exec_Speedup[0] = 1.0; /* 1 worker */
  583. Exec_Speedup[1] = 2.006;
  584. Exec_Speedup[2] = 2.814;
  585. Exec_Speedup[3] = 3.469;
  586. Exec_Speedup[4] = 4.029;
  587. Exec_Speedup[5] = 4.285;
  588. Exec_Speedup[6] = 4.646;
  589. Exec_Speedup[7] = 0.0;
  590. Exec_Latencies[0] = 28;
  591. Exec_Latencies[1] = 14;
  592. Exec_Latencies[2] = 10;
  593. Exec_Latencies[3] = 8;
  594. Exec_Latencies[4] = 7;
  595. Exec_Latencies[5] = 7;
  596. Exec_Latencies[6] = 6;
  597. Exec_Latencies[7] = 6;
  598. } else if (executed_app == FFT) {
  599. Exec_Speedup[0] = 1.0; /* 1 worker */
  600. Exec_Speedup[1] = 1.55;
  601. Exec_Speedup[2] = 0;
  602. Exec_Speedup[3] = 0;
  603. Exec_Speedup[4] = 0;
  604. Exec_Speedup[5] = 0;
  605. Exec_Speedup[6] = 0;
  606. Exec_Speedup[7] = 0;
  607. Exec_Latencies[0] = 772;
  608. Exec_Latencies[1] = 498;
  609. Exec_Latencies[2] = 0;
  610. Exec_Latencies[3] = 0;
  611. Exec_Latencies[4] = 0;
  612. Exec_Latencies[5] = 0;
  613. Exec_Latencies[6] = 0;
  614. Exec_Latencies[7] = 0;
  615. }
  616. }
  617. void app_init (char scen_directory[SCEN_DIR_SIZE], char scen_num[SCEN_NUM_SIZE]) {
  618. int i, j, pad_length = PADLENGTH;
  619. char buf[MAX_STR_NAME_SIZE], *buffer;
  620. FILE *matrix_input, *support_vectors_file, *coef_file, *test_vector_file, *umain_file, *umain2_file, *x_local_file;
  621. size_t bufsize = 32;
  622. if (executed_app == MATRIX_MUL) {
  623. cur_agent.array_size = MATRIX_ARRAY_SIZE;
  624. matrix = (int **) malloc(cur_agent.array_size * sizeof(int *));
  625. #ifdef PLAT_SCC
  626. strcpy(buf, "/shared/herc/");
  627. #else
  628. strcpy(buf, "../");
  629. #endif
  630. strcat(buf, scen_directory);
  631. strcat(buf, "/MATRIX-inputs/");
  632. strcat(buf, itoa(cur_agent.array_size));
  633. fprintf(log_file,"matrix file path = %s\n",buf);
  634. if ((matrix_input = fopen(buf, "r")) == NULL){
  635. printf("Cannot open input file with file path = %s ",buf);
  636. perror("open matrix_input");
  637. }
  638. for (i=0; i<cur_agent.array_size; i++) {
  639. matrix[i] = (int *) malloc(cur_agent.array_size * sizeof(int));
  640. for (j=0; j<cur_agent.array_size; j++)
  641. fscanf(matrix_input,"%d",&matrix[i][j]);
  642. }
  643. vector = (int *) malloc(cur_agent.array_size * sizeof(int));
  644. for (j=0; j<cur_agent.array_size; j++)
  645. fscanf(matrix_input,"%d",&vector[j]);
  646. fclose(matrix_input);
  647. } else if (executed_app == SVM) {
  648. #ifdef PLAT_SCC
  649. strcpy(buf, "/shared/herc/");
  650. #else
  651. strcpy(buf, "../");
  652. #endif
  653. strcat(buf,scen_directory);
  654. //strcat(buf,"/");
  655. //strcat(buf,scen_num);
  656. strcat(buf,"/SVM-inputs/support_vectors_N_sv_");
  657. strcat(buf,itoa(N_sv));
  658. strcat(buf,"_D_sv_");
  659. strcat(buf,itoa(D_sv));
  660. strcat(buf,".dat");
  661. fprintf(log_file,"svm file path = %s\n",buf);
  662. if ((support_vectors_file = fopen(buf,"r")) == NULL){
  663. printf("Cannot open input file with file path = %s ",buf);
  664. perror("open svm_input");
  665. }
  666. #ifdef PLAT_SCC
  667. strcpy(buf, "/shared/herc/");
  668. #else
  669. strcpy(buf, "../");
  670. #endif
  671. strcat(buf,scen_directory);
  672. //strcat(buf,"/");
  673. //strcat(buf,scen_num);
  674. strcat(buf,"/SVM-inputs/sv_coef_N_sv_");
  675. strcat(buf,itoa(N_sv));
  676. strcat(buf,"_D_sv_");
  677. strcat(buf,itoa(D_sv));
  678. strcat(buf,".dat");
  679. fprintf(log_file,"svm_coef file path = %s\n",buf);
  680. if ((coef_file = fopen(buf,"r")) == NULL){
  681. printf("Cannot open input file with file path = %s ",buf);
  682. perror("open svm_input");
  683. }
  684. #ifdef PLAT_SCC
  685. strcpy(buf, "/shared/herc/");
  686. #else
  687. strcpy(buf, "../");
  688. #endif
  689. strcat(buf,scen_directory);
  690. //strcat(buf,"/");
  691. //strcat(buf,scen_num);
  692. strcat(buf,"/SVM-inputs/test_vector_D_sv_");
  693. strcat(buf,itoa(D_sv));
  694. strcat(buf,".dat");
  695. fprintf(log_file,"test_vector file path = %s\n",buf);
  696. if ((test_vector_file = fopen(buf,"r")) == NULL){
  697. printf("Cannot open input file with file path = %s ",buf);
  698. perror("open svm_input");
  699. }
  700. svm_vectors = (float **)malloc(N_sv*sizeof(float *));
  701. if (svm_vectors == NULL){
  702. printf("--%d-- svm_vectors malloc fail!!\n", node_id);
  703. perror("malloc error");
  704. }
  705. svm_coef = (float *)malloc(N_sv*sizeof(float));
  706. if (svm_coef == NULL){
  707. printf("--%d-- svm_coef malloc fail!!\n", node_id);
  708. perror("malloc error");
  709. }
  710. buffer = (char *)malloc(bufsize * sizeof(char));
  711. for (i = 0; i < N_sv; i++) {
  712. svm_vectors[i] = (float *)malloc(D_sv*sizeof(float));
  713. if (svm_vectors[i] == NULL) {
  714. printf("--%d-- svm_vectors[%d] malloc fail!!\n", node_id, i);
  715. perror("malloc error");
  716. } else {
  717. for (j = 0; j < D_sv; j++) {
  718. /* Read support svm_vectors */
  719. if (j < D_sv){
  720. fscanf(support_vectors_file,"%f",&svm_vectors[i][j]);
  721. fgetc(support_vectors_file);
  722. }else{
  723. getline(&buffer,&bufsize,support_vectors_file);
  724. }
  725. }
  726. }
  727. }
  728. for (j = 0; j < N_sv; j++) {
  729. /* Read coefficients */
  730. fscanf(coef_file,"%f",&svm_coef[j]);
  731. fgetc(coef_file);
  732. }
  733. for (j = 0; j < D_sv; j++) {
  734. /* Read coefficients */
  735. fscanf(test_vector_file,"%f",&input_vector[j]);
  736. }
  737. cur_agent.array_size = -1;
  738. fclose(support_vectors_file);
  739. fclose(coef_file);
  740. fclose(test_vector_file);
  741. free(buffer);
  742. } else if (executed_app == FFT) {
  743. fprintf(log_file,"Initializing FFT application\n");
  744. x_local = (float *)malloc(2*(N+rootN*pad_length)*sizeof(float)+PAGE_SIZE);
  745. if (x_local == NULL){
  746. printf("Malloc error for x_local\n");
  747. perror("malloc error");
  748. exit(-1);
  749. }
  750. trans = (float *)malloc(2*(N+rootN*pad_length)*sizeof(float)+PAGE_SIZE);
  751. if (trans == NULL){
  752. printf("Malloc error for trans\n");
  753. perror("malloc error");
  754. exit(-1);
  755. }
  756. umain = (float *)malloc(2*rootN*sizeof(float));
  757. if (umain == NULL){
  758. printf("Malloc error for umain\n");
  759. perror("malloc error");
  760. exit(-1);
  761. }
  762. umain2 = (float *)malloc(2*(N+rootN*pad_length)*sizeof(float)+PAGE_SIZE);
  763. if (umain2 == NULL){
  764. printf("Malloc error for umain2\n");
  765. perror("malloc error");
  766. exit(-1);
  767. }
  768. upriv = (float *)malloc(2*(rootN-1)*sizeof(float));
  769. if (upriv == NULL){
  770. printf("--%d-- Malloc error for upriv\n", node_id);
  771. perror("malloc error");
  772. exit(-1);
  773. }
  774. #ifdef PLAT_SCC
  775. strcpy(buf, "/shared/herc/");
  776. #else
  777. strcpy(buf, "../");
  778. #endif
  779. strcat(buf,scen_directory);
  780. //strcat(buf,"/");
  781. //strcat(buf,scen_num);
  782. strcat(buf,"/FFT-inputs/umain_file");
  783. fprintf(log_file,"umain_file file path = %s\n",buf);
  784. if ((umain_file = fopen(buf,"r")) == NULL){
  785. printf("Cannot open input file with file path = %s ",buf);
  786. perror("open fft_input");
  787. }
  788. for (i=0; i<2*rootN; i++) {
  789. fscanf(umain_file,"%f",&umain[i]);
  790. }
  791. fclose(umain_file);
  792. #ifdef PLAT_SCC
  793. strcpy(buf, "/shared/herc/");
  794. #else
  795. strcpy(buf, "../");
  796. #endif
  797. strcat(buf,scen_directory);
  798. //strcat(buf,"/");
  799. //strcat(buf,scen_num);
  800. strcat(buf,"/FFT-inputs/umain2_file");
  801. fprintf(log_file,"umain2_file file path = %s\n",buf);
  802. if ((umain2_file = fopen(buf,"r")) == NULL){
  803. printf("Cannot open input file with file path = %s ",buf);
  804. perror("open umain_file");
  805. }
  806. //for (i=0; i<2*(N+rootN*pad_length)+PAGE_SIZE; i++) {
  807. for (i=0; i<2*(N+rootN*pad_length); i++) {
  808. fscanf(umain2_file,"%f",&umain2[i]);
  809. }
  810. fclose(umain2_file);
  811. #ifdef PLAT_SCC
  812. strcpy(buf, "/shared/herc/");
  813. #else
  814. strcpy(buf, "../");
  815. #endif
  816. strcat(buf,scen_directory);
  817. //strcat(buf,"/");
  818. //strcat(buf,scen_num);
  819. strcat(buf,"/FFT-inputs/x_local_file");
  820. fprintf(log_file,"x_local_file file path = %s\n",buf);
  821. if ((x_local_file = fopen(buf,"r")) == NULL){
  822. printf("Cannot open input file with file path = %s ",buf);
  823. perror("open x_local_file");
  824. }
  825. //for (i=0;i<2*(N+rootN*pad_length)+PAGE_SIZE;i++) {
  826. for (i=0; i<2*(N+rootN*pad_length); i++) {
  827. fscanf(x_local_file,"%f",&x_local[i]);
  828. }
  829. fclose(x_local_file);
  830. for (i = 0; i < 2*(rootN-1); i++){
  831. upriv[i] = umain[i];
  832. }
  833. }
  834. }
  835. int get_max_cores_count(app cur_app){
  836. /*if (cur_app.var < 1.0)
  837. return (int) ceilf(2.0*cur_app.A - 1);
  838. else
  839. return (int) ceilf(cur_app.A + cur_app.A*cur_app.var - cur_app.var);*/
  840. #ifdef SINGLE_WORKER
  841. return 2;
  842. #else
  843. if (executed_app == FFT) {
  844. return 3;
  845. } else {
  846. return MAX_WORKERS_COUNT;
  847. }
  848. #endif
  849. }
  850. float Speedup(app cur_app, int num_of_cores) {
  851. /*
  852. int type;
  853. if (cur_app.array_size == 1024) type = 0;
  854. else if (cur_app.array_size == 2048) type = 1;
  855. else if (cur_app.array_size == 4096) type = 2;
  856. else {
  857. fprintf(log_file, "Unknown array size = %d\n",cur_app.array_size);
  858. fflush(log_file);
  859. return 0.0;
  860. }
  861. return matr_speedup[type][num_of_cores-1];
  862. */
  863. if (num_of_cores > get_max_cores_count(cur_app)) {
  864. return 0;
  865. } else {
  866. return Exec_Speedup[num_of_cores-2];
  867. }
  868. }
  869. int get_times(app cur_app, int num_of_cores) {
  870. /*
  871. int type;
  872. if (cur_app.array_size == 1024) type = 0;
  873. else if (cur_app.array_size == 2048) type = 1;
  874. else if (cur_app.array_size == 4096) type = 2;
  875. else {
  876. fprintf(log_file, "Unknown array size = %d\n",cur_app.array_size);
  877. fflush(log_file);
  878. return 0.0;
  879. }
  880. return (cur_app.workld * matr_times[type][num_of_cores-2]);
  881. */
  882. return (cur_app.workld * Exec_Latencies[num_of_cores-2]);
  883. }
  884. void matrix_transpose(int n1, float *src, float *dest, int node_id, int myFirst, int myLast, int pad_length){
  885. int i;
  886. int j;
  887. int k;
  888. int l;
  889. int m;
  890. int blksize;
  891. int numblks;
  892. int firstfirst;
  893. int h_off;
  894. int v_off;
  895. int v;
  896. int h;
  897. int n1p;
  898. int row_count;
  899. //fprintf(log_file,"I am inside matrix_transpose-0 node_id is %d n1 %d\n",node_id,n1);
  900. blksize = myLast-myFirst;
  901. numblks = (2*blksize)/num_cache_lines;
  902. if (numblks * num_cache_lines != 2 * blksize) {
  903. numblks ++;
  904. }
  905. blksize = blksize / numblks;
  906. firstfirst = myFirst;
  907. row_count = n1/P;
  908. n1p = n1+pad_length;
  909. for (l=node_id+1;l<P;l++) {
  910. v_off = l*row_count;
  911. for (k=0; k<numblks; k++) {
  912. h_off = firstfirst;
  913. for (m=0; m<numblks; m++) {
  914. for (i=0; i<blksize; i++) {
  915. v = v_off + i;
  916. for (j=0; j<blksize; j++) {
  917. h = h_off + j;
  918. //fprintf(log_file,"Index dest is %d\n",2*(h*n1p+v));
  919. //fprintf(log_file,"Index src is %d\n",2*(v*n1p+h));
  920. //fprintf(log_file,"src = %f\n",src[2*(v*n1p+h)]);
  921. //fprintf(log_file,"src + 1 = %f\n",src[2*(v*n1p+h)+1]);
  922. //fprintf(log_file,"dest = %f\n",dest[2*(h*n1p+v)]);
  923. //fprintf(log_file,"dest + 1 = %f\n",dest[2*(h*n1p+v)+1]);
  924. //fflush(log_file);
  925. dest[2*(h*n1p+v)] = src[2*(v*n1p+h)];
  926. dest[2*(h*n1p+v)+1] = src[2*(v*n1p+h)+1];
  927. //fprintf(log_file,"yolo\n");
  928. }
  929. }
  930. h_off += blksize;
  931. }
  932. v_off+=blksize;
  933. }
  934. }
  935. //fprintf(log_file,"I am inside matrix_transpose-A\n");
  936. for (l=0;l<node_id;l++) {
  937. v_off = l*row_count;
  938. for (k=0; k<numblks; k++) {
  939. h_off = firstfirst;
  940. for (m=0; m<numblks; m++) {
  941. for (i=0; i<blksize; i++) {
  942. v = v_off + i;
  943. for (j=0; j<blksize; j++) {
  944. h = h_off + j;
  945. dest[2*(h*n1p+v)] = src[2*(v*n1p+h)];
  946. dest[2*(h*n1p+v)+1] = src[2*(v*n1p+h)+1];
  947. }
  948. }
  949. h_off += blksize;
  950. }
  951. v_off+=blksize;
  952. }
  953. }
  954. //fprintf(log_file,"I am inside matrix_transpose-B\n");
  955. v_off = node_id*row_count;
  956. for (k=0; k<numblks; k++) {
  957. h_off = firstfirst;
  958. for (m=0; m<numblks; m++) {
  959. for (i=0; i<blksize; i++) {
  960. v = v_off + i;
  961. for (j=0; j<blksize; j++) {
  962. h = h_off + j;
  963. dest[2*(h*n1p+v)] = src[2*(v*n1p+h)];
  964. dest[2*(h*n1p+v)+1] = src[2*(v*n1p+h)+1];
  965. }
  966. }
  967. h_off += blksize;
  968. }
  969. v_off+=blksize;
  970. }
  971. //fprintf(log_file,"I am inside matrix_transpose-C\n");
  972. }
  973. //FFT1D(1, M, N, x_local, trans, upriv, umain2, work_id, lower_bound, upper_bound, pad_length, P);
  974. void FFT1D(int direction, int M, int N, float *x, float *scratch, float *upriv, float *umain2, int node_id, int myFirst, int myLast, int pad_length, int P){
  975. int j, m1, n1;
  976. printf("I am %d and I am inside FFT1D\n",node_id);
  977. //fprintf(log_file,"I am inside FFT1D-A myFirst=%d myLast=%d\n",myFirst,myLast);
  978. m1 = M/2;
  979. n1 = 1 << m1;
  980. matrix_transpose(n1, x, scratch, node_id, myFirst, myLast, pad_length);
  981. //fprintf(log_file,"I am inside FFT1D-B\n");
  982. /* do n1 1D FFTs on columns */
  983. for (j = myFirst; j < myLast; j++){
  984. single_FFT1D(direction, m1, n1, upriv, &scratch[2*j*(n1+pad_length)]);
  985. twiddle_Col(direction, n1, N, j, umain2, &scratch[2*j*(n1+pad_length)],pad_length);
  986. }
  987. //fprintf(log_file,"I am inside FFT1D-C\n");
  988. matrix_transpose(n1, scratch, x, node_id, myFirst, myLast, pad_length);
  989. //fprintf(log_file,"I am inside FFT1D-D\n");
  990. /* do n1 1D FFTs on columns again */
  991. for (j = myFirst; j < myLast; j++) {
  992. single_FFT1D(direction, m1, n1, upriv, &x[2*j*(n1+pad_length)]);
  993. }
  994. //fprintf(log_file,"I am inside FFT1D-E\n");
  995. matrix_transpose(n1, x, scratch, node_id, myFirst, myLast, pad_length);
  996. //fprintf(log_file,"I am inside FFT1D-F\n");
  997. /*for (j = myFirst; j < myLast; j++){
  998. copyColumn(n1, &scratch[2*j*(n1+pad_length)], &x_shared[2*j*(n1+pad_length)]);
  999. }*/
  1000. return;
  1001. }
  1002. void copyColumn(int n1, float *src, float *dest){
  1003. int i;
  1004. for (i = 0; i < n1; i++) {
  1005. dest[2*i] = src[2*i];
  1006. dest[2*i+1] = src[2*i+1];
  1007. }
  1008. }
  1009. void single_FFT1D(int direction, int M, int N, float *u, float *x){
  1010. int j, k, q, L, r, Lstar;
  1011. float *u1, *x1, *x2;
  1012. float omega_r, omega_c, tau_r, tau_c, x_r, x_c;
  1013. reverse(N, M, x);
  1014. for (q=1; q<=M; q++) {
  1015. L = 1<<q; r = N/L; Lstar = L/2;
  1016. u1 = &u[2*(Lstar-1)];
  1017. for (k=0; k<r; k++) {
  1018. x1 = &x[2*(k*L)];
  1019. x2 = &x[2*(k*L+Lstar)];
  1020. for (j=0; j<Lstar; j++) {
  1021. omega_r = u1[2*j];
  1022. omega_c = direction*u1[2*j+1];
  1023. x_r = x2[2*j];
  1024. x_c = x2[2*j+1];
  1025. tau_r = omega_r*x_r - omega_c*x_c;
  1026. tau_c = omega_r*x_c + omega_c*x_r;
  1027. x_r = x1[2*j];
  1028. x_c = x1[2*j+1];
  1029. x2[2*j] = x_r - tau_r;
  1030. x2[2*j+1] = x_c - tau_c;
  1031. x1[2*j] = x_r + tau_r;
  1032. x1[2*j+1] = x_c + tau_c;
  1033. }
  1034. }
  1035. }
  1036. return;
  1037. }
  1038. void twiddle_Col(int direction, int n1, int N, int j, float *u, float *x, int pad_length){
  1039. int i;
  1040. float omega_c, omega_r, x_r, x_c;
  1041. for (i = 0; i < n1; i++) {
  1042. omega_r = u[2*(j*(n1+pad_length)+i)];
  1043. omega_c = direction*u[2*(j*(n1+pad_length)+i)+1];
  1044. x_r = x[2*i];
  1045. x_c = x[2*i+1];
  1046. x[2*i] = omega_r*x_r - omega_c*x_c;
  1047. x[2*i+1] = omega_r*x_c + omega_c*x_r;
  1048. }
  1049. return;
  1050. }
  1051. void reverse(int N, int M, float *x){
  1052. int j, k;
  1053. for (k = 0; k < N; k++){
  1054. j = reverse_bit(M, k);
  1055. if (j > k){
  1056. SWAP(x[2*j], x[2*k]);
  1057. SWAP(x[2*j+1], x[2*k+1]);
  1058. }
  1059. }
  1060. return;
  1061. }
  1062. int reverse_bit(int M, int k){
  1063. int i, j = 0, tmp = k;
  1064. for (i = 0; i < M; i++){
  1065. j = 2*j + (tmp&0x1);
  1066. tmp = tmp >> 1;
  1067. }
  1068. return j;
  1069. }