def test_global_max_g_is_max(self, backend, elems, x, g):
        k_code = """void k(unsigned int *x, unsigned int *g) {
          if ( *g < *x ) { *g = *x; }
        }"""

        k = op2.Kernel(k_code, 'k')

        g.data[0] = nelems * 2

        op2.par_loop(k, elems, x(op2.READ), g(op2.MAX))

        assert g.data[0] == nelems * 2
Verify that op2.MAX does not reduce a maximum value smaller than the
dataset/ETHPy150Open OP2/PyOP2/test/unit/test_direct_loop.py/TestDirectLoop.test_global_max_g_is_max