Octave

関連記事

• Matlab
• Sage-mathematics
• Mathematica

公式ウェブサイト より:

GNU Octave は主に数値解析を目的とした高水準インタプリタ言語です。線形ならびに非線形問題を数値的に解いたり、他の数値実験を行うことができます。データを可視化したり操作するためのグラフィック機能も備えています。通常、Octave は対話式のコマンドラインインターフェイスから使われますが、非対話式のプログラムを書くのにも使うことができます。Octave 言語は Matlab とよく似ているためプログラムを簡単に移植可能です。

インストール

octave パッケージをインストールしてください。

GUI インターフェイスは octave --gui で起動し、コマンドラインインターフェイスは octave-cli で起動します。

他のGUI インターフェイス

標準の GUI インターフェイスは octave パッケージに含まれています。それ以外の非公式 GUI フロントエンドも利用できます:

• Cantor — 数学的演算をいくつかのバックエンド (Scilab、Maxima、Octave など) の 1 つに委任するグラフィカルユーザーインターフェイス。

https://edu.kde.org/cantor/ || cantor

• JupyterLab — Octave を含む多くのプログラミングバックエンドをサポートするブラウザベースの対話型環境。

https://jupyterlab.readthedocs.io/en/stable/ || jupyterlab+jupyter-octave_kernel^AUR

パフォーマンス

Octave は、デフォルトで線形代数計算に blas パッケージを使用します。ただし、この実装では最新の CPU 命令を利用していません。パフォーマンスを高速化するために、openblas パッケージを blas のドロップイン置換としてインストールできます。Intel CPU の場合は intel-oneapi-mkl、NVIDIA GPU の場合は cuda など、利用可能なハードウェアに応じて他の BLAS 実装も使用できます。

この点を説明するために、次のコードを使用して、NxN 行列の乗算で実行される GFLOPS の推定値を取得できます:

N = 4096;
A = single(rand(N, N));
B = single(rand(N, N));
start = clock();
C = A * B;
elapsedTime = etime(clock(), start);
gFlops = 2 * N * N * N / (elapsedTime * 1e9)

Intel Core i7-9750H で次のコードを実行します。

octave ~/test_program.m

gFlops = 3.7222

openblas をインストールし、単一スレッドでプログラムを実行した後:

OMP_NUM_THREADS=1 octave ~/test_program.m

gFlops = 121.55

openblas に 9750H で利用可能な 12 スレッドをすべて使用させた後:

OMP_NUM_THREADS=12 octave ~/test_program.m

gFlops = 281.33

Octave-Forge

Octave は Matlab の Toolbox と同じようなパッケージのセットを Octave-Forge から提供しています。パッケージの完全なリストは こちら にあります。

パッケージは Octave から直接インストールすることも、AUR からインストールすることもできます。下を見て下さい。

Octave のインストーラーを使う

Octave のインストーラーを使うことでパッケージを管理できます。パッケージは ~/octave にインストールされます。-global オプションを使用した場合はシステムディレクトリにインストールされます。パッケージをインストールするには:

octave:1> pkg install -forge packagename

ノート: control など、Octave のパッケージの中には、コンパイルしてインストールするために Arch Linux のパッケージである gcc-fortran を必要とするものがあります。

パッケージをアンインストールするには:

octave:3> pkg uninstall packagename

Octave によって自動的にロードされるパッケージもありますが、ロードされない場合:

octave:4> pkg load packagename

または:

octave:5> pkg load all

パッケージがロードされているかどうかは pkg list を使って確認します。アスタリスクが付いているパッケージが既にロードされているパッケージです。

Octave の起動時に全てのパッケージがロードされるようにするには:

/usr/share/octave/site/m/startup/octaverc

## System-wide startup file for Octave.
##
## This file should contain any commands that should be executed each
## time Octave starts for every user at this site. 
 pkg load all

AUR を使う

いくつかのパッケージは AUR に存在します (パッケージを検索)。Octave-forge の新しいパッケージは Archlinux 用の Octave-forge ヘルパースクリプト を使うことで半自動的に Arch 用のパッケージを作成できます。

プロット

Qt がデフォルトのプロットバックエンドです:

>> available_graphics_toolkits
ans =
{
  [1,1] = fltk
  [1,2] = qt
}
>> graphics_toolkit
ans = qt

あるいは、FLTK または gnuplot をバックエンドに使うこともできます:

>> graphics_toolkit("gnuplot");

上記の変更を永続化させるには ~/.octaverc ファイルに追加してください。

Documentation タブのバグ

Octave バージョン 3.8.1 現在、公式 GUI の Documentation タブのリンクをクリックすると GUI がフリーズすることがあります。バグレポート #38305 によれば、Octave と Qt のシグナルハンドラが衝突するのが原因とされています。

解決法 としては次のコマンドで /usr/share/info/octave.info* にある Octave の Info ファイルを unzip することで解決します:

# gunzip /usr/share/info/octave.info*

もしくは Documentation タブは無視して別の Info リーダーを使って下さい。

Microsoft Excel スプレッドシートを読み込む

Octave で Microsoft Excel のファイルを読み込む方法は複数存在します。

オープンフォーマットに変換する

.xls ファイルを Octave で使う一番簡単な方法は LibreOffice の Calc (1024列まで制限) や Calligra Suite の Sheets (32768例まで制限) を使って .csv や .ods に変換することです。

変換したら Octave の csvread 関数を使って .csv ファイルを利用できます:

octave:1> csvread('myfile.csv');

.ods ファイルの場合、odsread 関数を含んでいる octave-io^AUR パッケージが必要です:

octave:1> odsread('myfile.ods');

.xlsx ファイルの場合 AUR の python-xlsx2csv^AUR パッケージを使います:

 xlsx2csv -t /path/to/save/location -x /path/to/myfile.xlsx 

Octave から xls ファイルを直接読み込む

XLS ファイルを使わざるをえず何らかの理由で CSV や ODS に変換できない場合、octave-io^AUR パッケージの xlsread 関数を使うことができます。

octave-io^AUR バージョン 1.2.5 から、'OCT' という名前のインターフェイスが追加されており、とくに依存パッケージを必要とせず .xlsx (.xls ではありません), .ods, .gnumeric の読み取りができます。ただし、Java ベースのインターフェイスは未だに存在します (特に .xls ファイルを読みこんだり、それらのファイル形式で書き出すときに必要)。