MAINT: integrate._qmc_quad: mark as private with preceding underscore

scipy/integrate/_quadrature.py

@@ -1049,11 +1049,11 @@ def newton_cotes(rn, equal=0):
 def _qmc_quad_iv(func, a, b, n_points, n_estimates, qrng, log):
 
  # lazy import to avoid issues with partially-initialized submodule
- if not hasattr(qmc_quad, 'qmc'):
+ if not hasattr(_qmc_quad, 'qmc'):
  from scipy import stats
- qmc_quad.stats = stats
+ _qmc_quad.stats = stats
  else:
- stats = qmc_quad.stats
+ stats = _qmc_quad.stats
 
  if not callable(func):
  message = "`func` must be callable."
@@ -1123,8 +1123,8 @@ def vfunc(x):
 QMCQuadResult = namedtuple('QMCQuadResult', ['integral', 'standard_error'])
 
 
-def qmc_quad(func, a, b, *, n_points=1024, n_estimates=8, qrng=None,
- log=False, args=None):
+def _qmc_quad(func, a, b, *, n_points=1024, n_estimates=8, qrng=None,
+  log=False, args=None):
  """
  Compute an integral in N-dimensions using Quasi-Monte Carlo quadrature.
 

scipy/integrate/tests/test_quadrature.py

@@ -7,7 +7,7 @@
 from scipy.integrate import (quadrature, romberg, romb, newton_cotes,
  cumulative_trapezoid, cumtrapz, trapz, trapezoid,
  quad, simpson, simps, fixed_quad, AccuracyWarning)
-from scipy.integrate._quadrature import qmc_quad
+from scipy.integrate._quadrature import _qmc_quad as qmc_quad
 from scipy import stats, special as sc