LAPACKにおいて逆行列を求める関数である"dgetri"のラッパーを作った。 using System; using System.Collections.Generic; using System.Linq; using System.Text; using System.Security; using System.Runtime.InteropServices; namespace mkl { public se…

MKLで逆行列を求めるには、行列を一度LU分解しておく必要がある。 ダウンロードしてきたMKLのC#のラッパーには、LAPACKのLU分解パッケージ"dgetrf"がなかったので、記載されていたものを参考に作成してみた。 using System; using System.Collections.Generi…

http://www.mlab.ice.uec.ac.jp/~ej-sib/numerical/numerical_lapack1.html関数に代入する数値の単精度、倍精度の違いは次の通り。 s...単精度の実数 d...倍精度の実数 c...単精度の複素数 z...倍精度の実数

C#が久々なので、文法をだいぶ忘れておる…。こことかこことか参考にさせていただいた。 using System; using System.Collections.Generic; using System.IO; using System.Linq; using System.Text; namespace ioTest { class ioTest { static void Main(str…

Intel Math Kernel Library (MKL)は、再配布が無償の数値計算用ライブラリである。 お値段は\53,550と良心的。Visual Studio 2010 Professionalを入れたので、C#でバリバリ開発しようとしていたら、MKLのソースはC++で書かれていた。なので、C#からMKLのDLL…

pythonで主成分分析をしてみた。 #! usr/bin/python # -*- coding:utf-8 -*- import numpy as np import scipy as sp from scipy import linalg import csv # imput data fileName = './input/pcaTest.csv' inputCSV = csv.reader(open(fileName,"rb")) #dat…

Rだと、read.table関数で一発でデータをデータフレームに変更してくれる。 pythonで似たようなことができないかなーと思い書いてみた。 #! usr/bin/python # -*- coding:utf-8 -*- import numpy as np import scipy as sp import csv filename = "./input/in…