code: def primes_from_2_to(n: int) -> np.ndarray:
      """Prime numbers from 2 to *n*.
  
      Parameters
      ----------
      n : int
          Sup bound with ``n >= 6``.
  
      Returns
      -------
      primes : list(int)
          Primes in ``2 <= p < n``.
  
      Notes
      -----
      Taken from [1]_ by P.T. Roy, written consent given on 23.04.2021
      by the original author, Bruno Astrolino, for free use in SciPy under
      the 3-clause BSD.
  
      References
      ----------
      .. [1] `StackOverflow <https://stackoverflow.com/questions/2068372>`_.
  
      """
      sieve = np.ones(n // 3 + (n % 6 == 2), dtype=bool)
      for i in range(1, int(n**0.5) // 3 + 1):
          k = 3 * i + 1 | 1
          sieve[k * k // 3 :: 2 * k] = False
          sieve[k * (k - 2 * (i & 1) + 4) // 3 :: 2 * k] = False
      return np.r_[2, 3, ((3 * np.nonzero(sieve)[0][1:] + 1) | 1)]
  
imports: from __future__ import annotations
  import os
  import copy
  import numbers
  from abc import ABC, abstractmethod
  import math
  from typing import ClassVar, List, Optional, overload, TYPE_CHECKING
  import warnings
  import numpy
  import scipy.stats
  from scipy._lib._util import rng_integers
  from scipy.stats._sobol import initialize_v, _cscramble, _fill_p_cumulative, _draw, _fast_forward, _categorize, initialize_direction_numbers, _MAXDIM, _MAXBIT
  from scipy.stats._qmc_cy import _cy_wrapper_centered_discrepancy, _cy_wrapper_wrap_around_discrepancy, _cy_wrapper_mixture_discrepancy, _cy_wrapper_l2_star_discrepancy, _cy_wrapper_update_discrepancy
label_desc: Use image node_id 7 for calling a global function with example usage: primes_from_2_to(n) and returns: np