有个问题drive me crazy - 未名空间MITBBS历史存档

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

Redian新闻

>未名空间

>Computation - 科学计算

有个问题drive me crazy

有个问题drive me crazy# Computation - 科学计算

w*n2004-05-12 07:05

1 楼

有谁能在Pentium 的CPU上得到这样的计算结果：
1.0/3.0=0.33333333333333333333 (3的个数多于16位）？
要是谁能写这么一个程序，我就佩服他了。当然是用最简单的办法。不论什么complior
和platform(Suse, redhat, or XP ...). 但是一定要使Pentium 的CPU.

t*e2004-05-12 07:05

2 楼

这是Pentium 4上得结果，正好16个3：
program main
real(8) :: a,b,c
c = 1.d0/3.d0
print'(F20.18)',c
end
0.333333333333333303

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

: 有谁能在Pentium 的CPU上得到这样的计算结果：
: 1.0/3.0=0.33333333333333333333 (3的个数多于16位）？
: 要是谁能写这么一个程序，我就佩服他了。当然是用最简单的办法。不论什么complior
: 和platform(Suse, redhat, or XP ...). 但是一定要使Pentium 的CPU.

w*n2004-05-12 07:05

3 楼

忘记说了，一定要用C++/C

不论什么complior

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

:
: 这是Pentium 4上得结果，正好16个3：
: program main
: real(8) :: a,b,c
: c = 1.d0/3.d0
: print'(F20.18)',c
: end
: 0.333333333333333303

i*w2004-05-12 07:05

4 楼

换UNIX系统。换REAL8. 好像能更多位吧

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

: 忘记说了，一定要用C++/C
:
: 不论什么complior

h*o2004-05-12 07:05

5 楼

Wed May 12 19:52:05 PDT 2004=====>[~/temp]
[email protected]: pts/0: 25 files 37Mb-> cat onethird.c
#include
main()
{
double a, b;
a = 1.e0;
b = 3.e0;
printf("%18.16lf\n", a/b);
}
Wed May 12 19:52:10 PDT 2004=====>[~/temp]
[email protected]: pts/0: 25 files 37Mb-> gcc -o onethird.exe onethird.c
Wed May 12 19:52:13 PDT 2004=====>[~/temp]
[email protected]: pts/0: 25 files 37Mb-> ./onethird.exe
0.3333333333333333
123456789ABCDEF0
十六个三了吧。
Wed May 12 19:52:17 PDT 2004=====>

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

: 忘记说了，一定要用C++/C
:
: 不论什么complior

a*s2004-05-12 07:05

6 楼

meaningless work

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

: 忘记说了，一定要用C++/C
:
: 不论什么complior

w*n2004-05-12 07:05

7 楼

what did u mean?

【在 a******s 的大作中提到】

: meaningless work

h*o2004-05-12 07:05

8 楼

[email protected]: pts/0: 25 files 37Mb-> cat /proc/cpuinfo
processor : 0
vendor_id : GenuineIntel
cpu family : 6
model : 9
model name : Intel(R) Pentium(R) M processor 1400MHz
stepping : 5
cpu MHz : 1395.881
cache size : 1024 KB
fdiv_bug : no
hlt_bug : no
f00f_bug : no
coma_bug : no
fpu : yes
fpu_exception : yes
cpuid level : 2
wp : yes
flags : fpu vme de pse tsc msr mce cx8 sep m

【在 h***o 的大作中提到】

: Wed May 12 19:52:05 PDT 2004=====>[~/temp]
: [email protected]: pts/0: 25 files 37Mb-> cat onethird.c
: #include
: main()
: {
: double a, b;
: a = 1.e0;
: b = 3.e0;
: printf("%18.16lf\n", a/b);
: }

w*n2004-05-12 07:05

9 楼

我是需要多于16个3。如果是double, 最多16个3。

【在 h***o 的大作中提到】

: Wed May 12 19:52:05 PDT 2004=====>[~/temp]
: [email protected]: pts/0: 25 files 37Mb-> cat onethird.c
: #include
: main()
: {
: double a, b;
: a = 1.e0;
: b = 3.e0;
: printf("%18.16lf\n", a/b);
: }

w*n2004-05-12 07:05

10 楼

what are you doing? u can only obtain 16 "3". :)

clflush dts acpi mmx fxsr sse sse2 tm pbe tm2 est

【在 h***o 的大作中提到】

: [email protected]: pts/0: 25 files 37Mb-> cat /proc/cpuinfo
: processor : 0
: vendor_id : GenuineIntel
: cpu family : 6
: model : 9
: model name : Intel(R) Pentium(R) M processor 1400MHz
: stepping : 5
: cpu MHz : 1395.881
: cache size : 1024 KB
: fdiv_bug : no

h*o2004-05-12 07:05

11 楼

那没戏，让(1.0+t == 1.0)为真的最大t值，在pentium机里差不多就是3e-17左右
所以没戏吧。

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

: 我是需要多于16个3。如果是double, 最多16个3。

h*o2004-05-12 07:05

12 楼

look at this one
[email protected]: pts/0: 25 files 37Mb-> gcc -o onethird.exe onethird.c
Wed May 12 20:00:26 PDT 2004=====>[~/temp]
[email protected]: pts/0: 25 files 37Mb-> ./onethird.exe
0.333333333333333314829616256247390992939472198486328125
I am pretty sure there are more than sixteen '3's....hoho

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

: what are you doing? u can only obtain 16 "3". :)
:
: clflush dts acpi mmx fxsr sse sse2 tm pbe tm2 est

t*e2004-05-12 07:05

13 楼

Why C?
I got more than 20 using fortran
program main
real :: c
c = 1./3.
print'(F22.20)',c
end
using ifc7.1
ifc -r16 -pc80 test.f90
./a.out
0.33333333333333333333(20 '3's)

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

: 忘记说了，一定要用C++/C
:
: 不论什么complior

w*n2004-05-12 07:05

14 楼

it depends on ur data type. GSL has extended precision IEEE ... but it's not
working in my PC/Suse.所以我想来这里问问有什么高人知道怎么干。

【在 h***o 的大作中提到】

: 那没戏，让(1.0+t == 1.0)为真的最大t值，在pentium机里差不多就是3e-17左右
: 所以没戏吧。

w*n2004-05-12 07:05

15 楼

这是我最想不通的事情。居然c/c++没有办法得到这样的double precision(extended
precision)？？？？？

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

: Why C?
: I got more than 20 using fortran
: program main
: real :: c
: c = 1./3.
: print'(F22.20)',c
: end
: using ifc7.1
: ifc -r16 -pc80 test.f90
: ./a.out

a*s2004-05-12 07:05

16 楼

-pc80好像没用吧？

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

: Why C?
: I got more than 20 using fortran
: program main
: real :: c
: c = 1./3.
: print'(F22.20)',c
: end
: using ifc7.1
: ifc -r16 -pc80 test.f90
: ./a.out

h*o2004-05-12 07:05

17 楼

16字节的real算起来巨慢，这个我有经验
我要不把那部分改了，再给我5年我也不可能毕业

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

: Why C?
: I got more than 20 using fortran
: program main
: real :: c
: c = 1./3.
: print'(F22.20)',c
: end
: using ifc7.1
: ifc -r16 -pc80 test.f90
: ./a.out

a*s2004-05-12 07:05

18 楼

你看看thinkme是怎么编译的，再看看hanzo那程序写法和编译法，就知道为什么了.:)

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

: 这是我最想不通的事情。居然c/c++没有办法得到这样的double precision(extended
: precision)？？？？？

w*n2004-05-12 07:05

19 楼

no. thinkme 用的是fortran. 和hanzo不是一回事情。

【在 a******s 的大作中提到】

: 你看看thinkme是怎么编译的，再看看hanzo那程序写法和编译法，就知道为什么了.:)

a*s2004-05-12 07:05

20 楼

我c和fortran还是能看出来的:),语言是区别，但是不是本质。

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

: no. thinkme 用的是fortran. 和hanzo不是一回事情。

w*n2004-05-12 07:05

21 楼

我就是找不到用gcc 编译的时候，如果我的data type在程序里是long double, 我是否要
在编译的时候加一些flag.

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

: Why C?
: I got more than 20 using fortran
: program main
: real :: c
: c = 1./3.
: print'(F22.20)',c
: end
: using ifc7.1
: ifc -r16 -pc80 test.f90
: ./a.out

w*n2004-05-12 07:05

22 楼

so, please tell me how to get it work in C, or it's just not going to work in
C. I will really appreicate. :)

.:)

【在 a******s 的大作中提到】

: 我c和fortran还是能看出来的:),语言是区别，但是不是本质。

a*s2004-05-12 07:05

23 楼

hanzo用的是gcc，是没有-r8 or -r16这样的选项的，他程序里面是double，这样就
相当于thinkme用-r8做出来的，而thinkme实际上用的是-r16，当然就比hanzo的高了，
但是-r16也是有限度的，
$ ./a.out
0.3333333333333333333333333333333333173000
解决的办法就是改gcc为icc，这样我估计就可以做到了，程序都是一样的，但是
declare的不同，当然结果就不同了。

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

: so, please tell me how to get it work in C, or it's just not going to work in
: C. I will really appreicate. :)
:
: .:)

t*e2004-05-12 07:05

24 楼

I cant make long double work properly using g++
here is the code
#include
int main()
{
long double a;
a = 1.0/3.0;
printf("%22.20lf\n",a);
return 0;
}
here is the output:
-0.00000000000000000000
which is obviously wrong.
If I use double, the output is:
0.33333333333333331483
still 16 '3's but not more than that

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

: so, please tell me how to get it work in C, or it's just not going to work in
: C. I will really appreicate. :)
:
: .:)

h*o2004-05-12 07:05

25 楼

gcc用long double能搞到19位3
[email protected]: pts/0: 26 files 37Mb-> ./onethird.exe
0.3333333333333333333423683514
加了-m128bit-long-double好象还是这么多

【在 a******s 的大作中提到】

: hanzo用的是gcc，是没有-r8 or -r16这样的选项的，他程序里面是double，这样就
: 相当于thinkme用-r8做出来的，而thinkme实际上用的是-r16，当然就比hanzo的高了，
: 但是-r16也是有限度的，
: $ ./a.out
: 0.3333333333333333333333333333333333173000
: 解决的办法就是改gcc为icc，这样我估计就可以做到了，程序都是一样的，但是
: declare的不同，当然结果就不同了。

h*o2004-05-12 07:05

26 楼

%22.20Lf

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

: I cant make long double work properly using g++
: here is the code
: #include
: int main()
: {
: long double a;
: a = 1.0/3.0;
: printf("%22.20lf\n",a);
: return 0;
: }

w*n2004-05-12 07:05

27 楼

hi, achillis, and thinkme and all:
thank you very much for helping out. it does sound like that i should try
icc. make sense. but what's the flag i shall use?
or any icc manual?

in

【在 a******s 的大作中提到】

w*n2004-05-12 07:05

28 楼

hi, hanzo,你是怎么搞到19位的？？？？with what flag??

，
work in

【在 h***o 的大作中提到】

: gcc用long double能搞到19位3
: [email protected]: pts/0: 26 files 37Mb-> ./onethird.exe
: 0.3333333333333333333423683514
: 加了-m128bit-long-double好象还是这么多

a*s2004-05-12 07:05

29 楼

install icc if you want to use c, ifc for fortran. You will find the document
in the doc directory. It should come with icc.

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

: hi, achillis, and thinkme and all:
: thank you very much for helping out. it does sound like that i should try
: icc. make sense. but what's the flag i shall use?
: or any icc manual?
:
: in

h*o2004-05-12 07:05

30 楼

man了一下gcc就明白了
-m128bit-long-double
These switches control the size of "long double" type. The i386
application binary interface specifies the size to be 96 bits, so
-m96bit-long-double is the default in 32 bit mode.
Modern architectures (Pentium and newer) would prefer "long double"
to be aligned to an 8 or 16 byte boundary. In arrays or structures
conforming to the ABI, this would not be possible. So specifying a
-m

【在 h***o 的大作中提到】

: gcc用long double能搞到19位3
: [email protected]: pts/0: 26 files 37Mb-> ./onethird.exe
: 0.3333333333333333333423683514
: 加了-m128bit-long-double好象还是这么多

a*s2004-05-12 07:05

31 楼

还是没我的多:)

【在 h***o 的大作中提到】

: gcc用long double能搞到19位3
: [email protected]: pts/0: 26 files 37Mb-> ./onethird.exe
: 0.3333333333333333333423683514
: 加了-m128bit-long-double好象还是这么多

h*o2004-05-12 07:05

32 楼

no flag needed....just use long double
printf use %Lf instead of %lf

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

: hi, hanzo,你是怎么搞到19位的？？？？with what flag??
:
: ，
: work in

w*n2004-05-12 07:05

33 楼

no, it's working here:
[email protected]:~/program/PRECISION> more long.c
#include
int main()
{
long double a;
a = 1.0/3.0;
printf("%22.20Lf\n",a);
return 0;
}
[email protected]cp57:~/program/PRECISION> gcc long.c
[email protected]:~/program/PRECISION> ./a.out
0.33333333333333331483

【在 h***o 的大作中提到】

: no flag needed....just use long double
: printf use %Lf instead of %lf

h*o2004-05-12 07:05

34 楼

你这个gcc太土了吧。hoho
[email protected]: pts/0: 26 files 37Mb-> gcc -v
Reading specs from /usr/lib/gcc-lib/i686-pc-linux-gnu/3.3.3/specs
Configured with: /var/tmp/portage/gcc-3.3.3-r3/work/gcc-3.3.3/configure --prefix=/usr --bindir=/usr/i686-pc-linux-gnu/gcc-bin/3.3 --includedir=/usr/lib/gcc-lib/i686-pc-linux-gnu/3.3.3/include --datadir=/usr/share/gcc-data/i686-pc-linux-gnu/3.3

w*n2004-05-12 07:05

35 楼

here it is:
[email protected]:~/program/PRECISION> gcc -v
Reading specs from /usr/lib/gcc-lib/i586-suse-linux/3.3.1/specs
Configured with: ../configure --enable-threads=posix --prefix=/usr

h*o2004-05-12 07:05

36 楼

要么是版本太低，要么是有些开关编译的时候没打开，比如

a*s2004-05-12 07:05

37 楼

这个是有可能的，编译器编译的时候有很多options的。
gentoo宰牛也没有icc牛，呵呵
我试了fedora的gcc
$ gcc a.c
$ ./a.out
0.33333333333333331483
#include
int main()
{
long double a;
a = 1.0/3.0;
printf("%22.20Lf\n",a);
return 0;
}

【在 h***o 的大作中提到】

: 要么是版本太低，要么是有些开关编译的时候没打开，比如

a*s2004-05-12 07:05

38 楼

$ gcc -v
Reading specs from /usr/lib/gcc-lib/i386-redhat-linux/3.3.2/specs
Configured with: ../configure --prefix=/usr --mandir=/usr/share/man --infodir=/usr/share/info --enable-shared --enable-threads=posix --disable-checking --with-system-zlib --enable-__cxa_atexit --host=i386-redhat-linux
Thread model: posix
gcc version 3.3.2 20031022 (Red Hat Linux 3.3.2-1)

【在 a******s 的大作中提到】

:
: 这个是有可能的，编译器编译的时候有很多options的。
: gentoo宰牛也没有icc牛，呵呵
: 我试了fedora的gcc
: $ gcc a.c
: $ ./a.out
: 0.33333333333333331483
: #include
: int main()
: {

w*n2004-05-12 07:05

39 楼

i just tried redhat gcc, it's not working. :(( hehe

【在 a******s 的大作中提到】

: $ gcc -v
: Reading specs from /usr/lib/gcc-lib/i386-redhat-linux/3.3.2/specs
: Configured with: ../configure --prefix=/usr --mandir=/usr/share/man --infodir=/usr/share/info --enable-shared --enable-threads=posix --disable-checking --with-system-zlib --enable-__cxa_atexit --host=i386-redhat-linux
: Thread model: posix
: gcc version 3.3.2 20031022 (Red Hat Linux 3.3.2-1)

a*s2004-05-12 07:05

40 楼

icc比gcc快多了，功能也多不少，尤其是-r选项，这么方便的东西为啥不用呢？
我都说了解决办法了，即使你重装机器，装gentoo的话，gentoo上的gcc也不比
icc快，何苦呢？:)

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

: i just tried redhat gcc, it's not working. :(( hehe

w*n2004-05-12 07:05

41 楼

hehe. achillis,这个问题我搞了很久乐，现在特别好奇.就想try try. i actually
like ur solution best.
我现在就是有点好奇，怎么hanzo的gcc 能work呢？其实20个3对我刚刚好。呵呵。

【在 a******s 的大作中提到】

: icc比gcc快多了，功能也多不少，尤其是-r选项，这么方便的东西为啥不用呢？
: 我都说了解决办法了，即使你重装机器，装gentoo的话，gentoo上的gcc也不比
: icc快，何苦呢？:)

a*s2004-05-12 07:05

42 楼

没什么好奇的，gentoo编译gcc的时候可能打开了那个option，kaka，如果你够坚定，
可以自己编个gcc和hanzo抗衡，kaka

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

: hehe. achillis,这个问题我搞了很久乐，现在特别好奇.就想try try. i actually
: like ur solution best.
: 我现在就是有点好奇，怎么hanzo的gcc 能work呢？其实20个3对我刚刚好。呵呵。

h*o2004-05-12 07:05

43 楼

你难道不知道现在amd的cpu比intel的卖得多？

【在 a******s 的大作中提到】

a*s2004-05-12 07:05

44 楼

做数值计算，我只用intel，kaka，管什么amd? amd现在和intel差多少钱了？

【在 h***o 的大作中提到】

: 你难道不知道现在amd的cpu比intel的卖得多？

a*s2004-05-12 07:05

45 楼

其实，你需不需要那么多位这才是我第一个回帖的意思。这个比icc还重要，kaka

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

f*r2004-05-12 07:05

46 楼

AMD is simply faster :)

【在 a******s 的大作中提到】

: 做数值计算，我只用intel，kaka，管什么amd? amd现在和intel差多少钱了？