21ic电子技术开发论坛 › 单片机与嵌入式系统 › 侃单片机论坛 › 最据山寨特色的《中颖快速浮点乘、除法程库》！！！ ...
123下一页
返回列表 发新帖我要提问

最据山寨特色的《中颖快速浮点乘、除法程库》！！！

[复制链接]
 楼主| Cortex-M0 发表于 2011-10-1 19:56 | 显示全部楼层 |阅读模式
ui, Value, 单片机, TE, 二进制
一、简介
浮点数的运算，一贯比定点数慢很多，在C51中，浮点数乘除法，由于缺少相适应指令，运算速度超级慢，比浮点数加减法，慢不是一个等级，极大的影响了51单片机的实时运算速度，本文中，依靠中颖SH79/88/89系列的特色----准16位乘除指令，对所有关联到的无符号、有符号int型，无符号、有符号long型 和 float型变量的乘除运算，均作了运算速度大幅度提升，为方便使用，以外挂库形式加入，实行对C51内库自动替换，无缝嵌入。

其中最具有山寨特色的是float型的除法运算，便用中颖SH79/88/89系列的准16位乘除指令，经过简单的组合运算，代替原Keil C51的加减法模拟除法程序，大幅度地提升了原float型的除法的运算速度，解决了原51浮点除法特慢的通病。


二、中颖专用乘除库清单
/*
extern uint16   imul(uint16 value1,uint16 value2);          // 16位*16位-->16位(无符号二进制乘法)   
extern uint16  uidiv(uint32 value1,uint16 value2);       // 16位/16位-->16位(无符号二进制除法)   
extern uint32   lmul(uint32 value1,uint32 value2);          // 32位*32位-->32位(无符号二进制乘法)   
extern uint32  uldiv(uint32 value1,uint32 value2);       // 32位/32位-->32位(无符号二进制除法)   
extern float   fpmul(float value1, float value2);           // 二进制浮点数快速乘法  
extern float   fpdiv(float value1, float value2);           // 二进制浮点数快速除法  
*/


三、便用说明
只需要在编译程序过程中国，加入《中颖专用乘除库》----MATH_SH88F51.LIB，其他什么事都不用干！在编译过程中实现自动替换原C51库函数，无缝嵌入。


四、友情提示
本《中颖专用乘除库》，仅适用于中颖SH79/88/89系列单片机，其他51单片机并不适用，由于Keil C51不支持中颖51单片机的软件仿真，所以必须要用中颖51单片机的硬件仿真或依靠中颖51单片机的实时运行，运算结果才能正确。
回复 赏

举报

相关帖子
 楼主| Cortex-M0 发表于 2011-10-1 19:57 | 显示全部楼层
《中颖专用乘除库》

本帖子中包含更多资源

您需要 登录 才可以下载或查看，没有账号？注册

×
回复 赏

举报

goodboy3021 发表于 2011-10-1 20:08 | 显示全部楼层
有用，占位支持！下载使用！
回复 赏

举报

 楼主| Cortex-M0 发表于 2011-10-1 20:16 | 显示全部楼层
本《中颖专用乘除库》，增加了若干非常实用的，无符号、有符号定点多字节乘法、除法子程序库，还有一个特色就是，增加了无符号、有符号定点多字节除法同时取商和余数子程序库，对同时需要用到这两个参数的场合，极大地提高了运行速度，使用时只需在C程序中，插入如下头本件 global.h 即可。

/***************************************************************************************
Model  : global.h
Description : Head file of defining global variable.
Author  : CLR
Create Time : 2011-09-29
Version ID  : 1.3
用  途      ：中颖SH79/88/89F51系列MCU
作  者      ：许意义
21ic  ID    ：LAOXU
中颖论坛    : bbs.21ic.com   
****************************************************************************************/
#ifndef  _GLOBAL_H
#define  _GLOBAL_H

#pragma  MODDP2     // 打开中颖单片机双DPTR功能，加速数据传送
//#pragma  NOMODDP2    // 关闭中颖单片机双DPTR功能

/*--------------------------------------------------------------------------------------*/
/*      Resource - math.asm            */
/*--------------------------------------------------------------------------------------*/
extern  uint8 ascii_hex(uint8 value1);     //  1位 ASCII        --> 1位 HEX(低半字节)
extern  uint8 hex_ascii(uint8 value1);     //  1位 HEX(低半字节)--> 1位 ASCII
extern  uint8 bcd8_hex(uint8 value1);         //  8位 BCD --> 8位 HEX
extern uint16  bcd16_hex(uint16 value1);        // 16位 BCD -->16位 HEX
extern uint32  bcd32_hex(uint32 value1);        // 32位 BCD -->32位 HEX      
extern  uint16 hex8_bcd(uint8 value1);         //  8位 HEX -->16位 BCD
extern uint32  hex16_bcd(uint16 value1);        // 16位 HEX -->32位 BCD
extern uint32  hex32_bcd(uint32 value1);        // 32位 HEX -->32位 BCD      
extern  uint32  mul16_8(uint16 value1,uint8 value2);  // 16位* 8位-->32位(无符号二进制乘法)  
extern  uint32  mul16(uint16 value1,uint16 value2);   // 16位*16位-->32位(无符号二进制乘法)
extern  uint32  mul32_8(uint32 value1,uint8 value2);  // 32位* 8位-->32位(无符号二进制乘法)  
extern  uint32  mul32_16(uint32 value1,uint16 value2);  // 32位*16位-->32位(无符号二进制乘法)  
extern  uint32  mul32(uint32 value1,uint32 value2);      // 32位*32位-->32位(无符号二进制乘法)  
extern uint16  div16_8(uint16 value1,uint8 value2);   // 16位/ 8位-->16位(无符号二进制除法)
extern  uint16  div16(uint16 value1,uint16 value2);   // 16位/16位-->16位(无符号二进制除法)
extern  uint32  div32_8(uint32 value1,uint8 value2);  // 32位/ 8位-->32位(无符号二进制除法)
extern  uint32  div32_16(uint32 value1,uint16 value2);  // 32位/16位-->32位(无符号二进制除法)
extern  uint32  div32(uint32 value1,uint32 value2);      // 32位/32位-->32位(无符号二进制除法)
extern uint8  mod16_8(uint16 value1,uint8 value2);  // 16位% 8位--> 8位(无符号二进制取余数)
extern uint16  mod16(uint16 value1,uint16 value2);   // 16位%16位-->16位(无符号二进制取余数)  
extern uint8   mod32_8(uint32 value1,uint8 value2);   // 32位% 8位--> 8位(无符号二进制取余数)
extern uint16 mod32_16(uint32 value1,uint16 value2);   // 32位%16位-->16位(无符号二进制取余数)  
extern uint32 mod32(uint32 value1,uint32 value2);      // 32位%32位-->32位(无符号二进制取余数)   
extern uint16  div16_8m(uint16 value1,uint8 value2,uint8 data *p);    // 16位/ 8位-->16位(无符号二进制除法及取余数，8位余数存*p)
extern uint16  div16m(uint16 value1,uint16 value2,uint16 data *p);    // 16位/16位-->16位(无符号二进制除法及取余数，16位余数存*p)
extern uint32  div32_8m(uint32 value1,uint8 value2,uint8 data *p);    // 32位/ 8位--> 8位(无符号二进制除法及取余数，8位余数存*p)
extern uint32 div32_16m(uint32 value1,uint16 value2,uint16 data *p); // 32位/16位-->16位(无符号二进制除法及取余数，16位余数存*p)
extern uint32 div32m(uint32 value1,uint32 value2,uint32 data *p);    // 32位/32位-->32位(无符号二进制除法及取余数，32位余数存*p)   
extern  int32  smul16_8(int16 value1,int8 value2);      // 16位* 8位-->32位(有符号二进制补码乘法)  
extern  int32  smul16(int16 value1,int16 value2);      // 16位*16位-->32位(有符号二进制补码乘法)
extern  int32  smul32_8(int32 value1,int8 value2);      // 32位* 8位-->32位(有符号二进制补码乘法)  
extern  int32  smul32_16(int32 value1,int16 value2);     // 32位*16位-->32位(有符号二进制补码乘法)  
extern  int32  smul32(int32 value1,int32 value2);      // 32位*32位-->32位(有符号二进制补码乘法)  
extern int16  sdiv16_8(int16 value1,int8 value2);       // 16位/ 8位-->16位(有符号二进制补码除法)
extern  int16  sdiv16(int16 value1,int16 value2);      // 16位/16位-->16位(有符号二进制补码除法)
extern  int32  sdiv32_8(int32 value1, int8 value2);      // 32位* 8位-->32位(有符号二进制补码乘法)   
extern  int32  sdiv32_16(int32 value1,int16 value2);     // 32位/16位-->32位(有符号二进制补码除法)
extern  int32  sdiv32(int32 value1,int32 value2);      // 32位/32位-->32位(有符号二进制补码除法)
extern int8   smod16_8(int16 value1,int8 value2);      // 16位% 8位--> 8位(有符号二进制补码取余数)
extern int16  smod16(int16 value1,int16 value2);      // 16位%16位-->16位(有符号二进制补码取余数)  
extern int8   smod32_8(int32 value1,int8 value2);       // 32位% 8位--> 8位(有符号二进制补码取余数)
extern int16  smod32_16(int32 value1,int16 value2);     // 32位%16位-->16位(有符号二进制补码取余数)  
extern int32  smod32(int32 value1,int32 value2);      // 32位%32位-->32位(有符号二进制补码取余数)   
extern int16  sdiv16_8m(int16 value1,int8 value2,int8 data *p); // 16位/ 8位-->16位(有符号二进制除法及取余数，8位余数存*p)
extern int16  sdiv16m(int16 value1,int16 value2,int16 data *p); // 16位/16位-->16位(有符号二进制除法及取余数，16位余数存*p)
extern int32  sdiv32_8m(int32 value1,int8 value2,int8 data *p);  // 32位/ 8位--> 8位(有符号二进制除法及取余数，8位余数存*p)
extern int32  sdiv32_16m(int32 value1,int16 value2,int16 data *p); // 32位/16位-->16位(有符号二进制除法及取余数，16位余数存*p)
extern int32  sdiv32m(int32 value1,int32 value2,int32 data *p); // 32位/32位-->32位(有符号二进制除法及取余数，32位余数存*p)   

extern  void  ld_pdata(int8 addr1, int8 addr2, int8 len);   //  @PDATA --> @IDATA（len字节PDATA读入片内IDATA RAM）
extern  void  st_pdata(int8 addr1, int8 addr2, int8 len);   //  @IDATA --> @PDATA（len字节片内IDATA RAM存入PDATA）

/*  中颖单片机C51定点库内部调用函数   
extern  int32 st_r03_r47(void);             //  R0R1R2R3 ---> R4R5R6R7
extern  int16 neg(void);                     //  -(R6R7)  ---> R6R7
*/

/*  中颖单片机C51定点库、浮点库自动调用专用运算库函数，合适有效地插入了中颖SH79/88/89系列的特色----准16位乘除指令，对所有
    关联到的无符号、有符号int型，无符号、有符号long型 和 float型变量的乘除运算，均作了运算速度大幅度提升，为方便使用，以外
挂库形式加入，实行对C51内库自动替换，无缝嵌入，其中最具有特色的是float型的除法运算，用中颖SH79/88/89系列的准16位乘除指
令，代替原Keil C51的加减法模拟除法程序，大幅度地提升了原float型的除法的运算速度，解决了原51浮点除法特慢的通病。

extern uint16   imul(uint16 value1,uint16 value2);          // 16位*16位-->16位(无符号二进制乘法)   
extern uint16  uidiv(uint32 value1,uint16 value2);       // 16位/16位-->16位(无符号二进制除法)   
extern uint32   lmul(uint32 value1,uint32 value2);          // 32位*32位-->32位(无符号二进制乘法)   
extern uint32  uldiv(uint32 value1,uint32 value2);       // 32位/32位-->32位(无符号二进制除法)   
extern float   fpmul(float value1, float value2);           // 二进制浮点数快速乘法  
extern float   fpdiv(float value1, float value2);           // 二进制浮点数快速除法  
*/

#endif


本帖子中包含更多资源

您需要 登录 才可以下载或查看，没有账号？注册

×
回复 赏

举报

 楼主| Cortex-M0 发表于 2011-10-1 20:19 | 显示全部楼层
欢迎各路盆友下载使用~~~
回复 赏

举报

highgear 发表于 2011-10-1 22:10 | 显示全部楼层
我顶，这对于算术运算的程序员简直是福音。
不做的人不知道，做过的人都知道。
回复 赏

举报

 楼主| Cortex-M0 发表于 2011-10-1 22:27 | 显示全部楼层
谢谢highgear老师支持~~~

这个程序库，最难的就是51的浮点数快速除法，今天花了大半天时间，来研究51的浮点数快速除法，充分使用了中颖SH79/88/89系列的特色----准16位乘除指令，对原Keil C51的加减法模拟除法库程序，进行改造，
经过简单的组合运算，成功替代原Keil C51的浮点除法库程序时，还是很兴奋的，程序编写其实很简单，关健是算法的从理论到实践的验证~~~

在8051上，要实现32位快速浮点数除法，并不容易，感谢中颖，感谢中颖51，感谢中颖51单片机的准16位乘除指令，让俺实现了用乘除组合运算，成功替代了原C51用加减法模拟除法的库函数，实际运算速度大幅度提升~~~
回复 赏

举报

 楼主| Cortex-M0 发表于 2011-10-1 22:36 | 显示全部楼层
将改写的中颖51快速浮点数乘法、除法哂一晒，供有情人参考~~~


中颖51快速浮点数乘法程序FPMUL.A51

复制
  1. $NOMOD51
    


  2. 

  3. NAME ?C?FPMUL
    


  4. 

  5. /**********************************************************************************************
    

  6. * float R4R5R6R7 fpmul(float R4R5R6R7, float R0R1R2R3);
    

  7. * result = value1 * value2
    

  8. * R4R5R6R7 = R4R5R6R7 * R0R1R2R3
    

  9. * 32位*32位-->32位(4字节浮点数快速乘法) 
    


  10. 

  11. 用 途 ：中颖SH79/88/89F51系列MCU
    

  12. 作 者 ：许意义
    

  13. 21ic ID ：LAOXU
    

  14. 中颖论坛 : bbs.21ic.com 
    

  15. 日 期 : 2011-09-29 
    


  16. 

  17. **********************************************************************************************/
    


  18. 

  19. ?PR?_FPMUL?FPMUL SEGMENT CODE 
    

  20. EXTRN CODE (_FPRET_1)
    

  21. EXTRN CODE (_FPRET_2)
    

  22. EXTRN CODE (_FPRET_3)
    

  23. EXTRN CODE (_FPRET_4)
    

  24. EXTRN CODE (_FPRET_5)
    

  25. EXTRN CODE (_FPRET_6)
    

  26. PUBLIC ?C?FPMUL
    

  27. RSEG ?PR?_FPMUL?FPMUL
    

  28. ?C?FPMUL: MOV A,R4
    

  29. ORL A,R5
    

  30. JZ FMUL01
    

  31. MOV A,R0
    

  32. ORL A,R1
    

  33. JNZ FMUL02
    

  34. MOV A,R5
    

  35. RLC A
    

  36. MOV A,R4
    

  37. RLC A
    

  38. INC A
    

  39. JZ FMUL03
    

  40. CLR A
    

  41. MOV R4,A
    

  42. MOV R7,A
    

  43. MOV R6,A
    

  44. MOV R5,A
    

  45. FMUL04: RET 
    

  46. FMUL01: MOV A,R1
    

  47. RLC A
    

  48. MOV A,R0
    

  49. RLC A
    

  50. INC A
    

  51. JNZ FMUL04
    

  52. FMUL03: LJMP _FPRET_4
    

  53. FMUL02: LCALL _FPRET_1
    

  54. ANL A,R0
    

  55. INC A
    

  56. JZ FMUL07
    

  57. CLR A
    

  58. XCH A,R4
    

  59. ADD A,#81H
    

  60. JNC FMUL08
    

  61. ADD A,R0
    

  62. JNC FMUL09
    

  63. FMUL07: LJMP _FPRET_6
    

  64. FMUL08: ADD A,R0
    

  65. JC FMUL09
    

  66. LJMP _FPRET_5
    

  67. FMUL09: PUSH ACC
    

  68. ORL 0x86,#0x0c
    

  69. MOV A,R3
    

  70. ORL A,R2
    

  71. JNZ FMUL12
    

  72. CJNE R1,#0x80,FMUL19
    

  73. FMUL15: POP ACC
    

  74. MOV R3,A
    

  75. LJMP _FPRET_3
    

  76. FMUL19: MOV A,R7
    

  77. ORL A,R6
    

  78. JNZ FMUL14
    

  79. CJNE R5,#0x80,FMUL20
    

  80. FMUL21: MOV A,R3
    

  81. MOV R7,A
    

  82. MOV A,R2
    

  83. MOV R6,A
    

  84. MOV A,R1
    

  85. MOV R5,A
    

  86. SJMP FMUL15
    

  87. FMUL20: MOV 0xF1,#0
    

  88. MOV A,R1
    

  89. MOV B,R5
    

  90. MUL AB
    

  91. MOV R6,A
    

  92. MOV A,B
    

  93. LJMP FMUL16
    

  94. FMUL18: MOV A,R1
    

  95. XCH A,R5
    

  96. MOV R1,A
    

  97. MOV A,R2
    

  98. MOV R6,A
    

  99. MOV A,R3
    

  100. MOV R7,A
    

  101. FMUL14: MOV 0xF1,R6
    

  102. MOV A,R7
    

  103. MOV B,R1
    

  104. MUL AB
    

  105. MOV R4,A
    

  106. MOV R7,B
    

  107. MOV A,0xF1
    

  108. XCH A,R5
    

  109. MOV 0xF1,#0
    

  110. MOV B,R1
    

  111. MUL AB
    

  112. ADD A,R5
    

  113. MOV R6,A
    

  114. CLR A
    

  115. ADDC A,B
    

  116. SJMP FMUL16
    

  117. FMUL12: MOV A,R7
    

  118. ORL A,R6
    

  119. JNZ FMUL17
    

  120. CJNE R5,#0x80,FMUL18
    

  121. SJMP FMUL21
    

  122. FMUL17: MOV 0xF1,R6
    

  123. MOV A,R7
    

  124. MOV B,R3
    

  125. MUL AB
    

  126. MOV R4,B
    

  127. MOV R0,0xF1
    

  128. MOV 0xF1,R6
    

  129. MOV A,R7
    

  130. MOV B,R2
    

  131. MUL AB
    

  132. ADD A,R4
    

  133. MOV A,B
    

  134. ADDC A,R0
    

  135. MOV R4,A
    

  136. CLR A
    

  137. ADDC A,0xF1
    

  138. MOV R0,A
    

  139. MOV 0xF1,R6
    

  140. MOV A,R7
    

  141. MOV B,R1
    

  142. MUL AB
    

  143. ADD A,R4
    

  144. MOV R4,A 
    

  145. MOV A,B
    

  146. ADDC A,R0
    

  147. MOV R7,A 
    

  148. CLR A
    

  149. ADDC A,0xF1
    

  150. MOV R0,A 
    

  151. MOV 0xF1,R2
    

  152. MOV A,R3
    

  153. MOV B,R5
    

  154. MUL AB
    

  155. ADD A,R4 
    

  156. MOV R4,A
    

  157. MOV A,B
    

  158. ADDC A,R7 
    

  159. MOV R7,A
    

  160. MOV A,0xF1
    

  161. ADDC A,R0
    

  162. MOV R6,A
    

  163. CLR A
    

  164. MOV 0xF1,A
    

  165. RLC A
    

  166. XCH A,R5
    

  167. MOV B,R1
    

  168. MUL AB
    

  169. ADD A,R6
    

  170. MOV R6,A
    

  171. MOV A,B
    

  172. ADDC A,R5
    

  173. FMUL16: MOV R5,A
    

  174. ANL 0x86,#0xf3
    

  175. RLC A
    

  176. POP ACC
    

  177. MOV R3,A
    

  178. JNC FMUL22
    

  179. INC R3
    

  180. CJNE R3,#0x00,FMUL23
    

  181. LJMP _FPRET_6
    

  182. FMUL22: MOV A,R4
    

  183. ADD A,R4
    

  184. MOV R4,A
    

  185. MOV A,R7
    

  186. RLC A
    

  187. MOV R7,A
    

  188. MOV A,R6
    

  189. RLC A
    

  190. MOV R6,A
    

  191. MOV A,R5
    

  192. RLC A
    

  193. MOV R5,A
    

  194. FMUL23: LJMP _FPRET_2
    


  195. 

  196. END
回复 赏

举报

 楼主| Cortex-M0 发表于 2011-10-1 22:39 | 显示全部楼层
中颖51快速浮点数除法程序FPDIV.A51

复制
  1. $NOMOD51
    


  2. 

  3. NAME ?C?FPDIV
    


  4. 

  5. /**********************************************************************************************
    

  6. * float R4R5R6R7 fpmul(float R4R5R6R7, float R0R1R2R3);
    

  7. * result = value1 / value2
    

  8. * R4R5R6R7 = R4R5R6R7 / R0R1R2R3
    

  9. * 32位/32位-->32位(4字节浮点数快速除法) 
    


  10. 

  11. 用 途 ：中颖SH79/88/89F51系列MCU
    

  12. 作 者 ：许意义
    

  13. 21ic ID ：LAOXU
    

  14. 中颖论坛 : bbs.21ic.com 
    

  15. 日 期 : 2011-09-29 
    


  16. 

  17. **********************************************************************************************/
    


  18. 

  19. ?PR?_FPDIV?FPDIV SEGMENT CODE 
    

  20. EXTRN CODE (_FPRET_1)
    

  21. EXTRN CODE (_FPRET_2)
    

  22. EXTRN CODE (_FPRET_3)
    

  23. EXTRN CODE (_FPRET_4)
    

  24. EXTRN CODE (_FPRET_5)
    

  25. EXTRN CODE (_FPRET_6)
    

  26. PUBLIC ?C?FPDIV
    

  27. DPL1 EQU 0x84
    

  28. DPH1 EQU 0x85
    

  29. RSEG ?PR?_FPDIV?FPDIV
    

  30. FPDIV_4: LJMP _FPRET_6
    

  31. ?C?FPDIV: MOV A,R4
    

  32. ANL A,R5
    

  33. INC A
    

  34. JZ FPDIV_1
    

  35. MOV A,R0
    

  36. ANL A,R1
    

  37. INC A
    

  38. JNZ FPDIV_2
    

  39. FPDIV_1: LJMP _FPRET_4
    

  40. FPDIV_2: LCALL _FPRET_1
    

  41. ANL A,R0
    

  42. INC A
    

  43. JZ FPDIV_1
    

  44. MOV A,R4
    

  45. ORL A,R0
    

  46. JZ FPDIV_1
    

  47. MOV A,R4
    

  48. JNZ FPDIV_3
    

  49. MOV R5,A
    

  50. MOV R6,A
    

  51. MOV R7,A
    

  52. RET 
    

  53. FPDIV_3: XCH A,R0
    

  54. JZ FPDIV_4
    

  55. ADD A,#0x81
    

  56. XCH A,R0
    

  57. JNC FPDIV_5
    

  58. CLR C
    

  59. SUBB A,R0
    

  60. JZ FPDIV_6
    

  61. JNC FPDIV_7
    

  62. FPDIV_6: LJMP _FPRET_5
    

  63. FPDIV_5: SUBB A,R0
    

  64. JNC FPDIV_4
    

  65. FPDIV_7: MOV DPL,A
    

  66. MOV A,R1
    

  67. ADD A,R1
    

  68. ORL A,R3
    

  69. ORL A,R2
    

  70. JNZ FPDIV_8
    

  71. MOV R3,DPL
    

  72. LJMP _FPRET_3
    

  73. FPDIV_8: CLR C
    

  74. MOV A,R7
    

  75. SUBB A,R3
    

  76. MOV A,R6
    

  77. SUBB A,R2
    

  78. MOV A,R5
    

  79. SUBB A,R1
    

  80. JC FPDIV_9
    

  81. CLR C
    

  82. MOV A,R5
    

  83. RRC A
    

  84. MOV R5,A
    

  85. MOV A,R6
    

  86. RRC A
    

  87. MOV R6,A
    

  88. MOV A,R7
    

  89. RRC A
    

  90. MOV R7,A
    

  91. INC DPL
    

  92. // H0(R0) = ffff / R1 
    

  93. FPDIV_9: ORL 0x86,#0x0c
    

  94. MOV A,#0xff 
    

  95. MOV 0xF1,A
    

  96. MOV B,R1
    

  97. DIV AB
    

  98. MOV R0,A 
    

  99. PUSH DPL
    

  100. // N0(R5R6R7R4) = R5R6R7<<8 + (R5R6R7*H0) (放大1字节） 
    

  101. MOV A,R7
    

  102. MOV 0xF1,R6
    

  103. MOV B,R0 
    

  104. MUL AB
    

  105. MOV R4,A
    

  106. MOV DPL1,A
    

  107. MOV A,R7
    

  108. ADD A,B
    

  109. MOV R7,A
    

  110. MOV DPH1,A
    

  111. MOV A,R6
    

  112. ADDC A,0xF1
    

  113. MOV R6,A
    

  114. MOV A,R5
    

  115. MOV 0xF1,#0
    

  116. MOV B,R0 
    

  117. MUL AB
    

  118. ADD A,R6
    

  119. MOV R6,A
    

  120. MOV DPL,A
    

  121. MOV A,B
    

  122. ADDC A,R5
    

  123. MOV R5,A
    

  124. MOV DPH,A // R5R6R7R4 ---> 双DPTR 
    

  125. // D0(R2R3) = R1R2R3 + (R1R2R3*H0>>8) 
    

  126. MOV A,R3
    

  127. MOV 0xF1,R2
    

  128. MOV B,R0 
    

  129. MUL AB
    

  130. MOV A,B
    

  131. ADD A,R3
    

  132. MOV R3,A
    

  133. MOV A,0xF1
    

  134. ADDC A,R2
    

  135. MOV R2,A
    

  136. MOV A,R1
    

  137. MOV B,R0 
    

  138. MUL AB
    

  139. ADD A,R2
    

  140. MOV R2,A
    

  141. // &(R2R3) = -D0 (2字节) 
    

  142. CLR C
    

  143. CLR A
    

  144. SUBB A,R3
    

  145. MOV R3,A
    

  146. CLR A
    

  147. SUBB A,R2
    

  148. MOV R2,A
    

  149. // Q = N0 + A + B + C
    

  150. // A = N0*& （丢3字节） R5R6R7R4(N0) * R2R3 ---> R5R6R7R4(A) | R5=0
    

  151. LCALL NFPMUL32
    

  152. LCALL NFPADD32
    

  153. // B = A*& （丢3字节） R5R6R7R4(A) * R2R3 ---> R5R6R7R4(B) | R5R6=0 
    

  154. LCALL NFPMUL32
    

  155. LCALL NFPADD32
    

  156. // C = B*& （丢3字节） R5R6R7R4(B) * R2R3 ---> R5R6R7R4(C) | R5R6R7=0 
    

  157. LCALL NFPMUL32
    

  158. LCALL NFPADD32
    

  159. ANL 0x86,#0xf3
    

  160. POP ACC
    

  161. DEC A
    

  162. MOV R3,A
    

  163. MOV R5,DPH
    

  164. MOV R6,DPL
    

  165. MOV R7,DPH1
    

  166. MOV R4,DPL1
    

  167. LJMP _FPRET_2
    

  168. // R5R6R7R4 * R2R3 ---> R5R6R7R4
    

  169. NFPMUL32: MOV A,R3 
    

  170. MOV 0xF1,R2
    

  171. MOV B,R4
    

  172. MUL AB
    

  173. MOV R1,B
    

  174. MOV R4,0xF1
    

  175. MOV A,R3 
    

  176. MOV 0xF1,R2
    

  177. MOV B,R7 
    

  178. MUL AB
    

  179. ADD A,R1
    

  180. MOV A,B 
    

  181. ADDC A,R4
    

  182. MOV R1,A 
    

  183. CLR A
    

  184. ADDC A,0xF1
    

  185. MOV R4,A
    

  186. MOV A,R5
    

  187. ORL A,R6
    

  188. JZ NFPML32_1
    

  189. MOV A,R3 
    

  190. MOV 0xF1,R2
    

  191. MOV B,R6
    

  192. MUL AB
    

  193. ADD A,R1
    

  194. MOV A,B 
    

  195. ADDC A,R4
    

  196. MOV R4,A 
    

  197. CLR A
    

  198. ADDC A,0xF1
    

  199. MOV R7,A
    

  200. MOV A,R5
    

  201. JZ NFPML32_2
    

  202. MOV B,A
    

  203. MOV A,R3
    

  204. MOV 0xF1,R2
    

  205. MUL AB
    

  206. ADD A,R4
    

  207. MOV R4,A 
    

  208. MOV A,B
    

  209. ADDC A,R7
    

  210. MOV R7,A
    

  211. CLR A
    

  212. MOV R5,A
    

  213. ADDC A,0xF1
    

  214. MOV R6,A
    

  215. RET
    

  216. NFPML32_1: MOV R7,#0
    

  217. NFPML32_2: MOV R6,#0
    

  218. MOV R5,#0
    

  219. RET
    

  220. // 双DPTR + R4R5R6R7 ---> 双DPTR 
    

  221. NFPADD32: MOV A,DPL1
    

  222. ADD A,R4
    

  223. MOV DPL1,A
    

  224. MOV A,DPH1
    

  225. ADDC A,R7
    

  226. MOV DPH1,A
    

  227. MOV A,DPL
    

  228. ADD A,R6
    

  229. MOV DPL,A
    

  230. MOV A,DPH
    

  231. ADDC A,R5
    

  232. MOV DPH,A
    

  233. RET
    


  234. 

  235. END
回复 赏

举报

 楼主| Cortex-M0 发表于 2011-10-1 22:42 | 显示全部楼层
中颖51快速浮点数乘除法通用返回处理程序FPRET.A51


复制
  1. $NOMOD51
    


  2. 

  3. NAME _FPRET
    


  4. 

  5. /**********************************************************************************************
    

  6. * float R4R5R6R7 fpret(float R4R5R6R7, float R0R1R2R3);
    

  7. * result = value1 *^/ value2
    

  8. * R4R5R6R7 = R4R5R6R7 *^/ R0R1R2R3
    

  9. * 32位*32位-->32位(4字节浮点数乘除法返回处理程序) 
    


  10. 

  11. 用 途 ：中颖SH79/88/89F51系列MCU
    

  12. 作 者 ：许意义
    

  13. 21ic ID ：LAOXU
    

  14. 中颖论坛 : bbs.21ic.com 
    

  15. 日 期 : 2011-09-29 
    


  16. 

  17. **********************************************************************************************/
    


  18. 

  19. ?PR?_FPRET?FPRET SEGMENT CODE 
    

  20. PUBLIC _FPRET_1
    

  21. PUBLIC _FPRET_2
    

  22. PUBLIC _FPRET_3
    

  23. PUBLIC _FPRET_4
    

  24. PUBLIC _FPRET_5
    

  25. PUBLIC _FPRET_6
    


  26. 

  27. RSEG ?PR?_FPRET?FPRET
    

  28. _FPRET_1: MOV A,R1
    

  29. SETB ACC.7
    

  30. XCH A,R1
    

  31. RLC A
    

  32. MOV A,R0
    

  33. RLC A
    

  34. MOV R0,A
    

  35. MOV PSW.5,C
    

  36. MOV A,R5
    

  37. SETB ACC.7
    

  38. XCH A,R5
    

  39. RLC A
    

  40. MOV A,R4
    

  41. RLC A
    

  42. MOV R4,A
    

  43. JNC FPT_01
    

  44. CPL PSW.5
    

  45. FPT_01: RET 
    

  46. _FPRET_2: MOV A,R4
    

  47. JNB ACC.7,_FPRET_3
    

  48. INC R7
    

  49. CJNE R7,#0x00,_FPRET_3
    

  50. INC R6
    

  51. CJNE R6,#0x00,_FPRET_3
    

  52. INC R5
    

  53. CJNE R5,#0x00,_FPRET_3
    

  54. INC R3
    

  55. MOV A,R3
    

  56. JZ _FPRET_6
    

  57. _FPRET_3: MOV C,PSW.5
    

  58. MOV A,R3
    

  59. RRC A
    

  60. MOV R4,A
    

  61. MOV A,R5
    

  62. MOV ACC.7,C
    

  63. MOV R5,A
    

  64. RET 
    

  65. _FPRET_4: MOV A,#0xFF
    

  66. FPT_02: MOV R4,A
    

  67. MOV R5,A
    

  68. FPT_03: MOV R6,A
    

  69. MOV R7,A
    

  70. RET 
    

  71. _FPRET_5: CLR A
    

  72. SJMP FPT_02
    

  73. _FPRET_6: MOV C,PSW.5
    

  74. MOV A,#0xFF
    

  75. RRC A
    

  76. MOV R4,A
    

  77. MOV R5,#0x80
    

  78. CLR A
    

  79. SJMP FPT_03
    


  80. 

  81. END
回复 赏

举报

zoomone 发表于 2011-10-2 12:56 | 显示全部楼层
有时间了就钻研钻研
回复 赏

举报

highgear 发表于 2011-10-5 08:23 | 显示全部楼层
顶出水面
回复 赏

举报

zxcscm 发表于 2011-10-18 08:59 | 显示全部楼层
来顶，顶~顶！
回复 赏

举报

xubx83 发表于 2011-10-18 16:38 | 显示全部楼层
顶个先
回复 赏

举报

zoomone 发表于 2011-10-20 12:57 | 显示全部楼层
能不能介绍一下怎么设置和使用的，我添加后编译老是出错
回复 赏

举报

 楼主| Cortex-M0 发表于 2011-10-20 16:42 | 显示全部楼层
关键这几点设置。







打包给个带工程的DEMO测试文件

本帖子中包含更多资源

您需要 登录 才可以下载或查看，没有账号？注册

×
回复 赏

举报

zoomone 发表于 2011-10-20 17:57 | 显示全部楼层
extern uint16   imul(uint16 value1,uint16 value2);          // 16位*16位-->16位(无符号二进制乘法)   
extern uint16  uidiv(uint32 value1,uint16 value2);       // 16位/16位-->16位(无符号二进制除法)   
extern uint32   lmul(uint32 value1,uint32 value2);          // 32位*32位-->32位(无符号二进制乘法)   
extern uint32  uldiv(uint32 value1,uint32 value2);       // 32位/32位-->32位(无符号二进制除法)   
extern float   fpmul(float value1, float value2);           // 二进制浮点数快速乘法  
extern float   fpdiv(float value1, float value2);           // 二进制浮点数快速除法  

为什么这几个函数不能用？
回复 赏

举报

 楼主| Cortex-M0 发表于 2011-10-20 18:58 | 显示全部楼层
不是这几个函数不能用，这几个函数和 Keil C51内部库函数重名，也就是说，只要在编译时加入 math_sh88f51.lib ,  将自动替代原 Keil C51内部库函数，实现无缝嵌入，在实际使用中，你根本不需要管这些库函数的存在。
例如:  int  a,b,c;
         float  x,y,z;
         c = a * b / c;
         z = x * y / z;
Keil C51将自动调用上述SH88F51专用库函数，替代原 Keil C51内部库函数，实现无缝嵌入。
回复 赏

举报

zoomone 发表于 2011-10-20 19:07 | 显示全部楼层
噢 看来不会的还很多:L
回复 赏

举报

 楼主| Cortex-M0 发表于 2011-10-20 19:52 | 显示全部楼层
友情提示:
假如你调用上述库函数，脱离SH88F51的真实运行环境，如用Keil C51软件仿真，或将程序装载到其他厂家的51单片机中运行，将出错，无法得到预期的结果~~~

本库只有在特定的环境中，即在SH88F51的硬件支持下，运行结果才会正确。
回复 赏

举报

下一页 »
123下一页
返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Cortex-M0
个人签名：学习chunyang和maychang两位老师，努力练好基本功。 学习HOT大叔，学好“Cortex-M0”，做“Cortex-M0”之菜鸟。

7

主题

4820

帖子

12

粉丝
快速回复 在线客服 返回列表 返回顶部