ldiag function

Left (i.e., pre) Multiply a matrix by a diagonal matrix