Scripts: dividelinks.py

# Insert N markers at each selected link.  Example Chimera command
#
#     runscript dividelinks.py 4
#
def divide_links(links, n):
    from VolumePath import Link
    for link in links:
        m1, m2 = link.marker1, link.marker2
        mprev = m1
        mset = m1.marker_set
        for i in range(n):
            f = float(i+1) / (n+1)
            xyz = tuple((1-f)*x1+f*x2 for x1, x2 in zip(m1.xyz(), m2.xyz()))
            m = mset.place_marker(xyz, link.rgba(), m1.radius())
            Link(mprev, m, link.rgba(), link.radius())
            mprev = m
        Link(mprev, m2, link.rgba(), link.radius())
        link.delete()

n = int(arguments[0])
from VolumePath import markerset
links = markerset.selected_links()
divide_links(links, n)