Previously, we have developed a methodology to quickly compare bulk structures of a solid electrolyte compound. We retrofit the symmetry checker from this work to be used in pyGlobOpt as a grouping algorithm. The grouping algorithm works by calculating an atomic distance matrix for each structure produced by pyGlobOpt . The distance matrix is converted in to a characteristic vector by taking the upper triangular of the distance matrix, converting it to a vector, and then sorting the values based on magnitude. The average absolute difference is calculated between all characteristic vectors and placed into a difference matrix of size N structures × N strutures . Next, the grouping process begins such that any structures with a difference score under an arbitrary difference threshold are grouped together. Once all structures are grouped the energy range of the clusters are checked, if this energy range is too high (by default anything over 0.1 eV is too high) then the difference threshold is halved and the structures are regrouped. This is repeated until all groups have a low-enough energy range.