Segfault when using sympy symbols as initial conditions.

This segfaults:

from diffeqpy import de
from sympy import Symbol

# Exponential decay
def f(u, p, t):
    return [-u[0] * p]

u0 = Symbol('u0') # or 1.0
p = Symbol('p') # or 1.0
t = (0.0, 1.0)

problem = de.ODEProblem(f, [u0], t, p)
solution = de.solve(problem)

Changing p and u0 to floats instead of symbolic variables fixes the problem:

`u0 = 1.0`	`u0 = Symbol('u0')`
`p = 1.0`	Correct answer	Segfault
`p = Symbol('p')`	Error	Segfault

Versions


Python 3.11.6 (main, Oct  2 2023, 20:46:14) [Clang 14.0.3 (clang-1403.0.22.14.1)] on darwin
diffeqpy 2.2.0
sympy 1.12
DifferentialEquations 7.11.0
Julia Version 1.9.3
Platform Info:
  OS: macOS (arm64-apple-darwin22.4.0)
  CPU: 8 × Apple M2
  WORD_SIZE: 64
  LIBM: libopenlibm
  LLVM: libLLVM-14.0.6 (ORCJIT, apple-m1)
  Threads: 1 on 4 virtual cores

SciML / diffeqpy

Segfault when using sympy symbols as initial conditions. #129