Fortran Coder

关于fortran调用matlab的问题!

查看数: 3911 | 评论数: 5 | 收藏 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2023-2-10 22:21

正文摘要:

本帖最后由 1033257950 于 2023-2-10 22:21 编辑 今天在尝试使用fortran调用matlab时,卡在了一个函数上:mxGetPr。程序写好后,运行时显示:forrtl: severe (157): Program Exception - access violation。逐行D ...

回复

辰囵1 发表于 2023-8-2 16:21:01
请问压缩包里面的文件您怎么调用呢,在哪儿设置的呢?
1033257950 发表于 2023-2-12 14:54:41
fcode 发表于 2023-2-12 11:39
不好说,我手上没有高版本的matlab,也没有资料。
从编程的角度,对64位系统来说,用于存放地址的变量,肯 ...

好的好的,感谢您的回复!
fcode 发表于 2023-2-12 11:39:28
不好说,我手上没有高版本的matlab,也没有资料。
从编程的角度,对64位系统来说,用于存放地址的变量,肯定是需要kind=8的。
其他的如果不确定,可以先定义成kind=8的,看看它的值会不会超过kind=4的上限(2147483648),如果超过了,那一般就需要kind=8
1033257950 发表于 2023-2-11 09:34:34
fcode 发表于 2023-2-10 23:38
你看的资料应该是针对32位的。而你用的64位,所以略有区别。一些变量需要定义为kind=8的才可以。

看注释

感谢您!问题已解决!我还想问一下,哪些变量需要定义成kind=8,我翻了翻matlab的帮助里,有一部分变量并没有提到必须定义成kind=8,有没有关于这个问题的参考资料?再次感谢!
fcode 发表于 2023-2-10 23:38:36
你看的资料应该是针对32位的。而你用的64位,所以略有区别。一些变量需要定义为kind=8的才可以。

看注释
[Fortran] 纯文本查看 复制代码
01MODULE MATLAB
02  IMPLICIT NONE
03  INTEGER(8)::EP !//kind=8
04  integer :: STATUS
05  INTEGER(8),EXTERNAL:: engOpen ,engPutVariable , engGetVariable ,engEvalString,engClose!//kind=8
06  INTEGER(8),EXTERNAL:: mxCreateDoubleMatrix, mxGetPr!//kind=8
07CONTAINS
08  SUBROUTINE STARTMATLAB()
09    WRITE(*,*) '正在打开matlab程序,请稍后...'
10    EP = engOpen('matlab')
11    IF(EP == 0) THEN
12      WRITE(*,*) '无法打开matlab'
13      STOP
14    ELSE
15      WRITE(*,*) '成功打开matlab'
16    ENDIF
17  END SUBROUTINE STARTMATLAB
18 
19  SUBROUTINE CLOSEMATLAB()
20    WRITE(*,*) '正在关闭matlab程序,请稍后...'
21    STATUS = engClose(EP)
22    IF(STATUS /= 0) THEN
23      WRITE(*,*) '无法关闭matlab'
24      STOP
25    ELSE
26      WRITE(*,*) '成功关闭matlab'
27    ENDIF
28  END SUBROUTINE CLOSEMATLAB
29 
30  SUBROUTINE MC(STRING)
31    CHARACTER(LEN=*)::STRING
32    PRINT * , '正在matlab中计算语句:     ', STRING
33 
34    IF(engEvalString(EP,STRING) /= 0) THEN
35      WRITE(*,*) 'engEvalString failed'
36      STOP
37    ENDIF
38  END SUBROUTINE MC
39 
40  SUBROUTINE F2M(FDATA,MSTRING,ROW,COL)
41    INTEGER :: ROW,COL
42    REAL(8),DIMENSION(ROW,COL) :: FDATA!//kind=8
43    CHARACTER(LEN=*) :: MSTRING
44    INTEGER(8) :: PTEMP!//kind=8
45    PTEMP = mxCreateDoubleMatrix(ROW,COL,0)
46 
47    IF (PTEMP == 0) THEN
48      WRITE(*,*) '无法申请到内存'
49      STOP
50    ENDIF
51    CALL mxCopyReal8ToPtr(fdata,mxGetPr(ptemp),ROW*COL)
52 
53    STATUS = engPutVariable(EP, MSTRING,PTEMP)
54    CALL mxDestroyArray(PTEMP)
55 
56    PRINT * , '正在matlab中生成矩阵:     ', MSTRING
57    IF(STATUS /= 0) THEN
58      WRITE(*,*) 'engPutVariable failed'
59      STOP
60    END IF
61  END SUBROUTINE F2M
62 
63  SUBROUTINE M2F(MSTRING,DDATA,ROW,COL)
64    INTEGER :: ROW,COL
65    REAL(8),DIMENSION(ROW,COL) :: DDATA!//kind=8
66    CHARACTER(LEN=*) :: MSTRING
67    INTEGER(8) :: PTEMP!//kind=8
68    PTEMP = engGetVariable(EP,MSTRING)
69    CALL mxCopyPtrToReal8( mxGetPr(ptemp),DDATA,ROW*COL)
70  END SUBROUTINE M2F
71 
72END MODULE MATLAB
73 
74PROGRAM MAIN
75  USE MATLAB
76  IMPLICIT NONE
77  INTEGER,PARAMETER::NDATA = 10
78  REAL(8),DIMENSION(NDATA)::X,Y,Z
79  INTEGER::M
80  DO M = 1,NDATA
81    X(M) = 3.14159*(M-1)/(NDATA-1)
82    Y(M) = SIN(X(M))
83  ENDDO
84  CALL STARTMATLAB()
85  CALL F2M(X,'X',1,NDATA)
86  CALL F2M(Y,'Y',1,NDATA)
87  CALL MC("plot(X,Y);title('y=sin(x)');xlabel('x');ylabel('y')")
88  CALL MC('Z=2*Y;')!//matlab区分大小写
89  CALL M2F('Z',z,1,NDATA)!//matlab区分大小写
90  PRINT * ,'      X      ' ,'      Y      ','      Z      ','   2*Y'
91  DO M = 1,NDATA
92    WRITE(*,'(G15.6,G15.6,G15.6,G15.6)') X(M),Y(M),Z(M),2*Y(M)
93  ENDDO
94 
95END PROGRAM MAIN

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

GMT+8, 2025-5-5 03:44

Powered by Discuz! X3.4

© 2013-2025 Comsenz Inc.

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