for a, b, c in cell_shift_indices:
                    new_coordinate_list.append(
                        [pos_in_origin_cell + np.array([a, b, c] * lattice_shift)][