最近はとりあえずPythonで実装してみて、よさげな結果が出そうならCで実装するという流れができつつある。 が、今回の場合はPython依存の関数を使いまくっているので、Cで実装するときどうするか、ということが問題だなぁ。 いよいよC++か？

Numpyの使い方などを勉強しつつIsomapを実装してみましょう、としてみたけど失敗。 固有値、固有ベクトルを求める辺りに問題がありそうななさそうな。 - 違った。 問題は固有値、固有ベクトルを求めた後に、固有値の大きいほうからn番目を使わないといけない…

そんなときはidentity()関数を使えばよい。 >>>from numpy import * >>>identity(2) >>>array([[1.,0.], [0.,1.]]) 同様のものにzeros()やones()がある。