Blender crashes when adding many objects with python?

Question

I am creating a program in blender that makes cubes based on the values of an array. For example, when one of the values is one it creates a cube at a specific point and then it goes on to print another cube until you have "printed" a numpy array. Here is my code:

import bpy
import numpy as np

def print_array(a, size=1):
    count = 0
    x = 0
    for i in a:
        y = 0
        for j in i:
            z = 0
            for k in j:
                if k == 1:
                    bpy.ops.mesh.primitive_cube_add(size=size, location=(x, y, z))
                    count += 1
                z += size
            y += size
        x += size
    return count

a = np.ones((3, 3, 3))

print_array(a)

While this does work for smaller numpy arrays blender crashes when I try to "print" larger arrays. Is there a less tasking way to do this without making the arrays smaller?