let bounded () =
    let init x y =
      (if y mod 2 = 0
      then [|x+1,y;   x,y+1; x-1,y+1; x-1,y; x-1,y-1;   x,y-1|]
      else [|x+1,y; x+1,y+1;   x,y+1; x-1,y;   x,y-1; x+1,y-1|])
      |> Array.to_list
      |> List.filter (fun (i,j) -> 0 <= i && i < width  &&  0 <= j && j < height)
      |> Array.of_list
      |> Array.map (fun (i,j) -> Femtolib.Array_simulate_2d.to_1d width height i j) in
    Femtolib.Array_simulate_2d.init width height init