Vdsp(bf561)中的浮点运算（2）：float的疑问

最新推荐文章于 2021-05-17 14:24:06 发布

iteye_18688

最新推荐文章于 2021-05-17 14:24:06 发布

阅读量105

点赞数

文章标签： F# Office Microsoft FP EXT

<?xml:namespace prefix = o ns = "urn:schemas-microsoft-com:office:office" />

快乐虾

http://blog.csdn.net/lights_joy/

lights@hb165.com

本文适用于

ADI BF561 DSP

优视BF561EVB开发板

Visual DSP++ 5.0 (update 5)

欢迎转载，但请保留作者信息

写一行很简单的C代码：

float a = 1234.56;

用vdsp编译后的汇编代码为：

R0 = 20972 ( X ) ;

R0.H = 17562 ;

[ FP + 0x10 ] = R0 ;

有点看不懂，呵呵，R0的值转换为十六进制就是0x<?xml:namespace prefix = st1 ns = "urn:schemas-microsoft-com:office:smarttags" /><chmetcnv w:st="on" tcsc="0" numbertype="1" negative="False" hasspace="False" sourcevalue="449" unitname="a">449A</chmetcnv>51EC。根据vdsp文档的说法，其单精度浮点数格式为：

<?xml:namespace prefix = v ns = "urn:schemas-microsoft-com:vml" /><shapetype id="_x0000_t75" stroked="f" filled="f" path="m@4@5l@4@11@9@11@9@5xe" o:preferrelative="t" o:spt="75" coordsize="21600,21600"><stroke joinstyle="miter"></stroke><formulas><f eqn="if lineDrawn pixelLineWidth 0"></f><f eqn="sum @0 1 0"></f><f eqn="sum 0 0 @1"></f><f eqn="prod @2 1 2"></f><f eqn="prod @3 21600 pixelWidth"></f><f eqn="prod @3 21600 pixelHeight"></f><f eqn="sum @0 0 1"></f><f eqn="prod @6 1 2"></f><f eqn="prod @7 21600 pixelWidth"></f><f eqn="sum @8 21600 0"></f><f eqn="prod @7 21600 pixelHeight"></f><f eqn="sum @10 21600 0"></f></formulas><path o:connecttype="rect" gradientshapeok="t" o:extrusionok="f"></path><lock aspectratio="t" v:ext="edit"></lock></shapetype><shape id="_x0000_i1025" style="WIDTH: 354pt; HEIGHT: 97.5pt" type="#_x0000_t75"><imagedata o:title="" src="file:///F:%5Cdevtemp%5Cmsohtml1%5C01%5Cclip_image001.png"></imagedata></shape>

浮点数计算公式：

将0x<chmetcnv w:st="on" tcsc="0" numbertype="1" negative="False" hasspace="False" sourcevalue="449" unitname="a">449A</chmetcnv>51EC转换为二进制：

0100 0100 1001 1010 0101 0001 1110 1100

可得：

Sign = 0

Mantissa = 001 1010 0101 0001 1110 1100

Exponent = 1000 1001

按照公式

将Exponent转换为十进制，其值为137。

将Mantissa转换为十进制，其值为：

2^-3 + 2^-4 + 2^-6 + 2^-9 + 2 ^-11 + 2^-15 + 2^-16 + 2^-17 + 2^-18 + 2^-20 + 2^{-21 =}

0.125 + 0.0625 + 0.015625 + 0.001953125 + 0.00048828125 + 0.000030517578125 + 0.0000152587890625 + 0.00000762939453125 + 0.000003814697265625 + 0.00000095367431640625 + 0.000000476837158203125 = 0.205625057220458984375

代入浮点数计算公式：

(-1)⁰ + 1. 205625057220458984375 * 2 ^(137-127)

= 1234.56005859375

那么编译器又是如何将1234.56转换为0x<chmetcnv w:st="on" tcsc="0" numbertype="1" negative="False" hasspace="False" sourcevalue="449" unitname="a">449A</chmetcnv>51EC的呢？

首先转换整数部分，用2除，取余数，其结果为：

1234 / 2 = 617 ….. 0

617 / 2 = 308 …….1

308 / 2 = 154 …….0

154 / 2 = 77 ……...0

77 / 2 = 38 ……….1

38 / 2 = 19 ……….0

19/ 2 = 9 ………….1

9 / 2 = 4………… .1

4/ 2 = 2 …………..0

2/2 = 1……………0

1/2 = 0……………1