2

我正在编写一个 MEX 文件，其中涉及从 Matlab 读取稀疏矩阵，但每次运行程序时，Matlab 都会崩溃。然后我写了一个测试文件，发现问题正好出在mxGetIr命令上。我不知道为什么会这样。请帮助我...非常感谢！

我在 Mac OS X 10.8.1 的 Matlab 2013a 下工作。

#include "fintrh.h"

  subroutine mexFunction(nlhs, plhs, nrhs, prhs)

  mwpointer plhs(*), prhs(*)
  integer nlhs, nrhs

  integer temp, nzmax
  integer, dimension(2) :: dimx

  real, dimension (:), allocatable :: sr
  integer, dimension (:), allocatable :: irs, jcs

  mwpointer mxGetPr
  mwpointer mxCreateDoubleMatrix
  mwsize  mxGetM, mxGetN, mxGetNzmax

  mwpointer mxGetIr, mxGetJc

  mwpointer temp_pr, spr

  dimx(1) = mxGetM(prhs(1))
  dimx(2) = mxGetN(prhs(1))
  nzmax = mxGetNzmax(prhs(1))
  allocate(sr(1:nzmax))
  allocate(irs(1:nzmax))
  allocate(jcs(1:(dimx(2)+1)))
  temp_pr = mxGetPr(prhs(1))
  call getreal(temp_pr,sr,nzmax)

  temp_pr = mxGetJc(prhs(1))

  plhs(1) = mxCreateDoubleMatrix(nzmax,1,0)
  temp_pr = mxGetPr(plhs(1))
  call putreal(sr, temp_pr, nzmax)

  end

  subroutine real8toreal(x, y, size)
  integer size
  real*8 x(size)
  real y(size)
  do 10 i=1,size
     y(i)= x(i)
  10   continue
  return
  end

  subroutine getreal(pr,x,size)
  mwpointer pr
  integer size
  real x(size)
  real*8, dimension (:), allocatable :: temp
  allocate(temp(1:size))
  call mxCopyPtrToReal8(pr,temp,size)
  call real8toreal(temp,x,size)
  deallocate(temp)      
  return
  end

  subroutine putreal(x,pr,size)
  mwpointer pr
  integer size
  real x(size)
  real*8, dimension (:), allocatable :: temp
  allocate(temp(1:size))
  call realtoreal8(x,temp,size)
  call mxCopyReal8ToPtr(temp,pr,size)
  deallocate(temp)      
  return
  end
4

1 回答 1

0

检查文档后，我认为您可能需要以不同的方式调用它。

从那里我猜是这样的：

mwPointer mxGetIr(pm)
mwPointer pm
于 2013-09-24T11:47:33.813 回答