<!-- Your title above should be a short description of what
was changed. Do not include the issue number in the title. -->

#### References to other Issues or PRs
<!-- If this pull request fixes an issue, write "Fixes #NNNN" in that exact
format, e.g. "Fixes #1234" (see
https://tinyurl.com/auto-closing for more information). Also, please
write a comment on that issue linking back to this pull request once it is
open. -->

Backport of https://github.com/sympy/sympy/pull/27065 and https://github.com/sympy/sympy/pull/27066 to the 1.13 release branch and bump versions ready to release 1.13.3

#### Brief description of what is fixed or changed


#### Other comments


#### Release Notes

<!-- Write the release notes for this release below between the BEGIN and END
statements. The basic format is a bulleted list with the name of the subpackage
and the release note for this PR. For example:

* solvers
  * Added a new solver for logarithmic equations.

* functions
  * Fixed a bug with log of integers. Formerly, `log(-x)` incorrectly gave `-log(x)`.

* physics.units
  * Corrected a semantical error in the conversion between volt and statvolt which
    reported the volt as being larger than the statvolt.

or if no release note(s) should be included use:

NO ENTRY

See https://github.com/sympy/sympy/wiki/Writing-Release-Notes for more
information on how to write release notes. The bot will check your release
notes automatically to see if they are formatted correctly. -->

<!-- BEGIN RELEASE NOTES -->
* matrices
   * Handle matrices of zero columns correctly in the `Matrix` constructor. In SymPy 1.12 these produced empty matrices of the wrong shape and in 1.13 they resulted in an error. Now `Matrix([[], []]).shape == (2,0)`.

* physics.wigner
   * NumPy numbers can be used with `clebsch_gordan`. Fixes a regression introduced in SymPy 1.13.0.
<!-- END RELEASE NOTES -->