浅谈sar指令在整数除法中的优化

最新推荐文章于 2021-05-04 10:52:14 发布
最新推荐文章于 2021-05-04 10:52:14 发布
阅读量911 收藏 1
点赞数 1
分类专栏： 常用
标 题:  【转载】浅谈sar指令在整数除法中的优化


初学《加密与解密3》看到第四章里的编译器的算法优化时的一点困惑，然后动手实践了一下得到一点心得体会，拿来跟大家分享！
   ps.本人汇编基础非常差，sar指令也是刚刚学会的，另外脑袋转的慢，别人用5分钟能解决的问题，我得用10分钟，因此这篇文章我花了6个多小时才完成。
如果有错误欢迎指出，以免误导他人！

浅谈Sar指令在整数除法中的优化

还不懂sar指令的同学马上去baidu啊，Google啊查去，我刚查完！
首先来看一段vc代码
代码: 
int get()
{
  return -8;
}

int main()
{
  int a = get()/2;
  printf("%d",a);
  return 0;
}
汇编代码（本文里反汇编码都是指Release版本下的OD返回编码）如下
代码: 
00401000  /$  B8 F8FFFFFF   mov     eax, -8   ; get子程序 总是返回-8
00401005  \.  C3            retn
00401006      90            nop
00401010  /$  E8 EBFFFFFF   call    00401000  ;得到返回值 eax = -8
00401015  |.  99            cdq                ;将eax的值符合扩展到edx
00401016  |.  2BC2          sub     eax, edx   ;等价于 eax = eax + 符号位（eax为正数或者0时符号位为0，负数时为1）
00401018  |.  D1F8          sar     eax, 1      ;带符合右移一位
0040101A  |.  50            push    eax
0040101B  |.  68 30704000   push    00407030      ;  ASCII "%d"
00401020  |.  E8 0B000000   call    00401030       ;printf函数
00401025  |.  83C4 08       add     esp, 8
00401028  |.  33C0          xor     eax, eax
0040102A  \.  C3            retn
0040102B      90            nop
简单的伪代码表示如下
代码: 
Mov eax, idata              ;idata表示一个整数
Add eax, eax的符号位
Sar eax, 1
上面是一个整数/2的情况 
如果一个整数/2^n(这里2^n表示2的n次方) 应该如何表示呢
我们将上面的VC代码修改这么一句
int a = get()/2; =>   int a = get()/16;
通过反汇编自己翻译成伪代码就是这样的
代码: 
Mov eax, idata              ;idata表示一个整数
Add eax, (eax的符号位)*(2^n - 1) ;这里n=4
Sar eax, n
上面这个就是一般的除数是2^n一般的公式了，正负整数都管用

粗略的用这样一个数学公式来表示：
Idata/(2^n) = YWn(idata+2^n - 1)
注意 YWn表示将其后面的数以sar方式移动n位 即 sar (数),n

再变换一下，将idata换成(idata – 2^n + 1)
(Idata – 2^n+1)/(2^n) = YWn(idata)             -----------最终公式

Ps.上面的公式是我OD反汇编总结出来的，缺少数学的严格推理证明。我不会证，谁来证明下啊？感激不尽！

下面来解释一个难点，也是写本文的初衷
看一段VC代码：
代码: 
int get()
{
  int i=rand()%2;
   if(i)
  return rand();
   else
   return rand()*(-1);
}

int main()
{
  int a = get()/9;
  printf("%d",a);
  return 0;
}
反汇编码：
代码: 
00401010  /$  E8 EBFFFFFF   call    00401000   ;调用get函数 eax=随即数
00401015  |.  8BC8          mov     ecx, eax  ; ecx = eax
00401017  |.  B8 398EE338   mov     eax, 38E38E39 ;这是一个编译器预先设定好的MagicNumber
0040101C  |.  F7E9          imul    ecx ;带符合的乘法，不多说
0040101E  |.  D1FA          sar     edx, 1 ;将乘法的高32位结果算数右移一位
00401020  |.  8BC2          mov     eax, edx ;eax = edx
00401022  |.  C1E8 1F       shr     eax, 1F   ;逻辑右移31位
00401025  |.  03D0          add     edx, eax ;以上这3句的意思就是  edx = edx+edx的符号位，本文最重要的目的就是解释他为什么加1
00401027  |.  52            push    edx
00401028  |.  68 30704000   push    00407030              ;  ASCII "%d"
0040102D  |.  E8 0E000000   call    00401040  ;printf edx
00401032  |.  83C4 08       add     esp, 8
00401035  |.  33C0          xor     eax, eax
00401037  \.  C3            retn
关于MagicNumber这里给你几个连接自己看，不多介绍
http://bbs.pediy.com/showthread.php?t=68849
http://bbs.pediy.com/showthread.php?t=81456
我们把上面的反汇编代码写成伪代码，如下
代码: 
idata 表示一个32位有符号数，这里一定是有符号的，否则逻辑错误
magic 表示立即数0x38E38E39  他是一个正数
FH(idata) 表示idata的最高位，即 FH(idtat) = (idata & 0x80000000)? 1:0
[edx,eax]表示 edx和eax组成的64位的寄存器

Mov [edx,eax], magic*idata
Sar [edx,eax],33      ;由于移位的操作此时eax里面的值已经没有意义了
Add edx + idata的符号位  ;因为magic是个正数，所以等价于magic*idata的符号位
而我们的本意是 idata/9
因此得到如下数学公式，我们就是要证明他的正确性
YW33(magic*idata) + FH(idata) = idata/9          -------------- (黄金公式)
①.  当FH(idata) == 0时
上述黄金公式即
Magic*idata/2^33 = idata/9
即
Magic*9 = 2^33
我们将magic = 0x0x38E38E39带入  然后用科学计算器验证了下结果如下
Magic*9 = 0x200000001 
2^33 = 0x200000000
这种优化的除法本身就是一种近似，这里不多讨论为什么存在一个1的差异
证毕。
②.  当FH(idata) == 1时
由 最终公式(Idata – 2^n+1)/(2^n) = YWn(idata)
得到
YW33(magic*idata) = (magic*idata – 2^33 + 1)/(2^33)
将右边带入黄金公式即：
(magic*idata – 2^33 + 1)/(2^33) + 1 = idata/9
即
Magic*idata/2^33 – 1 + 1 = idata/9          (;1/2^33 ≈ 0因此忽略掉)
即
Magic*idata/2^33 = idata/9
即
Magic*9 = 2^33

























浅谈Sar指令在整数除法中的优化

第一部分

首先来看一段vc代码

int get()

{

    return-8;

}

int main()

{

    inta = get()/2;

    printf("%d",a);

    return0;

}

汇编代码（本文里反汇编码都是指Release版本下的OD返回编码）如下

00401000  /$  B8 F8FFFFFF  mov     eax, -8   ; get子程序 总是返回-8

00401005  \. C3            retn

00401006      90            nop

00401010  /$  E8EBFFFFFF   call    00401000  ;得到返回值 eax = -8

00401015  |. 99            cdq              ;将eax的值符合扩展到edx

00401016  |. 2BC2          sub     eax, edx  ;等价于 eax = eax + 符号位（eax为正数或者0时符号位为0，负数时为1）

00401018  |.  D1F8         sar     eax, 1      ;带符合右移一位

0040101A  |. 50            push    eax

0040101B  |.  6830704000   push    00407030     ;  ASCII "%d"

00401020  |.  E80B000000   call    00401030       ;printf函数

00401025  |.  83C4 08      add     esp, 8

00401028  |.  33C0         xor     eax, eax

0040102A  \. C3            retn

0040102B      90            nop

 

简单的伪代码表示如下

moveax, idata              ;idata表示一个整数

addeax, eax的符号位

sareax, 1

 

上面是一个整数/2的情况

如果一个整数/2^n(这里2^n表示2的n次方) 应该如何表示呢

我们将上面的VC代码修改这么一句

inta = get()/2; => int a = get()/16;

通过反汇编自己翻译成伪代码就是这样的

moveax, idata              ;idata表示一个整数

addeax, (eax的符号位)*(2^n - 1) ;这里n=4

sareax, n

上面这个就是一般的除数是2^n一般的公式了，正负整数都管用

 

粗略的用这样一个数学公式来表示：

idata/(2^n) = YWn(idata+2^n - 1)

注意 YWn表示将其后面的数以sar方式移动n位即 sar (数),n

 

再变换一下，将idata换成(idata – 2^n + 1)

(idata – 2^n+1)/(2^n) = YWn(idata) -------  (最终公式)

 

ps.上面的公式是我OD反汇编总结出来的，缺少数学的严格推理证明。我不会证，谁来证明下啊？感激不尽！

注：(以上2个公式仅对负数管用，正数的公式是idata/(2^n)= YWn(idata) )

 

第二部分

再看一段VC代码：

int get()

{

    inti=rand()%2;

  if(i)

    returnrand();

  else

  return rand()*(-1);

}

int main()

{

    inta = get()/9;

    printf("%d",a);

    return0;

}

 

反汇编码：

00401010  /$  E8EBFFFFFF   call    00401000  ;调用get函数 eax=随即数

00401015  |. 8BC8          mov     ecx, eax ; ecx = eax

00401017  |.  B8398EE338   mov     eax,38E38E39 ;这是一个编译器预先设定好的MagicNumber

0040101C  |.  F7E9         imul    ecx       ;带符合的乘法，不多说

0040101E  |. D1FA          sar     edx, 1 ;将乘法的高32位结果算数右移一位

00401020  |. 8BC2          mov     eax, edx ;eax =edx

00401022  |. C1E8 1F       shr    eax, 1F   ;逻辑右移31位

00401025  |. 03D0          add     edx, eax ;以上这3句的意思就是  edx = edx+edx的符号位，本文最重要的目的就是解释他为什么加1

00401027  |. 52            push    edx

00401028  |.  6830704000   push    00407030  ;  ASCII"%d"

0040102D  |.  E80E000000   call    00401040 ;printf edx

00401032  |.  83C4 08       add    esp, 8

00401035  |.  33C0          xor     eax, eax

00401037  \. C3            retn

 

关于MagicNumber这里给你几个连接自己看，不多介绍

http://bbs.pediy.com/showthread.php?t=68849

http://bbs.pediy.com/showthread.php?t=81456

我们把上面的反汇编代码写成伪代码，如下

idata表示一个32位有符号数，这里一定是有符号的，否则逻辑错误

magic表示立即数0x38E38E39  他是一个正数

FH(idata)表示idata的最高位，即 FH(idtat) = (idata & 0x80000000)? 1:0

[edx,eax]表示 edx和eax组成的64位的寄存器

 

mov[edx,eax], magic*idata

sar[edx,eax],33 ; 移位的操作执行完本行指令后eax里面的值就没有意义了

addedx + idata的符号位 ;因为magic是个正数，所以等价于magic*idata的符号位

 

而我们的本意是 idata/9

因此得到如下数学公式，我们就是要证明他的正确性

YW33(magic*idata) + FH(idata) = idata/9 ---------- (黄金公式)

①.当FH(idata) == 0时

上述黄金公式即

magic*idata/2^33 = idata/9

即

magic*9 = 2^33

我们将magic = 0x0x38E38E39带入  然后用科学计算器验证了下结果如下

magic*9 = 0x200000001

2^33 = 0x200000000

这种优化的除法本身就是一种近似，这里不多讨论为什么存在一个1的差异

证毕。

②.当FH(idata) == 1时

由 最终公式 (Idata – 2^n+1)/(2^n) = YWn(idata)

得到

YW33(magic*idata) = (magic*idata – 2^33+ 1)/(2^33)

将右边带入黄金公式即：

(magic*idata – 2^33 + 1)/(2^33) + 1 =idata/9

即

Magic*idata/2^33 – 1 + 1 = idata/9       (;1/2^33 ≈ 0因此忽略掉)

即

magic*idata/2^33 = idata/9

即

magic*9 = 2^33

证毕。

 

 


对于公式Idata/(2^n)=YWn(idata+2^n-1)
有一点需要说明，即intel使用的有符号除法的商是采取向零圆整的方式(如：(-5)/2=-2……-1，5/2=2……1)，所以对于不能整除的数，移位的结果与商差1。这就是移位前加上(2^n-1)的原因。


假设负数idata的补码二进制表示为
bbbbb……b
    ______
      ↑
      n位


(2^n-1)二进制表示为
00001……1
    ______
      ↑
      n位


如果idata能够整除(2^n)，则其二进制补码末尾n位应该都是零，此时加上(2^n-1)再移位与直接移位结果相同，结果等于商
如果idata不能整除(2^n),则末尾n位中必有不是零的部分，此时加上(2^n-1)再移位相当于将移位结果加1，这样就与商相等了。


tianjun7
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
SAR雷达Yamaguchi四分量分解MATLAB代码（注释部分为极化角优化）
weixin_34161032的博客
04-12 706
function yamachuchi_matrix = yamaguchi(C) %对于给定的C矩阵做yamaguchi分解，返回结果是个各个散射的能量是一个row X col X 4的矩阵 %P_s（单词散射）,P_d（二面角）,P_v（体散射）,P_c（旋状体） row=size(C.C11,1); col=size(C.C11,2); i=sqrt(-1); yamachu...
数学形态学在SAR图像增强中的应用
01-19
与光学图像相比，SAR更容易辨别地面的伪装目标和隐藏于山林中的感兴趣目标，因此，它在国民经济和国防建设中有着及其重要的地位。 　自从SAR技术诞生以来，在军用和民用方面得到广泛的应用。 　合成孔径雷达由于...
sar命令使用详解
congjiu2607的博客
09-27 194
sar（System ActivityReporter系统活动情况报告）是目前Linux上最为全面的系统性能分析工具之一，可以从多方面对系统的活动进行报告，包括：文件的读写情况、系统调用的使用情况、磁盘I/O、CPU效率、内存使...
除法-cdq和移位(除数为2的幂)
qq_38184895的博客
01-25 255
c++ printf("nVarone / 8 = %d", nVarone / 8); 汇编代码 mov eax，dword ptr[ebp-4]] cdq and edx,7 //edx为高位 add eax.edx //eax+(2^n-1),n=3 sar eax,3 push... call...(printf函数
CS:APP二进制炸弹phase4
astrotycoon
06-14 6779
CS:APP二进制炸弹phase4
钱晓捷《32位汇编语言程序设计》第四章习题与代码解析
m0_46339652的博客
05-04 3611
目录4.15 EAX符号扩展到EDX4.17 输入输出0～9数字 4.15 EAX符号扩展到EDX IA-32处理器的指令CDQ将EAX符号扩展到EDX。假若没有该指令，编程实现该指令功能。 （1）按照符号扩展的含义编程，即：EAX最高为0，则EDX＝0；EAX最高为1，则EDX＝FFFFFFFFH。 程序代码： include io32.inc .data plmsg byte 'Input Number:' ,13,10,0 .code start: mov eax,offset
RFID技术中的数学形态学在SAR图像增强中的应用
11-06
在SAR图像增强中，数学形态学可以结合频域增强和空域增强方法，通过选择合适的结构元素和迭代次数，实现对噪声的有效滤除，同时保持或增强目标特征。这种非线性滤波方式特别适合处理SAR图像的相干斑噪声问题，因为它...
浅谈INFORMIX SQL优化.pdf
最新发布
01-04
"浅谈INFORMIX SQL优化" INFORMIX SQL优化是一种非常重要的数据库优化技术，在INFORMIX数据库管理系统中，它是数据库管理员、程序员和开发者的共同职责。INFORMIX SQL优化的主要目标是提高SQL语句的执行效率、降低...
Microchip推出全新SAR ADC系列产品，在恶劣环境中依然可实现高速、
01-20
今日宣布推出12款全新逐次逼近寄存器（SAR）模数转换器（ADC）以及专为新型SAR ADC产品系列设计的配套差分放大器，以满足应用市场对更高速度和更高分辨率模数转换的需求。MCP331x1（D）-xx系列产品专为高温和高电磁...
模拟技术中的关于SAR ADC的数字校准算法
10-22
在现有工艺水平下，由于受电容失配、系统失调以及噪声等因素的限制，采用电荷再分配结构的SAR ADC能够达到的最高精度被限制在12位左右。因此，高精度ADC设计必须依靠校准技术。一般校准技术有两类：模拟校准技术是在...
EAX、ECX、EDX、EBX寄存器的作用
01-15 4083
一般寄存器:AX、BX、CX、DX AX:累积暂存器，BX:基底暂存器，CX:计数暂存器，DX:资料暂存器 索引暂存器:SI、DI SI:来源索引暂存器，DI:目的索引暂存器 堆叠、基底暂存器:SP、BP SP:堆叠指标暂存器，BP:基底指标暂存器 EAX、ECX、EDX、EBX：為ax,bx,cx,dx的延伸，各為32位元 ESI、EDI、ESP、EBP：為si,di,
汇编常用指令的使用
zuoheizhu的博客
03-22 751
MOV A,123 就是直接给A赋值 这句话的意思，就是，把123放在A里面
(ZT)算术移位和逻辑移位的区别
谢绝留言与私信
10-04 1170
original url from : http://demon.tw/programming/assembly-sar-shr.html 汇编语言中SAR和SHR指令都是右移指令，SAR是算数右移指令（shift arithmetic right），而SHR是逻辑右移指令（shift logical right）。 两者的区别在于SAR右移时保留操作数的符号，即用符号位来
C语言的整数除法
不负长风
09-01 1万+
c语言的整数除法是向下取整的。
编译器整数除法的优化
qq_40710190的博客
08-16 1344
编译器整数除法的优化 整数的除法有2种情况 被除数是未知数 那么编译器会使用div指令，运行效率会下降 被除数/除数是常量 那么情况会比较复杂，如果除数是2n2^n2n，就会使用shr a,n。其中a是被除数，n是2的幂。如果除数不是2呢？那么情况就会更加复杂了，不同的编译器会有不同的算法。通常的思路是将除法转换成乘法来提高运行的效率，最常见的优化算法便是倒数相乘了 ​ ab=a∗(1b)\frac{a}{b}=a*(\frac{1}{b})ba​=a∗(b1​)
python基础--除法，取整，取模（取余）（/,//,%），以及int最大值，大数乘法
热门推荐
记录和分享程序人生的点点滴滴
10-16 15万+
1. 内容简介： 在算术运算中，加减乘比较容易理解，分别用+，-，*来表示；但是，经常容易混淆的是除法，取整，取模这三个操作，今天主要演示除法，取整，取模的用法，以及一些注意事项。
linux 性能监控命令7——sar 命令
changyanmanman的专栏
12-03 2919
sar（System Activity Reporter系统活动情况报告）是目前Linux 上最为全面的系统性能分析工具之一，可以从多方面对系统的活动进行报告，包括：文件的读写情况、系统调用的使用情况、磁盘I/O、CPU效率、内存使用状况、进程活动及IPC有关的活动等。本文主要以CentOS 6.3 x64系统为例，介绍sar命令。 sar命令常用格式 sar [options] [-
偏振SAR图像在海洋波测量中的应用
另一方面，α参数测量法依赖于SAR图像中H和V极化通道的相对相位差异，这是由于波浪对电磁波传播的影响导致的。这种方法能够揭示波浪在距离方向上的信息，从而补充方向角测量的数据。 2 测量方法 在方向角方向，测量...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值