Use sumpy expansion factory