commons-math3 解読 (0) ：体 Field, FieldElement ： math3 パッケージ

以前、「倭算数理研究所」で書いていた commons-math 解読シリーズの記事ですが明らかに書くブログを間違えたので、commons-math3 が出たこの期にこちらへ移転します。

今回は org.apache.commons.math3 パッケージ内に定義されている Field, FieldElement インターフェースについて見ていきます。

Field とは日本語で「体（タイ）」と呼ばれるもので、大雑把に言えば加減乗除の演算が定義された集合です。　具体的には有理数（分数）全体の集合、実数全体の集合など（通常の加減乗除の演算と合わせて）です。　整数は割り算ができないので体ではありません。　数学での定義は後に書いたので、興味のある方はそちらをどうぞ。　1つの体に関して Field オブジェクトは1つだけしか存在しない、つまり Field インターフェースを実装している具象クラスはシングルトンであることが望ましいようです。

FieldElement は上記の Field に含まれている要素です。　1つの有理数、1つの実数それぞれが FieldElement オブジェクトになります。

commons-math には代数系として「体」しか定義されていませんが、他にも「群 (Group)」・「環 (Ring)」などもあります。　それらを使いたい方は JScience の org.jscience.mathematics.structure パッケージに定義されているクラス群を使ってみるといいかもしれません。　ただし、そこに定義されている Field は、commons-math での FieldElement に対応するものなので注意が必要です。

では、以下で Field, FieldElement の定義と実装クラスを見ていきましょう。

Field<T> インターフェース

まずは Field<T> インターフェース。　型パラメータ T はこの Field の要素の型です：

package org.apache.commons.math3;

public interface Field<T>{
    T getZero();
    T getOne();
    Class<? extends FieldElement<T>> getRuntimeClass();
}

getZero(), getOne() メソッドはそれぞれ和演算 (+) と積演算 (*) の単位元を取得するメソッドです。　返り値はどちらも T になってますが、これらは後で見る FieldElement オブジェクトでないといけないはずなので、本当なら

public interface Field<T extends FieldElement<T>>{ ... }

となっているべきなんですが・・・　おそらく実装者が Generics に不慣れだったんでしょう。　もう変更不可能。　まぁ、そんなに困らないかと思いますが。　また、getRuntimeClass() メソッドはこの Field オブジェクトの要素となる FieldElement クラスの Class オブジェクトを返します。　まぁ、BigRealField オブジェクトに対して BigReal クラス、ComplexField オブジェクトに対して Complex クラスなどです。

commons-math 内に定義されている Field インターフェースの実装クラスには、以下のものがあります：

FractionField / BigFractionField (org.apache.commons.math3.fraction パッケージ)
BigRealField (org.apache.commons.math3.util パッケージ)
DfpField (org.apache.commons.math3.dfp パッケージ)
ComplexField (org.apache.commons.math3.complex パッケージ)

各クラスの詳細については後日に。

FieldElement<T> インターフェース

次は FieldElement<T> インターフェース。　実際にプログラミングで使うのはこのインターフェースを実装したクラスのオブジェクトです。

package org.apache.commons.math3;

public interface FieldElement<T>{

    // 体
    Field<T> getField();

    // 逆元
    T negate();
    T reciprocal();

    // 加減乗除
    T add(T a);
    T subtract(T a);
    T multiply(T a);
    T multiply(int n);
    T divide(T a);
}

この要素が属する「体」に対応する Field オブジェクトを取得するメソッド、和・積に関する逆元を返すメソッド、そして加減乗除の演算に対応するメソッドが定義されています。　この FieldElement インターフェースも、本当なら

public FieldElement<T extends FieldElement<T>>{ ... }

のように定義されるべきでしょうね。　まぁそれはともかく、commons-math3 内に定義されている FieldElement の実装クラスには

Fraction / BigFraction
BigReal
Dfp / DfpDec
Complex

などがあります。　パッケージは上記のそれぞれの体と同じです。　これらも後日に。

ちょっと数学を

代数学での「体」の定義は、以下のようになっています（x, y, z はこの体の要素とします）：

和演算 (+) に関して

閉じている

結合法則が成り立つ： (x+y)+z = x+(y+z)

交換法則が成り立つ： x+y = y+x

単位元 (0) が存在する： x+0 = 0+x = x

各元 x に対して逆元 -x が存在する： x+(-x) = (-x)+x = 0

積演算 (*) に関して

閉じている

結合法則が成り立つ： (x*y)*z = x*(y*z)

単位元 (1) が存在する： x*1 = 1*x = x

和演算の単位元 (0) 以外の各元 x に対して逆元 x^-1 が存在する： x*x^-1 = x^-1*x = 1

和演算・積演算に関して

分配法則が成り立つ： x*(y+z) = x*y + x*z, (y+z)*x = y*x + z*x

これに加えて積演算に関する交換法則 x*y = y*x が成り立つものを「可換体 (commutative field)」といいます。　commons-math で定義されている「体」は全て可換体です。

では、commons-math の Field, FieldElement が上記の体の定義を満たしているのかを見ていきましょう：

和演算 (FieldElement#add()) に関して
1. 閉じているか？・・・型パラメータ T にFieldElement のサブタイプを指定すれば OK
2. 結合法則・・・実装者への制約
3. 交換法則・・・実装者への制約
4. 単位元の存在・・・Field#getZero() で取得可
5. 逆元の存在・・・FieldElement#negate() で取得可
積演算 (FieldElement#multiply()) に関して
1. 閉じているか？・・・型パラメータ T にFieldElement のサブタイプを指定すれば OK
2. 結合法則・・・実装者への制約
3. 単位元の存在・・・Field#getOne() で取得可
4. 逆元の存在・・・FieldElement#reciprocal() で取得可
和演算・積演算に関して
1. 分配法則・・・実装者への制約