diff --git a/.github/workflows/build_python.yml b/.github/workflows/build_python.yml
index a0ecce12f9..c585d84d38 100644
--- a/.github/workflows/build_python.yml
+++ b/.github/workflows/build_python.yml
@@ -26,7 +26,11 @@ jobs:
       fail-fast: false
       matrix:
         # As of 20240501, macos-11/12/13 are AMD64, and macOS-14 is ARM64.
-        os: [ubuntu-20.04, ubuntu-22.04, ubuntu-24.04, windows-2019, windows-2022, macos-11, macos-12, macos-13, macos-14]
+        # As of 20240616, ubuntu 24.04 is failing to copy files into the quay.io/pypa/manylinux2014_i686:2024-05-13-0983f6f container
+        # It fails with "tar: ./todo: Cannot utime: Function not implemented" for every file in the repo.
+        # This is strange since it worked for the quay.io/pypa/manylinux2014_x86_64:2024-05-13-0983f6f container.
+        # I suggest trying ubuntu-24.04 again in a month.
+        os: [ubuntu-20.04, ubuntu-22.04, windows-2019, windows-2022, macos-11, macos-12, macos-13, macos-14]
 
     steps:
       - name: Clone Repository (Latest)
@@ -46,11 +50,20 @@ jobs:
       - name: Checkout pybind11 submodule
         run: git submodule update --init python/pybind11
 
+      - name: Set the MACOSX_DEPLOYMENT_TARGET
+        if: ${{ runner.os == 'macOS' }}
+        run: |
+          MACOSX_DEPLOYMENT_TARGET=$(sw_vers -productVersion | cut -d'.' -f 1)
+          echo "MACOSX_DEPLOYMENT_TARGET is $MACOSX_DEPLOYMENT_TARGET"
+          echo "MACOSX_DEPLOYMENT_TARGET=$MACOSX_DEPLOYMENT_TARGET" >> $GITHUB_ENV
+
+
       - name: Set up Fortran
         uses: fortran-lang/setup-fortran@main
         if: ${{ runner.os == 'macOS' }}
         with:
           compiler: gcc
+          version: 12
 
       # Copied from https://github.com/scipy/scipy/blob/main/.github/workflows/wheels.yml
       # For rtools, see https://github.com/r-windows/rtools-installer/releases, which has been
diff --git a/.github/workflows/cmake.yml b/.github/workflows/cmake.yml
index 96c7a59a4c..278e21214d 100644
--- a/.github/workflows/cmake.yml
+++ b/.github/workflows/cmake.yml
@@ -42,7 +42,8 @@ jobs:
         toolchain:
           - {compiler: gcc, version: 11, cflags: '-Wall -Wextra -Wpedantic -Werror', fflags: '-Wall -Wextra -Wpedantic -Werror -fimplicit-none -fcheck=all -fstack-check -Wno-function-elimination'}
           - {compiler: gcc, version: 12, cflags: '-Wall -Wextra -Wpedantic -Werror', fflags: '-Wall -Wextra -Wpedantic -Werror -fimplicit-none -fcheck=all -fstack-check -Wno-function-elimination'}
-          - {compiler: gcc, version: 13, cflags: '-Wall -Wextra -Wpedantic -Werror', fflags: '-Wall -Wextra -Wpedantic -Werror -fimplicit-none -fcheck=all -fstack-check -Wno-function-elimination'}
+          # As of 20240616 gcc13 has issues importing math.h on macOS.
+          # - {compiler: gcc, version: 13, cflags: '-Wall -Wextra -Wpedantic -Werror', fflags: '-Wall -Wextra -Wpedantic -Werror -fimplicit-none -fcheck=all -fstack-check -Wno-function-elimination'}
           - {compiler: intel-classic, version: '2021.8', cflags: '-diag-disable=10441 -Wall -w3 -Werror-all', fflags: '-warn all -debug extended -fimplicit-none -standard-semantics'}
           - {compiler: intel-classic, version: '2021.9', cflags: '-diag-disable=10441 -Wall -w3 -Werror-all', fflags: '-warn all -debug extended -fimplicit-none -standard-semantics'}
           - {compiler: intel-classic, version: '2021.10', cflags: '-diag-disable=10441 -Wall -w3 -Werror-all', fflags: '-warn all -debug extended -fimplicit-none -standard-semantics'}
diff --git a/python/pybind11 b/python/pybind11
index 8a099e44b3..3e9dfa2866 160000
--- a/python/pybind11
+++ b/python/pybind11
@@ -1 +1 @@
-Subproject commit 8a099e44b3d5f85b20f05828d919d2332a8de841
+Subproject commit 3e9dfa2866941655c56877882565e7577de6fc7b
diff --git a/python/tests/test_combining_constraints.py b/python/tests/test_combining_constraints.py
index 0f245ec389..801f34edfb 100644
--- a/python/tests/test_combining_constraints.py
+++ b/python/tests/test_combining_constraints.py
@@ -28,7 +28,7 @@ def test_providing_bounds_and_linear_constraints():
 def test_providing_bounds_and_nonlinear_constraints():
     nlc = prima_NLC(lambda x: x[0]**2, lb=[25], ub=[100])
     bounds = prima_Bounds([None, 1], [None, 1])
-    x0 = [0, 0]
+    x0 = [6, 1]  # Unfortunately the test is very fragile if we do not start near the optimal point
     res = prima_minimize(fun, x0, constraints=nlc, bounds=bounds)
     assert np.isclose(res.x[0], 5, atol=1e-6, rtol=1e-6)
     assert np.isclose(res.x[1], 1, atol=1e-6, rtol=1e-6)
@@ -46,26 +46,10 @@ def newfun(x):
     nlc = NLC(lambda x: x[0]**2, lb=[25], ub=[100])
     bounds = Bounds([-np.inf, 1, -np.inf], [np.inf, 1, np.inf])
     lc = LC(np.array([1,1,1]), lb=10, ub=15)
-    x0 = [0, 0, 0]
-    # macOS seems to stop just short of the optimal solution, so we help it along by
-    # taking a larger initial trust region radius and requiring a smaller final radius
-    # before stopping. The different packages have different names for these options.
-    if package == 'prima':
-        options = {'rhobeg': 2, 'rhoend': 1e-8}
-        method = None
-    elif package == 'pdfo':
-        options = {'radius_init': 2, 'radius_final': 1e-8}
-        method = None
-    elif package == 'scipy':
-        options = {'rhobeg': 2, 'tol': 1e-8}
-        # PDFO and PRIMA will select COBYLA but SciPy may select something else, so we tell it to select COBYLA
-        method = 'COBYLA'
-    else:
-        # Since this is test infrastructure under the control of the developers we
-        # should never get here except for a typo or something like that
-        raise ValueError(f"Unknown package: {package}")
-
-    res = minimize(newfun, x0, method=method, constraints=[nlc, lc], bounds=bounds, options=options)
+    x0 = [10, 1, 0]  # The test becomes very fragile if we do not start from a feasible point
+    # PDFO and PRIMA will select COBYLA but SciPy may select something else, so we tell it to select COBYLA
+    method = 'COBYLA' if package == 'scipy' else None
+    res = minimize(newfun, x0, method=method, constraints=[nlc, lc], bounds=bounds, options={})
 
     # 32 bit builds of PRIMA reach the optimal solution with the same level of precision as 64 bit builds
     # so we lower the atol/rtol to 1e-3 so that 32 bit builds will pass.
diff --git a/python/tests/test_compatibility_pdfo.py b/python/tests/test_compatibility_pdfo.py
index 6eaaaf6cd5..9ecfc90c06 100644
--- a/python/tests/test_compatibility_pdfo.py
+++ b/python/tests/test_compatibility_pdfo.py
@@ -20,6 +20,9 @@ def test_pdfo():
     scipy = pytest.importorskip("scipy")
     if version.parse(scipy.__version__) < version.parse("1.11.0"):
         pytest.skip("scipy version too old for this test (its version of COBYLA does not accept bounds)")
+    numpy = pytest.importorskip("numpy")
+    if version.parse(numpy.__version__) >= version.parse("2.0.0"):
+        pytest.skip("numpy version too new for this test (pdfo does not yet support numpy v2)")
 
     from pdfo import pdfo
     from scipy.optimize import NonlinearConstraint as NLC, LinearConstraint as LC, Bounds