| 查看: 2566 | 回复: 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朵）
» 猜你喜欢
材料调剂
已经有4人回复
-
面上模板改不了页边距吧？
已经有6人回复
-
307求调剂
已经有6人回复
-
304求调剂
已经有5人回复
-
317一志愿华南理工电气工程求调剂
已经有8人回复
-
272求调剂
已经有3人回复
-
化工专硕348，一志愿985求调剂
已经有6人回复
-
292求调剂
已经有3人回复
-
290求调剂
已经有6人回复
-
295求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助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人回复
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
dexter337
铁杆木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 7413.5
- 红花: 2
- 帖子: 77
- 在线: 118小时
- 虫号: 492054
- 注册: 2008-01-08
- 性别: GG
- 专业: 基础物理学
5楼2014-12-26 16:19:06
dexter337
铁杆木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 7413.5
- 红花: 2
- 帖子: 77
- 在线: 118小时
- 虫号: 492054
- 注册: 2008-01-08
- 性别: GG
- 专业: 基础物理学
6楼2014-12-26 16:30:11
dexter337
铁杆木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 7413.5
- 红花: 2
- 帖子: 77
- 在线: 118小时
- 虫号: 492054
- 注册: 2008-01-08
- 性别: GG
- 专业: 基础物理学
您都已经老师了,别跑出来捣乱哈.
7楼2014-12-26 16:31:10
8楼2016-01-17 16:20:34
问天阁
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 5173.8
- 红花: 5
- 帖子: 577
- 在线: 434小时
- 虫号: 1198420
- 注册: 2011-01-31
- 性别: GG
- 专业: 物理无机化学
9楼2016-01-18 12:38:06
