FORTRANでの計算方法について

Question

FORTRANで行列の計算を行いたいのですが、以下の式はどういう風に計算していけばいいでしょうか？

Z＝［W｛A（B+C)｝］・X

ABCは300×300の正方行列、Wは1×300の列ベクトル、Xは300×１の行ベクトルで、・Xは内積です。

Comments

[cateye]

前出の質問：https://teratail.com/questions/296386
・・・ソースは、どうなりました？

[puro]

自己解決しました。ありがとうございます。

Answer 1

式の通りに命令を並べれば大丈夫です。
z = dot_product(matmul(w, matmul(a, b + c)), x)

配列の形状は
integer, parameter :: n = 300
real :: a(n, n), b(n, n), c(n, n), w(n), x(n), z
となります。

行列積も 300x300 位なら Fortran の組み込み命令で十分です。何万回も繰り返すというのであれば BLAS 系のハードウェア最適化ルーチンを呼べばより一層効率的になりますが。