use larger output dtype
Merge remote-tracking branch 'esc/factorize' into factorize