Microsoft R Openは、デフォルトでMKLを使用する。使用スレッド数を1に制限。 MinicondaでPythonをインストールするとnumpy+MKLを使用するため速くなる。 最速はC++/Eigen/MKL。 Julia+MKLは未評価。MKLを使用するためにはソースからコンパイルする必要あり。 処理の大部分は行列(767x767)のコレスキー分解であり、MKLの効果が大きい。 サンプルデータ USPS handwritten digit data http://www.gaussianprocess.org/gpml/data/ ## データの入力 library("R.matlab") mat_data <- readMat('usps_resampled.mat') index3 = (mat_data$train.labels[3+1,] == 1) index5