CythonでFortran NumPy操作をNOGILで使用すると、NumPy multiplyに相当するFortranライブラリーは何ですか？

Question

私は別のSOユーザーを助けようとしています。その過程で、NumPy以外の単純な処理を行うCythonプログラムを作成することはできません。これは私にGILを強制します。そのため、OpenMP（マルチコア）を使用することは不可能です。興味深いポストを見つけましたが、SciPyからFortranのライブラリを直接、Cythonコード（BLAS、LAPACK）にインポートすることができます。私の場合はインテルMKL同等の機能でNumPyがインストールされています。私がしようとしているのは、1000x1000次元の2つのベクトルの単純なベクトル乗算で、もう1つは転置され、1000x1000行列になります。しかし、私はそのトリックを行う関連するFortranルーチン（NumPyの乗算に相当）を見つけることができません。すべてのルーチンは、代わりに行列の乗算を行うように見えます。したがって、SciPyの優れた機能は、import scipy.linalg.cython_blas as blasと cimport scipy.linalg.cython_lapack as lapackとなりました。理論的には、dgemmからFortranライブラリblas.dgemm(options)を呼び出して起動しましたが、単純に要素ごとの乗算ではなく行列積を行います。誰も1000x1000行列の結果、2つの1000x1ベクトル、単純な乗算を行うFortanモジュールを知っていますか？あなたが入力構文を追加することができればそれは素晴らしいでしょう。 Cの連続したメモリビューを関数[ie1] Cython NumPyベクトルに渡します。

Answer 1

説明しているのは、"broadcasting"という純粋なNumPy機能です。これらのブロードキャスト操作は、C（またはCython）コードを使用して行われます。 PyNumber_MultiplyのようなPython C APIを使っていつでもCythonでアクセスすることができます（ただしGILはリリースされません）。通常Cythonの通常の乗算​​はその関数に委譲する必要がありますので、通常は呼び出す必要はありません。それは直接です。

BLAS/LAPACKは、線形代数の多くで使用されています。この目的のためにそこに公開されている関数を「使用する」ことができたとしても、NumPyの使用法は同じではありません。