Search of (A-B)-(A-B)-(...) rings in ring_statistics about matscipy HOT 5 CLOSED

Perhaps generating all rings and then filtering based on type afterwards would be a simpler approach?

from matscipy.

Perhaps generating all rings and then filtering based on type afterwards would be a simpler approach?

Thanks for reply, i think the shortest path of A-A-A... rings and A-B-A-B... rings should be different, and a little modification of the shortest path algorithm forA-B-A-B rings should be better?

from matscipy.

Thinking about this a little more, I think you could actually do it with no modifications to the shortest path algorithm, just by passing a custom cutoff dictionary to neighbour_list. If you leave out A-A and B-B from the dictionary they will never considered to be neighbours.

from matscipy.

You can also explicitly remove A-A and B-B- neighbors from the neighbor list. Try the following (untested):

import _matscipy
from matscipy.neighbours import neighbour_list

i, j, r = neighbour_list('ijD', a, cutoff)
n = a.numbers
mask = a[i] != a[j]  # True for differing elements
i = i[mask]
j = j[mask]
r = r[mask]
d = _matscipy.distances_on_graph(i, j)
ring =_matscipy.find_sp_rings(i, j, r, d, maxlength)

from matscipy.

Thanks a lot, i checked the rings in my 816 atoms in AB2, and the rings i get are all A-B-A-B, it helps me a lot.
Best regards and keep healthy!

from matscipy.