Fortran Coder

查看: 10263|回复: 4
打印 上一主题 下一主题

[数学库] 调用IMSL的internal error问题

[复制链接]

2

帖子

1

主题

0

精华

新人

F 币
20 元
贡献
6 点
跳转到指定楼层
楼主
发表于 2015-2-5 14:28:39 | 显示全部楼层 回帖奖励 |倒序浏览 |阅读模式
64位系统 vs2010+ivf2011+imsl7.0
调用imsl函数编程序总是出现第一行错误,不知怎么解决,忘大神解答,感激不尽!!
错误1Compilation Aborted (code 1) c:\users\yzf\documents\visual studio 2010\Projects\Console2\Console2\Source1.f901

具体程序如下:
[Fortran] 纯文本查看 复制代码
program least_square
include 'link_fnl_shared.h' 
use imsl_libraries
integer, parameter :: m = 5
integer, parameter :: n = 2
integer :: i
double precision :: A(m,n), B(m), X(n)
double precision :: C(n, m), D(n,n), E(n)

data A / 2, 8, 2, 7, 4, -1, 4, 1, -1, 0 /
data B / 1, 0, 1, 8, 3 /

if( m > n ) then
C = transpose(a)
D = matmul( c, a )
E = matmul( c, b )
X = D .ix. E ! invert(D) * E
else if( m == n ) then
X = A .ix. B 
else
write(*,*) "No roots!"
end if
do i = 1, n
write(*,*) X(i)
end do
stop
end program least_square


分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

2

帖子

1

主题

0

精华

新人

F 币
20 元
贡献
6 点
沙发
 楼主| 发表于 2015-2-5 16:41:11 | 显示全部楼层
fcode 发表于 2015-2-5 15:40
在我这里,可以编译链接,得到结果:

0.792719919110212

编译后输出窗口是这样,还是看不明白啊,好像说是“if( m > n ) then”这里有问题,但我把if的过程删了,又显示另一处错误,您能帮我分析下吗?
------ Build started: Project: Console2, Configuration: Release|x64 ------

Compiling with Intel(R) Visual Fortran Compiler XE 12.0.5.221 [Intel(R) 64]...
Source1.f90
c:\users\yzf\documents\visual studio 2010\Projects\Console2\Console2\Source1.f90(13): internal error: Please visit 'http://www.intel.com/software/products/support' for assistance.
if( m > n ) then
^
[ Aborting due to internal error. ]
compilation aborted for c:\users\yzf\documents\visual studio 2010\Projects\Console2\Console2\Source1.f90 (code 1)

Build log written to  "file://c:\users\yzf\documents\visual studio 2010\Projects\Console2\Console2\x64\Release\BuildLog.htm"
Console2 - 1 error(s), 0 warning(s)


---------------------- Done ----------------------
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-5-20 13:14

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表