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 ファイルに追加してください。

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

.ods、.xls、および .xlsx ファイルは、odsread または xlsread 関数を使用して開くことができます。 octave-io^AUR パッケージが必要です。

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

CSV形式に変換する

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

変換が完了したら、.csv ファイルに対して組み込み Octave 関数 csvread を使用できます。

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

トラブルシューティング

Zsh デコード不可能なトークン

エラーが発生した場合

undecodable token: b(hex)[23m

grml-zsh-config をインストールして再ログインしてください。

vi モード解読不可能なトークン

たとえば、.inputrc が vi モード用に設定されているユーザー

~/.inputrc

$include /etc/inputrc
set editing-mode vi
$if mode=vi

set show-mode-in-prompt on
set vi-ins-mode-string \1\e[6 q\2
set vi-cmd-mode-string \1\e[2 q\2

set keymap vi-command
# these are for vi-command mode
Control-l: clear-screen
Control-a: beginning-of-line

set keymap vi-insert
# these are for vi-insert mode
Control-l: clear-screen
Control-a: beginning-of-line

$endif

Octave GUI プロンプトに q>> undecodable token: \001b(hex)[6\0020(hex) として表示される場合、破損している可能性があります。 これを解決するには、上記の .inputrc を次のように変更して、Octave の show-mode-in-prompt 設定を無効にします。

~/.inputrc

$include /etc/inputrc
set editing-mode vi
$if mode=vi

$if Octave
set show-mode-in-prompt off
$else
set show-mode-in-prompt on
set vi-ins-mode-string \1\e[6 q\2
set vi-cmd-mode-string \1\e[2 q\2

set keymap vi-command
# these are for vi-command mode
Control-l: clear-screen
Control-a: beginning-of-line

set keymap vi-insert
# these are for vi-insert mode
Control-l: clear-screen
Control-a: beginning-of-line
$endif

$endif

参照

• Wikipedia:ja:GNU Octave