code: def _np_bn_testing(x, scale, bias, rm, rv, momentum=0.1, e=1e-5):
      channel = x.shape[1]
      np.testing.assert_array_almost_equal(
          scale.shape, (1, channel, 1, 1)
      )
      np.testing.assert_array_almost_equal(
          bias.shape, (1, channel, 1, 1)
      )
      np.testing.assert_array_almost_equal(rm.shape, (1, channel, 1, 1))
      np.testing.assert_array_almost_equal(rv.shape, (1, channel, 1, 1))
      return scale * (x - rm) / np.sqrt(rv + e) + bias
  
imports: from __future__ import division
  import unittest
  import math
  import numpy
  from singa import singa_wrap
  from singa import tensor
  from cuda_helper import gpu_dev, cpu_dev
label_desc: Use image node_id 2 for calling a global function with example usage: _np_bn_testing(x, scale, bias, rm, rv, momentum, e) and returns: unknown