• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

php 的加法

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

无意间看到了php中关于加,减,乘,除 的计算方法

这里

http://lxr.php.net/source/xref/PHP-5.6/Zend/zend_operators.h#596

 

static zend_always_inline int fast_add_function(zval *result, zval *op1, zval *op2 TSRMLS_DC)
597{
598    if (EXPECTED(Z_TYPE_P(op1) == IS_LONG)) {
599        if (EXPECTED(Z_TYPE_P(op2) == IS_LONG)) {
600#if defined(__GNUC__) && defined(__i386__)
601        __asm__(
602            "movl    (%1), %%eax\n\t"
603            "addl   (%2), %%eax\n\t"
604            "jo     0f\n\t"
605            "movl   %%eax, (%0)\n\t"
606            "movb   %3, %c5(%0)\n\t"
607            "jmp    1f\n"
608            "0:\n\t"
609            "fildl    (%1)\n\t"
610            "fildl    (%2)\n\t"
611            "faddp    %%st, %%st(1)\n\t"
612            "movb   %4, %c5(%0)\n\t"
613            "fstpl    (%0)\n"
614            "1:"
615            :
616            : "r"(&result->value),
617              "r"(&op1->value),
618              "r"(&op2->value),
619              "n"(IS_LONG),
620              "n"(IS_DOUBLE),
621              "n"(ZVAL_OFFSETOF_TYPE)
622            : "eax","cc");
623#elif defined(__GNUC__) && defined(__x86_64__)
624        __asm__(
625            "movq    (%1), %%rax\n\t"
626            "addq   (%2), %%rax\n\t"
627            "jo     0f\n\t"
628            "movq   %%rax, (%0)\n\t"
629            "movb   %3, %c5(%0)\n\t"
630            "jmp    1f\n"
631            "0:\n\t"
632            "fildq    (%1)\n\t"
633            "fildq    (%2)\n\t"
634            "faddp    %%st, %%st(1)\n\t"
635            "movb   %4, %c5(%0)\n\t"
636            "fstpl    (%0)\n"
637            "1:"
638            :
639            : "r"(&result->value),
640              "r"(&op1->value),
641              "r"(&op2->value),
642              "n"(IS_LONG),
643              "n"(IS_DOUBLE),
644              "n"(ZVAL_OFFSETOF_TYPE)
645            : "rax","cc");
646#else
647            /*
648             * 'result' may alias with op1 or op2, so we need to
649             * ensure that 'result' is not updated until after we
650             * have read the values of op1 and op2.
651             */
652
653            if (UNEXPECTED((Z_LVAL_P(op1) & LONG_SIGN_MASK) == (Z_LVAL_P(op2) & LONG_SIGN_MASK)
654                && (Z_LVAL_P(op1) & LONG_SIGN_MASK) != ((Z_LVAL_P(op1) + Z_LVAL_P(op2)) & LONG_SIGN_MASK))) {
655                Z_DVAL_P(result) = (double) Z_LVAL_P(op1) + (double) Z_LVAL_P(op2);
656                Z_TYPE_P(result) = IS_DOUBLE;
657            } else {
658                Z_LVAL_P(result) = Z_LVAL_P(op1) + Z_LVAL_P(op2);
659                Z_TYPE_P(result) = IS_LONG;
660            }
661#endif
662            return SUCCESS;
663        } else if (EXPECTED(Z_TYPE_P(op2) == IS_DOUBLE)) {
664            Z_DVAL_P(result) = ((double)Z_LVAL_P(op1)) + Z_DVAL_P(op2);
665            Z_TYPE_P(result) = IS_DOUBLE;
666            return SUCCESS;
667        }
668    } else if (EXPECTED(Z_TYPE_P(op1) == IS_DOUBLE)) {
669        if (EXPECTED(Z_TYPE_P(op2) == IS_DOUBLE)) {
670            Z_DVAL_P(result) = Z_DVAL_P(op1) + Z_DVAL_P(op2);
671            Z_TYPE_P(result) = IS_DOUBLE;
672            return SUCCESS;
673        } else if (EXPECTED(Z_TYPE_P(op2) == IS_LONG)) {
674            Z_DVAL_P(result) = Z_DVAL_P(op1) + ((double)Z_LVAL_P(op2));
675            Z_TYPE_P(result) = IS_DOUBLE;
676            return SUCCESS;
677        }
678    }
679    return add_function(result, op1, op2 TSRMLS_CC);
680}

其中第653行中的宏 LONG_SIGN_MASK 定义为

#define LONG_SIGN_MASK (1L << (8*sizeof(long)-1))

 在64位机下,LONG_SIGN_MASK的值为 1L<< (8*8-1) = 1L<<63 = 2^63次方

 

 

Intel IA32 浮点运算

IA32处理器和很多其它一些处理器一样,有专门用于保存浮点数的寄存器,当在cpu中进行浮点数运算时,这些寄存器就用来保存输入输出及相关的中间结果。

但IA32有一个比较特别的地方,它的浮点数寄存器是80位的,而我们在程序中只用到32和64位两种类型,因此当把float,double放入到cpu中时,它们都会先被转换成了80位,然后以80位的方式进行运算,最后得到的结果再转换回来。这样特性使得浮点数的计算可以相对更精确些,但同时,一不小心很可能也会引出一些意想不到的问题。

你可能突然恍然大悟了,对的,我们最开始提到那个奇怪的问题就与此相关。

s/e得到结果是个80位的浮点数,由这个浮点数先转换成double再转成int,与直接就转换成int,结果很可能是不同的。

比如在我们的例子中,s/e ~ 29.999999....时,s/e转换成double使用round-to-even的方式,会得到也许是30.0000001,再转成整形时,得到30.

但如果直接由29.99999...转换成整型,得到却是29。

 

后来新出的系列Intel处理器,包括IA32及64位的处理器,提供了专门的硬件来直接处理浮点数,使得可以分开对待float型与double型,这些硬件特性在compiler的支持下,可以生相对高效的代码,同时也避免了我们上面所遇到的问题,有兴趣的读者可以google一下相关的关键字:sse。

 

 

 

 在64位机下
int           4个字节

long         8个字节

double        8个字节

float         4个字节

double long  8个字节

指针         8个字节

而在32位机下

int            4个字节

long         4个字节

double        8个字节

float         4个字节

double long  8个字节

指针                 4个字节

 


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
php与java通用AES加密解密算法发布时间:2022-07-10
下一篇:
PHP5与IIS的集成(Windows2003)发布时间:2022-07-10
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap