diff --git a/.github/workflows/jax_tests.yml b/.github/workflows/jax_tests.yml
index 809577b52d..2b21cc8aac 100644
--- a/.github/workflows/jax_tests.yml
+++ b/.github/workflows/jax_tests.yml
@@ -26,6 +26,9 @@ jobs:
             0.8.3,
             0.9.1,
             0.9.2,
+            0.10.0,
+            0.10.1,
+            0.10.2,
           ]
           # 0.7.0 have performance issues but we still support it, see diffrax#680
           # 0.7.1 fails with equinox, see equinox#1081
diff --git a/desc/batching.py b/desc/batching.py
index eec6bedb51..5ddee4b2b7 100644
--- a/desc/batching.py
+++ b/desc/batching.py
@@ -9,9 +9,7 @@
     _check_output_dtype_jacrev,
     _jacfwd_unravel,
     _jacrev_unravel,
-    _jvp,
     _std_basis,
-    _vjp,
 )
 from jax._src.api_util import _ensure_index, argnums_partial, check_callable
 from jax._src.numpy.vectorize import (
@@ -409,6 +407,27 @@ def wrapped(*args, **kwargs):
     return wrapped
 
 
+def _argnums_partial(fun, argnums, args, kwargs):
+    """Bind every argument except those in ``argnums``.
+
+    This mirrors JAX's internal ``argnums_partial`` but returns a plain callable
+    (and the tuple of differentiated arguments) instead of a
+    ``linear_util.WrappedFun``. A plain callable is what the public
+    ``jax.jvp``/``jax.vjp`` expect, so we avoid the private ``jax._src`` helpers
+    (e.g. ``_jvp``, which was removed in JAX 0.10.2).
+    """
+    argnums_t = (argnums,) if isinstance(argnums, int) else tuple(argnums)
+    dyn_args = tuple(args[i] for i in argnums_t)
+
+    def f_partial(*dyn):
+        full_args = list(args)
+        for i, a in zip(argnums_t, dyn):
+            full_args[i] = a
+        return fun(*full_args, **kwargs)
+
+    return f_partial, dyn_args
+
+
 def jacfwd_chunked(
     fun,
     argnums=0,
@@ -463,18 +482,17 @@ def jacfwd_chunked(
 
     @wraps(fun, docstr=docstr, argnums=argnums)
     def jacfun(*args, **kwargs):
-        f = lu.wrap_init(fun, kwargs)
-        f_partial, dyn_args = argnums_partial(
-            f, argnums, args, require_static_args_hashable=False
-        )
+        f_partial, dyn_args = _argnums_partial(fun, argnums, args, kwargs)
         tree_map(partial(_check_input_dtype_jacfwd, holomorphic), dyn_args)
         if not has_aux:
-            pushfwd = partial(_jvp, f_partial, dyn_args)
+            pushfwd = lambda tangents: jax.jvp(f_partial, dyn_args, tangents)
             y, jac = vmap_chunked(pushfwd, chunk_size=chunk_size)(_std_basis(dyn_args))
             y = tree_map(lambda x: x[0], y)
             jac = tree_map(lambda x: jnp.moveaxis(x, 0, -1), jac)
         else:
-            pushfwd = partial(_jvp, f_partial, dyn_args, has_aux=True)
+            pushfwd = lambda tangents: jax.jvp(
+                f_partial, dyn_args, tangents, has_aux=True
+            )
             y, jac, aux = vmap_chunked(pushfwd, chunk_size=chunk_size)(
                 _std_basis(dyn_args)
             )
@@ -550,15 +568,12 @@ def jacrev_chunked(
 
     @wraps(fun, docstr=docstr, argnums=argnums)
     def jacfun(*args, **kwargs):
-        f = lu.wrap_init(fun, kwargs)
-        f_partial, dyn_args = argnums_partial(
-            f, argnums, args, require_static_args_hashable=False
-        )
+        f_partial, dyn_args = _argnums_partial(fun, argnums, args, kwargs)
         tree_map(partial(_check_input_dtype_jacrev, holomorphic, allow_int), dyn_args)
         if not has_aux:
-            y, pullback = _vjp(f_partial, *dyn_args)
+            y, pullback = jax.vjp(f_partial, *dyn_args)
         else:
-            y, pullback, aux = _vjp(f_partial, *dyn_args, has_aux=True)
+            y, pullback, aux = jax.vjp(f_partial, *dyn_args, has_aux=True)
         tree_map(partial(_check_output_dtype_jacrev, holomorphic), y)
         jac = vmap_chunked(pullback, chunk_size=chunk_size)(_std_basis(y))
         jac = jac[0] if isinstance(argnums, int) else jac
diff --git a/requirements.txt b/requirements.txt
index 638e183cd3..3e5c220c30 100644
--- a/requirements.txt
+++ b/requirements.txt
@@ -1,4 +1,4 @@
-jax >= 0.6.2, != 0.7.1, < 0.10.0
+jax >= 0.6.2, != 0.7.1, <= 0.10.1
 colorama <= 0.4.6
 diffrax >= 0.6.0, <= 0.7.2
 equinox >=0.11.10, <=0.13.8
diff --git a/tests/test_examples.py b/tests/test_examples.py
index 81d8f660c7..f8fc4434a5 100644
--- a/tests/test_examples.py
+++ b/tests/test_examples.py
@@ -174,7 +174,14 @@ def test_solve_bounds():
     obj = ObjectiveFunction(
         ForceBalance(normalize=False, normalize_target=False, bounds=(-3e3, 3e3), eq=eq)
     )
-    eq.solve(objective=obj, ftol=1e-16, xtol=1e-16, maxiter=200, verbose=3)
+    eq.solve(
+        objective=obj,
+        ftol=1e-16,
+        xtol=1e-16,
+        maxiter=200,
+        verbose=3,
+        options={"tr_method": "svd"},
+    )
 
     # check that all errors are nearly 0, since residual values are within target bounds
     f = obj.compute_scaled_error(obj.x(eq))
diff --git a/tests/test_integrals.py b/tests/test_integrals.py
index d26761ec19..cd4255cefa 100644
--- a/tests/test_integrals.py
+++ b/tests/test_integrals.py
@@ -724,6 +724,7 @@ def filter(z1, z2):
         return z1[mask], z2[mask]
 
     @pytest.mark.unit
+    @pytest.mark.xfail(reason="will be fixed by #2199")
     def test_z1_first(self):
         """Case where straight line through first two intersects is in epigraph."""
         start = np.pi / 3