外部プログラムをWolfram言語から呼び出す
外部プログラム中に定義してある関数をWolfram言語から呼び出すためには,関数に引数を渡したり,関数の結果を戻すためのWolfram Symbolic Transfer Protocol (WSTP)コードを外部プログラムに追加する必要がある.
:Begin:
:Function: f
:Pattern: f[x_Integer, y_Integer]
:Arguments: {x, y}
:ArgumentTypes: {Integer, Integer}
:ReturnType: Integer
:End:
:Begin: | 特定の関数に関するテンプレートを開始する |
:Function: | 外部プログラム中の関数名 |
:Pattern: | 関数を呼び出すように定義されたパターン |
:Arguments: | 関数に渡す引数 |
:ArgumentTypes: | 関数に渡す引数の型 |
:ReturnType: | 関数の戻り値の型 |
:End: | ひとつの関数に関するテンプレートを終了する |
:Evaluate: | 関数をインストールする際に評価されるWolfram言語入力 |
外部関数に対応するWSTPテンプレートは関数のソースコードと結合させる必要がある.ソースコードがC言語で書かれていると仮定すると,それは単にWSTPのヘッダファイルを含むコード1行と小さいmainプログラムをソースコードに追加するだけで済む.
#include "wstp.h"
int f(int x, int y) {
return x+y;
}
int main(int argc, char *argv[]) {
return WSMain(argc, argv);
}
WSTPテンプレートは通常 file.tmという名前のファイルにする.そのファイルにはさまざまな関数のテンプレートに混ざって,Cのソースコードも挿入することができる.
Unixを例に取ると,WSTPデベロッパキットにはwsccと呼ばれるプログラムがあり,このプログラムは名前が.tmで終るWSTPテンプレートファイルを前処理してCのソースコードを生成し,そのソースコードをccでコンパイルする.wsccはコマンドラインオプションや他のファイルは直接ccに渡す.
wscc -o f.exe f.tm f.c
Windowsでは,WSTPデベロッパーキットには,mprepプログラムが入っていて,これを直接呼び出して,前処理したい.tmファイルを処理する.mprepはCのソースコードを出力するので,それをCコンパイラにかける.
Install["prog"] | 外部プログラムをインストールする |
Uninstall[link] | 外部プログラムをアンインストールする |
Links["prog"] | "prog"に関連するアクティブなリンクを表示する |
Links[] | アクティブなリンクをすべて表示する |
LinkPatterns[link] | あるリンク上で評価されるパターンを表示する |
WSTPテンプレートが処理されるとき,基本的に2つのことが行われる.まず,:Pattern:との仕様に従い,WSTP経由で外部プログラムを呼び出すWolfram言語の定義を生成する.次に,:Function:,との仕様に従って外部プログラム中の関数を呼び出すCのソースコードが生成される.
:Begin:
:Function: prog_add
:Pattern: SkewAdd[x_Integer, y_Integer:1]
:Arguments: {x, If[x > 1, y, y + x - 2]}
:ArgumentTypes: {Integer, Integer}
:ReturnType: Integer
:End:
WSTPテンプレート中の:Pattern:とには,どんなWSTPの式を与えても構わない.に与えたものはすべて,その外部関数の呼出しごとに毎回評価される.評価の結果は引数のリストとして関数に渡される.
その場合,WSTPテンプレート中の:Evaluate:を利用する.:Evaluate:の後ろに与える式は何行になっても構わない.空行,あるいはスペース文字やタブ文字で始まらない行がきたら,式が終ったとみなされる.
:Evaluate: SkewAdd::usage = "SkewAdd[x, y] performs
a skew addition in an external program."
外部プログラムがインストールされるときは,WSTPテンプレートファイルに記述された順に仕様が読み込まれる.:Begin:より前の:Evaluate:に与えられた式は外部関数の定義が読み込まれる前に評価される.
:Evaluate: BeginPackage["XPack`"]
:Evaluate: XF1::usage = "XF1[x, y]は1つの外部関数です."
:Evaluate: XF2::usage = "XF2[x]は他の外部関数です."
:Evaluate: Begin["`Private`"]
:Begin:
:Function: f
:Pattern: XF1[x_Integer, y_Integer]
:Arguments: {x, y}
:ArgumentTypes: {Integer, Integer}
:ReturnType: Integer
:End:
:Begin:
:Function: g
:Pattern: XF2[x_?NumberQ]
:Arguments: {x}
:ArgumentTypes: {Real}
:ReturnType: Real
:End:
:Evaluate: End[ ]
:Evaluate: EndPackage[ ]
int f(int i, int j) {
return i + j;
}
double g(double x) {
return x*x;
}
:Evaluate:を利用して,外部プログラムが最初にインストールされるときに評価されるWolfram言語の式を指定することができる.また,外部プログラムのmain()中のWSMain()呼出しの前に挿入したコードも,外部プログラムのインストール前に実行される.外部関数が実行される前に,外部プログラムを初期化する場合等に使うとよいだろう.
WSEvaluateString(stdlink,"string") | 文字列をWolfram言語の入力として評価する |
int diff(int i, int j) {
if (i < j) WSEvaluateString(stdlink, "Print[\"negative\"]");
return i - j;
}
WSEvaluateString()の評価結果は,どんなものであろうと単に無視されることに注意すること.戻り値の利用にはWolfram言語と外部プログラムとの間で完全な双方向の通信が必要になる.これに関しては「外部プログラムとの双方向通信」で述べる.