test_pool2d_op.py 4.5 KB
Newer Older
C
chengduoZH 已提交
1 2 3 4 5
import unittest
import numpy as np
from op_test import OpTest


C
chengduoZH 已提交
6
def max_pool2D_forward_naive(x, ksize, strides, paddings=[0, 0], global_pool=0):
C
chengduoZH 已提交
7 8

    N, C, H, W = x.shape
C
chengduoZH 已提交
9 10
    if global_pool == 1:
        ksize = [H, W]
C
chengduoZH 已提交
11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
    H_out = (H - ksize[0] + 2 * paddings[0]) / strides[0] + 1
    W_out = (W - ksize[1] + 2 * paddings[1]) / strides[1] + 1
    out = np.zeros((N, C, H_out, W_out))
    for i in xrange(H_out):
        for j in xrange(W_out):
            r_start = np.max((i * strides[0] - paddings[0], 0))
            r_end = np.min((i * strides[0] + ksize[0] - paddings[0], H))
            c_start = np.max((j * strides[1] - paddings[1], 0))
            c_end = np.min((j * strides[1] + ksize[1] - paddings[1], W))
            x_masked = x[:, :, r_start:r_end, c_start:c_end]

            out[:, :, i, j] = np.max(x_masked, axis=(2, 3))
    return out


C
chengduoZH 已提交
26
def avg_pool2D_forward_naive(x, ksize, strides, paddings=[0, 0], global_pool=0):
C
chengduoZH 已提交
27 28

    N, C, H, W = x.shape
C
chengduoZH 已提交
29 30
    if global_pool == 1:
        ksize = [H, W]
C
chengduoZH 已提交
31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51
    H_out = (H - ksize[0] + 2 * paddings[0]) / strides[0] + 1
    W_out = (W - ksize[1] + 2 * paddings[1]) / strides[1] + 1
    out = np.zeros((N, C, H_out, W_out))
    for i in xrange(H_out):
        for j in xrange(W_out):
            r_start = np.max((i * strides[0] - paddings[0], 0))
            r_end = np.min((i * strides[0] + ksize[0] - paddings[0], H))
            c_start = np.max((j * strides[1] - paddings[1], 0))
            c_end = np.min((j * strides[1] + ksize[1] - paddings[1], W))
            x_masked = x[:, :, r_start:r_end, c_start:c_end]

            out[:, :, i, j] = np.sum(x_masked, axis=(2, 3)) / (
                (r_end - r_start) * (c_end - c_start))
    return out


class TestPool2d_Op(OpTest):
    def setUp(self):
        self.initTestCase()
        input = np.random.random(self.shape).astype("float32")
        output = self.pool2D_forward_naive(input, self.ksize, self.strides,
C
chengduoZH 已提交
52
                                           self.paddings, self.global_pool)
53
        self.inputs = {'X': input}
C
chengduoZH 已提交
54 55 56 57 58

        self.attrs = {
            'strides': self.strides,
            'paddings': self.paddings,
            'ksize': self.ksize,
59
            'poolingType': self.pool_type,
C
chengduoZH 已提交
60
            'globalPooling': self.global_pool,
C
chengduoZH 已提交
61 62
        }

63
        self.outputs = {'Out': output}
C
chengduoZH 已提交
64 65 66 67 68

    def test_check_output(self):
        self.check_output()

    def test_check_grad(self):
C
chengduoZH 已提交
69
        if self.pool_type != "max":
70
            self.check_grad(set(['X']), 'Out', max_relative_error=0.07)
C
chengduoZH 已提交
71 72

    def initTestCase(self):
C
chengduoZH 已提交
73 74
        self.global_pool = True
        self.op_type = "pool2d"
C
chengduoZH 已提交
75 76
        self.pool_type = "avg"
        self.pool2D_forward_naive = avg_pool2D_forward_naive
C
chengduoZH 已提交
77 78 79 80 81 82
        self.shape = [2, 3, 5, 5]
        self.ksize = [3, 3]
        self.strides = [1, 1]
        self.paddings = [0, 0]


C
chengduoZH 已提交
83 84
class TestCase1(TestPool2d_Op):
    def initTestCase(self):
C
chengduoZH 已提交
85
        self.global_pool = False
C
chengduoZH 已提交
86
        self.op_type = "pool2d"
C
chengduoZH 已提交
87 88
        self.pool_type = "avg"
        self.pool2D_forward_naive = avg_pool2D_forward_naive
C
chengduoZH 已提交
89
        self.shape = [2, 3, 7, 7]
C
chengduoZH 已提交
90 91
        self.ksize = [3, 3]
        self.strides = [1, 1]
C
chengduoZH 已提交
92
        self.paddings = [0, 0]
C
chengduoZH 已提交
93 94


C
chengduoZH 已提交
95 96
class TestCase2(TestPool2d_Op):
    def initTestCase(self):
C
chengduoZH 已提交
97
        self.global_pool = False
C
chengduoZH 已提交
98
        self.op_type = "pool2d"
C
chengduoZH 已提交
99 100
        self.pool_type = "avg"
        self.pool2D_forward_naive = avg_pool2D_forward_naive
C
chengduoZH 已提交
101
        self.shape = [2, 3, 7, 7]
C
chengduoZH 已提交
102 103 104 105 106
        self.ksize = [3, 3]
        self.strides = [1, 1]
        self.paddings = [1, 1]


C
chengduoZH 已提交
107 108
class TestCase3(TestPool2d_Op):
    def initTestCase(self):
C
chengduoZH 已提交
109
        self.global_pool = True
C
chengduoZH 已提交
110 111 112 113 114 115
        self.op_type = "pool2d"
        self.pool_type = "max"
        self.pool2D_forward_naive = max_pool2D_forward_naive
        self.shape = [2, 3, 5, 5]
        self.ksize = [3, 3]
        self.strides = [1, 1]
C
chengduoZH 已提交
116
        self.paddings = [0, 0]
C
chengduoZH 已提交
117 118


C
chengduoZH 已提交
119
class TestCase3(TestPool2d_Op):
C
chengduoZH 已提交
120
    def initTestCase(self):
C
chengduoZH 已提交
121
        self.global_pool = False
C
chengduoZH 已提交
122 123 124
        self.op_type = "pool2d"
        self.pool_type = "max"
        self.pool2D_forward_naive = max_pool2D_forward_naive
C
chengduoZH 已提交
125 126 127 128 129 130 131 132 133 134 135 136 137
        self.shape = [2, 3, 7, 7]
        self.ksize = [3, 3]
        self.strides = [1, 1]
        self.paddings = [0, 0]


class TestCase3(TestPool2d_Op):
    def initTestCase(self):
        self.global_pool = False
        self.op_type = "pool2d"
        self.pool_type = "max"
        self.pool2D_forward_naive = max_pool2D_forward_naive
        self.shape = [2, 3, 7, 7]
C
chengduoZH 已提交
138 139 140 141
        self.ksize = [3, 3]
        self.strides = [1, 1]
        self.paddings = [1, 1]

C
chengduoZH 已提交
142 143 144

if __name__ == '__main__':
    unittest.main()