duty 发表于 2014-4-7 11:02:30

调用MKL_DSS无法保存结果到文件


错误提示如图
代码如下
    PROGRAM Test
      INCLUDE 'mkl_dss.f77'

      INTEGER tIU(:), tJU(:), Idum(1), ITER
      DOUBLE PRECISION tUM(:), tXDIFF(:), tX(:)
      ALLOCATABLE tIU, tJU, tUM, tXDIFF, tX
      
      INTEGER nNonZeros,nRow
      
      INTEGER*8 handle
      INTEGER error

      OPEN (UNIT=10, FILE='testInfo.dat')
      READ (10,7004) nNonZeros, nRow
      CLOSE(10)
   
      ALLOCATE(tIU(nRow+1))
      ALLOCATE(tXDIFF(nRow))
      ALLOCATE(tX(nRow))
      ALLOCATE(tJU(nNonZeros))
      ALLOCATE(tUM(nNonZeros))

      
      OPEN (UNIT=10, FILE='columns.dat')
      READ (10,7002) (tJU(II),II=1,nNonZeros)
      CLOSE(10)
      
      OPEN (UNIT=10, FILE='rowIndex.dat')
      READ (10,7002) (tIU(II),II=1,nRow+1)
      CLOSE(10)
      
      OPEN (UNIT=10, FILE='values.dat')
      READ (10,7001) (tUM(II),II=1,nNonZeros)
      CLOSE(10)
      
      OPEN (UNIT=10, FILE='loadVector.dat')
      READ (10,7001) (tXDIFF(II),II=1,nRow)
      CLOSE(10)   
      


      ITER=1
C---------------------------------------------------------------------------
C Initialize the solver
C---------------------------------------------------------------------------
      error = dss_create(handle, MKL_DSS_DEFAULTS+MKL_DSS_OOC_VARIABLE)
      IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
C---------------------------------------------------------------------------

C---------------------------------------------------------------------------
C Define the non-zero structure of the matrix
C---------------------------------------------------------------------------
      error = dss_define_structure( handle, MKL_DSS_NON_SYMMETRIC,
   & tIU, nRow, nRow, tJU, nNonZeros )
      IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
C---------------------------------------------------------------------------

C---------------------------------------------------------------------------
C Reorder the matrix
C---------------------------------------------------------------------------

      IF (ITER.EQ.1) THEN
          error = dss_reorder( handle, MKL_DSS_GET_ORDER, Idum)
          IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
      ELSE
          error = dss_reorder( handle, MKL_DSS_MY_ORDER, Idum)
          IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
      ENDIF

C---------------------------------------------------------------------------
C Factor the matrix
C---------------------------------------------------------------------------
      error = dss_factor_real( handle, MKL_DSS_DEFAULTS, tUM)
      IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
      
C---------------------------------------------------------------------------
C Get the solution vector
C---------------------------------------------------------------------------
      error = dss_solve_real( handle, MKL_DSS_DEFAULTS,
   & tXDIFF, 1, tX)
      IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
C---------------------------------------------------------------------------
   
C---------------------------------------------------------------------------
C Deallocate solver storage
C---------------------------------------------------------------------------
      error = dss_delete( handle, MKL_DSS_DEFAULTS )
      IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
C---------------------------------------------------------------------------

      OPEN (UNIT=10, FILE='OutValues.dat')
      WRITE (10,7005) (tX(II),II=1,nRow)
      CLOSE(10)
      GOTO 1999
      
7001FORMAT(1X,F)   
7002FORMAT(1X,I5)
7003FORMAT(1X,I,1X,E)
7004FORMAT(1X,I,1X,I)
7005FORMAT(1X,D)
      
999   WRITE(*,*) "Solver returned error code ", error
      
1999DEALLOCATE(tIU,tJU,tUM,tXDIFF,tX)
      
END


附件为所需输入文件~

fcode 发表于 2014-4-7 21:50:13

这确实是一个很诡异的情况。

IVF 的错误列表上这样说：
severe (8):Internal consistencycheckfailure
FOR$IOS_BUG_CHECK.   Internalerror.   Pleasecheckthatthe program is correct .Recompileifan error existedintheprogram.   Ifthis error persists, submit a problem report .

我想这个错误应该是 IVF 内部的缺陷导致的。而不是任何语法错误或调用方式错误。

通过更底层的一些逆向分析，得到一个比较奇特的“不正式”的解决方法：

（其中 !// ****** 是我修改后的行，不要问为什么，我也不知道，呵呵）

      PROGRAM Test
      INCLUDE 'mkl_dss.f77'

      INTEGER tIU(:), tJU(:), Idum(1), ITER
      DOUBLE PRECISION tUM(:), tXDIFF(:), tX(:)
      ALLOCATABLE tIU, tJU, tUM, tXDIFF, tX
      real*8 :: av(3) = 58.0
      INTEGER nNonZeros,nRow
      
      INTEGER*8 handle
      INTEGER error

      OPEN (UNIT=10, FILE='testInfo.dat')
      READ (10,7004) nNonZeros, nRow
      CLOSE(10)
   
      ALLOCATE(tIU(nRow+1))
      ALLOCATE(tXDIFF(nRow))
      ALLOCATE(tX(nRow))
      ALLOCATE(tJU(nNonZeros))
      ALLOCATE(tUM(nNonZeros))
      OPEN (UNIT=11, FILE='OutValues.dat') !// ******
      write(11,*)!// ******
      write(*,*) 'Program Begin'!// ******
      OPEN (UNIT=10, FILE='columns.dat')
      READ (10,7002) (tJU(II),II=1,nNonZeros)
      CLOSE(10)
      
      OPEN (UNIT=10, FILE='rowIndex.dat')
      READ (10,7002) (tIU(II),II=1,nRow+1)
      CLOSE(10)
      
      OPEN (UNIT=10, FILE='values.dat')
      READ (10,7001) (tUM(II),II=1,nNonZeros)
      CLOSE(10)
      
      OPEN (UNIT=10, FILE='loadVector.dat')
      READ (10,7001) (tXDIFF(II),II=1,nRow)
      CLOSE(10)   
      


      ITER=1
C---------------------------------------------------------------------------
C Initialize the solver
C---------------------------------------------------------------------------
      error = dss_create(handle, MKL_DSS_DEFAULTS+MKL_DSS_OOC_VARIABLE)
      IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
C---------------------------------------------------------------------------

C---------------------------------------------------------------------------
C Define the non-zero structure of the matrix
C---------------------------------------------------------------------------
      error = dss_define_structure( handle, MKL_DSS_NON_SYMMETRIC,
   & tIU, nRow, nRow, tJU, nNonZeros )
      IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
C---------------------------------------------------------------------------

C---------------------------------------------------------------------------
C Reorder the matrix
C---------------------------------------------------------------------------

      IF (ITER.EQ.1) THEN
          error = dss_reorder( handle, MKL_DSS_GET_ORDER, Idum)
          IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
      ELSE
          error = dss_reorder( handle, MKL_DSS_MY_ORDER, Idum)
          IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
      ENDIF

C---------------------------------------------------------------------------
C Factor the matrix
C---------------------------------------------------------------------------
      error = dss_factor_real( handle, MKL_DSS_DEFAULTS, tUM)
      IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
      
C---------------------------------------------------------------------------
C Get the solution vector
C---------------------------------------------------------------------------
      error = dss_solve_real( handle, MKL_DSS_DEFAULTS,
   & tXDIFF, 1, tX)
      IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
C---------------------------------------------------------------------------
   
C---------------------------------------------------------------------------
C Deallocate solver storage
C---------------------------------------------------------------------------
      
C---------------------------------------------------------------------------

      Rewind(11)!// ******
      WRITE (11,7005) (tX(II),II=1,nRow) !// ******
      CLOSE(11) !// ******
      error = dss_delete( handle, MKL_DSS_DEFAULTS )
      IF (error .NE. MKL_DSS_SUCCESS ) GOTO 999
      GOTO 1999
      
7001FORMAT(1X,F)   
7002FORMAT(1X,I5)
7003FORMAT(1X,I,1X,E)
7004FORMAT(1X,I,1X,I)
7005FORMAT(1X,D)
      
999   WRITE(*,*) "Solver returned error code ", error
      
1999DEALLOCATE(tIU,tJU,tUM,tXDIFF,tX)
      write(*,*) 'Program End' !// ******
      END

duty 发表于 2014-4-7 22:51:16

fcode 发表于 2014-4-7 21:50
这确实是一个很诡异的情况。

IVF 的错误列表上这样说：


谢谢管理员啊~确实可以有结果输出啦~
可是程序为啥终止不了啊~

duty 发表于 2014-4-10 04:41:02

程序还是不能终止啊~这是内存crash么？

楚香饭 发表于 2014-4-10 17:17:46

这是IVF的BUG吧

aliouying 发表于 2014-4-17 09:30:47

为什么不用pardiso呢？
页: [1]
查看完整版本: 调用MKL_DSS无法保存结果到文件