from caffe2.python import core
import caffe2.python.hypothesis_test_util as hu
import caffe2.python.serialized_test.serialized_test_util as serial
import hypothesis.strategies as st
import itertools as it
import numpy as np
class TestMomentsOp(serial.SerializedTestCase):
def run_moments_test(self, X, axes, keepdims, gc, dc):
if axes is None:
op = core.CreateOperator(
"Moments",
["X"],
["mean", "variance"],
keepdims=keepdims,
)
else:
op = core.CreateOperator(
"Moments",
["X"],
["mean", "variance"],
axes=axes,
keepdims=keepdims,
)
def ref(X):
mean = np.mean(X, axis=None if axes is None else tuple(
axes), keepdims=keepdims)
variance = np.var(X, axis=None if axes is None else tuple(
axes), keepdims=keepdims)
return [mean, variance]
self.assertReferenceChecks(gc, op, [X], ref)
self.assertDeviceChecks(dc, op, [X], [0, 1])
self.assertGradientChecks(gc, op, [X], 0, [0, 1])
@serial.given(X=hu.tensor(dtype=np.float32), keepdims=st.booleans(),
num_axes=st.integers(1, 4), **hu.gcs)
def test_moments(self, X, keepdims, num_axes, gc, dc):
self.run_moments_test(X, None, keepdims, gc, dc)
num_dims = len(X.shape)
if num_dims < num_axes:
self.run_moments_test(X, range(num_dims), keepdims, gc, dc)
else:
for axes in it.combinations(range(num_dims), num_axes):
self.run_moments_test(X, axes, keepdims, gc, dc)