问个double和long double的问题 - 未名空间MITBBS历史存档

国际科技财经博客移民网络热点娱乐民生时事公众号

Redian新闻

>未名空间

>Programming - 葵花宝典

问个double和long double的问题

问个double和long double的问题# Programming - 葵花宝典

l*i2008-02-01 08:02

1 楼

偶有个做计算的c写的函数库，定义的时候所有的变量类型都是double，但是
因为计算精度的问题，现在想把所有的变量都变成long double，除了挨个改源文件
之外，有没有什么简单点的方法呢（比如compile的时候有没有什么特殊的flags
之类能干这个的）？
系统是64位的linux, compiler有gcc 4.2和icc 10.1(不知道是什么，貌似
gcc比icc慢不少，所以最后可能会用icc)

y*w2008-02-01 08:02

2 楼

用MACRO就可以了
如果一次性的改，用sed也可以

t*t2008-02-01 08:02

3 楼

x86吗? 注意x86上, 计算一律是80位的, double和long double的区别应该就是变量保
留多少

【在 l*********i 的大作中提到】

: 偶有个做计算的c写的函数库，定义的时候所有的变量类型都是double，但是
: 因为计算精度的问题，现在想把所有的变量都变成long double，除了挨个改源文件
: 之外，有没有什么简单点的方法呢（比如compile的时候有没有什么特殊的flags
: 之类能干这个的）？
: 系统是64位的linux, compiler有gcc 4.2和icc 10.1(不知道是什么，貌似
: gcc比icc慢不少，所以最后可能会用icc)

l*i2008-02-01 08:02

4 楼

是x86_64,Intel Core 2 duo，计算一律是80位是啥意思啊，如果一个变
量i (double或者long double)经过很多步算术运算，是说仅仅是最后一步的
的数值的保留位数不同么？

【在 t****t 的大作中提到】

: x86吗? 注意x86上, 计算一律是80位的, double和long double的区别应该就是变量保
: 留多少

l*i2008-02-01 08:02

5 楼

啊，太感谢了，果然，用#define double long double好像就搞定了。

【在 y*w 的大作中提到】

: 用MACRO就可以了
: 如果一次性的改，用sed也可以

t*t2008-02-01 08:02

6 楼

就是说,x86/x86_64的CPU,浮点运算单元就是80位的
你变量可以是96位或者128位,但是后面的啥都没有

【在 l*********i 的大作中提到】

: 啊，太感谢了，果然，用#define double long double好像就搞定了。

c*e2008-02-01 08:02

7 楼

用了这个宏就不用修改函数了？比方说sqrt不用改成sqrtl？
long double sqrt(long double x); [C++ only]
long double sqrtl(long double x); [optional]
C++里面还是用sqrt，但是一般的C里面还要改吗？谢谢！

【在 l*********i 的大作中提到】

: 啊，太感谢了，果然，用#define double long double好像就搞定了。

l*i2008-02-01 08:02

8 楼

哦，原来是这样，长知识了，呵呵。

【在 t****t 的大作中提到】

: 就是说,x86/x86_64的CPU,浮点运算单元就是80位的
: 你变量可以是96位或者128位,但是后面的啥都没有

l*i2008-02-01 08:02

9 楼

不好意思，再请教一下，我看到gcc在x86_64下面缺省的-mfpmath的flag的设置
是sse，在网上google到这么一段:
"sse- take a benefit of SSE2 (SSE3) instruction set (Streaming SIMD
Extension, SIMD - Single Instruction Multiple Data), a default mode for x86-
64 processors, to use this option one should also enable SSE and SSE2 (SSE3)
instruction set to be used in machine code, the FP computations are done
simultaniously on a 16 new 128bit registers (xmm0-xmm15), each of xmm
register makes the same operation on 2 doubles (or 4 floats), the re

【在 t****t 的大作中提到】

: 就是说,x86/x86_64的CPU,浮点运算单元就是80位的
: 你变量可以是96位或者128位,但是后面的啥都没有

t*t2008-02-01 08:02

10 楼

你再仔细读一遍，我觉得这个说得很清楚
SSEx做的是64b或32b的浮点运算，比x87的80位更少
所谓128位是说可以两个64b或者4个32b一起运算，当然会快一些。

x86-
SSE3)
is

【在 l*********i 的大作中提到】

: 不好意思，再请教一下，我看到gcc在x86_64下面缺省的-mfpmath的flag的设置
: 是sse，在网上google到这么一段:
: "sse- take a benefit of SSE2 (SSE3) instruction set (Streaming SIMD
: Extension, SIMD - Single Instruction Multiple Data), a default mode for x86-
: 64 processors, to use this option one should also enable SSE and SSE2 (SSE3)
: instruction set to be used in machine code, the FP computations are done
: simultaniously on a 16 new 128bit registers (xmm0-xmm15), each of xmm
: register makes the same operation on 2 doubles (or 4 floats), the re

l*i2008-02-01 08:02

11 楼

多谢多谢，也就是说如果我比较在意精度，又不想折腾去搞那个无限精度
库的话，最好的选择是用long double再选择x87的FPU？

【在 t****t 的大作中提到】

: 你再仔细读一遍，我觉得这个说得很清楚
: SSEx做的是64b或32b的浮点运算，比x87的80位更少
: 所谓128位是说可以两个64b或者4个32b一起运算，当然会快一些。
:
: x86-
: SSE3)
: is