OK, updated and that test added. I had to reorganize things a bit unfortunately. Because the result type should treat scalars differently, I cannot move the `asarray` to `braodcast_arrays`.
Just noticed, that this function has no support whatsoever of structured types. But the old one didn't have either, so I think this is for another day.