s = complete_sandpile(4)
d = SandpileDivisor(s,[0,1,-2,3])
d.r_of_D()

from sage.interfaces.four_ti_2 import four_ti_2
four_ti_2.graver([1,2,3]) # optional - 4ti2