The DoubleCoset function in GAP has not been wrapped yet, so this is currently "not implemented". But, you can always ask GAP to compute (double) cosets from within Sage:
sage: gap.eval('g:=Group((1,2,3,4),(1,2));;')