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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 3G开发 > 查看源码
ReedArith.bsv
资源名称：mit-ofdm-wifi.rar　[点击查看]
上传用户：aoptech
上传日期：2014-09-22
资源大小：784k
文件大小：15k
源码类别：

3G开发

开发平台：

Others

ReedArith.bsv：源码内容
  1. //----------------------------------------------------------------------//
  2. // The MIT License 
  3. // 
  4. // Copyright (c) 2007 Alfred Man Cheuk Ng, mcn02@mit.edu 
  5. // 
  6. // Permission is hereby granted, free of charge, to any person 
  7. // obtaining a copy of this software and associated documentation 
  8. // files (the "Software"), to deal in the Software without 
  9. // restriction, including without limitation the rights to use,
  10. // copy, modify, merge, publish, distribute, sublicense, and/or sell
  11. // copies of the Software, and to permit persons to whom the
  12. // Software is furnished to do so, subject to the following conditions:
  13. // 
  14. // The above copyright notice and this permission notice shall be
  15. // included in all copies or substantial portions of the Software.
  16. // 
  17. // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
  18. // EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES
  19. // OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
  20. // NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
  21. // HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY,
  22. // WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
  23. // FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
  24. // OTHER DEALINGS IN THE SOFTWARE.
  25. //----------------------------------------------------------------------//
  27. //**********************************************************************
  28. // Galois field arithmetic
  29. //----------------------------------------------------------------------
  30. // $Id: Arith.bsv
  31. //
  33. import ReedTypes::*;  
  34. import Vector::*;
  36. // -----------------------------------------------------------
  37. // The primitive polynomial defines the Galois field in which 
  38. // Reed-Solomon decoder operates, and all the following   
  39. // arithmetic operations are defined under.  Changing this    
  40. // value cause the whole Reed-Solomon decoder to operate   
  41. // under the new primitive polynomial.
  42. // primitive_polynomial[i] = Coefficient of x**i for i = 0:7
  44. // -----------------------------------------------------------
  45. Byte    n_param              = 8'd255;
  46. Byte    t_param              = 8'd16;
  48. // -----------------------------------------------------------
  49. //(* noinline *) 
  50. function Byte gf_mult (Polynomial primitive_polynomial, Byte left, Byte right);
  51.    
  52.    Bit#(15) first  = 15'b0;
  53.    Bit#(15) result = 15'b0;
  54.    
  55.    // this function bring back higher degree values back to the field
  56.    function Byte getAddSymbol(Bool ctrl, Integer shift);
  57.       Byte shiftPoly = primitive_polynomial << shift; 
  58.       Byte addSymbol = ctrl ? shiftPoly : 0; 
  59.       return addSymbol;
  60.    endfunction
  61.    
  62.   for (Integer i = 0; i < 8; i = i + 1)
  63.      for (Integer j = 0; j < 8 ; j = j + 1)
  64.         begin
  65.            if (first[i+j] == 0) // initialize result[i+j]
  66.               result[i+j] = (left[i] & right[j]);
  67.            else                 // accumulate
  68.               result[i+j] = result[i+j] ^ (left[i] & right[j]);
  69.            first[i+j] = 1; // only initialize each signal once 
  70.         end
  71.    
  72.   Vector#(7,Bool) ctrlVec = unpack(result[14:8]);
  73.   Byte finalResult = fold( ^ ,cons(result[7:0],zipWith(getAddSymbol,ctrlVec,genVector)));
  74.    
  75.   return finalResult;
  77. //    result [0] = result [0] ^ (left [0] & right [0]);
  78. //    result [1] = result [1] ^ (left [1] & right [0]) ^ (left [0] & right [1]);
  79. //    result [2] = result [2] ^ (left [2] & right [0]) ^ (left [1] & right [1]) ^ (left [0] & right [2]);
  80. //    result [3] = result [3] ^ (left [3] & right [0]) ^ (left [2] & right [1]) ^ (left [1] & right [2]) ^ (left [0] & right [3]);
  81. //    result [4] = result [4] ^ (left [4] & right [0]) ^ (left [3] & right [1]) ^ (left [2] & right [2]) ^ (left [1] & right [3]) ^ (left [0] & right [4]);
  82. //    result [5] = result [5] ^ (left [5] & right [0]) ^ (left [4] & right [1]) ^ (left [3] & right [2]) ^ (left [2] & right [3]) ^ (left [1] & right [4]) ^ (left [0] & right [5]);
  83. //    result [6] = result [6] ^ (left [6] & right [0]) ^ (left [5] & right [1]) ^ (left [4] & right [2]) ^ (left [3] & right [3]) ^ (left [2] & right [4]) ^ (left [1] & right [5]) ^ (left [0] & right [6]);
  84. //    result [7] = result [7] ^ (left [7] & right [0]) ^ (left [6] & right [1]) ^ (left [5] & right [2]) ^ (left [4] & right [3]) ^ (left [3] & right [4]) ^ (left [2] & right [5]) ^ (left [1] & right [6]) ^ (left [0] & right [7]);
  85. //    result [8] = result [8] ^ (left [7] & right [1]) ^ (left [6] & right [2]) ^ (left [5] & right [3]) ^ (left [4] & right [4]) ^ (left [3] & right [5]) ^ (left [2] & right [6]) ^ (left [1] & right [7]);
  86. //    result [9] = result [9] ^ (left [7] & right [2]) ^ (left [6] & right [3]) ^ (left [5] & right [4]) ^ (left [4] & right [5]) ^ (left [3] & right [6]) ^ (left [2] & right [7]);
  87. //    result [10] = result [10] ^ (left [7] & right [3]) ^ (left [6] & right [4]) ^ (left [5] & right [5]) ^ (left [4] & right [6]) ^ (left [3] & right [7]);
  88. //    result [11] = result [11] ^ (left [7] & right [4]) ^ (left [6] & right [5]) ^ (left [5] & right [6]) ^ (left [4] & right [7]);
  89. //    result [12] = result [12] ^ (left [7] & right [5]) ^ (left [6] & right [6]) ^ (left [5] & right [7]);
  90. //    result [13] = result [13] ^ (left [7] & right [6]) ^ (left [6] & right [7]);
  91. //    result [14] = result [14] ^ (left [7] & right [7]);
  93. //    if (result [14] == 1'b1)
  94. //       result = result ^ ((zeroExtend (primitive_polynomial)) << (6));
  95. //    if (result [13] == 1'b1)
  96. //       result = result ^ ((zeroExtend (primitive_polynomial)) << (5));
  97. //    if (result [12] == 1'b1)
  98. //       result = result ^ ((zeroExtend (primitive_polynomial)) << (4));
  99. //    if (result [11] == 1'b1)
  100. //       result = result ^ ((zeroExtend (primitive_polynomial)) << (3));
  101. //    if (result [10] == 1'b1)
  102. //       result = result ^ ((zeroExtend (primitive_polynomial)) << (2));
  103. //    if (result [9] == 1'b1)
  104. //       result = result ^ ((zeroExtend (primitive_polynomial)) << (1));
  105. //    if (result [8] == 1'b1)
  106. //       result = result ^ ((zeroExtend (primitive_polynomial)) << (0));
  108. //   return (result [7:0]);
  110. endfunction
  112. // -----------------------------------------------------------
  113. function Byte gf_add (Byte left, Byte right);
  114.    return (left ^ right);
  115. endfunction
  117. // -----------------------------------------------------------
  118. //(* noinline *) 
  119. function Byte alpha_n (Polynomial primitive_polynomial, Byte n);
  120.    return times_alpha_n(primitive_polynomial, 1,n);
  121. endfunction
  123. // -----------------------------------------------------------
  124. //(* noinline *) 
  125. function Byte times_alpha_n (Polynomial primitive_polynomial, Byte a, Byte n);
  126. //    Byte b=a;
  127. //    for (Byte i = 0; i < n; i = i + 1)
  128. //       b=times_alpha(primitive_polynomial, b);
  129. //    return b;
  131.    Byte multVal = 1 << n;
  132.    return gf_mult(primitive_polynomial,a,multVal);
  133. endfunction
  135. // -----------------------------------------------------------
  136. //(* noinline *) 
  137. function Byte times_alpha (Polynomial primitive_polynomial, Byte a);
  138.    return gf_mult(primitive_polynomial, a, 2);
  139. //   return (a<<1)^({a[7],a[7],a[7],a[7],a[7],a[7],a[7],a[7]} & primitive_polynomial);
  140. endfunction
  142. // // -----------------------------------------------------------
  143. // function Byte gf_inv (Byte a);
  145. //    case (a) matches
  146. //         0 : return         2;
  147. //         1 : return         1;
  148. //         2 : return       142;
  149. //         3 : return       244;
  150. //         4 : return        71;
  151. //         5 : return       167;
  152. //         6 : return       122;
  153. //         7 : return       186;
  154. //         8 : return       173;
  155. //         9 : return       157;
  156. //        10 : return       221;
  157. //        11 : return       152;
  158. //        12 : return        61;
  159. //        13 : return       170;
  160. //        14 : return        93;
  161. //        15 : return       150;
  162. //        16 : return       216;
  163. //        17 : return       114;
  164. //        18 : return       192;
  165. //        19 : return        88;
  166. //        20 : return       224;
  167. //        21 : return        62;
  168. //        22 : return        76;
  169. //        23 : return       102;
  170. //        24 : return       144;
  171. //        25 : return       222;
  172. //        26 : return        85;
  173. //        27 : return       128;
  174. //        28 : return       160;
  175. //        29 : return       131;
  176. //        30 : return        75;
  177. //        31 : return        42;
  178. //        32 : return       108;
  179. //        33 : return       237;
  180. //        34 : return        57;
  181. //        35 : return        81;
  182. //        36 : return        96;
  183. //        37 : return        86;
  184. //        38 : return        44;
  185. //        39 : return       138;
  186. //        40 : return       112;
  187. //        41 : return       208;
  188. //        42 : return        31;
  189. //        43 : return        74;
  190. //        44 : return        38;
  191. //        45 : return       139;
  192. //        46 : return        51;
  193. //        47 : return       110;
  194. //        48 : return        72;
  195. //        49 : return       137;
  196. //        50 : return       111;
  197. //        51 : return        46;
  198. //        52 : return       164;
  199. //        53 : return       195;
  200. //        54 : return        64;
  201. //        55 : return        94;
  202. //        56 : return        80;
  203. //        57 : return        34;
  204. //        58 : return       207;
  205. //        59 : return       169;
  206. //        60 : return       171;
  207. //        61 : return        12;
  208. //        62 : return        21;
  209. //        63 : return       225;
  210. //        64 : return        54;
  211. //        65 : return        95;
  212. //        66 : return       248;
  213. //        67 : return       213;
  214. //        68 : return       146;
  215. //        69 : return        78;
  216. //        70 : return       166;
  217. //        71 : return         4;
  218. //        72 : return        48;
  219. //        73 : return       136;
  220. //        74 : return        43;
  221. //        75 : return        30;
  222. //        76 : return        22;
  223. //        77 : return       103;
  224. //        78 : return        69;
  225. //        79 : return       147;
  226. //        80 : return        56;
  227. //        81 : return        35;
  228. //        82 : return       104;
  229. //        83 : return       140;
  230. //        84 : return       129;
  231. //        85 : return        26;
  232. //        86 : return        37;
  233. //        87 : return        97;
  234. //        88 : return        19;
  235. //        89 : return       193;
  236. //        90 : return       203;
  237. //        91 : return        99;
  238. //        92 : return       151;
  239. //        93 : return        14;
  240. //        94 : return        55;
  241. //        95 : return        65;
  242. //        96 : return        36;
  243. //        97 : return        87;
  244. //        98 : return       202;
  245. //        99 : return        91;
  246. //       100 : return       185;
  247. //       101 : return       196;
  248. //       102 : return        23;
  249. //       103 : return        77;
  250. //       104 : return        82;
  251. //       105 : return       141;
  252. //       106 : return       239;
  253. //       107 : return       179;
  254. //       108 : return        32;
  255. //       109 : return       236;
  256. //       110 : return        47;
  257. //       111 : return        50;
  258. //       112 : return        40;
  259. //       113 : return       209;
  260. //       114 : return        17;
  261. //       115 : return       217;
  262. //       116 : return       233;
  263. //       117 : return       251;
  264. //       118 : return       218;
  265. //       119 : return       121;
  266. //       120 : return       219;
  267. //       121 : return       119;
  268. //       122 : return         6;
  269. //       123 : return       187;
  270. //       124 : return       132;
  271. //       125 : return       205;
  272. //       126 : return       254;
  273. //       127 : return       252;
  274. //       128 : return        27;
  275. //       129 : return        84;
  276. //       130 : return       161;
  277. //       131 : return        29;
  278. //       132 : return       124;
  279. //       133 : return       204;
  280. //       134 : return       228;
  281. //       135 : return       176;
  282. //       136 : return        73;
  283. //       137 : return        49;
  284. //       138 : return        39;
  285. //       139 : return        45;
  286. //       140 : return        83;
  287. //       141 : return       105;
  288. //       142 : return         2;
  289. //       143 : return       245;
  290. //       144 : return        24;
  291. //       145 : return       223;
  292. //       146 : return        68;
  293. //       147 : return        79;
  294. //       148 : return       155;
  295. //       149 : return       188;
  296. //       150 : return        15;
  297. //       151 : return        92;
  298. //       152 : return        11;
  299. //       153 : return       220;
  300. //       154 : return       189;
  301. //       155 : return       148;
  302. //       156 : return       172;
  303. //       157 : return         9;
  304. //       158 : return       199;
  305. //       159 : return       162;
  306. //       160 : return        28;
  307. //       161 : return       130;
  308. //       162 : return       159;
  309. //       163 : return       198;
  310. //       164 : return        52;
  311. //       165 : return       194;
  312. //       166 : return        70;
  313. //       167 : return         5;
  314. //       168 : return       206;
  315. //       169 : return        59;
  316. //       170 : return        13;
  317. //       171 : return        60;
  318. //       172 : return       156;
  319. //       173 : return         8;
  320. //       174 : return       190;
  321. //       175 : return       183;
  322. //       176 : return       135;
  323. //       177 : return       229;
  324. //       178 : return       238;
  325. //       179 : return       107;
  326. //       180 : return       235;
  327. //       181 : return       242;
  328. //       182 : return       191;
  329. //       183 : return       175;
  330. //       184 : return       197;
  331. //       185 : return       100;
  332. //       186 : return         7;
  333. //       187 : return       123;
  334. //       188 : return       149;
  335. //       189 : return       154;
  336. //       190 : return       174;
  337. //       191 : return       182;
  338. //       192 : return        18;
  339. //       193 : return        89;
  340. //       194 : return       165;
  341. //       195 : return        53;
  342. //       196 : return       101;
  343. //       197 : return       184;
  344. //       198 : return       163;
  345. //       199 : return       158;
  346. //       200 : return       210;
  347. //       201 : return       247;
  348. //       202 : return        98;
  349. //       203 : return        90;
  350. //       204 : return       133;
  351. //       205 : return       125;
  352. //       206 : return       168;
  353. //       207 : return        58;
  354. //       208 : return        41;
  355. //       209 : return       113;
  356. //       210 : return       200;
  357. //       211 : return       246;
  358. //       212 : return       249;
  359. //       213 : return        67;
  360. //       214 : return       215;
  361. //       215 : return       214;
  362. //       216 : return        16;
  363. //       217 : return       115;
  364. //       218 : return       118;
  365. //       219 : return       120;
  366. //       220 : return       153;
  367. //       221 : return        10;
  368. //       222 : return        25;
  369. //       223 : return       145;
  370. //       224 : return        20;
  371. //       225 : return        63;
  372. //       226 : return       230;
  373. //       227 : return       240;
  374. //       228 : return       134;
  375. //       229 : return       177;
  376. //       230 : return       226;
  377. //       231 : return       241;
  378. //       232 : return       250;
  379. //       233 : return       116;
  380. //       234 : return       243;
  381. //       235 : return       180;
  382. //       236 : return       109;
  383. //       237 : return        33;
  384. //       238 : return       178;
  385. //       239 : return       106;
  386. //       240 : return       227;
  387. //       241 : return       231;
  388. //       242 : return       181;
  389. //       243 : return       234;
  390. //       244 : return         3;
  391. //       245 : return       143;
  392. //       246 : return       211;
  393. //       247 : return       201;
  394. //       248 : return        66;
  395. //       249 : return       212;
  396. //       250 : return       232;
  397. //       251 : return       117;
  398. //       252 : return       127;
  399. //       253 : return       255;
  400. //       254 : return       126;
  401. //       255 : return       253;
  402. //    endcase
  403. // endfunction