* Previous implementation would silently ignore explicit axis for dimensions that were not broadcastable.
* Added note about difference in behavior between Aesara and Numpy functions

Closes #830
Co-authored-by: Brandon T. Willard <971601+brandonwillard@users.noreply.github.com>