import damask
import numpy as np
size = np.ones(3)*1e-5
cells = [64,64,64]
N_grains = 200
seeds = damask.seeds.from_random(size,N_grains,cells)
grid = damask.GeomGrid.from_Voronoi_tessellation(cells,size,seeds)
grid.save(f'Polycystal_{N_grains}_{cells[0]}x{cells[1]}x{cells[2]}')
grid
cells: 64 × 64 × 64 size: 1e-05 × 1e-05 × 1e-05 m³ origin: 0.0 0.0 0.0 m # materials: 200