Remove constraints before comparison

Fixes a bug where leaving constraints on throws an error within the ase comparison algorithm

pynta/utils.py

@@ -31,6 +31,7 @@ def get_unique_sym(geoms):
     for geom in geoms:
         adsorbate_atom_obj = read(geom)
         adsorbate_atom_obj.pbc = True
+        adsorbate_atom_obj.set_constraint() # Reset constraints for comparison
         comparision = comparator.compare(
             adsorbate_atom_obj, good_adsorbates_atom_obj_list)
 
@@ -56,6 +57,7 @@ def get_unique_sym_indices(geoms):
     for geom in geoms:
         adsorbate_atom_obj = read(geom)
         adsorbate_atom_obj.pbc = True
+        adsorbate_atom_obj.set_constraint() # Reset constraints before comparison
         comparision = comparator.compare(
             adsorbate_atom_obj, good_adsorbates_atom_obj_list)
 
@@ -86,6 +88,7 @@ def get_unique_sym_structs(geoms):
     for i,geom in enumerate(geoms_copy):
         adsorbate_atom_obj = geom
         adsorbate_atom_obj.pbc = True
+        adsorbate_atom_obj.set_constraint() # Reset constraints before comparison
         comparision = comparator.compare(
             adsorbate_atom_obj, good_adsorbates_atom_obj_list)
 
@@ -114,6 +117,7 @@ def get_unique_sym_struct_indices(geoms):
     for i,geom in enumerate(geoms_copy):
         adsorbate_atom_obj = geom
         adsorbate_atom_obj.pbc = True
+        adsorbate_atom_obj.set_constraint() # Reset constraints before comparison
         comparision = comparator.compare(
             adsorbate_atom_obj, good_adsorbates_atom_obj_list)
 
@@ -142,6 +146,7 @@ def get_unique_sym_struct_index_clusters(geoms):
     for i,geom in enumerate(geoms_copy):
         adsorbate_atom_obj = geom
         adsorbate_atom_obj.pbc = True
+        adsorbate_atom_obj.set_constraint() # Reset constraints before comparison
         comparison = None
         for j,adlist in enumerate(good_adsorbates_atom_obj_list):
             comparison = comparator.compare(adsorbate_atom_obj, [adlist[0]])
@@ -176,6 +181,7 @@ def filter_nonunique_TS_guess_indices(geoms,Es):
     for j,geom in enumerate(geoms):
         adsorbate_atom_obj = read(geom)
         adsorbate_atom_obj.pbc = True
+        adsorbate_atom_obj.set_constraint() # Reset constraints before comparison
         for i,good_adsorbate in enumerate(good_adsorbates_atom_obj_list):
             comparison = comparator.compare(adsorbate_atom_obj,good_adsorbate)
             if comparison and Es[j] < Esout[i]:
@@ -249,4 +255,4 @@ def construct_constraint(d):
     constraint_dict = copy.deepcopy(d)
     constructor = getattr(ase.constraints,constraint_dict["type"])
     del constraint_dict["type"]
-    return constructor(**constraint_dict)
+    return constructor(**constraint_dict)