def get_bonds(input_group):
out_list = []<EOL>for i in range(len(input_group.bond_order_list)):<EOL><INDENT>out_list.append((input_group.bond_atom_list[i * <NUM_LIT:2>], input_group.bond_atom_list[i * <NUM_LIT:2> + <NUM_LIT:1>],))<EOL><DEDENT>return out_list<EOL>
Utility function to get indices (in pairs) of the bonds.