24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3207  |  回复: 10

fay_L

铁虫 (初入文坛)

[求助] FORTRAN中比双精度(real(kind=8))更高的精度的数据类型是什么?已有4人参与

最近在编程解决问题,但是求矩阵的行列式值时总遇到infinity的情况,发现是精度不够的问题,已换成双精度的还是不行,请问还有比双精度更高精度的数据类型吗?如何定义呢?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

卡尔俾斯

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
绝对不是精度问题,双精度足够用了

[ 发自小木虫客户端 ]
4楼2014-12-20 08:48:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinjinliu

银虫 (著名写手)

引用回帖:
4楼: Originally posted by 卡尔俾斯 at 2014-12-20 08:48:45
绝对不是精度问题,双精度足够用了

我也觉得是,一般都是这样的,人家语言解决了多少问题啊,多少大数字问题,你这就行不通,可定是程序有问题,或者解法本身存在问题。
5楼2014-12-20 09:19:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

引用回帖:
6楼: Originally posted by chairman_gcd at 2014-12-20 20:51:20
FM Fortran 是编译器吗?...

FM Fortran 不是编译器,是一module,可以定义任意长度的实型数据及其运算。
8楼2014-12-21 00:51:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

感谢参与,应助指数 +1
双精度己经是VF最好的精度了。个别编译器有real*16。如果你非要坚持更高精度,可以试试FM Fortran, 那里你可规定任意位,比如200位有效数字

[ 发自小木虫客户端 ]
2楼2014-12-19 14:07:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qchem

铁杆木虫 (著名写手)


【答案】应助回帖

感谢参与,应助指数 +1
需要64位的编译器
gfortran就可以更高

[ 发自小木虫客户端 ]
3楼2014-12-19 16:57:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chairman_gcd

禁虫 (小有名气)

本帖内容被屏蔽

6楼2014-12-20 20:51:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

引用回帖:
4楼: Originally posted by 卡尔俾斯 at 2014-12-20 08:48:45
绝对不是精度问题,双精度足够用了

"绝对不是精度问题",这话就太绝对了。我曾经做过一个时间递推问题。由于问题不稳定,每一时间步大概损失近一位数字的精度,而我需要推进100多步。假如只有双精度,大概经过20步,所有精度损失殆尽,然后就发散了。但我用了FM Fortran, 定义了120位有效数字,就解决了原本是不稳定的数值推进问题。可见,提高有效位数还是有作用的。

[ 发自小木虫客户端 ]
7楼2014-12-21 00:32:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

--八云--

专家顾问 (知名作家)

博览众家之长,终集百家之短

【答案】应助回帖

感谢参与,应助指数 +1
如果你知道了。。。。告诉我一声,我一直以为双精度是fortran最高的精度,除了自定义以外。。。
人原本是不死的,后来他创造了死亡,故而人们歌颂他。
9楼2014-12-21 07:49:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qchem

铁杆木虫 (著名写手)


引用回帖:
8楼: Originally posted by pippi6 at 2014-12-21 00:51:27
FM Fortran 不是编译器,是一module,可以定义任意长度的实型数据及其运算。...

能给个FM Fortran的链接吗?谢谢
10楼2017-02-10 11:15:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fay_L 的主题更新
信息提示
请填处理意见