@@ -1133,7 +1133,7 @@ def test_interpolate_1_plus_2_cycloaddition(self):
11331133 self .assertEqual (len (ts_xyz ['symbols' ]), 10 )
11341134 self .assertFalse (colliding_atoms (ts_xyz ),
11351135 msg = f'Collision in 1+2_Cycloaddition TS:\n { xyz_to_str (ts_xyz )} ' )
1136- expected_ts = """C 1.59999925 -0.11618654 -0.14166302
1136+ expected_ts_1 = """C 1.59999925 -0.11618654 -0.14166302
11371137C 0.29517860 -0.02143486 -0.02613492
11381138C -1.15821797 -1.12490772 0.14486040
11391139C -0.81238032 0.84414025 0.04444949
@@ -1143,7 +1143,18 @@ def test_interpolate_1_plus_2_cycloaddition(self):
11431143H -1.52801447 -1.64655150 -0.72678867
11441144H -0.94547237 1.40230195 0.96062403
11451145H -1.11212744 1.33826544 -0.86912905"""
1146- self .assertTrue (any (almost_equal_coords (ts , str_to_xyz (expected_ts )) for ts in ts_xyzs ))
1146+ expected_ts_2 = """C 1.59999925 -0.11618654 -0.14166302
1147+ C 0.29517860 -0.02143486 -0.02613492
1148+ C -0.92013120 -0.71833111 0.10894610
1149+ C -0.99229728 1.28107087 0.04554500
1150+ H 2.21797993 0.77036923 -0.22897655
1151+ H 2.09015362 -1.08321135 -0.15246324
1152+ H -1.12327237 -1.17593811 1.06705013
1153+ H -1.28992770 -1.23997489 -0.76270297
1154+ H -1.12538933 1.83923257 0.96171954
1155+ H -1.29204440 1.77519606 -0.86803354"""
1156+ self .assertTrue (any (almost_equal_coords (ts , str_to_xyz (expected_ts_1 )) for ts in ts_xyzs ) or
1157+ any (almost_equal_coords (ts , str_to_xyz (expected_ts_2 )) for ts in ts_xyzs ))
11471158 # The TS should have extended forming bonds compared to the product.
11481159 # Get atom map to find the carbene atom in the product.
11491160 atom_map = map_rxn (rxn = rxn )
0 commit comments