设为首页收藏本站如何选择编译器？在线F77转F90繁体浏览
 切换到窄版 开启辅助访问

Fortran Coder

Fortran Coder»论坛 › 研讨区 › 编程工具交流 › 调用MKL_DSS无法保存结果到文件
返回列表 发新帖
查看: 10013|回复: 5
打印 上一主题 下一主题

[数学库] 调用MKL_DSS无法保存结果到文件

[复制链接]
duty

21

帖子

6

主题

0

精华

入门

F 币
83 元
贡献
50 点
跳转到指定楼层
楼主
发表于 2014-4-7 11:02:30 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

错误提示如图
代码如下
[Fortran] 纯文本查看 复制代码
    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
      
7001  FORMAT(1X,F)     
7002  FORMAT(1X,I5)
7003  FORMAT(1X,I,1X,E)
7004  FORMAT(1X,I,1X,I)
7005  FORMAT(1X,D)
      
999   WRITE(*,*) "Solver returned error code ", error
      
1999  DEALLOCATE(tIU,tJU,tUM,tXDIFF,tX)
      
 END



附件为所需输入文件~

TestData.zip

14.23 KB, 下载次数: 4

测试输入文件
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩
回复

使用道具 举报

fcode实名认证

2022

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1598 元
贡献
689 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章
沙发
发表于 2014-4-7 21:50:13 | 只看该作者
这确实是一个很诡异的情况。

IVF 的错误列表上这样说：
severe (8):  Internal consistency  check  failure
FOR$IOS_BUG_CHECK.   Internal  error.   Please  check  that  the program is correct .  Recompile  if  an error existed  in  the  program.   If  this error persists, submit a problem report .


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

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

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

[Fortran] 纯文本查看 复制代码
      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
      
7001  FORMAT(1X,F)     
7002  FORMAT(1X,I5)
7003  FORMAT(1X,I,1X,E)
7004  FORMAT(1X,I,1X,I)
7005  FORMAT(1X,D)
      
999   WRITE(*,*) "Solver returned error code ", error
      
1999  DEALLOCATE(tIU,tJU,tUM,tXDIFF,tX)
      write(*,*) 'Program End' !// ******
      END
回复 支持 反对

使用道具 举报

duty

21

帖子

6

主题

0

精华

入门

F 币
83 元
贡献
50 点
板凳
 楼主| 发表于 2014-4-7 22:51:16 | 只看该作者
fcode 发表于 2014-4-7 21:50
这确实是一个很诡异的情况。

IVF 的错误列表上这样说：

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

使用道具 举报

duty

21

帖子

6

主题

0

精华

入门

F 币
83 元
贡献
50 点
地板
 楼主| 发表于 2014-4-10 04:41:02 | 只看该作者
程序还是不能终止啊~这是内存crash么？
回复 支持 反对

使用道具 举报

楚香饭

725

帖子

4

主题

0

精华

大师

农村外出务工人员

F 币
657 元
贡献
337 点

新人勋章爱心勋章水王勋章元老勋章热心勋章
5#
发表于 2014-4-10 17:17:46 | 只看该作者
这是IVF的BUG吧
回复 支持 反对

使用道具 举报

aliouying

136

帖子

3

主题

0

精华

版主

F 币
1964 元
贡献
1677 点

帅哥勋章管理勋章爱心勋章新人勋章热心勋章元老勋章
6#
发表于 2014-4-17 09:30:47 | 只看该作者
为什么不用pardiso呢？
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

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

GMT+8, 2024-11-24 02:01

Powered by Tencent X3.4

© 2013-2024 Tencent

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