Fortran Coder

查看: 76|回复: 4

[数学库] IMSL报错

[复制链接]

17

帖子

2

主题

0

精华

入门

F 币
84 元
贡献
56 点
发表于 2018-5-25 22:09:39 | 显示全部楼层 |阅读模式
[Fortran] 纯文本查看 复制代码
    Program  sacling_factor  
    include 'link_fnl_shared.h' 
    use imsl_libraries
    implicit none
    integer,parameter ::N=180          
    integer :: i 
    real*8 ::X(N,1),Y(N,1)
    real*8 ::XT(1,N)
    real*8 XTX(1,1),XTY(1,1)
    open(11,file='1.txt')
   ! open(14,file='1.txt',status='old')      
    open(12,file='2.txt',status='old')          
    open(13,file='The_result.txt',status='replace')
      do i=1,N    
            read(11,'(10X,D15.3)') X(i,1)
            write(*,*) X(i,1)
            read(12,'(10X,D15.3)') Y(i,1)
      enddo
  ! Y=XA
  ! we want to compulate the parameter A
  ! XTY=XTXA====>A=(XTY)/(XTX)

        XT=transpose(X)
        XTX=matmul(XT,X)
        XTY=matmul(XT,Y)
     write(13,*) (XTY)/(XTX)

    end program sacling_factor

报错信息为 [ Aborting due to internal error. ]

回复

使用道具 举报

537

帖子

0

主题

0

精华

大师

F 币
629 元
贡献
288 点

规矩勋章元老勋章新人勋章水王勋章

QQ
发表于 2018-5-25 23:28:52 | 显示全部楼层
你这代码里哪儿用到 IMSL 了?

17

帖子

2

主题

0

精华

入门

F 币
84 元
贡献
56 点
 楼主| 发表于 2018-5-27 16:16:40 | 显示全部楼层
vvt 发表于 2018-5-25 23:28
你这代码里哪儿用到 IMSL 了?

       XT=transpose(X)
        XTX=matmul(XT,X)
        XTY=matmul(XT,Y)

17

帖子

2

主题

0

精华

入门

F 币
84 元
贡献
56 点
 楼主| 发表于 2018-5-27 16:18:53 | 显示全部楼层
vvt 发表于 2018-5-25 23:28
你这代码里哪儿用到 IMSL 了?

发现把前面include及use删掉竟然出结果了,之前的电脑上编译竟然可以通过正常运行的

537

帖子

0

主题

0

精华

大师

F 币
629 元
贡献
288 点

规矩勋章元老勋章新人勋章水王勋章

QQ
发表于 2018-5-27 19:53:46 | 显示全部楼层
jeydragon 发表于 2018-5-27 16:16
XT=transpose(X)
        XTX=matmul(XT,X)
        XTY=matmul(XT,Y)

这俩函数都是标准语法的函数,并不是IMSL提供的。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2018-6-25 12:23

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

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