theanoのガンマ関数でヘッセ行列を得ることは可能ですか？

Question

b=T.fvector() 
theano.gradient.hessian(T.gammaln(b).sum(), [b]) 

は、NotImplementedErrorとなります。同等のコードは、例えば、 T.logの代わりにT.gammaln私はtheanoはちょうどgammalnの2番目の派生を知らないと思いますか？

Answer 1

gammalnのgradはpsiを使用します。しかし、私たちはpsiのgradを実装していませんでした。

クラスpsiのgrad（）メソッドは、ファイルtheano/scalar/basic_scipy.pyに実装できます。これにより、このケースが機能します。

ドキュメントのこのページでは、（Theanoに新しいオペアンプを追加して大学院生を実装する方法を説明する方法を説明します）方法：http://deeplearning.net/software/theano/extending/extending_theano.html

あなたはこれを実装するヘルプが必要な場合は、theano-devメーリングリストを使用すると、より適しています継続的な議論のためのstackoverflow。

これを行うと、Theanoに大きな貢献をします。