Save result of construct_zone_axes
It would be very convenient to be able to save the atomic planes found using construct_zone_axes()
in an HDF5 file. Currently, a workaround is to save each set of parallel planes using numpy:
sublattice.construct_zone_axes()
zone_vector0 = sublattice.zones_axis_average_distances[0]
atom_planes0 = sublattice.atom_planes_by_zone_vector[zone_vector0]
np.save('atom_planes0.npy', atom_planes0)
but with many atoms the recursion limit is quickly reached.