Try to add work-around for sympy/sympy#26645

symengine/lib/symengine_wrapper.in.pyx

@@ -31,6 +31,13 @@ try:
 except ImportError:
     have_numpy = False
 
+try:
+    import flint as _flint
+    have_flint_py = True
+except ImportError:
+    _flint = None
+    have_flint_py = False
+
 
 class SympifyError(Exception):
     pass
@@ -499,6 +506,17 @@ def sympy2symengine(a, raise_error=False):
         elif isinstance(a, sympy.ConditionSet):
             return conditionset(*(a.args))
 
+    if have_flint_py:
+        if isinstance(a, _flint.types.nmod.nmod):
+            # Work around for sympy/sympy#26645
+            class _modular_integer(sympy.polys.domains.modularinteger.ModularInteger):
+                mod = int(a.modulus())
+                dom = sympy.polys.domains.ZZ
+                sym = True
+
+            b = _modular_integer(int(a))
+            return PyNumber(b, sympy_module)
+
     if raise_error:
         raise SympifyError(("sympy2symengine: Cannot convert '%r' (of type %s)"
                             " to a symengine type.") % (a, type(a)))

symengine/tests/test_sympy_conv.py

@@ -778,6 +778,7 @@ def test_pynumber():
     assert isinstance(b, PyNumber)
     assert b == a                  # Check equality via SymEngine
     assert a == b                  # Check equality via SymPy
+    assert (b-a).simplify() == 0
     assert str(a) == str(b)
 
     a = 1 - a