LA: restrict modified sparse*diagonal matrix product to square diagonal matrices only