jeydragon 发表于 2018-5-25 22:09:39

IMSL报错

    Programsacling_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. ]

vvt 发表于 2018-5-25 23:28:52

你这代码里哪儿用到 IMSL 了？

jeydragon 发表于 2018-5-27 16:16:40

vvt 发表于 2018-5-25 23:28
你这代码里哪儿用到 IMSL 了？

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

jeydragon 发表于 2018-5-27 16:18:53

vvt 发表于 2018-5-25 23:28
你这代码里哪儿用到 IMSL 了？

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

vvt 发表于 2018-5-27 19:53:46

jeydragon 发表于 2018-5-27 16:16
XT=transpose(X)
      XTX=matmul(XT,X)
      XTY=matmul(XT,Y)

这俩函数都是标准语法的函数，并不是IMSL提供的。
页: [1]
查看完整版本: IMSL报错