Question

我正在编写这个F90程序来计算fortran中的函数，该函数从.mat文件中获取输入并将结果保存在另一个.mat文件中。

我按照this回答来编译并正确链接代码。这是我的makefile命令：

gfortran -g -fcheck=all binhkorn_mat.F90 -I/usr/local/MATLAB/R2015b/extern/include/ -L/usr/local/MATLAB/R2015b/bin/glnxa64 -cpp -o binhkorn_mat -lmat -lmx -Wl,-rpath /usr/local/MATLAB/R2015b/bin/glnxa64/

输出文件显然已正确编译，但是一旦我运行程序，就会出现以下SF（我正在使用LINUX Ubuntu 14.04 LTS）：

    Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0  0x7FD650063F27
#1  0x7FD6500644F4
#2  0x7FD64FCBCD3F
#3  0x7FD64FDD7AF6
#4  0x400A3F in binhkorn_mat at binhkorn_mat.F90:17 (discriminator 2)
 ./binhkorn_mat: Segmentation fault

我无法弄清楚编译器是否存在错误，或者我是否对指针/函数定义做了错误。这是代码（binhkorn_mat.F90）：

 #include "fintrf.h"
PROGRAM binhkorn_mat
IMPLICIT NONE
mwPointer matOpen, matGetVariable, matPutVariable
mwPointer mpin, mpX, mpout, mpcf
INTEGER :: i,j
REAL*8, DIMENSION(2) :: x
REAL*8, DIMENSION(4) :: cf

!input/output through .mat f
mpin = matOpen('X.mat', 'u')
mpX = matGetVariable(mpin, 'X')
CALL mxCopyPtrToReal8(mpX, x, 2)
CALL matClose(mpin)

!fitness functions
 cf(1) = ((x(1)-2)**2 + (x(2)-1)**2 + 2)
 cf(2) = (9*x(1) + (x(2)-1)**2)

!constraints
 cf(3) = x(1)*x(1) + x(2)*x(2) - 225
 cf(4) = x(1) - 3*x(2) + 10

!output file created
CALL mxCopyReal8ToPtr(cf, mpcf, 4)
mpout = matOpen('cf.mat', 'w')
mpcf = matPutVariable(mpout, 'cf', mpcf)
CALL matClose(mpout)

END PROGRAM

X.mat文件由外部Matlab脚本正确创建，并包含一个名为X的变量，它是一个2元素的行向量。

Answer 1

我基本上误解了如何使用很多功能。我作为输入提供给它们中的一些指针不是正确的。我在这里发布了工作解决方案：

#include "fintrf.h"
PROGRAM binhkorn_mat
IMPLICIT NONE
mwPointer matOpen, matGetVariable!, matPutVariable
mwPointer mxGetData, mxGetNumberOfElements, mxCreateNumericArray
mwPointer mpin, mpX, mpout, mpcf
mwSize ndim
mwSize dims(2)
INTEGER :: s
INTEGER*4 mxClassIDFromClassName
 CHARACTER (LEN = 6) :: classname
REAL*8, DIMENSION(2) :: x
REAL*8, DIMENSION(4) :: cf

!input/output through .mat f
mpin = matOpen('X.mat', 'r')
mpX = matGetVariable(mpin, 'X')
CALL mxCopyPtrToReal8(mxGetData(mpX), x, mxGetNumberOfElements(mpX))
!CALL matClose(mpin)

!fitness functions
 cf(1) = ((x(1)-2)**2 + (x(2)-1)**2 + 2)
 cf(2) = (9*x(1) + (x(2)-1)**2)

!constraints
 cf(3) = x(1)*x(1) + x(2)*x(2) - 225
 cf(4) = x(1) - 3*x(2) + 10

!output .mat file created and filled
s = size(cf)
ndim = 2
 classname = 'double'
dims(1) = 1
dims(2) = s
mpcf = mxCreateNumericArray(ndim, dims, mxClassIDFromClassName(classname), 0)
CALL mxCopyReal8ToPtr(cf, mxGetData(mpcf), mxGetNumberOfElements(mpcf))
mpout = matOpen('cf.mat', 'w')
CALL matPutVariable(mpout, 'cf', mpcf)
!CALL matClose(mpout)

END PROGRAM

在Fortran中导入.mat文件 - 分段错误错误

1 个答案: