Fixed (hid?) bug in calculateBondSymmetryNumber().

The function Structure.split() returns a list of structures that result from
the splitting. In calculateBondSymmetryNumber(), this function was assumed to
return exactly two structures, which it should if things are going normally.
However, this behavior is not guaranteed, so a check was added to ensure that
exactly two structures are returned from split(). If not, the function returns
1, which is as if it was not called for that bond.

source/rmg/structure.py

@@ -1021,9 +1021,10 @@ def calculateBondSymmetryNumber(self, bond):
 				else:
 					structure = Structure(self.atoms(), self.bonds())
 					structure.removeBond(bond)
-					structure1, structure2 = structure.split()
-
+					structures = structure.split()
+					if len(structures) != 2: return symmetryNumber
 
+					structure1, structure2 = structures
 					if bond.atoms[0] in structure1.atoms(): structure1.removeAtom(bond.atoms[0])
 					if bond.atoms[1] in structure1.atoms(): structure1.removeAtom(bond.atoms[1])
 					if bond.atoms[0] in structure2.atoms(): structure2.removeAtom(bond.atoms[0])