Game of life with scipy convolution

gol.py

import numpy as np
from scipy.ndimage import convolve

dim = 10
board = np.zeros((dim, dim), dtype=np.uint8)
kernel = np.array([[1,1,1],
                   [1,0,1],
                   [1,1,1]], dtype=np.uint8)

def step(board, kernel=kernel):
    convolved = convolve(board, kernel, mode="wrap")
    return (((board == 1) & (convolved > 1) & (convolved < 4))
           | ((board == 0) & (convolved == 3))).astype(np.uint8)