私はあなたの質問が正しい場合、基本的に4Dテンソルになるミニバッチの畳み込みレイヤーを行う方法を尋ねています。
簡単に言えば、入力を一括して処理し、それぞれに畳み込みを適用するとします。ループを使用してベクトル化せずにコード化するのはかなり簡単です。
ベクトル化の実装は、多くの場合、im2col techniqueに基づいており、基本的に4-D入力テンソルを巨大行列に変換し、行列乗算を実行します。ここではPythonでnumpy.lib.stride_tricks
を使用して往路の実装があります:それはnumpy
に実装されている線形代数ライブラリ、いくつかの非自明な機能を使用していますが、あなたのライブラリにないかもしれ
import numpy as np
def conv_forward(x, w, b, stride, pad):
N, C, H, W = x.shape
F, _, HH, WW = w.shape
# Check dimensions
assert (W + 2 * pad - WW) % stride == 0, 'width does not work'
assert (H + 2 * pad - HH) % stride == 0, 'height does not work'
# Pad the input
p = pad
x_padded = np.pad(x, ((0, 0), (0, 0), (p, p), (p, p)), mode='constant')
# Figure out output dimensions
H += 2 * pad
W += 2 * pad
out_h = (H - HH)/stride + 1
out_w = (W - WW)/stride + 1
# Perform an im2col operation by picking clever strides
shape = (C, HH, WW, N, out_h, out_w)
strides = (H * W, W, 1, C * H * W, stride * W, stride)
strides = x.itemsize * np.array(strides)
x_stride = np.lib.stride_tricks.as_strided(x_padded,
shape=shape, strides=strides)
x_cols = np.ascontiguousarray(x_stride)
x_cols.shape = (C * HH * WW, N * out_h * out_w)
# Now all our convolutions are a big matrix multiply
res = w.reshape(F, -1).dot(x_cols) + b.reshape(-1, 1)
# Reshape the output
res.shape = (F, N, out_h, out_w)
out = res.transpose(1, 0, 2, 3)
out = np.ascontiguousarray(out)
return out
注こと。
ところで、一般的に、データセット全体を1つのバッチとしてプッシュしたくない場合は、複数のバッチに分割します。
これにはライブラリを使用していますか?一般的なチャンネルでは、CNNが非常に扱いやすいものです。 –
@JonasAdler私は線形代数を処理するためにC++にarmadilloライブラリを使用しています。具体的には、私は2Dと3Dの行列それぞれに対して行列とキューブのデータ構造を使用します。 –