Extend the `linalg.reduce` kernel to support reductions across multiple axes

This issue has no description.

The issue names the linalg.reduce kernel but gives no file, test, or entry point. Start by locating the kernel implementation and existing axis-handling tests, then determine the intended behavior for reductions across multiple axes. Done means the supported behavior is implemented and covered by appropriate tests.

Written by the indexing model from the issue text.