ENH: support more integration algorithms

.cspell.json

@@ -99,6 +99,7 @@
     "fromdict",
     "gellmann",
     "genindex",
+    "getattr",
     "getsource",
     "gitpod",
     "gridspec",
@@ -115,7 +116,9 @@
     "isfunction",
     "isidentifier",
     "isinstance",
+    "isnan",
     "jpsi",
+    "jslink",
     "jupyterlab",
     "katex",
     "kernelspec",
@@ -173,12 +176,15 @@
     "pypi",
     "pyproject",
     "pytestconfig",
+    "quadcc",
+    "quadgk",
     "rankdir",
     "relim",
     "repr",
     "richman",
     "rightarrow",
     "risch",
+    "rombergts",
     "rtfd",
     "savefig",
     "sdist",
@@ -260,6 +266,7 @@
     "isospin",
     "itertools",
     "jupyter",
+    "Kronrod",
     "Källén",
     "lambdification",
     "lambdified",
@@ -278,8 +285,10 @@
     "pyplot",
     "pytest",
     "qrules",
+    "quadax",
     "rsplit",
     "Schwarz",
+    "SciPy",
     "setuptools",
     "spflueger",
     "struct",
@@ -292,6 +301,7 @@
     "unevaluatable",
     "unitarity",
     "unphysical",
+    "vectorizing",
     "venv",
     "weisskopf",
     "yticklabels",

.pre-commit-config.yaml

@@ -13,7 +13,7 @@ repos:
       - id: check-useless-excludes
 
   - repo: https://github.com/ComPWA/policy
-    rev: b122464266f519d41284e4a51673fb1232d1fd0b
+    rev: 0.8.3
     hooks:
       - id: check-dev-files
         args:
@@ -58,7 +58,7 @@ repos:
             metadata.vscode
 
   - repo: https://github.com/astral-sh/ruff-pre-commit
-    rev: v0.15.0
+    rev: v0.15.1
     hooks:
       - id: ruff-check
         args: [--fix]
@@ -112,7 +112,7 @@ repos:
           - --in-place
 
   - repo: https://github.com/python-jsonschema/check-jsonschema
-    rev: 0.36.1
+    rev: 0.36.2
     hooks:
       - id: check-jsonschema
         name: Check CITATION.cff
@@ -152,6 +152,6 @@ repos:
         types_or: [python, pyi, jupyter]
 
   - repo: https://github.com/astral-sh/uv-pre-commit
-    rev: 0.10.0
+    rev: 0.10.2
     hooks:
       - id: uv-lock

docs/bibliography.bib

@@ -67,7 +67,7 @@ @techreport{Chung:2014-SpinFormalismsUpdated
   url = {https://cds.cern.ch/record/186421/files/BNL-76975-2006-IR.pdf}
 }
 
-@techreport{chungFormulasAngularMomentumBarrier2015,
+@techreport{Chung:2015-FormulasAngularMomentumBarrier,
   title = {Formulas for {{Angular-Momentum Barrier Factors}} ({{Version II}})},
   author = {Chung, Suh-Urk},
   year = 2015,

docs/conf.py

@@ -290,10 +290,12 @@ def _get_scipy_url() -> str:
     "compwa-report": ("https://compwa.github.io/report", None),
     "graphviz": ("https://graphviz.readthedocs.io/en/stable", None),
     "ipywidgets": (f"https://ipywidgets.readthedocs.io/en/{pin('ipywidgets')}", None),
+    "jax": ("https://docs.jax.dev/en/latest", None),
     "matplotlib": (f"https://matplotlib.org/{pin('matplotlib')}", None),
     "numpy": (f"https://numpy.org/doc/{pin_minor('numpy')}", None),
     "python": ("https://docs.python.org/3", None),
     "qrules": (f"https://qrules.readthedocs.io/{pin('qrules')}", None),
+    "quadax": ("https://quadax.readthedocs.io/en/stable", None),
     "scipy": (_get_scipy_url(), None),
     "spb": (
         f"https://sympy-plot-backends.readthedocs.io/en/v{pin('sympy-plot-backends')}",

docs/usage/amplitude.ipynb

@@ -1011,7 +1011,7 @@
     "            continue\n",
     "        if isinstance(node, (float, sp.Float)):\n",
     "            expression = expression.subs(node, round(node, n_decimals))\n",
-    "        if isinstance(node, sp.Pow) and node.args[1] == 1 / 2:\n",
+    "        if isinstance(node, sp.Pow) and node.args[1] == 1 / 2:  # noqa: RUF069\n",
     "            expression = expression.subs(node, round(node.n(), n_decimals))\n",
     "    return expression"
    ]

docs/usage/dynamics.ipynb

@@ -32,6 +32,7 @@
     "dynamics/analytic-continuation\n",
     "dynamics/k-matrix\n",
     "dynamics/riemann-sheets\n",
+    "dynamics/integration-algorithms\n",
     "```\n",
     "\n",
     "```{autolink-skip}\n",

docs/usage/dynamics/analytic-continuation.ipynb

@@ -530,7 +530,7 @@
    "outputs": [],
    "source": [
     "L = sp.Symbol(\"L\", integer=True, nonnegative=True)\n",
-    "integral_expr = ChewMandelstamIntegral(s, m1, m2, L)\n",
+    "integral_expr = ChewMandelstamIntegral(s, m1, m2, L, configuration={\"epsrel\": 1e-3})\n",
     "integral_expr.doit(deep=False)"
    ]
   },
@@ -543,8 +543,7 @@
     "integral_s_wave_func = sp.lambdify(\n",
     "    [s, m1, m2, integral_expr.epsilon],\n",
     "    integral_expr.subs(L, 0).doit(),\n",
-    ")\n",
-    "integral_s_wave_func = np.vectorize(integral_s_wave_func)"
+    ")"
    ]
   },
   {
@@ -566,8 +565,16 @@
     "integral_p_wave_func = sp.lambdify(\n",
     "    [s, m1, m2, integral_expr.epsilon],\n",
     "    integral_expr.subs(L, 1).doit(),\n",
-    ")\n",
-    "integral_p_wave_func = np.vectorize(integral_p_wave_func)"
+    ")"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    ":::{warning}\n",
+    "In some cases, if the `epsabs` or `epsrel` arguments are set too low, {func}`~scipy.integrate.quad_vec` returns an array that [consists entirely of `nan` values](https://github.com/ComPWA/ampform/issues/488#issuecomment-3868515010).\n",
+    ":::"
    ]
   },
   {
@@ -587,6 +594,9 @@
     "s_wave_values = integral_s_wave_func(s_values, m1_val, m2_val, epsilon=1e-5)\n",
     "p_wave_values = integral_p_wave_func(s_values, m1_val, m2_val, epsilon=1e-5)\n",
     "\n",
+    "assert not np.all(np.isnan(s_wave_values)), \"All values NaN\"\n",
+    "assert not np.all(np.isnan(p_wave_values)), \"All values NaN\"\n",
+    "\n",
     "l_val = [0, 1]\n",
     "fig, axes = plt.subplots(figsize=(6, 7), nrows=2, sharex=True)\n",
     "fig.patch.set_facecolor(\"none\")\n",