module mkConvolution3x3#(parameter WIDTH = 256, HEIGHT = 256)(
    input bit [WIDTH-1:0][HEIGHT-1:0] inputFeatureMap,
    input bit [8:0] kernel,  // 3x3カーネル、フラット化された9ビットの値
    output bit [WIDTH-2:0][HEIGHT-2:0] outputFeatureMap // 出力サイズはストライドとパディングに依存
);
    for (Integer row = 1; row < HEIGHT-1; row = row + 1) {
        for (Integer col = 1; col < WIDTH-1; col = col + 1) {
            Integer sum = 0;
            for (Integer kr = -1; kr <= 1; kr = kr + 1) {
                for (Integer kc = -1; kc <= 1; kc = kc + 1) {
                    // 入力特徴マップとカーネルの要素積の和を計算
                    sum += inputFeatureMap[row+kr][col+kc] * kernel[(kr+1)*3 + (kc+1)];
                }
            }
            // 出力特徴マップに畳み込み結果を代入
            outputFeatureMap[row-1][col-1] = sum;
        }
    }
endmodule