| 查看: 2959 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
dexter337铁杆木虫 (小有名气)
|
[交流]
最近碰到的一个Intel visual Fortran 计算精度的问题.大家的机器上是否也是这样? 已有5人参与
|
||
|
我 以前都是用compaq visual Fortran 6.6的imsl库进行计算的.但是近期需要加大计算的矩阵维数.因此需要改用Intel visual Fortran的mkl函数库.我的程序中用到了zgesv这个函数求解一个复系数的方程组. 原来用compaq visual Fortran的时候用的是imsl中的dlsacg精度很好.但是在Intel visual Fortran上跑的结果误差很严重.于是开始查找出错原因. 后来我发现一个奇怪的现象.就是用imsl手册中的一个计算例子,在Intel 上出的结果精度好像不够.希望大家提提意见,多多讨论,小虫不胜感激. Intel visual Fortran中程序例子如下: program main implicit none integer(8),parameter::n=3 complex(8)::a(n,n),b(n) integer(8)::info,ipiv(n) a=reshape((/(3.0d0,-2.0d0),(1.0d0,1.0d0),(4.0d0,0.0d0),(2.0d0,4.0d0),(2.0d0,-6.0d0),& (-5.0d0,1.0d0),(0.0d0,-3.0d0),(1.0d0,2.0d0),(3.0d0,-2.0d0)/),(/n,n/)) b=(/(10.0d0,5.0d0),(6.0d0,-7.0d0),(-1.0d0,2.0d0)/) call zgesv(n,1,a,n,ipiv,b,n,info) open(2,file='x.txt') write(2,5)b 5 format(2f20.16) end program main 这个例子求解方程a*x=b这样的问题.结果是很整齐的解析结果是x=(1-i,2+i,3i).用compaq visual Fortran的结果是. 1.0000000000000000 -1.0000000000000000 2.0000000000000000 1.0000000000000000 0.0000000000000000 3.0000000000000000 基本上看不出误差.相应的dlsacg函数在我要计算的模型中工作的精度都很好. 可是问题来了,在Intel中(也就是上面的程序.用mkl库给出的结果是. 1.0000000000000024 -1.0000000000000004 1.9999999999999987 1.0000000000000000 -0.0000000000000039 2.9999999999999991 这个结果有些出人意料,希望大家不要笑话我.因为我接下来用到这个函数的地方矩阵的性质比例子中的矩阵差的非常多.计算的结果精度非常不理想. 想问问论坛里的虫友们,大家有没有遇到过这样的问题.有没有什么调节或者校准mkl或者Intel visual Fortran的方法或者设置呢? 我使用的是 Intel visual Fortran 11.1.065专业版.自带的mkl是10.2版的. 在这里先谢谢大家了,祝虫友们圣诞快乐! |
回复此楼» 本帖已获得的红花（最新10朵）
» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
国自然评审意见
已经有17人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助fortran的计算精度
已经有9人回复
- intel visual fortran中MKL如何使用，是用use module还是include？
已经有4人回复
- 对于Fortran计算结果，结果发散？应如何处理？
已经有5人回复
- 【2013版】Intel.Visual.Fortran.Composer.XE.2013.0.089
已经有147人回复
- 求助：intel visual fortran的格式问题
已经有14人回复
- 请热心人能提供一个带imsl库的intel visual fortran编译器下载链接
已经有6人回复
- Intel Visual Fortran应用程序开发
已经有259人回复
- 如何改变fortran的浮点精度
已经有28人回复
- Intel_Visual_Fortran_Compiler_v10_1_019_keygen
已经有121人回复
- 请教一下Visual fortran的project的head files 和 resources files的用处，谢谢
已经有5人回复
- intel visual fortran 多核并行计算
已经有14人回复
- 【分享】PGI Visual Fortran 2010 v11.1.with.VS2010.Shell.X64
已经有5人回复
- 【求助】请问下FORTRAN的计算精度能不能达到MATLAB的计算精度高呢？
已经有3人回复
- 【求助】我下载的INTEL VISUAL FORTRAN FOR WIN中没有.for模板？【已完结】
已经有5人回复
- 【求助】Intel Visual Fortran与Compaq Visual Fortran的区别？
已经有18人回复
- 【求助】visual fortran6.5 的使用
已经有5人回复
dexter337
铁杆木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 7413.5
- 红花: 2
- 帖子: 77
- 在线: 118小时
- 虫号: 492054
- 注册: 2008-01-08
- 性别: GG
- 专业: 基础物理学
5楼2014-12-26 16:19:06
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
2楼2014-12-26 04:28:54
送红花一朵
 |
3楼2014-12-26 08:43:39
4楼2014-12-26 09:17:54
30