[发明专利]应用于嵌入式监测设备的IAPWS-IF97物性幂次计算方法在审
| 申请号： | 201410837601.8 | 申请日： | 2014-12-29 |
| 公开（公告）号： | CN104461450A | 公开（公告）日： | 2015-03-25 |
| 发明（设计）人： | 刘争光;程懋华;孙振业 | 申请（专利权）人： | 东南大学 |
| 主分类号： | G06F7/552 | 分类号： | G06F7/552 |
| 代理公司： | 南京苏高专利商标事务所(普通合伙) 32204 | 代理人： | 李昊 |
| 地址： | 210096 *** | 国省代码： | 江苏;32 |
| 权利要求书： | 查看更多 | 说明书： | 查看更多 |
| 摘要： | 本发明公开一种应用于嵌入式监测设备的IAPWS-IF97物性幂次计算方法，其特征在于，对于IAPWS-IF97公式中的任意常数x的xp幂次计算，调用如下方法实现：1)定义数组Xi(i＝1～p正整数)，并赋值X1＝x；2)当P＝1时，直接返回X1的计算值；当P＝2时，计算X2＝X1*X1，然后返回X2的计算值；当P＝3时，依次计算X2＝X1*X1、X3＝X2*X1，然后返回X3的计算值；当P＝4时，依次计算X2＝X1*X1、X4＝X2*X2，然后返回X4的计算值。 | ||
| 搜索关键词： | 应用于 嵌入式 监测 设备 iapws if97 物性 计算方法 | ||
【主权项】：
一种应用于嵌入式监测设备的IAPWS‑IF97物性幂次计算方法，其特征在于，对于IAPWS‑IF97公式中的任意常数x的xp幂次计算，调用如下方法实现：1)定义数组Xi(i＝1～p正整数)，并赋值X1＝x；2)当P＝1时，直接返回X1的计算值；当P＝2时，计算X2＝X1*X1，然后返回X2的计算值；当P＝3时，依次计算X2＝X1*X1、X3＝X2*X1，然后返回X3的计算值；当P＝4时，依次计算X2＝X1*X1、X4＝X2*X2，然后返回X4的计算值；当P＝5时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2，然后返回X5的计算值；当P＝6时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3，然后返回X6的计算值；当P＝7时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X7＝X5*X2，然后返回X7的计算值；当P＝8时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4，然后返回X8的计算值；当P＝9时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1，然后返回X9的计算值；当P＝10时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X10＝X5*X5，然后返回X10的计算值；当P＝11时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X10＝X5*X5、X11＝X10*X1，然后返回X11的计算值；当P＝12时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6，然后返回X12的计算值；当P＝13时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1、X13＝X9*X4，然后返回X13的计算值；当P＝14时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X7＝X5*X2、X14＝X7*X7，然后返回X14的计算值；当P＝15时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X15＝X12*X3，然后返回X15的计算值；当P＝16时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X16＝X8*X8，然后返回X16的计算值；当P＝17时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1、X17＝X9*X8，然后返回X17的计算值；当P＝18时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X16＝X8*X8、X18＝X16*X2，然后返回X18的计算值；当P＝19时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X16＝X8*X8、X18＝X16*X2、X19＝X18*X1，然后返回X19的计算值；当P＝20时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X10＝X5*X5、X20＝X10*X10，然后返回X20的计算值；当P＝21时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X15＝X12*X3、X21＝X15*X6，然后返回X21的计算值；当P＝22时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X10＝X5*X5、X11＝X10*X1、X22＝X11*X11，然后返回X22的计算值；当P＝23时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X10＝X5*X5、X20＝X10*X10、X23＝X20*X3，然后返回X23的计算值；当P＝24时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X24＝X12*X12，然后返回X24的计算值；当P＝25时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1、X17＝X9*X8、X25＝X17*X8，然后返回X25的计算值；当P＝26时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1、X13＝X9*X4、X26＝X13*X13，然后返回X26的计算值；当P＝27时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X24＝X12*X12、X27＝X24*X3，然后返回X27的计算值；当P＝28时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X7＝X5*X2、X14＝X7*X7、X28＝X14*X14，然后返回X28的计算值；当P＝29时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1、X17＝X9*X8、X25＝X17*X8、X29＝X25*X4，然后返回X29的计算值；当P＝30时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X15＝X12*X3、X30＝X15*X15，然后返回X30的计算值；当P＝31时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X7＝X5*X2、X14＝X7*X7、X28＝X14*X14、X31＝X28*X3，然后返回X31的计算值；当P＝32时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X16＝X8*X8、X32＝X16*X16，然后返回X32的计算值；当P＝33时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X16＝X8*X8、X32＝X16*X16、X33＝X32*X1，然后返回X33的计算值；当P＝34时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1、X17＝X9*X8、X34＝X17*X17，然后返回X34的计算值；当P＝35时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1、X13＝X9*X4、X26＝X13*X13、X35＝X26*X9，然后返回X35的计算值；当P＝36时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X16＝X8*X8、X18＝X16*X2、X36＝X18*X18，然后返回X36的计算值；当P＝37时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X16＝X8*X8、X18＝X16*X2、X36＝X18*X18、X37＝X36*X1，然后返回X37的计算值；当P＝38时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X16＝X8*X8、X18＝X16*X2、X19＝X18*X1、X38＝X19*X19，然后返回X38的计算值；当P＝39时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X24＝X12*X12、X27＝X24*X3、X39＝X27*X12，然后返回X39的计算值；当P＝40时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X10＝X5*X5、X20＝X10*X10、X40＝X20*X20，然后返回X40的计算值；当P＝41时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X10＝X5*X5、X20＝X10*X10、X40＝X40*X40、X41＝X40*X1，然后返回X41的计算值；当P＝42时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X15＝X12*X3、X21＝X15*X6、X42＝X21*X21，然后返回X42的计算值；当P＝43时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1、X17＝X9*X8、X34＝X17*X17、X43＝X34*X9，然后返回X43的计算值；当P＝44时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X10＝X5*X5、X11＝X10*X1、X22＝X11*X11、X44＝X22*X22，然后返回X44的计算值；当P＝45时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X15＝X12*X3、X30＝X15*X15、X45＝X30*X15，然后返回X45的计算值；当P＝46时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X10＝X5*X5、X20＝X10*X10、X23＝X20*X3、X46＝X23*X23，然后返回X46的计算值；当P＝47时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X10＝X5*X5、X20＝X10*X10、X23＝X20*X3、X46＝X23*X23，X47＝X46*X1然后返回X47的计算值；当P＝48时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X24＝X12*X12、X48＝X24*X24，然后返回X48的计算值；当P＝49时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X16＝X8*X8、X32＝X16*X16、X33＝X32*X1、X49＝X33*X16，然后返回X49的计算值；当P＝50时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1、X17＝X9*X8、X25＝X17*X8、X50＝X25*X25，然后返回X50的计算值；当P＝51时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X24＝X12*X12、X48＝X24*X24、X51＝X48*X3，然后返回X51的计算值；当P＝52时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1、X13＝X9*X4、X26＝X13*X13、X52＝X26*X26，然后返回X52的计算值；当P＝53时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X16＝X8*X8、X18＝X16*X2、X36＝X18*X18、X37＝X36*X1，X53＝X37*X16然后返回X53的计算值；当P＝54时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X24＝X12*X12、X27＝X24*X3、X54＝X27*X27，然后返回X54的计算值；当P＝55时，依次计算X2＝X1*X1、X3＝X2*X1、X6＝X3*X3、X12＝X6*X6、X24＝X12*X12、X27＝X24*X3、X54＝X27*X27，X55＝X54*X1然后返回X55的计算值；当P＝56时，依次计算X2＝X1*X1、X3＝X2*X1、X5＝X3*X2、X7＝X5*X2、X14＝X7*X7、X28＝X14*X14、X56＝X28*X28然后返回X56的计算值；当P＝57时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X16＝X8*X8、X32＝X16*X16、X33＝X32*X1、X49＝X33*X16，X57＝X49*X8然后返回X57的计算值；当P＝58时，依次计算X2＝X1*X1、X4＝X2*X2、X8＝X4*X4、X9＝X8*X1、X17＝X9*X8、X25＝X17*X8、X29＝X25*X4，X58＝X29*X29然后返回X58的计算值。
下载完整专利技术内容需要扣除积分，VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息，商用须获得专利权人授权。该专利全部权利属于东南大学，未经东南大学许可，擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作，请联系【客服】
本文链接：http://www.vipzhuanli.com/patent/201410837601.8/，转载请声明来源钻瓜专利网。
- 上一篇：用于微处理器的动态重设方法
- 下一篇：显示控制的方法、装置及移动终端
- 同类专利
- 用于执行浮点平方根运算的装置和方法-201510423618.3
- 内尔·伯吉斯;大卫·雷蒙德·鲁茨 - ARM有限公司
- 2015-07-17 - 2019-07-12 - G06F7/552
- 本发明涉及用于执行浮点平方根运算的装置和方法。该数据处理装置具有用于对被开方数值R执行浮点平方根运算以生成结果值的处理电路。处理电路具有用于处理并非是二的精确乘方的被开方数值R的第一平方根处理电路，以及用于处理作为2的精确乘方的被开方数值的第二平方根处理电路。二的乘方检测电路检测被开方数值是否是二的精确乘方并适当地选择第一平方根处理电路的输出或第二平方根处理电路的输出。这使得当被开方数是2的乘方时能够以较少的处理循环生成结果。
- 一种便于硬件实现的开立方算法-201811009126.X
- 旷旭 - 成都锐成芯微科技股份有限公司
- 2018-08-31 - 2019-01-11 - G06F7/552
- 本发明公开了一种便于硬件实现的开立方算法，涉及集成电路技术领域。该算法包括以下步骤：S1、输入被开立方数至寄存器；S2、根据被开立方数的二进制数值以得到被开立方数的有效位最高六位，并与预设常数进行比较，根据比较结果以确定立方根的有效位最高两位，并根据比较结果设置相应的临时立方根；S3、将临时立方根进行立方运算以得到临时立方数；S4、将被开立方数与临时立方数进行比较，以确定立方根的下一有效位，并根据被开立方数与临时立方数的比较结果以设置新的临时立方根；S5、判断是否已确定立方根的所有有效位，若已确定，则输出被开立方数的立方根；若未确定，则返回步骤S3，直至确定立方根。本发明减少了计算中迭代步骤、节约了资源。
- 便于硬件实现的开平方算法-201810056951.9
- 旷旭 - 成都锐成芯微科技股份有限公司
- 2018-01-22 - 2018-07-24 - G06F7/552
- 本发明公开了一种便于硬件实现的开平方算法，涉及集成电路技术领域。该算法包括以下步骤：S1、输入被开方数至寄存器；S2、根据被开方数的二进制数值进行初始化以得到初始被减数，并预设初始临时减数；S3、将初始被减数和初始临时减数相减；S4、根据计算得到的差及被开方数生成新的被减数和新的临时减数，并将新的被减数与新的临时减数相减以得到新的差；S5、判断是否已求出平方根，若已求出平方根，则输出所有的差以得到被开方数的平方根；若未求出平方根，则返回步骤S4，直至求出被开方数的平方根。本发明技术方案使得开方运算计算速度增加，同时使用的硬件面积更小、便于硬件实现。
- 一种应用于GPU中的幂指数运算装置和运算方法-201210019688.9
- 魏继增;李炳超;常轶松;郭炜;孙济洲 - 天津大学
- 2012-01-21 - 2017-09-29 - G06F7/552
- 本发明涉及集成电路领域，特别是公开了一种应用于GPU中的幂指数运算硬件装置和运算方法第一步，计算log2cosθ第二步，计算xlog2cosθ；第三步，计算2mul_result；第四步，得到2mul_result之后，就可以用IEEE754单精度浮点数表示cosxθ的近似值。与现有技术相比，本发明能够在满足嵌入式系统因为空间大小和功耗而受限制的计算能力的同时，快速而又比较准确的图形处理单元完成Phong模型中的高光计算，为集成电路的设计提供了具有更好的图形处理能力的嵌入式系统设计。
- 一种编译器前端乘幂运算的方法-201410541113.2
- 不公告发明人 - 镇江鼎拓科技信息有限公司
- 2014-10-14 - 2016-05-11 - G06F7/552
- 一种编译器前端乘幂运算的方法就是利用乘幂运算进行，它是最基本的数学表达式之一，并且被大多数编译器支持，在高性能计和密码算法中应用十分频繁。
- 一种平方根计算电路-201420438128.1
- 李文忠 - 李文忠
- 2014-08-01 - 2015-04-22 - G06F7/552
- 本专利公开了一种平方根计算电路，其由减法器、移位寄存器、计数器、比较器、锁存器及其它一些电路组成。
- 应用于嵌入式监测设备的IAPWS-IF97物性幂次计算方法-201410837601.8
- 刘争光;程懋华;孙振业 - 东南大学
- 2014-12-29 - 2015-03-25 - G06F7/552
- 本发明公开一种应用于嵌入式监测设备的IAPWS-IF97物性幂次计算方法，其特征在于，对于IAPWS-IF97公式中的任意常数x的xp幂次计算，调用如下方法实现：1)定义数组Xi(i＝1～p正整数)，并赋值X1＝x；2)当P＝1时，直接返回X1的计算值；当P＝2时，计算X2＝X1*X1，然后返回X2的计算值；当P＝3时，依次计算X2＝X1*X1、X3＝X2*X1，然后返回X3的计算值；当P＝4时，依次计算X2＝X1*X1、X4＝X2*X2，然后返回X4的计算值。
- 用于在密码系统中抗故障求幂的方法以及设备-201410087819.6
- M.乔耶 - 汤姆逊许可公司
- 2014-03-11 - 2014-09-17 - G06F7/552
- 在设备（100）中的处理器（120）实施抗故障求幂，其使用输入x以及秘密指数d通过如下方法得到结果S：使用先验选择的整数r以及选择的随机元素α∈{0，...，r-1}来形成扩展的基数
使得
一般地，关于先验选择的、与模数N互质的整数t=br2（其中b是整数），处理器（110）进行模数的求逆iN=N-1modt。处理器（120）通过计算
来产生（S2）扩展的基数，然后计算（S3）扩展的模数
计算（S4）
检验（S5）是否Sr≡1+dar(modr2)，并且当且仅当其如此，经由接口（110）返回（S6）结果S=SrmodN。
- 图形芯片设计中求幂算法的实现-201110202083.9
- 林苍松;陈宝民;陈怒兴 - 长沙景嘉微电子有限公司
- 2011-07-19 - 2011-11-09 - G06F7/552
- 本发明公开了一种图形芯片设计中的求幂算法的实现方式，它分为三个步骤：(1)用查找表法将底数求基于2的对数；(2)将对数与指数相乘；(3)对乘数结果求基于2的指数，本发明中将第一步又划分为三个步骤：(1)浮点数定点化；(2)将基于2的对数查找表多级分解成若干小表；(3)查表结果处理，由于在光照算法中将会不可避免的出现大量求幂计算，而直接求幂计算将给硬件带来难以接受的资源压力，这种实现方式既能保证求幂算法的计算精度，又能大幅度减少芯片的资源使用。
- 使用缩减宽度数据进行计算的方法和装置-200910179405.5
- T·唐;高派·科里 - 马维尔国际贸易有限公司
- 2004-11-17 - 2010-03-24 - G06F7/552
- 公开了使用缩减宽度数据进行计算的方法的装置。具体来说，示例性的方法确定与生成和估算函数相关联的缩减宽度数据值。缩减宽度数据值中的一些在编译阶段被存入指令存储器中的指令中，并且在运行时阶段从指令存储器中取回。
- 用于算术处理的装置、方法和程序-200910159027.4
- 茂木幸彦 - 索尼株式会社
- 2009-07-29 - 2010-02-03 - G06F7/552
- 本发明公开了用于算术处理的装置、方法和程序。尾数/指数分割器将输入值X＝(1+X1/223)x(2^X2)分割为尾数X1和指数X2。插值处理器利用尾数X1来参考尾数/指数分割器并且通过插值处理确定幂值(log2(1+X1/223))。对数计算器根据指数X2和来自插值处理器的幂值来确定对数值Z＝log2XY＝Y(X2+log2(1+X1/223))。整数/分数分割器将对数值Z分割为整数Zint和分数Zamari。插值处理器响应于分数Zamari而参考分数幂次表存储单元并且通过插值处理确定幂值(2^Zamari)。幂计算器确定XY＝2^Z＝(2^Zamari)x(2^Zint)，从而得到输入值X的Y次幂。
- 专利分类
