Mathematica 使用C语言编写库函数教程 带示例

其实帮助文档 LibraryLink/tutorial/Overview（或者在线的）写的很详细了，这里只说一下重点

• 用 MArgument_getXXXX 获得函数参数
• 用 MArgument_setXXXX 设置返回值
• 编译链接 WolframRTL.dll

比如

f2.c

#include <math.h>

#include "WolframImageLibrary.h"
#include "WolframLibrary.h"

DLLEXPORT int WolframLibrary_getVersion() {
    return WolframLibraryVersion;
}
DLLEXPORT int WolframLibrary_initialize(WolframLibraryData libData) {
    return LIBRARY_NO_ERROR;
}
DLLEXPORT void WolframLibrary_uninitialize(WolframLibraryData libData) {
    return;
}

// 这是一个正常的C函数 唯一要注意的是加static
static double f1(double z) {
    double sum = 0;
    for (int m = 1; m <= 20; m++) {
        sum += pow(9, -m) * pow(z * z * z, m) * tgamma(2 / 3.) / (tgamma(2 / 3. + m) * tgamma(1 + m));
    }
    return sum;
}

// f2 = LibraryFunctionLoad["f2.dll", "f2", {Real}, Real]
DLLEXPORT int f2(WolframLibraryData libData, mint Argc, MArgument *Args, MArgument Res) {
    MArgument_setReal(Res, f1(MArgument_getReal(Args[0])));
    return LIBRARY_NO_ERROR;
}

Makefile

MMA_DIR=D:\Program Files\Wolfram Research\Mathematica\12.0

CC=gcc
CFLAGS=-O3 -s -DNDEBUG

MMA_INC=$(MMA_DIR)\SystemFiles\IncludeFiles\C
MMA_DLL=$(MMA_DIR)\SystemFiles\Libraries\Windows-x86-64\WolframRTL.dll
MMA_STR=-I "$(MMA_INC)" "$(MMA_DLL)"

f2.dll: f2.c
	$(CC) f2.c $(CFLAGS) $(MMA_STR) -shared -o f2.dll

然后

f2 = LibraryFunctionLoad["f2.dll", "f2", {Real}, Real]

VS的话去 SystemFiles\Libraries\Windows-x86-64 那个目录里找 WolframRTL.lib 链接上就是了

顺便吐槽一句：Mathematica 的 ffi 慢的一匹，上面的库函数方法好多了，和调用自己的函数速度一样。

再吐槽一句：Mathematica 的 ffi 和 LibraryFunctionLoad 参数顺序特么是反的？？