{ "__type": "IngestedDoc", "__tag": 4010, "_content": { "Notes": { "__type": "Section", "__tag": 4015, "children": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "For valid results, the integral must converge; behavior for divergent integrals is not guaranteed." } ] }, { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Strong", "__tag": 4048, "children": [ { "__type": "Text", "__tag": 4046, "value": "Details of QUADPACK level routines" } ] } ] }, { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "CrossRef", "__tag": 4002, "value": "quad", "reference": { "__type": "LocalRef", "__tag": 4022, "kind": "module", "path": "scipy.integrate._quadpack_py:quad" }, "kind": "module" }, { "__type": "Text", "__tag": 4046, "value": " calls routines from the FORTRAN library QUADPACK. This section provides details on the conditions for each routine to be called and a short description of each routine. For each level of integration, " }, { "__type": "InlineCode", "__tag": 4051, "value": "qagse" }, { "__type": "Text", "__tag": 4046, "value": " is used for finite limits or " }, { "__type": "InlineCode", "__tag": 4051, "value": "qagie" }, { "__type": "Text", "__tag": 4046, "value": " is used if either limit (or both!) are infinite. The following provides a short description from " }, { "__type": "FootnoteReference", "__tag": 4066, "label": "1" }, { "__type": "Text", "__tag": 4046, "value": " for each routine." } ] }, { "__type": "DefList", "__tag": 4033, "children": [ { "__type": "DefListItem", "__tag": 4037, "dt": { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "qagse" } ] }, "dd": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "is an integrator based on globally adaptive interval subdivision in connection with extrapolation, which will eliminate the effects of integrand singularities of several types. The integration is is performed using a 21-point Gauss-Kronrod quadrature within each subinterval." } ] } ] }, { "__type": "DefListItem", "__tag": 4037, "dt": { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "qagie" } ] }, "dd": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "handles integration over infinite intervals. The infinite range is mapped onto a finite interval and subsequently the same strategy as in " }, { "__type": "InlineCode", "__tag": 4051, "value": "QAGS" }, { "__type": "Text", "__tag": 4046, "value": " is applied." } ] } ] } ] }, { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Strong", "__tag": 4048, "children": [ { "__type": "Text", "__tag": 4046, "value": "Array API Standard Support" } ] } ] }, { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "CrossRef", "__tag": 4002, "value": "dblquad", "reference": { "__type": "LocalRef", "__tag": 4022, "kind": "module", "path": "scipy.integrate._quadpack_py:dblquad" }, "kind": "module" }, { "__type": "Text", "__tag": 4046, "value": " has experimental support for Python Array API Standard compatible backends in addition to NumPy. Please consider testing these features by setting an environment variable " }, { "__type": "InlineCode", "__tag": 4051, "value": "SCIPY_ARRAY_API=1" }, { "__type": "Text", "__tag": 4046, "value": " and providing CuPy, PyTorch, JAX, or Dask arrays as array arguments. The following combinations of backend and device (or other capability) are supported." } ] }, { "__type": "Code", "__tag": 4050, "value": "==================== ==================== ====================\nLibrary CPU GPU\n==================== ==================== ====================\nNumPy ✅ n/a \nCuPy n/a ⛔ \nPyTorch ⛔ ⛔ \nJAX ⛔ ⛔ \nDask ⛔ n/a \n==================== ==================== ====================", "execution_status": null }, { "__type": "Blockquote", "__tag": 4059, "children": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "See " }, { "__type": "InlineRole", "__tag": 4003, "value": "dev-arrayapi", "domain": null, "role": "ref", "inventory": null }, { "__type": "Text", "__tag": 4046, "value": " for more information." } ] } ] } ], "title": [], "level": 0, "target": null }, "Warns": { "__type": "Section", "__tag": 4015, "children": [], "title": [], "level": 0, "target": null }, "Raises": { "__type": "Section", "__tag": 4015, "children": [], "title": [], "level": 0, "target": null }, "Yields": { "__type": "Section", "__tag": 4015, "children": [], "title": [], "level": 0, "target": null }, "Methods": { "__type": "Section", "__tag": 4015, "children": [], "title": [], "level": 0, "target": null }, "Returns": { "__type": "Section", "__tag": 4015, "children": [ { "__type": "Parameters", "__tag": 4026, "children": [ { "__type": "DocParam", "__tag": 4016, "name": "y", "annotation": "float", "desc": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "The resultant integral." } ] } ] }, { "__type": "DocParam", "__tag": 4016, "name": "abserr", "annotation": "float", "desc": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "An estimate of the error." } ] } ] } ] } ], "title": [], "level": 0, "target": null }, "Summary": { "__type": "Section", "__tag": 4015, "children": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "Compute a double integral." } ] } ], "title": [], "level": 0, "target": null }, "Receives": { "__type": "Section", "__tag": 4015, "children": [], "title": [], "level": 0, "target": null }, "Warnings": { "__type": "Section", "__tag": 4015, "children": [], "title": [], "level": 0, "target": null }, "Attributes": { "__type": "Section", "__tag": 4015, "children": [], "title": [], "level": 0, "target": null }, "Parameters": { "__type": "Section", "__tag": 4015, "children": [ { "__type": "Parameters", "__tag": 4026, "children": [ { "__type": "DocParam", "__tag": 4016, "name": "func", "annotation": "callable", "desc": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "A Python function or method of at least two variables: y must be the first argument and x the second argument." } ] } ] }, { "__type": "DocParam", "__tag": 4016, "name": "a, b", "annotation": "float", "desc": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "The limits of integration in x: " }, { "__type": "ParamRef", "__tag": 4071, "name": "a" }, { "__type": "Text", "__tag": 4046, "value": " < " }, { "__type": "ParamRef", "__tag": 4071, "name": "b" } ] } ] }, { "__type": "DocParam", "__tag": 4016, "name": "gfun", "annotation": "callable or float", "desc": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "The lower boundary curve in y which is a function taking a single floating point argument (x) and returning a floating point result or a float indicating a constant boundary curve." } ] } ] }, { "__type": "DocParam", "__tag": 4016, "name": "hfun", "annotation": "callable or float", "desc": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "The upper boundary curve in y (same requirements as " }, { "__type": "ParamRef", "__tag": 4071, "name": "gfun" }, { "__type": "Text", "__tag": 4046, "value": ")." } ] } ] }, { "__type": "DocParam", "__tag": 4016, "name": "args", "annotation": "sequence, optional", "desc": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "Extra arguments to pass to " }, { "__type": "ParamRef", "__tag": 4071, "name": "func" }, { "__type": "Text", "__tag": 4046, "value": "." } ] } ] }, { "__type": "DocParam", "__tag": 4016, "name": "epsabs", "annotation": "float, optional", "desc": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "Absolute tolerance passed directly to the inner 1-D quadrature integration. Default is 1.49e-8. " }, { "__type": "InlineCode", "__tag": 4051, "value": "dblquad" }, { "__type": "Text", "__tag": 4046, "value": " tries to obtain an accuracy of " }, { "__type": "InlineCode", "__tag": 4051, "value": "abs(i-result) <= max(epsabs, epsrel*abs(i))" }, { "__type": "Text", "__tag": 4046, "value": " where " }, { "__type": "InlineCode", "__tag": 4051, "value": "i" }, { "__type": "Text", "__tag": 4046, "value": " = inner integral of " }, { "__type": "InlineCode", "__tag": 4051, "value": "func(y, x)" }, { "__type": "Text", "__tag": 4046, "value": " from " }, { "__type": "InlineCode", "__tag": 4051, "value": "gfun(x)" }, { "__type": "Text", "__tag": 4046, "value": " to " }, { "__type": "InlineCode", "__tag": 4051, "value": "hfun(x)" }, { "__type": "Text", "__tag": 4046, "value": ", and " }, { "__type": "InlineCode", "__tag": 4051, "value": "result" }, { "__type": "Text", "__tag": 4046, "value": " is the numerical approximation. See " }, { "__type": "ParamRef", "__tag": 4071, "name": "epsrel" }, { "__type": "Text", "__tag": 4046, "value": " below." } ] } ] }, { "__type": "DocParam", "__tag": 4016, "name": "epsrel", "annotation": "float, optional", "desc": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "Relative tolerance of the inner 1-D integrals. Default is 1.49e-8. If " }, { "__type": "InlineCode", "__tag": 4051, "value": "epsabs <= 0" }, { "__type": "Text", "__tag": 4046, "value": ", " }, { "__type": "ParamRef", "__tag": 4071, "name": "epsrel" }, { "__type": "Text", "__tag": 4046, "value": " must be greater than both 5e-29 and " }, { "__type": "InlineCode", "__tag": 4051, "value": "50 * (machine epsilon)" }, { "__type": "Text", "__tag": 4046, "value": ". See " }, { "__type": "ParamRef", "__tag": 4071, "name": "epsabs" }, { "__type": "Text", "__tag": 4046, "value": " above." } ] } ] } ] } ], "title": [], "level": 0, "target": null }, "Extended Summary": { "__type": "Section", "__tag": 4015, "children": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "Return the double (definite) integral of " }, { "__type": "InlineCode", "__tag": 4051, "value": "func(y, x)" }, { "__type": "Text", "__tag": 4046, "value": " from " }, { "__type": "InlineCode", "__tag": 4051, "value": "x = a..b" }, { "__type": "Text", "__tag": 4046, "value": " and " }, { "__type": "InlineCode", "__tag": 4051, "value": "y = gfun(x)..hfun(x)" }, { "__type": "Text", "__tag": 4046, "value": "." } ] } ], "title": [], "level": 0, "target": null }, "Other Parameters": { "__type": "Section", "__tag": 4015, "children": [], "title": [], "level": 0, "target": null } }, "_ordered_sections": [ "Summary", "Extended Summary", "Parameters", "Attributes", "Methods", "Returns", "Yields", "Receives", "Other Parameters", "Raises", "Warns", "Warnings", "Notes" ], "item_file": "/scipy/integrate/_quadpack_py.py", "item_line": 698, "item_type": "function", "aliases": [ "scipy.integrate.dblquad" ], "example_section_data": { "__type": "Section", "__tag": 4015, "children": [ { "__type": "Text", "__tag": 4046, "value": "Compute the double integral of ``x * y**2`` over the box\n``x`` ranging from 0 to 2 and ``y`` ranging from 0 to 1.\nThat is, :math:`\\int^{x=2}_{x=0} \\int^{y=1}_{y=0} x y^2 \\,dy \\,dx`.\n\n" }, { "__type": "Code", "__tag": 4050, "value": "import numpy as np\nfrom scipy import integrate\nf = lambda y, x: x*y**2\n", "execution_status": "success" }, { "__type": "Code", "__tag": 4050, "value": "integrate.dblquad(f, 0, 2, 0, 1)\n", "execution_status": "failure" }, { "__type": "Text", "__tag": 4046, "value": "\nCalculate :math:`\\int^{x=\\pi/4}_{x=0} \\int^{y=\\cos(x)}_{y=\\sin(x)} 1\n\\,dy \\,dx`.\n\n" }, { "__type": "Code", "__tag": 4050, "value": "f = lambda y, x: 1\n", "execution_status": "success" }, { "__type": "Code", "__tag": 4050, "value": "integrate.dblquad(f, 0, np.pi/4, np.sin, np.cos)\n", "execution_status": "failure" }, { "__type": "Text", "__tag": 4046, "value": "\nCalculate :math:`\\int^{x=1}_{x=0} \\int^{y=2-x}_{y=x} a x y \\,dy \\,dx`\nfor :math:`a=1, 3`.\n\n" }, { "__type": "Code", "__tag": 4050, "value": "f = lambda y, x, a: a*x*y\n", "execution_status": "success" }, { "__type": "Code", "__tag": 4050, "value": "integrate.dblquad(f, 0, 1, lambda x: x, lambda x: 2-x, args=(1,))\nintegrate.dblquad(f, 0, 1, lambda x: x, lambda x: 2-x, args=(3,))\n", "execution_status": "failure" }, { "__type": "Text", "__tag": 4046, "value": "\nCompute the two-dimensional Gaussian Integral, which is the integral of the\nGaussian function :math:`f(x,y) = e^{-(x^{2} + y^{2})}`, over\n:math:`(-\\infty,+\\infty)`. That is, compute the integral\n:math:`\\iint^{+\\infty}_{-\\infty} e^{-(x^{2} + y^{2})} \\,dy\\,dx`.\n\n" }, { "__type": "Code", "__tag": 4050, "value": "f = lambda x, y: np.exp(-(x ** 2 + y ** 2))\n", "execution_status": "success" }, { "__type": "Code", "__tag": 4050, "value": "integrate.dblquad(f, -np.inf, np.inf, -np.inf, np.inf)\n", "execution_status": "failure" } ], "title": [], "level": 0, "target": null }, "see_also": [ { "__type": "SeeAlsoItem", "__tag": 4028, "name": { "__type": "CrossRef", "__tag": 4002, "value": "fixed_quad", "reference": { "__type": "LocalRef", "__tag": 4022, "kind": "module", "path": "scipy.integrate._quadrature:fixed_quad" }, "kind": "module" }, "descriptions": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "fixed-order Gaussian quadrature" } ] } ], "type": "func" }, { "__type": "SeeAlsoItem", "__tag": 4028, "name": { "__type": "CrossRef", "__tag": 4002, "value": "nquad", "reference": { "__type": "LocalRef", "__tag": 4022, "kind": "module", "path": "scipy.integrate._quadpack_py:nquad" }, "kind": "module" }, "descriptions": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "N-dimensional integrals" } ] } ], "type": "func" }, { "__type": "SeeAlsoItem", "__tag": 4028, "name": { "__type": "CrossRef", "__tag": 4002, "value": "quad", "reference": { "__type": "LocalRef", "__tag": 4022, "kind": "module", "path": "scipy.integrate._quadpack_py:quad" }, "kind": "module" }, "descriptions": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "single integral" } ] } ], "type": "func" }, { "__type": "SeeAlsoItem", "__tag": 4028, "name": { "__type": "CrossRef", "__tag": 4002, "value": "romb", "reference": { "__type": "LocalRef", "__tag": 4022, "kind": "module", "path": "scipy.integrate._quadrature:romb" }, "kind": "module" }, "descriptions": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "integrator for sampled data" } ] } ], "type": "func" }, { "__type": "SeeAlsoItem", "__tag": 4028, "name": { "__type": "CrossRef", "__tag": 4002, "value": "scipy.special", "reference": { "__type": "LocalRef", "__tag": 4022, "kind": "module", "path": "scipy.special" }, "kind": "module" }, "descriptions": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "for coefficients and roots of orthogonal polynomials" } ] } ], "type": "func" }, { "__type": "SeeAlsoItem", "__tag": 4028, "name": { "__type": "CrossRef", "__tag": 4002, "value": "simpson", "reference": { "__type": "LocalRef", "__tag": 4022, "kind": "module", "path": "scipy.integrate._quadrature:simpson" }, "kind": "module" }, "descriptions": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "integrator for sampled data" } ] } ], "type": "func" }, { "__type": "SeeAlsoItem", "__tag": 4028, "name": { "__type": "CrossRef", "__tag": 4002, "value": "tplquad", "reference": { "__type": "LocalRef", "__tag": 4022, "kind": "module", "path": "scipy.integrate._quadpack_py:tplquad" }, "kind": "module" }, "descriptions": [ { "__type": "Paragraph", "__tag": 4045, "children": [ { "__type": "Text", "__tag": 4046, "value": "triple integral" } ] } ], "type": "func" } ], "signature": { "__type": "SignatureNode", "__tag": 4029, "kind": "function", "parameters": [ { "__type": "SigParam", "__tag": 4030, "name": "func", "annotation": { "__type": "Empty", "__tag": 4031 }, "kind": "POSITIONAL_OR_KEYWORD", "default": { "__type": "Empty", "__tag": 4031 } }, { "__type": "SigParam", "__tag": 4030, "name": "a", "annotation": { "__type": "Empty", "__tag": 4031 }, "kind": "POSITIONAL_OR_KEYWORD", "default": { "__type": "Empty", "__tag": 4031 } }, { "__type": "SigParam", "__tag": 4030, "name": "b", "annotation": { "__type": "Empty", "__tag": 4031 }, "kind": "POSITIONAL_OR_KEYWORD", "default": { "__type": "Empty", "__tag": 4031 } }, { "__type": "SigParam", "__tag": 4030, "name": "gfun", "annotation": { "__type": "Empty", "__tag": 4031 }, "kind": "POSITIONAL_OR_KEYWORD", "default": { "__type": "Empty", "__tag": 4031 } }, { "__type": "SigParam", "__tag": 4030, "name": "hfun", "annotation": { "__type": "Empty", "__tag": 4031 }, "kind": "POSITIONAL_OR_KEYWORD", "default": { "__type": "Empty", "__tag": 4031 } }, { "__type": "SigParam", "__tag": 4030, "name": "args", "annotation": { "__type": "Empty", "__tag": 4031 }, "kind": "POSITIONAL_OR_KEYWORD", "default": "()" }, { "__type": "SigParam", "__tag": 4030, "name": "epsabs", "annotation": { "__type": "Empty", "__tag": 4031 }, "kind": "POSITIONAL_OR_KEYWORD", "default": "1.49e-08" }, { "__type": "SigParam", "__tag": 4030, "name": "epsrel", "annotation": { "__type": "Empty", "__tag": 4031 }, "kind": "POSITIONAL_OR_KEYWORD", "default": "1.49e-08" } ], "return_annotation": { "__type": "Empty", "__tag": 4031 }, "target_name": "dblquad" }, "references": [ ".. [1] Piessens, Robert; de Doncker-Kapenga, Elise;", " Überhuber, Christoph W.; Kahaner, David (1983).", " QUADPACK: A subroutine package for automatic integration.", " Springer-Verlag.", " ISBN 978-3-540-12553-2." ], "qa": "scipy.integrate._quadpack_py:dblquad", "arbitrary": [], "local_refs": [ "a", "abserr", "args", "b", "epsabs", "epsrel", "func", "gfun", "hfun", "y" ] }