ENH: add function for mulm