let set width height array x y elt =
(*   assert (0 <= x && x < width  &&  0 <= y && y < height); *)
  array.(to_1d width height x y) <- elt