diff --git a/__init__.py b/__init__.py
index e13647ef8ea294a02b8556cdf510998ff13d96ed..3b6abc9dd38428828a23bb51c16e4f1accba78d6 100644
--- a/__init__.py
+++ b/__init__.py
@@ -1,3 +1,4 @@
+import pystencils.sympy_gmpy_bug_workaround
 from pystencils.field import Field, FieldType
 from pystencils.data_types import TypedSymbol
 from pystencils.slicing import make_slice
diff --git a/sympy_gmpy_bug_workaround.py b/sympy_gmpy_bug_workaround.py
new file mode 100644
index 0000000000000000000000000000000000000000..ce2bfa3c8123c2263291210fc2d98ff1bc67de5c
--- /dev/null
+++ b/sympy_gmpy_bug_workaround.py
@@ -0,0 +1,16 @@
+# Disable gmpy backend until this bug is resolved if joblib serialize
+# See https://github.com/sympy/sympy/pull/13530
+import os
+import warnings
+
+os.environ['MPMATH_NOGMPY'] = '1'
+try:
+    import mpmath.libmp
+    # In case the user has imported sympy first, then pystencils
+    if mpmath.libmp.BACKEND == 'gmpy':
+        warnings.warn("You are using the gmpy backend. You might encounter an error 'argument is not an mpz sympy'."
+                      "This is due to a known bug in sympy/gmpy library. "
+                      "To prevent this, import pystencils first then sympy or set the environment variable "
+                      "MPMATH_NOGMPY=1")
+except ImportError:
+    pass