关于fortran调用matlab的问题!
本帖最后由 1033257950 于 2023-2-10 22:21 编辑今天在尝试使用fortran调用matlab时,卡在了一个函数上:mxGetPr。程序写好后,运行时显示:forrtl: severe (157): Program Exception - access violation。逐行Debug后发现,问题出现在程序第55行里:
CALL mxCopyPtrToReal8( mxGetPr(ptemp),DDATA,ROW*COL)
在55行设置断点后,发现问题如图2所示,上网找了好久也没发现类似问题,因此求助论坛各位大佬。代码、所需库文件、头文件一并附上,感谢各位!
你看的资料应该是针对32位的。而你用的64位,所以略有区别。一些变量需要定义为kind=8的才可以。
看注释
MODULE MATLAB
IMPLICIT NONE
INTEGER(8)::EP !//kind=8
integer :: STATUS
INTEGER(8),EXTERNAL:: engOpen ,engPutVariable , engGetVariable ,engEvalString,engClose!//kind=8
INTEGER(8),EXTERNAL:: mxCreateDoubleMatrix, mxGetPr!//kind=8
CONTAINS
SUBROUTINE STARTMATLAB()
WRITE(*,*) '正在打开matlab程序,请稍后...'
EP = engOpen('matlab')
IF(EP == 0) THEN
WRITE(*,*) '无法打开matlab'
STOP
ELSE
WRITE(*,*) '成功打开matlab'
ENDIF
END SUBROUTINE STARTMATLAB
SUBROUTINE CLOSEMATLAB()
WRITE(*,*) '正在关闭matlab程序,请稍后...'
STATUS = engClose(EP)
IF(STATUS /= 0) THEN
WRITE(*,*) '无法关闭matlab'
STOP
ELSE
WRITE(*,*) '成功关闭matlab'
ENDIF
END SUBROUTINE CLOSEMATLAB
SUBROUTINE MC(STRING)
CHARACTER(LEN=*)::STRING
PRINT * , '正在matlab中计算语句: ', STRING
IF(engEvalString(EP,STRING) /= 0) THEN
WRITE(*,*) 'engEvalString failed'
STOP
ENDIF
END SUBROUTINE MC
SUBROUTINE F2M(FDATA,MSTRING,ROW,COL)
INTEGER :: ROW,COL
REAL(8),DIMENSION(ROW,COL) :: FDATA!//kind=8
CHARACTER(LEN=*) :: MSTRING
INTEGER(8) :: PTEMP!//kind=8
PTEMP = mxCreateDoubleMatrix(ROW,COL,0)
IF (PTEMP == 0) THEN
WRITE(*,*) '无法申请到内存'
STOP
ENDIF
CALL mxCopyReal8ToPtr(fdata,mxGetPr(ptemp),ROW*COL)
STATUS = engPutVariable(EP, MSTRING,PTEMP)
CALL mxDestroyArray(PTEMP)
PRINT * , '正在matlab中生成矩阵: ', MSTRING
IF(STATUS /= 0) THEN
WRITE(*,*) 'engPutVariable failed'
STOP
END IF
END SUBROUTINE F2M
SUBROUTINE M2F(MSTRING,DDATA,ROW,COL)
INTEGER :: ROW,COL
REAL(8),DIMENSION(ROW,COL) :: DDATA!//kind=8
CHARACTER(LEN=*) :: MSTRING
INTEGER(8) :: PTEMP!//kind=8
PTEMP = engGetVariable(EP,MSTRING)
CALL mxCopyPtrToReal8( mxGetPr(ptemp),DDATA,ROW*COL)
END SUBROUTINE M2F
END MODULE MATLAB
PROGRAM MAIN
USE MATLAB
IMPLICIT NONE
INTEGER,PARAMETER::NDATA = 10
REAL(8),DIMENSION(NDATA)::X,Y,Z
INTEGER::M
DO M = 1,NDATA
X(M) = 3.14159*(M-1)/(NDATA-1)
Y(M) = SIN(X(M))
ENDDO
CALL STARTMATLAB()
CALL F2M(X,'X',1,NDATA)
CALL F2M(Y,'Y',1,NDATA)
CALL MC("plot(X,Y);title('y=sin(x)');xlabel('x');ylabel('y')")
CALL MC('Z=2*Y;')!//matlab区分大小写
CALL M2F('Z',z,1,NDATA)!//matlab区分大小写
PRINT * ,' X ' ,' Y ',' Z ',' 2*Y'
DO M = 1,NDATA
WRITE(*,'(G15.6,G15.6,G15.6,G15.6)') X(M),Y(M),Z(M),2*Y(M)
ENDDO
END PROGRAM MAIN fcode 发表于 2023-2-10 23:38
你看的资料应该是针对32位的。而你用的64位,所以略有区别。一些变量需要定义为kind=8的才可以。
看注释
感谢您!问题已解决!我还想问一下,哪些变量需要定义成kind=8,我翻了翻matlab的帮助里,有一部分变量并没有提到必须定义成kind=8,有没有关于这个问题的参考资料?再次感谢! 不好说,我手上没有高版本的matlab,也没有资料。
从编程的角度,对64位系统来说,用于存放地址的变量,肯定是需要kind=8的。
其他的如果不确定,可以先定义成kind=8的,看看它的值会不会超过kind=4的上限(2147483648),如果超过了,那一般就需要kind=8 fcode 发表于 2023-2-12 11:39
不好说,我手上没有高版本的matlab,也没有资料。
从编程的角度,对64位系统来说,用于存放地址的变量,肯 ...
好的好的,感谢您的回复! 请问压缩包里面的文件您怎么调用呢,在哪儿设置的呢?
页:
[1]