On Mon, 6 Nov 2023 at 03:22, T-Rex wrote:
>
> Thanks @oscar! I didn't think I would get a reply from the blog's owner :-)
>
> Re setting of the environment variable SYMPY_GROUND_TYPES=flint: (1) Would
> that be in conflict with other parts of sympy/numpy/scipy? (2) After it is
> set to
The problem is, as I noted, you are using NumPy arrays to store SymPy
expressions instead of using SymPy matrices. NumPy arrays do not
define * as matrix multiplication like sympy matrix does, but rather
as elementwise multiplication with broadcasting. It's generally a good
idea to use @ for