24小时热门版块排行榜    

查看: 3601  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 CyTWnZiM0cS6 2026-08-02 3/150 2026-08-02 20:51 by 3VGOm6paSgqv
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 CyTWnZiM0cS6 2026-08-02 3/150 2026-08-02 20:44 by 3VGOm6paSgqv
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 7/350 2026-08-02 20:24 by 3VGOm6paSgqv
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 20:24 by 3VGOm6paSgqv
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 AS6fX1tr1on4 2026-08-02 6/300 2026-08-02 20:21 by 3VGOm6paSgqv
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 12/600 2026-08-02 20:21 by drbart
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 AS6fX1tr1on4 2026-08-02 9/450 2026-08-02 20:18 by 3VGOm6paSgqv
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 AS6fX1tr1on4 2026-08-02 14/700 2026-08-02 20:11 by 3VGOm6paSgqv
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 AS6fX1tr1on4 2026-08-02 11/550 2026-08-02 20:11 by 3VGOm6paSgqv
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 7/350 2026-08-02 20:11 by 3VGOm6paSgqv
[基金申请] 面上提前没消息,有中的吗 +11 archvillain 2026-08-02 13/650 2026-08-02 20:09 by symmetry
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 cBGhfJFuNMca 2026-08-02 13/650 2026-08-02 20:04 by 3VGOm6paSgqv
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 cBGhfJFuNMca 2026-08-02 17/850 2026-08-02 19:58 by 3VGOm6paSgqv
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 cBGhfJFuNMca 2026-08-02 21/1050 2026-08-02 19:52 by 3VGOm6paSgqv
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
请填处理意见