ComplexLibrary.bsv
上传用户:aoptech
上传日期:2014-09-22
资源大小:784k
文件大小:3k
源码类别:

3G开发

开发平台:

Others

  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. //----------------------------------------------------------------------//
  26. //////////////////////////////////////////////////////////
  27. // Some useful functions for Complex Type
  28. // Author: Alfred Man C Ng 
  29. // Email: mcn02@mit.edu
  30. // Data: 9-29-2006
  31. /////////////////////////////////////////////////////////
  32. import Complex::*;
  33. // complex conjugate
  34. function Complex#(a) cmplxConj(Complex#(a) x) provisos (Arith#(a));
  35.     return cmplx(x.rel, negate(x.img));
  36. endfunction // Complex
  37. instance Bounded#(Complex#(a)) provisos(Bounded#(a));
  38.     function Complex#(a) minBound();
  39.         a min = minBound;
  40.         return cmplx(min,min);
  41.     endfunction // Complex
  42.     function Complex#(a) maxBound();
  43.         a max = maxBound;
  44.         return cmplx(max,max);
  45.     endfunction // Complex
  46. endinstance
  47. // for complex single bit multiply
  48. function Complex#(Bit#(rsz)) cmplxSignExtend(Complex#(Bit#(asz)) a)
  49.   provisos (Add#(xxA,asz,rsz));
  50.       let rel = signExtend(a.rel);
  51.       let img = signExtend(a.img);
  52.       return cmplx(rel, img);
  53. endfunction // Complex
  54. // for complex modulus = rel^2 + img^2, ri = 2ai + 1, rf = 2af
  55. function Bit#(ri)  cmplxModSq(Complex#(Bit#(ai)) a)
  56.   provisos (Add#(ai,ai,ci), Add#(1,ci,ri), Add#(xxA,ai,ri));
  57.       return ((signExtend(a.rel) * signExtend(a.rel))  + (signExtend(a.img) * signExtend(a.img)));
  58. endfunction // FixedPoint
  59.