GNU Scientific Library

C言語で記述された科学技術計算関数ライブラリ

GNU Scientific Library (GSL) は、C言語で記述された科学技術計算関数のライブラリである。オープンソースであり、GNU General Public Licenseのもとで配布されている。

GNU Scientific Library
最新版
2.7 / 2021年6月1日 (3年前) (2021-06-01)
リポジトリ ウィキデータを編集
対応言語 C言語
種別 ライブラリ
ライセンス GNU General Public License
公式サイト www.gnu.org/software/gsl/
テンプレートを表示

このプロジェクトは1996年にロスアラモス国立研究所のDr. M. GalassiとDr. J. Theilerの着想に始まり、計算物理の専門家集団(Dr G. Jungman、Dr B. Gough、Dr J. Davies、R. Priedhorsky、Dr M. Booth、Dr F. Rossi、Dr D. Eddelbuettelら)を中心に作成された。

線形計算についてはBLASをサポートしており、CBLAS インターフェイスを実装している。

Linuxをはじめ、Unix系OSを中心にサポートしている。Microsoft Visual Studio用のバイナリもある。

参考論文のレファレンス、文献等も充実している (リファレンス・マニュアルには日本語訳がある)。リファレンス・マニュアルにはサンプル・コードも多数収録されている。PSPPPerl Data Language (en)、O2scl[1]などのフリーウェア・プロジェクトでも利用されている。

複素数型やベクトル/行列型などは ANSI C で規定されている構造体で実装されており、C++ のクラスではない。そのためたとえば、複素数オブジェクト同士の加算が + 演算子で行えるようになっている訳ではなく、加算のための関数 (この場合 gsl_complex_add) を、二つの複素数オブジェクトを引数として呼ばねばならない。なお複素数型について、gsl-2.7 以降で C11に対応するCコンパイラを用い、 ヘッダファイルを適切にインクルードすれば、gsl_complex 型を C11標準の複素数型と同様に扱えるようになった[2]

拡張倍精度以上の精度における計算は、変数の内部表現が言語仕様で標準化されておらず、さらに環境に依存して精度が大きく変化するために対応していない。

開発

編集

GSL の開発チームは、GSL が GNUプロジェクトのパッケージであることを明示しており、GSL のコードはすべて誰もが自由に、どんな用途にでも使えることを保証している。そのために、プロプライエタリなコードや (たとえば 書籍 Numerical Recipes ISBN 4874085601 などの) GNU の定義するフリーソフトウェアに該当しないコードとは対立した開発姿勢をとっている。

2011年現在、年に1〜2回のメンテナンスリリースによるバグ修正対応が基本になっている。一方でブロックやスライスといったデータ構造の有用性や C++ 対応の是非についての議論も ML 上で行われており、もし議論が収束して開発陣での合意が形成されれば将来のバージョンで反映される可能性があるが、具体的なスケジュールを考慮するような段階ではない。

C++ サポート

編集

GSL は C 言語ライブラリであるため、C++のクラスから利用できる。しかしメンバー関数へのポインタは、その型が関数へのポインタとは異なる[3]ため、利用できない。関数へのポインタは、静的に定義された関数 (C言語における一般的な関数定義によるもの) に対して利用する必要がある。

C++から GSL を利用するためのラッパーも複数あるが、いずれも不完全であり、どの開発もあまり活発ではないか、停止している。GSL のヘルプ・メイリング・リストの議論[4]では、ラッパーを使わなくても普通に関数を呼ぶのに支障はない、ベクトル用にラッパーを作って使っている、線形代数が目的ならEigen[5]がある、などの情報が寄せられている。

なお、Microsoft Visual Studio 2008用にまとめられたパッケージが公開されている[6]

他の言語とのバインディング

編集

FORTRAN から GSL の各関数を利用するためのインターフェイスが開発中であり、β版がリリースされている[7]。またラッパーを自作したい場合の情報も公開されている[8]

GNU Octave から GSL の特殊関数を利用できるパッケージがリリースされている[9]

GSL のホームページ[10] に紹介されているもののうち、現在も開発が続いている主な言語バインディングを以下に挙げる。

プログラム例

編集

ベッセル関数の値を計算する C プログラムの例を以下に示す[11]

#include <stdio.h>
#include <gsl/gsl_sf_bessel.h>

int main(void)
{
  double x = 5.0;
  double y = gsl_sf_bessel_J0(x);
  printf("J0(%g) = %.18e\n", x, y);
  return 0;
}

GNU Make を使って上のプログラムをコンパイルし、GSL とリンクしようとする場合、そのコマンドは Makefile ファイル中では以下のようになる。

gcc $(gsl-config --cflags) example.c $(gsl-config --libs)

上のコマンドで生成された実行ファイルを実行すると、以下のように出力する。計算値の精度は倍精度実数である。

J0(5) = -1.775967713143382920e-01

提供する機能

編集

関連項目

編集

脚注

編集
  1. ^ Object-oriented Scientific Computing Library — o2scl 0.926 documentation”. neutronstars.utk.edu. 2021年12月18日閲覧。
  2. ^ Complex Numbers — GSL 2.7 documentation”. www.gnu.org. 2021年12月18日閲覧。
  3. ^ pointer to member function (英語)
  4. ^ C++ラッパーの有無に関する質問と返答 (英語)
  5. ^ Eigen ホームページ (英語)
  6. ^ GSL Visual Studio 移植版
  7. ^ http://www.lrz-muenchen.de/services/software/mathematik/gsl/fortran/index.html FGSL: A Fortran interface to the GNU Scientific Library (英語)]FGSL 原著論文 (英語)
  8. ^ Fortranでgsl (日本語)
  9. ^ GNU Octave の `gsl' パッケージ
  10. ^ GSL のホームページ (英語)
  11. ^ オンライン・マニュアル中のサンプル・プログラム

関連書籍

編集
  • Mark Galassi, Jim Davies, James Theiler, Brian Gough, Gerard Jungman, Michael Booth, Fabrice Rossi (2001). Gnu Scientific Library Reference Manual, (English), Network Theory Ltd. ISBN 978-0-9541617-0-5 .

外部リンク

編集