def METHOD_NAME(language):
    @types( 'float32[:]', 'float32', 'int')
    @types( 'double[:]' , 'double' , 'int')
    def array_real_1d_scalar_add( x, a, x_len ):
        for i in range(x_len):
            x[i] += a
    f1 = array_real_1d_scalar_add
    f2 = epyccel(f1, language=language)
    for t in float_types:
        size = randint(1, 30)
        x1 = uniform(np.finfo(t).max / 2, size=size)
        x2 = np.copy(x1)
        a = uniform(np.finfo(t).max / 2)
        f1(x1, a, size)
        f2(x2, a, size)
        assert np.array_equal( x1, x2 )