FM / operator

/usr/local/lib64/python3.6/site-packages/caffe2/python/operator_test/

/usr/local/lib64/python3.6/site-packages/caffe2/python/operator_test

Name	Size	Mode	Actions
__pycache__/	-	0755	rm
activation_ops_test.py	9691	0644	edit dl rm
adadelta_test.py	7932	0644	edit dl rm
adagrad_test.py	7586	0644	edit dl rm
adagrad_test_helper.py	5181	0644	edit dl rm
adam_test.py	21559	0644	edit dl rm
affine_channel_op_test.py	3784	0644	edit dl rm
alias_with_name_test.py	938	0644	edit dl rm
apmeter_test.py	2738	0644	edit dl rm
arg_ops_test.py	1917	0644	edit dl rm
assert_test.py	797	0644	edit dl rm
async_net_barrier_test.py	946	0644	edit dl rm
atomic_ops_test.py	4104	0644	edit dl rm
basic_rnn_test.py	4720	0644	edit dl rm
batch_box_cox_test.py	5080	0644	edit dl rm
batch_bucketize_op_test.py	3730	0644	edit dl rm
batch_moments_op_test.py	2795	0644	edit dl rm
batch_sparse_to_dense_op_test.py	4199	0644	edit dl rm
bbox_transform_test.py	12258	0644	edit dl rm
bisect_percentile_op_test.py	6227	0644	edit dl rm
blobs_queue_db_test.py	3240	0644	edit dl rm
boolean_mask_test.py	16389	0644	edit dl rm
boolean_unmask_test.py	1711	0644	edit dl rm
box_with_nms_limit_op_test.py	8750	0644	edit dl rm
bucketize_op_test.py	930	0644	edit dl rm
cast_op_test.py	1600	0644	edit dl rm
ceil_op_test.py	888	0644	edit dl rm
channel_backprop_stats_op_test.py	2131	0644	edit dl rm
channel_shuffle_test.py	1794	0644	edit dl rm
channel_stats_op_test.py	2639	0644	edit dl rm
checkpoint_test.py	1500	0644	edit dl rm
clip_op_test.py	1984	0644	edit dl rm
clip_tensor_op_test.py	2076	0644	edit dl rm
collect_and_distribute_fpn_rpn_proposals_op_test.py	11269	0644	edit dl rm
concat_op_cost_test.py	2858	0644	edit dl rm
concat_split_op_test.py	7266	0644	edit dl rm
conditional_test.py	995	0644	edit dl rm
conftest.py	1446	0644	edit dl rm
conv_test.py	32473	0644	edit dl rm
conv_transpose_test.py	15945	0644	edit dl rm
copy_ops_test.py	7374	0644	edit dl rm
copy_rows_to_tensor_op_test.py	2526	0644	edit dl rm
cosine_embedding_criterion_op_test.py	1953	0644	edit dl rm
counter_ops_test.py	3348	0644	edit dl rm
crf_test.py	5315	0644	edit dl rm
cross_entropy_ops_test.py	10085	0644	edit dl rm
ctc_beam_search_decoder_op_test.py	5197	0644	edit dl rm
ctc_greedy_decoder_op_test.py	4743	0644	edit dl rm
cudnn_recurrent_test.py	5817	0644	edit dl rm
dataset_ops_test.py	23847	0644	edit dl rm
data_couple_op_test.py	858	0644	edit dl rm
decay_adagrad_test.py	2694	0644	edit dl rm
deform_conv_test.py	19276	0644	edit dl rm
dense_vector_to_id_list_op_test.py	2044	0644	edit dl rm
depthwise_3x3_conv_test.py	1863	0644	edit dl rm
detectron_keypoints.py	7973	0644	edit dl rm
distance_op_test.py	4351	0644	edit dl rm
dropout_op_test.py	2971	0644	edit dl rm
duplicate_operands_test.py	734	0644	edit dl rm
elementwise_linear_op_test.py	1382	0644	edit dl rm
elementwise_logical_ops_test.py	4617	0644	edit dl rm
elementwise_ops_test.py	33354	0644	edit dl rm
elementwise_op_broadcast_test.py	17466	0644	edit dl rm
emptysample_ops_test.py	1977	0644	edit dl rm
enforce_finite_op_test.py	1286	0644	edit dl rm
ensure_clipped_test.py	1505	0644	edit dl rm
ensure_cpu_output_op_test.py	1244	0644	edit dl rm
erf_op_test.py	749	0644	edit dl rm
expand_op_test.py	2109	0644	edit dl rm
fc_operator_test.py	3720	0644	edit dl rm
feature_maps_ops_test.py	21492	0644	edit dl rm
filler_ops_test.py	8476	0644	edit dl rm
find_op_test.py	1316	0644	edit dl rm
flatten_op_test.py	922	0644	edit dl rm
flexible_top_k_test.py	2609	0644	edit dl rm
floor_op_test.py	894	0644	edit dl rm
fused_nbit_rowwise_conversion_ops_test.py	14077	0644	edit dl rm
fused_nbit_rowwise_test_helper.py	2693	0644	edit dl rm
gather_ops_test.py	9216	0644	edit dl rm
gather_ranges_op_test.py	9125	0644	edit dl rm
given_tensor_byte_string_to_uint8_fill_op_test.py	1392	0644	edit dl rm
given_tensor_fill_op_test.py	1503	0644	edit dl rm
glu_op_test.py	1212	0644	edit dl rm
group_conv_test.py	2870	0644	edit dl rm
group_norm_op_test.py	5252	0644	edit dl rm
gru_test.py	12932	0644	edit dl rm
heatmap_max_keypoint_op_test.py	4770	0644	edit dl rm
histogram_test.py	3097	0644	edit dl rm
hsm_test.py	9456	0644	edit dl rm
hyperbolic_ops_test.py	1472	0644	edit dl rm
im2col_col2im_test.py	4311	0644	edit dl rm
image_input_op_test.py	17345	0644	edit dl rm
index_hash_ops_test.py	2885	0644	edit dl rm
index_ops_test.py	4597	0644	edit dl rm
instance_norm_test.py	9917	0644	edit dl rm
integral_image_ops_test.py	3419	0644	edit dl rm
jsd_ops_test.py	1044	0644	edit dl rm
key_split_ops_test.py	1289	0644	edit dl rm
lars_test.py	1354	0644	edit dl rm
layer_norm_op_test.py	14983	0644	edit dl rm
leaky_relu_test.py	5639	0644	edit dl rm
learning_rate_adaption_op_test.py	2837	0644	edit dl rm
learning_rate_op_test.py	8652	0644	edit dl rm
lengths_pad_op_test.py	1625	0644	edit dl rm
lengths_reducer_fused_nbit_rowwise_ops_test.py	15495	0644	edit dl rm
lengths_tile_op_test.py	1332	0644	edit dl rm
lengths_top_k_ops_test.py	2371	0644	edit dl rm
length_split_op_test.py	4868	0644	edit dl rm
listwise_l2r_operator_test.py	8740	0644	edit dl rm
load_save_test.py	33241	0644	edit dl rm
locally_connected_op_test.py	7761	0644	edit dl rm
loss_ops_test.py	902	0644	edit dl rm
lpnorm_op_test.py	2725	0644	edit dl rm
map_ops_test.py	2249	0644	edit dl rm
margin_ranking_criterion_op_test.py	1816	0644	edit dl rm
math_ops_test.py	1603	0644	edit dl rm
matmul_op_test.py	10096	0644	edit dl rm
mean_op_test.py	1469	0644	edit dl rm
merge_id_lists_op_test.py	2989	0644	edit dl rm
mkl_conv_op_test.py	1547	0644	edit dl rm
mkl_packed_fc_op_test.py	2647	0644	edit dl rm
mod_op_test.py	1459	0644	edit dl rm
moments_op_test.py	1722	0644	edit dl rm
momentum_sgd_test.py	6480	0644	edit dl rm
mpi_test.py	8154	0644	edit dl rm
mul_gradient_benchmark.py	1509	0644	edit dl rm
negate_gradient_op_test.py	1518	0644	edit dl rm
ngram_ops_test.py	2327	0644	edit dl rm
normalize_op_test.py	1679	0644	edit dl rm
numpy_tile_op_test.py	1924	0644	edit dl rm
one_hot_ops_test.py	7478	0644	edit dl rm
onnx_while_test.py	3070	0644	edit dl rm
order_switch_test.py	1306	0644	edit dl rm
pack_ops_test.py	12634	0644	edit dl rm
pack_rnn_sequence_op_test.py	2891	0644	edit dl rm
pad_test.py	1377	0644	edit dl rm
partition_ops_test.py	6838	0644	edit dl rm
percentile_op_test.py	4427	0644	edit dl rm
piecewise_linear_transform_test.py	6187	0644	edit dl rm
pooling_test.py	16508	0644	edit dl rm
prepend_dim_test.py	1505	0644	edit dl rm
python_op_test.py	1312	0644	edit dl rm
quantile_test.py	3276	0644	edit dl rm
rand_quantization_op_speed_test.py	3128	0644	edit dl rm
rank_loss_operator_test.py	5752	0644	edit dl rm
rebatching_queue_test.py	9047	0644	edit dl rm
record_queue_test.py	3125	0644	edit dl rm
recurrent_network_test.py	14048	0644	edit dl rm
recurrent_net_executor_test.py	10922	0644	edit dl rm
reduce_ops_test.py	17341	0644	edit dl rm
reduction_ops_test.py	4664	0644	edit dl rm
reshape_ops_test.py	8211	0644	edit dl rm
resize_op_test.py	9417	0644	edit dl rm
rmac_regions_op_test.py	3178	0644	edit dl rm
rms_norm_op_test.py	1325	0644	edit dl rm
rnn_cell_test.py	59707	0644	edit dl rm
roi_align_rotated_op_test.py	7567	0644	edit dl rm
rowwise_counter_test.py	2205	0644	edit dl rm
scale_op_test.py	2177	0644	edit dl rm
segment_ops_test.py	25745	0644	edit dl rm
self_binning_histogram_test.py	12915	0644	edit dl rm
selu_op_test.py	3232	0644	edit dl rm
sequence_ops_test.py	16000	0644	edit dl rm
shape_inference_test.py	25708	0644	edit dl rm
sinusoid_position_encoding_op_test.py	2308	0644	edit dl rm
softmax_ops_test.py	23685	0644	edit dl rm
softplus_op_test.py	516	0644	edit dl rm
sparse_dropout_with_replacement_op_test.py	2885	0644	edit dl rm
sparse_gradient_checker_test.py	1294	0644	edit dl rm
sparse_itemwise_dropout_with_replacement_op_test.py	2913	0644	edit dl rm
sparse_lengths_sum_benchmark.py	4159	0644	edit dl rm
sparse_lp_regularizer_test.py	2553	0644	edit dl rm
sparse_normalize_test.py	3136	0644	edit dl rm
sparse_ops_test.py	3469	0644	edit dl rm
sparse_to_dense_mask_op_test.py	3693	0644	edit dl rm
spatial_bn_op_test.py	20182	0644	edit dl rm
specialized_segment_ops_test.py	11775	0644	edit dl rm
split_op_cost_test.py	8645	0644	edit dl rm
square_root_divide_op_test.py	2179	0644	edit dl rm
stats_ops_test.py	1789	0644	edit dl rm
stats_put_ops_test.py	6596	0644	edit dl rm
storm_test.py	6507	0644	edit dl rm
string_ops_test.py	4154	0644	edit dl rm
text_file_reader_test.py	2517	0644	edit dl rm
thresholded_relu_op_test.py	2323	0644	edit dl rm
tile_op_test.py	3887	0644	edit dl rm
top_k_test.py	9113	0644	edit dl rm
torch_integration_test.py	39941	0644	edit dl rm
transpose_op_test.py	2722	0644	edit dl rm
trigonometric_op_test.py	1715	0644	edit dl rm
unique_ops_test.py	2255	0644	edit dl rm
unique_uniform_fill_op_test.py	1335	0644	edit dl rm
unsafe_coalesce_test.py	2940	0644	edit dl rm
upsample_op_test.py	7308	0644	edit dl rm
utility_ops_test.py	15054	0644	edit dl rm
video_input_op_test.py	10503	0644	edit dl rm
weighted_multi_sample_test.py	1997	0644	edit dl rm
weighted_sample_test.py	2739	0644	edit dl rm
weighted_sum_test.py	3052	0644	edit dl rm
weight_scale_test.py	2057	0644	edit dl rm
wngrad_test.py	8279	0644	edit dl rm
__init__.py	0	0644	edit dl rm

Edit: /usr/local/lib64/python3.6/site-packages/caffe2/python/operator_test/segment_ops_test.py (25745B)

from functools import partial
from hypothesis import given, settings

import numpy as np
import unittest
import hypothesis.strategies as st

from caffe2.python import core, workspace
import caffe2.python.hypothesis_test_util as hu
import caffe2.python.serialized_test.serialized_test_util as serial

def sparse_lengths_sum_ref(D, I, L, normalize_by_lengths=False):
    R = np.zeros(shape=(L.size,) + D.shape[1:], dtype=np.float32)
    line = 0
    for g in range(L.size):
        for _ in range(L[g]):
            if len(D.shape) > 1:
                R[g, :] += D[I[line], :]
            else:
                R[g] += D[I[line]]
            line += 1

if normalize_by_lengths and L[g] > 1:
            if len(D.shape) > 1:
                R[g, :] = R[g, :] / L[g]
            else:
                R[g] = R[g] / L[g]

return [R]

def sparse_lengths_mean_ref(D, I, L):
    return sparse_lengths_sum_ref(D, I, L, normalize_by_lengths=True)

class TesterBase:
    def segment_reduce_op(self, data, segment_ids, reducer, indices=None):
        segments = self.split(data, segment_ids, indices)
        output = np.zeros((len(segments), ) + data.shape[1:])
        for i, segment in enumerate(segments):
            if len(segment) > 0:
                output[i] = reducer(segment)
            else:
                output[i] = 0.0
        return output

def segment_reduce_grad_op(
        self,
        data,
        segment_ids,
        reducer_grad,
        grad_out,
        output,
        indices=None
    ):
        segments = self.split(data, segment_ids, indices)
        segment_grads = [
            reducer_grad(grad_out[i], [output[i]], [segment])
            for i, segment in enumerate(segments)
        ]
        return self.unsplit(data.shape[1:], segment_grads, segment_ids)

def _test(self, prefix, input_strategy, refs, gpu=False, **kwargs):
        tester = self
        operator_args = kwargs.pop('operator_args', {})
        threshold = kwargs.pop('threshold', 1e-4)
        grad_check = kwargs.pop('grad_check', True)

@given(X=input_strategy, **hu.gcs)
        def test_segment_ops(self, X, gc, dc):
            if not gpu and gc.device_type > 0:
                return
            for op_name, ref, grad_ref in refs:
                inputs = ['input%d' % i for i in range(0, len(X))]
                op = core.CreateOperator(
                    prefix + op_name, inputs, ['output'], **operator_args
                )
                print('Operator %s, ' % op.type, gc.device_type)

def seg_reduce(data, *args):
                    indices, segments = (
                        args if len(args) == 2 else (None, args[0])
                    )
                    out = tester.segment_reduce_op(
                        data=data,
                        segment_ids=segments,
                        indices=indices,
                        reducer=ref
                    )
                    return (out, )

def seg_reduce_grad(grad_out, outputs, inputs):
                    data = inputs[0]
                    args = inputs[1:]
                    indices, segments = (
                        args if len(args) == 2 else (None, args[0])
                    )
                    # grad r.t. data
                    grad_val = tester.segment_reduce_grad_op(
                        data, segments, grad_ref, grad_out, outputs[0], indices
                    )
                    # if sparse, include indices along with data gradient
                    data_grad_slice = (
                        (grad_val, indices) if indices is not None else grad_val
                    )
                    # other inputs don't have gradient
                    return (data_grad_slice, ) + (None, ) * (len(inputs) - 1)

kwargs = {}
                if grad_check:
                    kwargs['output_to_grad'] = 'output'
                    kwargs['grad_reference'] = seg_reduce_grad
                self.assertReferenceChecks(
                    device_option=gc,
                    op=op,
                    inputs=X,
                    reference=seg_reduce,
                    threshold=threshold,
                    **kwargs
                )
        return test_segment_ops

class SegmentsTester(TesterBase):
    def split(self, data, segment_ids, indices=None):
        """
        Given:
          data[M1 x M2 x ... x Md]
                          the input data
          indices[N]      the index of each entry of segment_ids into data,
                          where 0 <= index[i] < M1,
                          with default indices=[0,1,...N]
          segment_ids[N]  the segment_id for each entry of indices,

returns K outputs, each one containing data entries corresponding
        to one of the segments present in `segment_ids`.
        """
        if segment_ids.size == 0:
            return []
        K = max(segment_ids) + 1
        outputs = [
            np.zeros(
                (np.count_nonzero(segment_ids == seg_id), ) + data.shape[1:],
                dtype=data.dtype
            ) for seg_id in range(0, K)
        ]
        counts = np.zeros(K, dtype=int)
        for i, seg_id in enumerate(segment_ids):
            data_idx = i if indices is None else indices[i]
            outputs[seg_id][counts[seg_id]] = data[data_idx]
            counts[seg_id] += 1
        return outputs

def unsplit(self, extra_shape, inputs, segment_ids):
        """ Inverse operation to `split`, with indices=None """
        output = np.zeros((len(segment_ids), ) + extra_shape)
        if len(segment_ids) == 0:
            return output
        K = max(segment_ids) + 1
        counts = np.zeros(K, dtype=int)
        for i, seg_id in enumerate(segment_ids):
            output[i] = inputs[seg_id][counts[seg_id]]
            counts[seg_id] += 1
        return output

class LengthsTester(TesterBase):
    def split(self, data, lengths, indices=None):
        K = len(lengths)
        outputs = [
            np.zeros((lengths[seg_id], ) + data.shape[1:],
                     dtype=data.dtype) for seg_id in range(0, K)
        ]
        start = 0
        for i in range(0, K):
            for j in range(0, lengths[i]):
                data_index = start + j
                if indices is not None:
                    data_index = indices[data_index]
                outputs[i][j] = data[data_index]
            start += lengths[i]
        return outputs

def unsplit(self, extra_shape, inputs, lengths):
        N = sum(lengths)
        output = np.zeros((N, ) + extra_shape)
        K = len(lengths)
        assert len(inputs) == K
        current = 0
        for i in range(0, K):
            for j in range(0, lengths[i]):
                output[current] = inputs[i][j]
                current += 1
        return output

def sum_grad(grad_out, outputs, inputs):
    return np.repeat(
        np.expand_dims(grad_out, axis=0),
        inputs[0].shape[0],
        axis=0
    )

def logsumexp(x):
    return np.log(np.sum(np.exp(x), axis=0))

def logsumexp_grad(grad_out, outputs, inputs):
    sum_exps = np.sum(np.exp(inputs[0]), axis=0)
    return np.repeat(
        np.expand_dims(grad_out / sum_exps, 0),
        inputs[0].shape[0],
        axis=0
    ) * np.exp(inputs[0])

def logmeanexp(x):
    return np.log(np.mean(np.exp(x), axis=0))

def mean(x):
    return np.mean(x, axis=0)

def mean_grad(grad_out, outputs, inputs):
    return np.repeat(
        np.expand_dims(grad_out / inputs[0].shape[0], 0),
        inputs[0].shape[0],
        axis=0
    )

def max_fwd(x):
    return np.amax(x, axis=0)

def max_grad(grad_out, outputs, inputs):
    flat_inputs = inputs[0].flatten()
    flat_outputs = np.array(outputs[0]).flatten()
    flat_grad_in = np.zeros(flat_inputs.shape)
    flat_grad_out = np.array(grad_out).flatten()
    blocks = inputs[0].shape[0]
    if blocks == 0:
        return np.zeros(inputs[0].shape)
    block_size = flat_inputs.shape[0] // blocks

for i in range(block_size):
        out_grad = flat_grad_out[i]
        out = flat_outputs[i]
        for j in range(blocks):
            idx = j * block_size + i
            # we can produce multiple outputs for max
            if out == flat_inputs[idx]:
                flat_grad_in[idx] = out_grad

return np.resize(flat_grad_in, inputs[0].shape)

REFERENCES_ALL = [
    ('Sum', partial(np.sum, axis=0), sum_grad),
    ('Mean', partial(np.mean, axis=0), mean_grad),
]

REFERENCES_SORTED = [
    ('RangeSum', partial(np.sum, axis=0), sum_grad),
    ('RangeLogSumExp', logsumexp, logsumexp_grad),
    # gradient is the same as sum
    ('RangeLogMeanExp', logmeanexp, logsumexp_grad),
    ('RangeMean', mean, mean_grad),
    ('RangeMax', max_fwd, max_grad),
]

REFERENCES_LENGTHS_ONLY = [
    ('Max', partial(np.amax, axis=0), max_grad),
]

def sparse_lengths_weighted_sum_ref(D, W, I, L):
    R = np.zeros(shape=(len(L), ) + D.shape[1:], dtype=D.dtype)
    line = 0
    for g in range(len(L)):
        for _ in range(L[g]):
            if len(D.shape) > 1:
                R[g, :] += W[line] * D[I[line], :]
            else:
                R[g] += W[line] * D[I[line]]
            line += 1
    return [R]

def sparse_lengths_weighted_sum_grad_ref(
        GO, fwd_out, fwd_in, grad_on_weights=False):
    D, W, I, L = fwd_in
    GI = np.zeros(shape=(len(I), ) + D.shape[1:], dtype=D.dtype)
    GW = np.zeros(shape=W.shape, dtype=W.dtype) if grad_on_weights else None
    line = 0
    for g in range(len(L)):
        for _ in range(L[g]):
            if len(GO.shape) > 1:
                GI[line, :] = W[line] * GO[g, :]
            else:
                GI[line] = W[line] * GO[g]
            if GW is not None:
                if len(GO.shape) > 1:
                    GW[line] = np.dot(GO[g].flatten(), D[I[line], :].flatten())
                else:
                    GW[line] = np.dot(GO[g].flatten(), D[I[line]].flatten())
            line += 1
    print(GW)
    return [(GI, I), GW, None, None]

class TestSegmentOps(hu.HypothesisTestCase):
    def test_sorted_segment_ops(self):
        SegmentsTester()._test(
            'SortedSegment',
            hu.segmented_tensor(
                dtype=np.float32,
                is_sorted=True,
                allow_empty=True
            ),
            REFERENCES_ALL + REFERENCES_SORTED
        )(self)

def test_unsorted_segment_ops(self):
        SegmentsTester()._test(
            'UnsortedSegment',
            hu.segmented_tensor(
                dtype=np.float32,
                is_sorted=False,
                allow_empty=True
            ),
            REFERENCES_ALL,
        )(self)

def test_unsorted_segment_ops_gpu(self):
        SegmentsTester()._test(
            'UnsortedSegment',
            hu.segmented_tensor(
                dtype=np.float32,
                is_sorted=False,
                allow_empty=True,
            ),
            REFERENCES_ALL,
            gpu=workspace.has_gpu_support,
            grad_check=False,
        )(self)

def test_sparse_sorted_segment_ops(self):
        SegmentsTester()._test(
            'SparseSortedSegment',
            hu.sparse_segmented_tensor(
                dtype=np.float32,
                is_sorted=True,
                allow_empty=True
            ),
            REFERENCES_ALL
        )(self)

def test_sparse_unsorted_segment_ops(self):
        SegmentsTester()._test(
            'SparseUnsortedSegment',
            hu.sparse_segmented_tensor(
                dtype=np.float32,
                is_sorted=False,
                allow_empty=True
            ),
            REFERENCES_ALL
        )(self)

def test_lengths_ops(self):
        LengthsTester()._test(
            'Lengths',
            hu.lengths_tensor(
                dtype=np.float32,
                min_value=1,
                max_value=5,
                allow_empty=True
            ),
            REFERENCES_ALL + REFERENCES_LENGTHS_ONLY,
        )(self)

def test_sparse_lengths_ops(self):
        for itype in [np.int32, np.int64]:
            LengthsTester()._test(
                'SparseLengths',
                hu.sparse_lengths_tensor(
                    dtype=np.float32,
                    min_value=1,
                    max_value=5,
                    allow_empty=True,
                    itype=itype,
                ),
                REFERENCES_ALL,
            )(self)

@unittest.skipIf(not workspace.has_gpu_support, "No gpu support")
    @given(**hu.gcs)
    def test_unsorted_sums_large(self, gc, dc):
        X = np.random.rand(10000, 32, 12).astype(np.float32)
        segments = np.random.randint(0, 10000, size=10000).astype(np.int32)
        op = core.CreateOperator("UnsortedSegmentSum", ["X", "segments"], "out")
        self.assertDeviceChecks(dc, op, [X, segments], [0])

@unittest.skipIf(not workspace.has_gpu_support, "No gpu support")
    @given(**hu.gcs)
    def test_sorted_segment_range_mean(self, gc, dc):
        X = np.random.rand(6, 32, 12).astype(np.float32)
        segments = np.array([0, 0, 1, 1, 2, 3]).astype(np.int32)
        op = core.CreateOperator(
            "SortedSegmentRangeMean",
            ["X", "segments"],
            "out"
        )
        self.assertDeviceChecks(dc, op, [X, segments], [0])
        self.assertGradientChecks(gc, op, [X, segments], 0, [0])

@unittest.skipIf(not workspace.has_gpu_support, "No gpu support")
    @given(**hu.gcs)
    def test_sorted_segment_range_log_mean_exp(self, gc, dc):
        X = np.random.rand(7, 32, 12).astype(np.float32)
        segments = np.array([0, 0, 1, 1, 2, 2, 3]).astype(np.int32)
        op = core.CreateOperator(
            "SortedSegmentRangeLogMeanExp",
            ["X", "segments"],
            "out"
        )
        self.assertDeviceChecks(dc, op, [X, segments], [0])
        self.assertGradientChecks(gc, op, [X, segments], 0, [0])

@unittest.skipIf(not workspace.has_gpu_support, "No gpu support")
    @given(**hu.gcs)
    def test_unsorted_means_large(self, gc, dc):
        X = np.random.rand(10000, 31, 19).astype(np.float32)
        segments = np.random.randint(0, 10000, size=10000).astype(np.int32)
        op = core.CreateOperator("UnsortedSegmentMean", ["X", "segments"], "out")
        self.assertDeviceChecks(dc, op, [X, segments], [0])

@serial.given(
        inputs=hu.lengths_tensor(
            dtype=np.float32,
            min_value=1,
            max_value=5,
            allow_empty=True,
        ),
        **hu.gcs
    )
    def test_lengths_sum(self, inputs, gc, dc):
        X, Y = inputs
        op = core.CreateOperator("LengthsSum", ["X", "Y"], "out")

def ref(D, L):
            R = np.zeros(shape=(L.size, ) + D.shape[1:], dtype=D.dtype)
            line = 0
            for g in range(L.size):
                for _ in range(L[g]):
                    if len(D.shape) > 1:
                        R[g, :] += D[line, :]
                    else:
                        R[g] += D[line]
                    line += 1
            return [R]

self.assertReferenceChecks(gc, op, [X, Y], ref)
        self.assertDeviceChecks(dc, op, [X, Y], [0])
        self.assertGradientChecks(gc, op, [X, Y], 0, [0])

@serial.given(
        inputs=hu.sparse_lengths_tensor(
            dtype=np.float32,
            min_value=1,
            max_value=5,
            allow_empty=True
        ),
        **hu.gcs
    )
    def test_sparse_lengths_sum(self, inputs, gc, dc):
        X, Y, Z = inputs
        op = core.CreateOperator("SparseLengthsSum", ["X", "Y", "Z"], "out")

def ref(D, I, L):
            R = np.zeros(shape=(L.size, ) + D.shape[1:], dtype=D.dtype)
            line = 0
            for g in range(L.size):
                for _ in range(L[g]):
                    if len(D.shape) > 1:
                        R[g, :] += D[I[line], :]
                    else:
                        R[g] += D[I[line]]
                    line += 1
            return [R]

self.assertReferenceChecks(gc, op, [X, Y, Z], ref)
        self.assertDeviceChecks(dc, op, [X, Y, Z], [0])
        self.assertGradientChecks(gc, op, [X, Y, Z], 0, [0])

@serial.given(
        inputs=hu.lengths_tensor(
            dtype=np.float32,
            min_value=1,
            max_value=5,
            allow_empty=True,
        ),
        **hu.gcs
    )
    def test_lengths_mean(self, inputs, gc, dc):
        X, Y = inputs
        op = core.CreateOperator("LengthsMean", ["X", "Y"], "out")

return [R]

self.assertReferenceChecks(gc, op, [X, Y], ref)
        self.assertDeviceChecks(dc, op, [X, Y], [0])
        self.assertGradientChecks(gc, op, [X, Y], 0, [0])

@serial.given(
        inputs=hu.sparse_lengths_tensor(
            dtype=np.float32,
            min_value=1,
            max_value=5,
            allow_empty=True
        ),
        **hu.gcs
    )
    def test_sparse_lengths_mean(self, inputs, gc, dc):
        X, Y, Z = inputs
        op = core.CreateOperator("SparseLengthsMean", ["X", "Y", "Z"], "out")

if L[g] > 1:
                    if len(D.shape) > 1:
                        R[g, :] = R[g, :] / L[g]
                    else:
                        R[g] = R[g] / L[g]

return [R]

self.assertReferenceChecks(gc, op, [X, Y, Z], ref)
        self.assertDeviceChecks(dc, op, [X, Y, Z], [0])
        self.assertGradientChecks(gc, op, [X, Y, Z], 0, [0])

@serial.given(
        grad_on_weights=st.booleans(),
        inputs=hu.sparse_lengths_tensor(
            dtype=np.float32,
            min_value=1,
            max_value=5,
            allow_empty=True
        ),
        seed=st.integers(min_value=0, max_value=100),
        **hu.gcs
    )
    def test_sparse_lengths_weighted_sum(
            self, grad_on_weights, inputs, seed, gc, dc):
        D, I, L = inputs

np.random.seed(int(seed))

W = np.random.rand(I.size).astype(np.float32)
        op = core.CreateOperator(
            "SparseLengthsWeightedSum",
            ["D", "W", "I", "L"],
            "out",
            grad_on_weights=grad_on_weights)
        self.assertDeviceChecks(dc, op, [D, W, I, L], [0])
        self.assertReferenceChecks(
            device_option=gc,
            op=op,
            inputs=[D, W, I, L],
            reference=sparse_lengths_weighted_sum_ref,
            threshold=1e-4,
            output_to_grad='out',
            grad_reference=partial(
                sparse_lengths_weighted_sum_grad_ref,
                grad_on_weights=grad_on_weights),
        )
        self.assertGradientChecks(gc, op, [D, W, I, L], 0, [0])
        if grad_on_weights:
            self.assertGradientChecks(gc, op, [D, W, I, L], 1, [0])

@given(**hu.gcs)
    def test_sparse_lengths_indices_in_gradient_sum_gpu(self, gc, dc):
        X = np.random.rand(3, 3, 4, 5).astype(np.float32)
        Y = np.asarray([3, 3, 2]).astype(np.int32)
        Z = np.random.randint(0, 50, size=8).astype(np.int64)
        op = core.CreateOperator(
            "SparseLengthsIndicesInGradientSumGradient", ["X", "Y", "Z"], "out"
        )
        self.assertDeviceChecks(dc, op, [X, Y, Z], [0])

@given(**hu.gcs)
    def test_sparse_lengths_indices_in_gradient_mean_gpu(self, gc, dc):
        X = np.random.rand(3, 3, 4, 5).astype(np.float32)
        Y = np.asarray([3, 3, 2]).astype(np.int32)
        Z = np.random.randint(0, 50, size=8).astype(np.int64)
        op = core.CreateOperator(
            "SparseLengthsIndicesInGradientMeanGradient", ["X", "Y", "Z"], "out"
        )
        self.assertDeviceChecks(dc, op, [X, Y, Z], [0])

@given(**hu.gcs_cpu_only)
    def test_legacy_sparse_and_lengths_sum_gradient(self, gc, dc):
        X = np.random.rand(3, 64).astype(np.float32)
        Y = np.asarray([20, 20, 10]).astype(np.int32)
        workspace.FeedBlob("X", X)
        workspace.FeedBlob("Y", Y)
        test_net = core.Net("test_net")
        test_net.SparseLengthsSumGradient(["X", "Y"], "out1")
        test_net.LengthsSumGradient(["X", "Y"], "out2")
        workspace.RunNetOnce(test_net)
        out1 = workspace.FetchBlob("out1")
        out2 = workspace.FetchBlob("out2")
        self.assertTrue((out1 == out2).all())

@given(**hu.gcs)
    @settings(deadline=10000)
    def test_sparse_lengths_sum_invalid_index(self, gc, dc):
        D = np.random.rand(50, 3, 4, 5).astype(np.float32)
        I = (np.random.randint(0, 10000, size=10) + 10000).astype(np.int64)
        L = np.asarray([4, 4, 2]).astype(np.int32)
        op = core.CreateOperator(
            "SparseLengthsSum",
            ["D", "I", "L"],
            "out")
        workspace.FeedBlob('D', D)
        workspace.FeedBlob('I', I)
        workspace.FeedBlob('L', L)
        with self.assertRaises(RuntimeError):
            workspace.RunOperatorOnce(op)

@serial.given(**hu.gcs_cpu_only)
    def test_sparse_lengths_positional_weighted_sum(
            self, gc, dc):
        D = np.random.rand(50, 3, 4, 5).astype(np.float32)
        W = np.random.rand(50).astype(np.float32)
        indices = np.random.randint(0, 50, size=10).astype(np.int64)
        L = np.asarray([4, 4, 2]).astype(np.int32)
        op = core.CreateOperator(
            "SparseLengthsPositionalWeightedSum",
            ["D", "W", "indices", "L"],
            "out")

def ref_sparse(D, W, indices, L):
            workspace.FeedBlob("L", L)
            lengths_range_fill_op = core.CreateOperator(
                "LengthsRangeFill", ["L"], ["L_pos_seq"])
            workspace.RunOperatorOnce(lengths_range_fill_op)

workspace.FeedBlob("W", W)
            gather_op = core.CreateOperator(
                "Gather", ["W", "L_pos_seq"], ["W_gathered"])
            workspace.RunOperatorOnce(gather_op)

workspace.FeedBlob("D", D)
            workspace.FeedBlob("indices", indices)
            sparse_op = core.CreateOperator(
                "SparseLengthsWeightedSum",
                ["D", "W_gathered", "indices", "L"],
                "out_ref")
            workspace.RunOperatorOnce(sparse_op)

return (workspace.FetchBlob("out_ref"),)

self.assertReferenceChecks(
            gc, op, [D, W, indices, L], ref_sparse)

@unittest.skipIf(not workspace.has_gpu_support, "No GPU support")
    @given(
        input=hu.tensor(min_dim=2, max_dim=2, max_value=20, dtype=np.float16),
        data_strategy=st.data(),
        is_mean=st.booleans(),
        **hu.gcs
    )
    @settings(deadline=None)
    def test_sparse_lengths_fp16(self, input, data_strategy, is_mean, gc, dc):
        m = input.shape[0]

lengths = data_strategy.draw(
            hu.tensor(
                max_dim=1,
                max_value=input.shape[0],
                dtype=np.int32,
                elements=st.integers(min_value=0, max_value=27),
            )
        )
        lengths_sum = int(np.sum(lengths).item())

indices = data_strategy.draw(
            hu.arrays(
                [lengths_sum], dtype=np.int64, elements=st.sampled_from(np.arange(m))
            )
        )
        if is_mean:
            op = core.CreateOperator(
                "SparseLengthsMean", ["input", "indices", "lengths"], "out"
            )
            self.assertReferenceChecks(gc, op, [input, indices, lengths], sparse_lengths_mean_ref)

else:
            op = core.CreateOperator(
                "SparseLengthsSum", ["input", "indices", "lengths"], "out"
            )
            self.assertReferenceChecks(gc, op, [input, indices, lengths], sparse_lengths_sum_ref)

# @given(
   #     inputs=hu.lengths_tensor(
   #         dtype=np.float32,
   #         min_value=1,
   #         max_value=5,
   #         min_dim=1,
   #         max_dim=1,
   #         allow_empty=False,
   #     ),
   #     **hu.gcs
   # )
   # def test_lengths_max_gpu(self, inputs, gc, dc):
   #     def lengths_max_ref(I, L):
   #         R = np.zeros(shape=(len(L)), dtype=I.dtype)
   #         line = 0
   #         for g in range(len(L)):
   #             for i in range(L[g]):
   #                 if i == 0:
   #                     R[g] = I[line]
   #                 else:
   #                     R[g] = max(R[g], I[line])
   #                 line += 1
   #         return [R]

#     X, lengths = inputs
   #     op = core.CreateOperator("LengthsMax", ["X", "lengths"], "out")
   #     self.assertDeviceChecks(dc, op, [X, lengths], [0])
   #     self.assertReferenceChecks(
   #         device_option=gc,
   #         op=op,
   #         inputs=[X, lengths],
   #         reference=lengths_max_ref,
   #         threshold=1e-4,
   #         output_to_grad='out',
   #     )

if __name__ == "__main__":
    import unittest
    unittest.main()

cmd: