我有一个工作代码，可以通过使用 Petsc 的 KSPSolve 解决某个问题（使用 gmres 求解器和 ilu 预处理器）。我想从 MATLAB 调用这段代码，所以我编写了一个 mex 文件。但是，当调用 KSPSolve 时，它会使我的 mex 文件和 MATLAB 崩溃。
我使用了 valgrind，似乎第一个错误出现在某些 sse2 代码中。如果我删除对 KSPSolve 的调用，一切正常。如果我在 MATLAB 之外调用整个代码，它也可以工作。
有人知道 Petsc 和 MATLAB 之间可能存在的兼容性问题吗？这可以解决吗？谢谢，