整数开方运算的DSP实现

开方运算的DSP实现

(2014-04-07 10:59:36)

转载�

标签： c dsp

分类：DSP

1. //===============================================

3. //函数名：VSqrt3

4. //功能： 实现对32位定点数的开方

5. //性能： 60M主频28015硬件下运行时间小于10us

6. //转自：http://read.pudn.com/downloads180/sourcecode/mpi/840129/sqrt.c__.htm

8. unsigned longVSqrt3(unsigned longx )

9. {

10. unsigned longx1;

11. ints=1;

12. unsigned longg0,g1;

14. if(x<=1) returnx;

16. x1=x-1;

17. if(x1>65535)

18. {

19. s+=8;

20. x1>>=16;

21. }

22. if(x1>255)

23. {

24. s+=4;

25. x1>>=8;

26. }

27. if(x1>15)

28. {

29. s+=2;

30. x1>>=4;

31. }

32. if(x1>3)

33. {

34. s+=1;

35. }

36. g0=1;

37. g0=g0<<s;

38. g1 =(g0 +(x>>s))>>1;

40. while(g1

41. {

42. g0=g1;

43. g1=(g0+x/g0)>>1;

44. }

45. returng0;

46. }

48. //===============================================

49.