Generated by Cython 3.0.11

Yellow lines hint at Python interaction.
Click on a line that starts with a "+" to see the C code that Cython generated for it.

Raw output: fragments_helpers.c

+001: #define NPY_NO_DEPRECATED_API NPY_1_7_API_VERSION
  __pyx_t_7 = __Pyx_PyDict_NewPresized(0); if (unlikely(!__pyx_t_7)) __PYX_ERR(0, 1, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_7);
  if (PyDict_SetItem(__pyx_d, __pyx_n_s_test, __pyx_t_7) < 0) __PYX_ERR(0, 1, __pyx_L1_error)
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
 002: #cython: language_level=3
 003: cimport cython
+004: import numpy as np
  __pyx_t_7 = __Pyx_ImportDottedModule(__pyx_n_s_numpy, NULL); if (unlikely(!__pyx_t_7)) __PYX_ERR(0, 4, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_7);
  if (PyDict_SetItem(__pyx_d, __pyx_n_s_np, __pyx_t_7) < 0) __PYX_ERR(0, 4, __pyx_L1_error)
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
 005: cimport numpy as np
+006: np.import_array()
  __pyx_t_9 = __pyx_f_5numpy_import_array(); if (unlikely(__pyx_t_9 == ((int)-1))) __PYX_ERR(0, 6, __pyx_L1_error)
 007: 
+008: INT64 = np.int64
  __Pyx_GetModuleGlobalName(__pyx_t_7, __pyx_n_s_np); if (unlikely(!__pyx_t_7)) __PYX_ERR(0, 8, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_7);
  __pyx_t_4 = __Pyx_PyObject_GetAttrStr(__pyx_t_7, __pyx_n_s_int64); if (unlikely(!__pyx_t_4)) __PYX_ERR(0, 8, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_4);
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
  if (PyDict_SetItem(__pyx_d, __pyx_n_s_INT64, __pyx_t_4) < 0) __PYX_ERR(0, 8, __pyx_L1_error)
  __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
+009: ctypedef np.int64_t INT64_t
typedef __pyx_t_5numpy_int64_t __pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t;
+010: FLOAT64 = np.float64
  __Pyx_GetModuleGlobalName(__pyx_t_4, __pyx_n_s_np); if (unlikely(!__pyx_t_4)) __PYX_ERR(0, 10, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_4);
  __pyx_t_7 = __Pyx_PyObject_GetAttrStr(__pyx_t_4, __pyx_n_s_float64); if (unlikely(!__pyx_t_7)) __PYX_ERR(0, 10, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_7);
  __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
  if (PyDict_SetItem(__pyx_d, __pyx_n_s_FLOAT64, __pyx_t_7) < 0) __PYX_ERR(0, 10, __pyx_L1_error)
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
 011: ctypedef np.float64_t FLOAT64_t
 012: 
+013: @cython.boundscheck(False) # turn off bounds-checking for entire function
/* Python wrapper */
static PyObject *__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_1extract_fragments(PyObject *__pyx_self, 
#if CYTHON_METH_FASTCALL
PyObject *const *__pyx_args, Py_ssize_t __pyx_nargs, PyObject *__pyx_kwds
#else
PyObject *__pyx_args, PyObject *__pyx_kwds
#endif
); /*proto*/
static PyMethodDef __pyx_mdef_11chromatinhd_7loaders_17fragments_helpers_1extract_fragments = {"extract_fragments", (PyCFunction)(void*)(__Pyx_PyCFunction_FastCallWithKeywords)__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_1extract_fragments, __Pyx_METH_FASTCALL|METH_KEYWORDS, 0};
static PyObject *__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_1extract_fragments(PyObject *__pyx_self, 
#if CYTHON_METH_FASTCALL
PyObject *const *__pyx_args, Py_ssize_t __pyx_nargs, PyObject *__pyx_kwds
#else
PyObject *__pyx_args, PyObject *__pyx_kwds
#endif
) {
  __Pyx_memviewslice __pyx_v_cellxgene_oi = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_cellxgene_indptr = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_coordinates = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_genemapping = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_out_coordinates = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_out_genemapping = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_out_local_cellxgene_ix = { 0, 0, { 0 }, { 0 }, { 0 } };
  #if !CYTHON_METH_FASTCALL
  CYTHON_UNUSED Py_ssize_t __pyx_nargs;
  #endif
  CYTHON_UNUSED PyObject *const *__pyx_kwvalues;
  PyObject *__pyx_r = 0;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("extract_fragments (wrapper)", 0);
  #if !CYTHON_METH_FASTCALL
  #if CYTHON_ASSUME_SAFE_MACROS
  __pyx_nargs = PyTuple_GET_SIZE(__pyx_args);
  #else
  __pyx_nargs = PyTuple_Size(__pyx_args); if (unlikely(__pyx_nargs < 0)) return NULL;
  #endif
  #endif
  __pyx_kwvalues = __Pyx_KwValues_FASTCALL(__pyx_args, __pyx_nargs);
  {
    PyObject **__pyx_pyargnames[] = {&__pyx_n_s_cellxgene_oi,&__pyx_n_s_cellxgene_indptr,&__pyx_n_s_coordinates,&__pyx_n_s_genemapping,&__pyx_n_s_out_coordinates,&__pyx_n_s_out_genemapping,&__pyx_n_s_out_local_cellxgene_ix,0};
  PyObject* values[7] = {0,0,0,0,0,0,0};
    if (__pyx_kwds) {
      Py_ssize_t kw_args;
      switch (__pyx_nargs) {
        case  7: values[6] = __Pyx_Arg_FASTCALL(__pyx_args, 6);
        CYTHON_FALLTHROUGH;
        case  6: values[5] = __Pyx_Arg_FASTCALL(__pyx_args, 5);
        CYTHON_FALLTHROUGH;
        case  5: values[4] = __Pyx_Arg_FASTCALL(__pyx_args, 4);
        CYTHON_FALLTHROUGH;
        case  4: values[3] = __Pyx_Arg_FASTCALL(__pyx_args, 3);
        CYTHON_FALLTHROUGH;
        case  3: values[2] = __Pyx_Arg_FASTCALL(__pyx_args, 2);
        CYTHON_FALLTHROUGH;
        case  2: values[1] = __Pyx_Arg_FASTCALL(__pyx_args, 1);
        CYTHON_FALLTHROUGH;
        case  1: values[0] = __Pyx_Arg_FASTCALL(__pyx_args, 0);
        CYTHON_FALLTHROUGH;
        case  0: break;
        default: goto __pyx_L5_argtuple_error;
      }
      kw_args = __Pyx_NumKwargs_FASTCALL(__pyx_kwds);
      switch (__pyx_nargs) {
        case  0:
        if (likely((values[0] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_cellxgene_oi)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[0]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 13, __pyx_L3_error)
        else goto __pyx_L5_argtuple_error;
        CYTHON_FALLTHROUGH;
        case  1:
        if (likely((values[1] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_cellxgene_indptr)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[1]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 13, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("extract_fragments", 1, 7, 7, 1); __PYX_ERR(0, 13, __pyx_L3_error)
        }
        CYTHON_FALLTHROUGH;
        case  2:
        if (likely((values[2] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_coordinates)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[2]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 13, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("extract_fragments", 1, 7, 7, 2); __PYX_ERR(0, 13, __pyx_L3_error)
        }
        CYTHON_FALLTHROUGH;
        case  3:
        if (likely((values[3] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_genemapping)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[3]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 13, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("extract_fragments", 1, 7, 7, 3); __PYX_ERR(0, 13, __pyx_L3_error)
        }
        CYTHON_FALLTHROUGH;
        case  4:
        if (likely((values[4] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_out_coordinates)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[4]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 13, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("extract_fragments", 1, 7, 7, 4); __PYX_ERR(0, 13, __pyx_L3_error)
        }
        CYTHON_FALLTHROUGH;
        case  5:
        if (likely((values[5] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_out_genemapping)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[5]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 13, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("extract_fragments", 1, 7, 7, 5); __PYX_ERR(0, 13, __pyx_L3_error)
        }
        CYTHON_FALLTHROUGH;
        case  6:
        if (likely((values[6] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_out_local_cellxgene_ix)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[6]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 13, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("extract_fragments", 1, 7, 7, 6); __PYX_ERR(0, 13, __pyx_L3_error)
        }
      }
      if (unlikely(kw_args > 0)) {
        const Py_ssize_t kwd_pos_args = __pyx_nargs;
        if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_kwvalues, __pyx_pyargnames, 0, values + 0, kwd_pos_args, "extract_fragments") < 0)) __PYX_ERR(0, 13, __pyx_L3_error)
      }
    } else if (unlikely(__pyx_nargs != 7)) {
      goto __pyx_L5_argtuple_error;
    } else {
      values[0] = __Pyx_Arg_FASTCALL(__pyx_args, 0);
      values[1] = __Pyx_Arg_FASTCALL(__pyx_args, 1);
      values[2] = __Pyx_Arg_FASTCALL(__pyx_args, 2);
      values[3] = __Pyx_Arg_FASTCALL(__pyx_args, 3);
      values[4] = __Pyx_Arg_FASTCALL(__pyx_args, 4);
      values[5] = __Pyx_Arg_FASTCALL(__pyx_args, 5);
      values[6] = __Pyx_Arg_FASTCALL(__pyx_args, 6);
    }
    __pyx_v_cellxgene_oi = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[0], PyBUF_WRITABLE); if (unlikely(!__pyx_v_cellxgene_oi.memview)) __PYX_ERR(0, 17, __pyx_L3_error)
    __pyx_v_cellxgene_indptr = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[1], PyBUF_WRITABLE); if (unlikely(!__pyx_v_cellxgene_indptr.memview)) __PYX_ERR(0, 18, __pyx_L3_error)
    __pyx_v_coordinates = __Pyx_PyObject_to_MemoryviewSlice_d_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[2], PyBUF_WRITABLE); if (unlikely(!__pyx_v_coordinates.memview)) __PYX_ERR(0, 19, __pyx_L3_error)
    __pyx_v_genemapping = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[3], PyBUF_WRITABLE); if (unlikely(!__pyx_v_genemapping.memview)) __PYX_ERR(0, 20, __pyx_L3_error)
    __pyx_v_out_coordinates = __Pyx_PyObject_to_MemoryviewSlice_d_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[4], PyBUF_WRITABLE); if (unlikely(!__pyx_v_out_coordinates.memview)) __PYX_ERR(0, 21, __pyx_L3_error)
    __pyx_v_out_genemapping = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[5], PyBUF_WRITABLE); if (unlikely(!__pyx_v_out_genemapping.memview)) __PYX_ERR(0, 22, __pyx_L3_error)
    __pyx_v_out_local_cellxgene_ix = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[6], PyBUF_WRITABLE); if (unlikely(!__pyx_v_out_local_cellxgene_ix.memview)) __PYX_ERR(0, 23, __pyx_L3_error)
  }
  goto __pyx_L6_skip;
  __pyx_L5_argtuple_error:;
  __Pyx_RaiseArgtupleInvalid("extract_fragments", 1, 7, 7, __pyx_nargs); __PYX_ERR(0, 13, __pyx_L3_error)
  __pyx_L6_skip:;
  goto __pyx_L4_argument_unpacking_done;
  __pyx_L3_error:;
  {
    Py_ssize_t __pyx_temp;
    for (__pyx_temp=0; __pyx_temp < (Py_ssize_t)(sizeof(values)/sizeof(values[0])); ++__pyx_temp) {
      __Pyx_Arg_XDECREF_FASTCALL(values[__pyx_temp]);
    }
  }
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_cellxgene_oi, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_cellxgene_indptr, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_coordinates, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_genemapping, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_out_coordinates, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_out_genemapping, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_out_local_cellxgene_ix, 1);
  __Pyx_AddTraceback("chromatinhd.loaders.fragments_helpers.extract_fragments", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __Pyx_RefNannyFinishContext();
  return NULL;
  __pyx_L4_argument_unpacking_done:;
  __pyx_r = __pyx_pf_11chromatinhd_7loaders_17fragments_helpers_extract_fragments(__pyx_self, __pyx_v_cellxgene_oi, __pyx_v_cellxgene_indptr, __pyx_v_coordinates, __pyx_v_genemapping, __pyx_v_out_coordinates, __pyx_v_out_genemapping, __pyx_v_out_local_cellxgene_ix);
  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;
  int __pyx_clineno = 0;

  /* function exit code */
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_cellxgene_oi, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_cellxgene_indptr, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_coordinates, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_genemapping, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_out_coordinates, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_out_genemapping, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_out_local_cellxgene_ix, 1);
  {
    Py_ssize_t __pyx_temp;
    for (__pyx_temp=0; __pyx_temp < (Py_ssize_t)(sizeof(values)/sizeof(values[0])); ++__pyx_temp) {
      __Pyx_Arg_XDECREF_FASTCALL(values[__pyx_temp]);
    }
  }
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

static PyObject *__pyx_pf_11chromatinhd_7loaders_17fragments_helpers_extract_fragments(CYTHON_UNUSED PyObject *__pyx_self, __Pyx_memviewslice __pyx_v_cellxgene_oi, __Pyx_memviewslice __pyx_v_cellxgene_indptr, __Pyx_memviewslice __pyx_v_coordinates, __Pyx_memviewslice __pyx_v_genemapping, __Pyx_memviewslice __pyx_v_out_coordinates, __Pyx_memviewslice __pyx_v_out_genemapping, __Pyx_memviewslice __pyx_v_out_local_cellxgene_ix) {
  __pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t __pyx_v_out_ix;
  __pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t __pyx_v_local_cellxgene_ix;
  __pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t __pyx_v_cellxgene_ix;
  __pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t __pyx_v_position;
  PyObject *__pyx_r = NULL;
/* … */
  /* function exit code */
  __pyx_L1_error:;
  __PYX_XCLEAR_MEMVIEW(&__pyx_t_8, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_t_9, 1);
  __Pyx_XDECREF(__pyx_t_12);
  __Pyx_AddTraceback("chromatinhd.loaders.fragments_helpers.extract_fragments", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = NULL;
  __pyx_L0:;
  __Pyx_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}
/* … */
  __pyx_tuple__22 = PyTuple_Pack(11, __pyx_n_s_cellxgene_oi, __pyx_n_s_cellxgene_indptr, __pyx_n_s_coordinates, __pyx_n_s_genemapping, __pyx_n_s_out_coordinates, __pyx_n_s_out_genemapping, __pyx_n_s_out_local_cellxgene_ix, __pyx_n_s_out_ix, __pyx_n_s_local_cellxgene_ix, __pyx_n_s_cellxgene_ix, __pyx_n_s_position); if (unlikely(!__pyx_tuple__22)) __PYX_ERR(0, 13, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_tuple__22);
  __Pyx_GIVEREF(__pyx_tuple__22);
/* … */
  __pyx_t_7 = __Pyx_CyFunction_New(&__pyx_mdef_11chromatinhd_7loaders_17fragments_helpers_1extract_fragments, 0, __pyx_n_s_extract_fragments, NULL, __pyx_n_s_chromatinhd_loaders_fragments_he, __pyx_d, ((PyObject *)__pyx_codeobj__23)); if (unlikely(!__pyx_t_7)) __PYX_ERR(0, 13, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_7);
  if (PyDict_SetItem(__pyx_d, __pyx_n_s_extract_fragments, __pyx_t_7) < 0) __PYX_ERR(0, 13, __pyx_L1_error)
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
  __pyx_codeobj__23 = (PyObject*)__Pyx_PyCode_New(7, 0, 0, 11, 0, CO_OPTIMIZED|CO_NEWLOCALS, __pyx_empty_bytes, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_tuple__22, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_kp_s_src_chromatinhd_loaders_fragment, __pyx_n_s_extract_fragments, 13, __pyx_empty_bytes); if (unlikely(!__pyx_codeobj__23)) __PYX_ERR(0, 13, __pyx_L1_error)
 014: @cython.wraparound(False)  # turn off negative index wrapping for entire function
 015: @cython.cdivision
 016: def extract_fragments(
 017:     INT64_t [::1] cellxgene_oi,
 018:     INT64_t [::1] cellxgene_indptr,
 019:     INT64_t [:,::1] coordinates,
 020:     INT64_t [::1] genemapping,
 021:     INT64_t [:,::1] out_coordinates,
 022:     INT64_t [::1] out_genemapping,
 023:     INT64_t [::1] out_local_cellxgene_ix,
 024: ):
 025:     cdef INT64_t out_ix, local_cellxgene_ix, cellxgene_ix, position
+026:     out_ix = 0 # will store where in the output array we are currently
  __pyx_v_out_ix = 0;
+027:     local_cellxgene_ix = 0 # will store the current fragment counting from 0
  __pyx_v_local_cellxgene_ix = 0;
 028: 
+029:     with nogil:
  {
      #ifdef WITH_THREAD
      PyThreadState *_save;
      _save = NULL;
      Py_UNBLOCK_THREADS
      __Pyx_FastGIL_Remember();
      #endif
      /*try:*/ {
/* … */
      /*finally:*/ {
        /*normal exit:*/{
          #ifdef WITH_THREAD
          __Pyx_FastGIL_Forget();
          Py_BLOCK_THREADS
          #endif
          goto __pyx_L5;
        }
        __pyx_L4_error: {
          #ifdef WITH_THREAD
          __Pyx_FastGIL_Forget();
          Py_BLOCK_THREADS
          #endif
          goto __pyx_L1_error;
        }
        __pyx_L5:;
      }
  }
+030:         for local_cellxgene_ix in range(cellxgene_oi.shape[0]):
        __pyx_t_1 = (__pyx_v_cellxgene_oi.shape[0]);
        __pyx_t_2 = __pyx_t_1;
        for (__pyx_t_3 = 0; __pyx_t_3 < __pyx_t_2; __pyx_t_3+=1) {
          __pyx_v_local_cellxgene_ix = __pyx_t_3;
+031:             cellxgene_ix = cellxgene_oi[local_cellxgene_ix]
          __pyx_t_4 = __pyx_v_local_cellxgene_ix;
          __pyx_v_cellxgene_ix = (*((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_cellxgene_oi.data) + __pyx_t_4)) )));
+032:             for position in range(cellxgene_indptr[cellxgene_ix], cellxgene_indptr[cellxgene_ix+1]):
          __pyx_t_4 = (__pyx_v_cellxgene_ix + 1);
          __pyx_t_5 = (*((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_cellxgene_indptr.data) + __pyx_t_4)) )));
          __pyx_t_4 = __pyx_v_cellxgene_ix;
          __pyx_t_6 = __pyx_t_5;
          for (__pyx_t_7 = (*((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_cellxgene_indptr.data) + __pyx_t_4)) ))); __pyx_t_7 < __pyx_t_6; __pyx_t_7+=1) {
            __pyx_v_position = __pyx_t_7;
+033:                 out_coordinates[out_ix] = coordinates[position]
            __pyx_t_8.data = __pyx_v_coordinates.data;
            __pyx_t_8.memview = __pyx_v_coordinates.memview;
            __PYX_INC_MEMVIEW(&__pyx_t_8, 0);
            {
    Py_ssize_t __pyx_tmp_idx = __pyx_v_position;
    Py_ssize_t __pyx_tmp_stride = __pyx_v_coordinates.strides[0];
        __pyx_t_8.data += __pyx_tmp_idx * __pyx_tmp_stride;
}

__pyx_t_8.shape[0] = __pyx_v_coordinates.shape[1];
__pyx_t_8.strides[0] = __pyx_v_coordinates.strides[1];
    __pyx_t_8.suboffsets[0] = -1;

__pyx_t_9.data = __pyx_v_out_coordinates.data;
            __pyx_t_9.memview = __pyx_v_out_coordinates.memview;
            __PYX_INC_MEMVIEW(&__pyx_t_9, 0);
            {
    Py_ssize_t __pyx_tmp_idx = __pyx_v_out_ix;
    Py_ssize_t __pyx_tmp_stride = __pyx_v_out_coordinates.strides[0];
        __pyx_t_9.data += __pyx_tmp_idx * __pyx_tmp_stride;
}

__pyx_t_9.shape[0] = __pyx_v_out_coordinates.shape[1];
__pyx_t_9.strides[0] = __pyx_v_out_coordinates.strides[1];
    __pyx_t_9.suboffsets[0] = -1;

if (unlikely((__pyx_memoryview_copy_contents(__pyx_t_8, __pyx_t_9, 1, 1, 0) < 0))) __PYX_ERR(0, 33, __pyx_L4_error)
            __PYX_XCLEAR_MEMVIEW(&__pyx_t_9, 0);
            __pyx_t_9.memview = NULL; __pyx_t_9.data = NULL;
            __PYX_XCLEAR_MEMVIEW(&__pyx_t_8, 0);
            __pyx_t_8.memview = NULL; __pyx_t_8.data = NULL;
+034:                 out_genemapping[out_ix] = genemapping[position]
            __pyx_t_10 = __pyx_v_position;
            __pyx_t_11 = __pyx_v_out_ix;
            *((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_out_genemapping.data) + __pyx_t_11)) )) = (*((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_genemapping.data) + __pyx_t_10)) )));
+035:                 out_local_cellxgene_ix[out_ix] = local_cellxgene_ix
            __pyx_t_10 = __pyx_v_out_ix;
            *((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_out_local_cellxgene_ix.data) + __pyx_t_10)) )) = __pyx_v_local_cellxgene_ix;
 036: 
+037:                 out_ix += 1
            __pyx_v_out_ix = (__pyx_v_out_ix + 1);
          }
        }
      }
 038: 
+039:     return out_ix
  __Pyx_XDECREF(__pyx_r);
  __pyx_t_12 = __Pyx_PyInt_From_npy_int64(__pyx_v_out_ix); if (unlikely(!__pyx_t_12)) __PYX_ERR(0, 39, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_12);
  __pyx_r = __pyx_t_12;
  __pyx_t_12 = 0;
  goto __pyx_L0;
 040: 
 041: 
+042: @cython.boundscheck(False) # turn off bounds-checking for entire function
/* Python wrapper */
static PyObject *__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_3multiple_arange(PyObject *__pyx_self, 
#if CYTHON_METH_FASTCALL
PyObject *const *__pyx_args, Py_ssize_t __pyx_nargs, PyObject *__pyx_kwds
#else
PyObject *__pyx_args, PyObject *__pyx_kwds
#endif
); /*proto*/
static PyMethodDef __pyx_mdef_11chromatinhd_7loaders_17fragments_helpers_3multiple_arange = {"multiple_arange", (PyCFunction)(void*)(__Pyx_PyCFunction_FastCallWithKeywords)__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_3multiple_arange, __Pyx_METH_FASTCALL|METH_KEYWORDS, 0};
static PyObject *__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_3multiple_arange(PyObject *__pyx_self, 
#if CYTHON_METH_FASTCALL
PyObject *const *__pyx_args, Py_ssize_t __pyx_nargs, PyObject *__pyx_kwds
#else
PyObject *__pyx_args, PyObject *__pyx_kwds
#endif
) {
  __Pyx_memviewslice __pyx_v_a = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_b = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_ixs = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_local_cellxregion_ix = { 0, 0, { 0 }, { 0 }, { 0 } };
  #if !CYTHON_METH_FASTCALL
  CYTHON_UNUSED Py_ssize_t __pyx_nargs;
  #endif
  CYTHON_UNUSED PyObject *const *__pyx_kwvalues;
  PyObject *__pyx_r = 0;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("multiple_arange (wrapper)", 0);
  #if !CYTHON_METH_FASTCALL
  #if CYTHON_ASSUME_SAFE_MACROS
  __pyx_nargs = PyTuple_GET_SIZE(__pyx_args);
  #else
  __pyx_nargs = PyTuple_Size(__pyx_args); if (unlikely(__pyx_nargs < 0)) return NULL;
  #endif
  #endif
  __pyx_kwvalues = __Pyx_KwValues_FASTCALL(__pyx_args, __pyx_nargs);
  {
    PyObject **__pyx_pyargnames[] = {&__pyx_n_s_a,&__pyx_n_s_b,&__pyx_n_s_ixs,&__pyx_n_s_local_cellxregion_ix,0};
  PyObject* values[4] = {0,0,0,0};
    if (__pyx_kwds) {
      Py_ssize_t kw_args;
      switch (__pyx_nargs) {
        case  4: values[3] = __Pyx_Arg_FASTCALL(__pyx_args, 3);
        CYTHON_FALLTHROUGH;
        case  3: values[2] = __Pyx_Arg_FASTCALL(__pyx_args, 2);
        CYTHON_FALLTHROUGH;
        case  2: values[1] = __Pyx_Arg_FASTCALL(__pyx_args, 1);
        CYTHON_FALLTHROUGH;
        case  1: values[0] = __Pyx_Arg_FASTCALL(__pyx_args, 0);
        CYTHON_FALLTHROUGH;
        case  0: break;
        default: goto __pyx_L5_argtuple_error;
      }
      kw_args = __Pyx_NumKwargs_FASTCALL(__pyx_kwds);
      switch (__pyx_nargs) {
        case  0:
        if (likely((values[0] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_a)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[0]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 42, __pyx_L3_error)
        else goto __pyx_L5_argtuple_error;
        CYTHON_FALLTHROUGH;
        case  1:
        if (likely((values[1] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_b)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[1]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 42, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("multiple_arange", 1, 4, 4, 1); __PYX_ERR(0, 42, __pyx_L3_error)
        }
        CYTHON_FALLTHROUGH;
        case  2:
        if (likely((values[2] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_ixs)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[2]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 42, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("multiple_arange", 1, 4, 4, 2); __PYX_ERR(0, 42, __pyx_L3_error)
        }
        CYTHON_FALLTHROUGH;
        case  3:
        if (likely((values[3] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_local_cellxregion_ix)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[3]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 42, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("multiple_arange", 1, 4, 4, 3); __PYX_ERR(0, 42, __pyx_L3_error)
        }
      }
      if (unlikely(kw_args > 0)) {
        const Py_ssize_t kwd_pos_args = __pyx_nargs;
        if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_kwvalues, __pyx_pyargnames, 0, values + 0, kwd_pos_args, "multiple_arange") < 0)) __PYX_ERR(0, 42, __pyx_L3_error)
      }
    } else if (unlikely(__pyx_nargs != 4)) {
      goto __pyx_L5_argtuple_error;
    } else {
      values[0] = __Pyx_Arg_FASTCALL(__pyx_args, 0);
      values[1] = __Pyx_Arg_FASTCALL(__pyx_args, 1);
      values[2] = __Pyx_Arg_FASTCALL(__pyx_args, 2);
      values[3] = __Pyx_Arg_FASTCALL(__pyx_args, 3);
    }
    __pyx_v_a = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[0], PyBUF_WRITABLE); if (unlikely(!__pyx_v_a.memview)) __PYX_ERR(0, 46, __pyx_L3_error)
    __pyx_v_b = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[1], PyBUF_WRITABLE); if (unlikely(!__pyx_v_b.memview)) __PYX_ERR(0, 47, __pyx_L3_error)
    __pyx_v_ixs = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[2], PyBUF_WRITABLE); if (unlikely(!__pyx_v_ixs.memview)) __PYX_ERR(0, 48, __pyx_L3_error)
    __pyx_v_local_cellxregion_ix = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[3], PyBUF_WRITABLE); if (unlikely(!__pyx_v_local_cellxregion_ix.memview)) __PYX_ERR(0, 49, __pyx_L3_error)
  }
  goto __pyx_L6_skip;
  __pyx_L5_argtuple_error:;
  __Pyx_RaiseArgtupleInvalid("multiple_arange", 1, 4, 4, __pyx_nargs); __PYX_ERR(0, 42, __pyx_L3_error)
  __pyx_L6_skip:;
  goto __pyx_L4_argument_unpacking_done;
  __pyx_L3_error:;
  {
    Py_ssize_t __pyx_temp;
    for (__pyx_temp=0; __pyx_temp < (Py_ssize_t)(sizeof(values)/sizeof(values[0])); ++__pyx_temp) {
      __Pyx_Arg_XDECREF_FASTCALL(values[__pyx_temp]);
    }
  }
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_a, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_b, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_ixs, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_local_cellxregion_ix, 1);
  __Pyx_AddTraceback("chromatinhd.loaders.fragments_helpers.multiple_arange", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __Pyx_RefNannyFinishContext();
  return NULL;
  __pyx_L4_argument_unpacking_done:;
  __pyx_r = __pyx_pf_11chromatinhd_7loaders_17fragments_helpers_2multiple_arange(__pyx_self, __pyx_v_a, __pyx_v_b, __pyx_v_ixs, __pyx_v_local_cellxregion_ix);
  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;
  int __pyx_clineno = 0;

  /* function exit code */
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_a, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_b, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_ixs, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_local_cellxregion_ix, 1);
  {
    Py_ssize_t __pyx_temp;
    for (__pyx_temp=0; __pyx_temp < (Py_ssize_t)(sizeof(values)/sizeof(values[0])); ++__pyx_temp) {
      __Pyx_Arg_XDECREF_FASTCALL(values[__pyx_temp]);
    }
  }
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

static PyObject *__pyx_pf_11chromatinhd_7loaders_17fragments_helpers_2multiple_arange(CYTHON_UNUSED PyObject *__pyx_self, __Pyx_memviewslice __pyx_v_a, __Pyx_memviewslice __pyx_v_b, __Pyx_memviewslice __pyx_v_ixs, __Pyx_memviewslice __pyx_v_local_cellxregion_ix) {
  __pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t __pyx_v_out_ix;
  __pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t __pyx_v_pair_ix;
  __pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t __pyx_v_position;
  PyObject *__pyx_r = NULL;
/* … */
  /* function exit code */
  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_9);
  __Pyx_AddTraceback("chromatinhd.loaders.fragments_helpers.multiple_arange", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = NULL;
  __pyx_L0:;
  __Pyx_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}
/* … */
  __pyx_tuple__24 = PyTuple_Pack(7, __pyx_n_s_a, __pyx_n_s_b, __pyx_n_s_ixs, __pyx_n_s_local_cellxregion_ix, __pyx_n_s_out_ix, __pyx_n_s_pair_ix, __pyx_n_s_position); if (unlikely(!__pyx_tuple__24)) __PYX_ERR(0, 42, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_tuple__24);
  __Pyx_GIVEREF(__pyx_tuple__24);
/* … */
  __pyx_t_7 = __Pyx_CyFunction_New(&__pyx_mdef_11chromatinhd_7loaders_17fragments_helpers_3multiple_arange, 0, __pyx_n_s_multiple_arange, NULL, __pyx_n_s_chromatinhd_loaders_fragments_he, __pyx_d, ((PyObject *)__pyx_codeobj__25)); if (unlikely(!__pyx_t_7)) __PYX_ERR(0, 42, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_7);
  if (PyDict_SetItem(__pyx_d, __pyx_n_s_multiple_arange, __pyx_t_7) < 0) __PYX_ERR(0, 42, __pyx_L1_error)
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
  __pyx_codeobj__25 = (PyObject*)__Pyx_PyCode_New(4, 0, 0, 7, 0, CO_OPTIMIZED|CO_NEWLOCALS, __pyx_empty_bytes, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_tuple__24, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_kp_s_src_chromatinhd_loaders_fragment, __pyx_n_s_multiple_arange, 42, __pyx_empty_bytes); if (unlikely(!__pyx_codeobj__25)) __PYX_ERR(0, 42, __pyx_L1_error)
 043: @cython.wraparound(False)  # turn off negative index wrapping for entire function
 044: @cython.cdivision
 045: def multiple_arange(
 046:     INT64_t [::1] a,
 047:     INT64_t [::1] b,
 048:     INT64_t [::1] ixs,
 049:     INT64_t [::1] local_cellxregion_ix,
 050: ):
 051:     cdef INT64_t out_ix, pair_ix, position
+052:     out_ix = 0 # will store where in the output array we are currently
  __pyx_v_out_ix = 0;
+053:     pair_ix = 0 # will store the current a, b pair index
  __pyx_v_pair_ix = 0;
 054: 
+055:     with nogil:
  {
      #ifdef WITH_THREAD
      PyThreadState *_save;
      _save = NULL;
      Py_UNBLOCK_THREADS
      __Pyx_FastGIL_Remember();
      #endif
      /*try:*/ {
/* … */
      /*finally:*/ {
        /*normal exit:*/{
          #ifdef WITH_THREAD
          __Pyx_FastGIL_Forget();
          Py_BLOCK_THREADS
          #endif
          goto __pyx_L5;
        }
        __pyx_L5:;
      }
  }
+056:         for pair_ix in range(a.shape[0]):
        __pyx_t_1 = (__pyx_v_a.shape[0]);
        __pyx_t_2 = __pyx_t_1;
        for (__pyx_t_3 = 0; __pyx_t_3 < __pyx_t_2; __pyx_t_3+=1) {
          __pyx_v_pair_ix = __pyx_t_3;
+057:             for position in range(a[pair_ix], b[pair_ix]):
          __pyx_t_4 = __pyx_v_pair_ix;
          __pyx_t_5 = (*((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_b.data) + __pyx_t_4)) )));
          __pyx_t_4 = __pyx_v_pair_ix;
          __pyx_t_6 = __pyx_t_5;
          for (__pyx_t_7 = (*((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_a.data) + __pyx_t_4)) ))); __pyx_t_7 < __pyx_t_6; __pyx_t_7+=1) {
            __pyx_v_position = __pyx_t_7;
+058:                 ixs[out_ix] = position
            __pyx_t_8 = __pyx_v_out_ix;
            *((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_ixs.data) + __pyx_t_8)) )) = __pyx_v_position;
+059:                 local_cellxregion_ix[out_ix] = pair_ix
            __pyx_t_8 = __pyx_v_out_ix;
            *((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_local_cellxregion_ix.data) + __pyx_t_8)) )) = __pyx_v_pair_ix;
+060:                 out_ix += 1
            __pyx_v_out_ix = (__pyx_v_out_ix + 1);
          }
+061:             pair_ix += 1
          __pyx_v_pair_ix = (__pyx_v_pair_ix + 1);
        }
      }
 062: 
+063:     return out_ix
  __Pyx_XDECREF(__pyx_r);
  __pyx_t_9 = __Pyx_PyInt_From_npy_int64(__pyx_v_out_ix); if (unlikely(!__pyx_t_9)) __PYX_ERR(0, 63, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_r = __pyx_t_9;
  __pyx_t_9 = 0;
  goto __pyx_L0;
 064: 
 065: 
 066: 
+067: @cython.boundscheck(False) # turn off bounds-checking for entire function
/* Python wrapper */
static PyObject *__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_5multiple_arange(PyObject *__pyx_self, 
#if CYTHON_METH_FASTCALL
PyObject *const *__pyx_args, Py_ssize_t __pyx_nargs, PyObject *__pyx_kwds
#else
PyObject *__pyx_args, PyObject *__pyx_kwds
#endif
); /*proto*/
static PyMethodDef __pyx_mdef_11chromatinhd_7loaders_17fragments_helpers_5multiple_arange = {"multiple_arange", (PyCFunction)(void*)(__Pyx_PyCFunction_FastCallWithKeywords)__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_5multiple_arange, __Pyx_METH_FASTCALL|METH_KEYWORDS, 0};
static PyObject *__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_5multiple_arange(PyObject *__pyx_self, 
#if CYTHON_METH_FASTCALL
PyObject *const *__pyx_args, Py_ssize_t __pyx_nargs, PyObject *__pyx_kwds
#else
PyObject *__pyx_args, PyObject *__pyx_kwds
#endif
) {
  __Pyx_memviewslice __pyx_v_a = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_b = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_ixs = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_local_cellxregion_ix = { 0, 0, { 0 }, { 0 }, { 0 } };
  #if !CYTHON_METH_FASTCALL
  CYTHON_UNUSED Py_ssize_t __pyx_nargs;
  #endif
  CYTHON_UNUSED PyObject *const *__pyx_kwvalues;
  PyObject *__pyx_r = 0;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("multiple_arange (wrapper)", 0);
  #if !CYTHON_METH_FASTCALL
  #if CYTHON_ASSUME_SAFE_MACROS
  __pyx_nargs = PyTuple_GET_SIZE(__pyx_args);
  #else
  __pyx_nargs = PyTuple_Size(__pyx_args); if (unlikely(__pyx_nargs < 0)) return NULL;
  #endif
  #endif
  __pyx_kwvalues = __Pyx_KwValues_FASTCALL(__pyx_args, __pyx_nargs);
  {
    PyObject **__pyx_pyargnames[] = {&__pyx_n_s_a,&__pyx_n_s_b,&__pyx_n_s_ixs,&__pyx_n_s_local_cellxregion_ix,0};
  PyObject* values[4] = {0,0,0,0};
    if (__pyx_kwds) {
      Py_ssize_t kw_args;
      switch (__pyx_nargs) {
        case  4: values[3] = __Pyx_Arg_FASTCALL(__pyx_args, 3);
        CYTHON_FALLTHROUGH;
        case  3: values[2] = __Pyx_Arg_FASTCALL(__pyx_args, 2);
        CYTHON_FALLTHROUGH;
        case  2: values[1] = __Pyx_Arg_FASTCALL(__pyx_args, 1);
        CYTHON_FALLTHROUGH;
        case  1: values[0] = __Pyx_Arg_FASTCALL(__pyx_args, 0);
        CYTHON_FALLTHROUGH;
        case  0: break;
        default: goto __pyx_L5_argtuple_error;
      }
      kw_args = __Pyx_NumKwargs_FASTCALL(__pyx_kwds);
      switch (__pyx_nargs) {
        case  0:
        if (likely((values[0] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_a)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[0]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 67, __pyx_L3_error)
        else goto __pyx_L5_argtuple_error;
        CYTHON_FALLTHROUGH;
        case  1:
        if (likely((values[1] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_b)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[1]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 67, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("multiple_arange", 1, 4, 4, 1); __PYX_ERR(0, 67, __pyx_L3_error)
        }
        CYTHON_FALLTHROUGH;
        case  2:
        if (likely((values[2] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_ixs)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[2]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 67, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("multiple_arange", 1, 4, 4, 2); __PYX_ERR(0, 67, __pyx_L3_error)
        }
        CYTHON_FALLTHROUGH;
        case  3:
        if (likely((values[3] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_local_cellxregion_ix)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[3]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 67, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("multiple_arange", 1, 4, 4, 3); __PYX_ERR(0, 67, __pyx_L3_error)
        }
      }
      if (unlikely(kw_args > 0)) {
        const Py_ssize_t kwd_pos_args = __pyx_nargs;
        if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_kwvalues, __pyx_pyargnames, 0, values + 0, kwd_pos_args, "multiple_arange") < 0)) __PYX_ERR(0, 67, __pyx_L3_error)
      }
    } else if (unlikely(__pyx_nargs != 4)) {
      goto __pyx_L5_argtuple_error;
    } else {
      values[0] = __Pyx_Arg_FASTCALL(__pyx_args, 0);
      values[1] = __Pyx_Arg_FASTCALL(__pyx_args, 1);
      values[2] = __Pyx_Arg_FASTCALL(__pyx_args, 2);
      values[3] = __Pyx_Arg_FASTCALL(__pyx_args, 3);
    }
    __pyx_v_a = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[0], PyBUF_WRITABLE); if (unlikely(!__pyx_v_a.memview)) __PYX_ERR(0, 71, __pyx_L3_error)
    __pyx_v_b = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[1], PyBUF_WRITABLE); if (unlikely(!__pyx_v_b.memview)) __PYX_ERR(0, 72, __pyx_L3_error)
    __pyx_v_ixs = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[2], PyBUF_WRITABLE); if (unlikely(!__pyx_v_ixs.memview)) __PYX_ERR(0, 73, __pyx_L3_error)
    __pyx_v_local_cellxregion_ix = __Pyx_PyObject_to_MemoryviewSlice_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[3], PyBUF_WRITABLE); if (unlikely(!__pyx_v_local_cellxregion_ix.memview)) __PYX_ERR(0, 74, __pyx_L3_error)
  }
  goto __pyx_L6_skip;
  __pyx_L5_argtuple_error:;
  __Pyx_RaiseArgtupleInvalid("multiple_arange", 1, 4, 4, __pyx_nargs); __PYX_ERR(0, 67, __pyx_L3_error)
  __pyx_L6_skip:;
  goto __pyx_L4_argument_unpacking_done;
  __pyx_L3_error:;
  {
    Py_ssize_t __pyx_temp;
    for (__pyx_temp=0; __pyx_temp < (Py_ssize_t)(sizeof(values)/sizeof(values[0])); ++__pyx_temp) {
      __Pyx_Arg_XDECREF_FASTCALL(values[__pyx_temp]);
    }
  }
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_a, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_b, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_ixs, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_local_cellxregion_ix, 1);
  __Pyx_AddTraceback("chromatinhd.loaders.fragments_helpers.multiple_arange", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __Pyx_RefNannyFinishContext();
  return NULL;
  __pyx_L4_argument_unpacking_done:;
  __pyx_r = __pyx_pf_11chromatinhd_7loaders_17fragments_helpers_4multiple_arange(__pyx_self, __pyx_v_a, __pyx_v_b, __pyx_v_ixs, __pyx_v_local_cellxregion_ix);
  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;
  int __pyx_clineno = 0;

  /* function exit code */
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_a, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_b, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_ixs, 1);
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_local_cellxregion_ix, 1);
  {
    Py_ssize_t __pyx_temp;
    for (__pyx_temp=0; __pyx_temp < (Py_ssize_t)(sizeof(values)/sizeof(values[0])); ++__pyx_temp) {
      __Pyx_Arg_XDECREF_FASTCALL(values[__pyx_temp]);
    }
  }
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

static PyObject *__pyx_pf_11chromatinhd_7loaders_17fragments_helpers_4multiple_arange(CYTHON_UNUSED PyObject *__pyx_self, __Pyx_memviewslice __pyx_v_a, __Pyx_memviewslice __pyx_v_b, __Pyx_memviewslice __pyx_v_ixs, __Pyx_memviewslice __pyx_v_local_cellxregion_ix) {
  __pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t __pyx_v_out_ix;
  __pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t __pyx_v_pair_ix;
  __pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t __pyx_v_position;
  PyObject *__pyx_r = NULL;
/* … */
  /* function exit code */
  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_9);
  __Pyx_AddTraceback("chromatinhd.loaders.fragments_helpers.multiple_arange", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = NULL;
  __pyx_L0:;
  __Pyx_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}
/* … */
  __pyx_t_7 = __Pyx_CyFunction_New(&__pyx_mdef_11chromatinhd_7loaders_17fragments_helpers_5multiple_arange, 0, __pyx_n_s_multiple_arange, NULL, __pyx_n_s_chromatinhd_loaders_fragments_he, __pyx_d, ((PyObject *)__pyx_codeobj__26)); if (unlikely(!__pyx_t_7)) __PYX_ERR(0, 67, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_7);
  if (PyDict_SetItem(__pyx_d, __pyx_n_s_multiple_arange, __pyx_t_7) < 0) __PYX_ERR(0, 67, __pyx_L1_error)
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
/* … */
  __pyx_codeobj__26 = (PyObject*)__Pyx_PyCode_New(4, 0, 0, 7, 0, CO_OPTIMIZED|CO_NEWLOCALS, __pyx_empty_bytes, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_tuple__24, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_kp_s_src_chromatinhd_loaders_fragment, __pyx_n_s_multiple_arange, 67, __pyx_empty_bytes); if (unlikely(!__pyx_codeobj__26)) __PYX_ERR(0, 67, __pyx_L1_error)
 068: @cython.wraparound(False)  # turn off negative index wrapping for entire function
 069: @cython.cdivision
 070: def multiple_arange(
 071:     INT64_t [::1] a,
 072:     INT64_t [::1] b,
 073:     INT64_t [::1] ixs,
 074:     INT64_t [::1] local_cellxregion_ix,
 075: ):
 076:     cdef INT64_t out_ix, pair_ix, position
+077:     out_ix = 0 # will store where in the output array we are currently
  __pyx_v_out_ix = 0;
+078:     pair_ix = 0 # will store the current a, b pair index
  __pyx_v_pair_ix = 0;
 079: 
+080:     with nogil:
  {
      #ifdef WITH_THREAD
      PyThreadState *_save;
      _save = NULL;
      Py_UNBLOCK_THREADS
      __Pyx_FastGIL_Remember();
      #endif
      /*try:*/ {
/* … */
      /*finally:*/ {
        /*normal exit:*/{
          #ifdef WITH_THREAD
          __Pyx_FastGIL_Forget();
          Py_BLOCK_THREADS
          #endif
          goto __pyx_L5;
        }
        __pyx_L5:;
      }
  }
+081:         for pair_ix in range(a.shape[0]):
        __pyx_t_1 = (__pyx_v_a.shape[0]);
        __pyx_t_2 = __pyx_t_1;
        for (__pyx_t_3 = 0; __pyx_t_3 < __pyx_t_2; __pyx_t_3+=1) {
          __pyx_v_pair_ix = __pyx_t_3;
+082:             for position in range(a[pair_ix], b[pair_ix]):
          __pyx_t_4 = __pyx_v_pair_ix;
          __pyx_t_5 = (*((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_b.data) + __pyx_t_4)) )));
          __pyx_t_4 = __pyx_v_pair_ix;
          __pyx_t_6 = __pyx_t_5;
          for (__pyx_t_7 = (*((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_a.data) + __pyx_t_4)) ))); __pyx_t_7 < __pyx_t_6; __pyx_t_7+=1) {
            __pyx_v_position = __pyx_t_7;
+083:                 ixs[out_ix] = position
            __pyx_t_8 = __pyx_v_out_ix;
            *((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_ixs.data) + __pyx_t_8)) )) = __pyx_v_position;
+084:                 local_cellxregion_ix[out_ix] = pair_ix
            __pyx_t_8 = __pyx_v_out_ix;
            *((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) __pyx_v_local_cellxregion_ix.data) + __pyx_t_8)) )) = __pyx_v_pair_ix;
+085:                 out_ix += 1
            __pyx_v_out_ix = (__pyx_v_out_ix + 1);
          }
+086:             pair_ix += 1
          __pyx_v_pair_ix = (__pyx_v_pair_ix + 1);
        }
      }
 087: 
+088:     return out_ix
  __Pyx_XDECREF(__pyx_r);
  __pyx_t_9 = __Pyx_PyInt_From_npy_int64(__pyx_v_out_ix); if (unlikely(!__pyx_t_9)) __PYX_ERR(0, 88, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_r = __pyx_t_9;
  __pyx_t_9 = 0;
  goto __pyx_L0;
 089: 
+090: @cython.boundscheck(False) # turn off bounds-checking for entire function
/* Python wrapper */
static PyObject *__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_7count_A_cython(PyObject *__pyx_self, 
#if CYTHON_METH_FASTCALL
PyObject *const *__pyx_args, Py_ssize_t __pyx_nargs, PyObject *__pyx_kwds
#else
PyObject *__pyx_args, PyObject *__pyx_kwds
#endif
); /*proto*/
static PyMethodDef __pyx_mdef_11chromatinhd_7loaders_17fragments_helpers_7count_A_cython = {"count_A_cython", (PyCFunction)(void*)(__Pyx_PyCFunction_FastCallWithKeywords)__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_7count_A_cython, __Pyx_METH_FASTCALL|METH_KEYWORDS, 0};
static PyObject *__pyx_pw_11chromatinhd_7loaders_17fragments_helpers_7count_A_cython(PyObject *__pyx_self, 
#if CYTHON_METH_FASTCALL
PyObject *const *__pyx_args, Py_ssize_t __pyx_nargs, PyObject *__pyx_kwds
#else
PyObject *__pyx_args, PyObject *__pyx_kwds
#endif
) {
  PyObject *__pyx_v_s = 0;
  __Pyx_memviewslice __pyx_v_out_onehot = { 0, 0, { 0 }, { 0 }, { 0 } };
  #if !CYTHON_METH_FASTCALL
  CYTHON_UNUSED Py_ssize_t __pyx_nargs;
  #endif
  CYTHON_UNUSED PyObject *const *__pyx_kwvalues;
  PyObject *__pyx_r = 0;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("count_A_cython (wrapper)", 0);
  #if !CYTHON_METH_FASTCALL
  #if CYTHON_ASSUME_SAFE_MACROS
  __pyx_nargs = PyTuple_GET_SIZE(__pyx_args);
  #else
  __pyx_nargs = PyTuple_Size(__pyx_args); if (unlikely(__pyx_nargs < 0)) return NULL;
  #endif
  #endif
  __pyx_kwvalues = __Pyx_KwValues_FASTCALL(__pyx_args, __pyx_nargs);
  {
    PyObject **__pyx_pyargnames[] = {&__pyx_n_s_s,&__pyx_n_s_out_onehot,0};
  PyObject* values[2] = {0,0};
    if (__pyx_kwds) {
      Py_ssize_t kw_args;
      switch (__pyx_nargs) {
        case  2: values[1] = __Pyx_Arg_FASTCALL(__pyx_args, 1);
        CYTHON_FALLTHROUGH;
        case  1: values[0] = __Pyx_Arg_FASTCALL(__pyx_args, 0);
        CYTHON_FALLTHROUGH;
        case  0: break;
        default: goto __pyx_L5_argtuple_error;
      }
      kw_args = __Pyx_NumKwargs_FASTCALL(__pyx_kwds);
      switch (__pyx_nargs) {
        case  0:
        if (likely((values[0] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_s)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[0]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 90, __pyx_L3_error)
        else goto __pyx_L5_argtuple_error;
        CYTHON_FALLTHROUGH;
        case  1:
        if (likely((values[1] = __Pyx_GetKwValue_FASTCALL(__pyx_kwds, __pyx_kwvalues, __pyx_n_s_out_onehot)) != 0)) {
          (void)__Pyx_Arg_NewRef_FASTCALL(values[1]);
          kw_args--;
        }
        else if (unlikely(PyErr_Occurred())) __PYX_ERR(0, 90, __pyx_L3_error)
        else {
          __Pyx_RaiseArgtupleInvalid("count_A_cython", 1, 2, 2, 1); __PYX_ERR(0, 90, __pyx_L3_error)
        }
      }
      if (unlikely(kw_args > 0)) {
        const Py_ssize_t kwd_pos_args = __pyx_nargs;
        if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_kwvalues, __pyx_pyargnames, 0, values + 0, kwd_pos_args, "count_A_cython") < 0)) __PYX_ERR(0, 90, __pyx_L3_error)
      }
    } else if (unlikely(__pyx_nargs != 2)) {
      goto __pyx_L5_argtuple_error;
    } else {
      values[0] = __Pyx_Arg_FASTCALL(__pyx_args, 0);
      values[1] = __Pyx_Arg_FASTCALL(__pyx_args, 1);
    }
    __pyx_v_s = ((PyObject*)values[0]);
    __pyx_v_out_onehot = __Pyx_PyObject_to_MemoryviewSlice_d_dc_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t(values[1], PyBUF_WRITABLE); if (unlikely(!__pyx_v_out_onehot.memview)) __PYX_ERR(0, 93, __pyx_L3_error)
  }
  goto __pyx_L6_skip;
  __pyx_L5_argtuple_error:;
  __Pyx_RaiseArgtupleInvalid("count_A_cython", 1, 2, 2, __pyx_nargs); __PYX_ERR(0, 90, __pyx_L3_error)
  __pyx_L6_skip:;
  goto __pyx_L4_argument_unpacking_done;
  __pyx_L3_error:;
  {
    Py_ssize_t __pyx_temp;
    for (__pyx_temp=0; __pyx_temp < (Py_ssize_t)(sizeof(values)/sizeof(values[0])); ++__pyx_temp) {
      __Pyx_Arg_XDECREF_FASTCALL(values[__pyx_temp]);
    }
  }
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_out_onehot, 1);
  __Pyx_AddTraceback("chromatinhd.loaders.fragments_helpers.count_A_cython", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __Pyx_RefNannyFinishContext();
  return NULL;
  __pyx_L4_argument_unpacking_done:;
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_s), (&PyBytes_Type), 1, "s", 1))) __PYX_ERR(0, 93, __pyx_L1_error)
  __pyx_r = __pyx_pf_11chromatinhd_7loaders_17fragments_helpers_6count_A_cython(__pyx_self, __pyx_v_s, __pyx_v_out_onehot);
  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;
  int __pyx_clineno = 0;

  /* function exit code */
  goto __pyx_L0;
  __pyx_L1_error:;
  __pyx_r = NULL;
  __pyx_L0:;
  __PYX_XCLEAR_MEMVIEW(&__pyx_v_out_onehot, 1);
  {
    Py_ssize_t __pyx_temp;
    for (__pyx_temp=0; __pyx_temp < (Py_ssize_t)(sizeof(values)/sizeof(values[0])); ++__pyx_temp) {
      __Pyx_Arg_XDECREF_FASTCALL(values[__pyx_temp]);
    }
  }
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

static PyObject *__pyx_pf_11chromatinhd_7loaders_17fragments_helpers_6count_A_cython(CYTHON_UNUSED PyObject *__pyx_self, PyObject *__pyx_v_s, __Pyx_memviewslice __pyx_v_out_onehot) {
  char *__pyx_v_cstr;
  int __pyx_v_i;
  int __pyx_v_length;
  PyObject *__pyx_r = NULL;
/* … */
  /* function exit code */
  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_9);
  __Pyx_AddTraceback("chromatinhd.loaders.fragments_helpers.count_A_cython", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = NULL;
  __pyx_L0:;
  __Pyx_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}
/* … */
  __pyx_tuple__27 = PyTuple_Pack(5, __pyx_n_s_s, __pyx_n_s_out_onehot, __pyx_n_s_cstr, __pyx_n_s_i, __pyx_n_s_length); if (unlikely(!__pyx_tuple__27)) __PYX_ERR(0, 90, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_tuple__27);
  __Pyx_GIVEREF(__pyx_tuple__27);
/* … */
  __pyx_t_7 = __Pyx_CyFunction_New(&__pyx_mdef_11chromatinhd_7loaders_17fragments_helpers_7count_A_cython, 0, __pyx_n_s_count_A_cython, NULL, __pyx_n_s_chromatinhd_loaders_fragments_he, __pyx_d, ((PyObject *)__pyx_codeobj__28)); if (unlikely(!__pyx_t_7)) __PYX_ERR(0, 90, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_7);
  if (PyDict_SetItem(__pyx_d, __pyx_n_s_count_A_cython, __pyx_t_7) < 0) __PYX_ERR(0, 90, __pyx_L1_error)
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
 091: @cython.wraparound(False)  # turn off negative index wrapping for entire function
 092: @cython.cdivision
 093: def count_A_cython(bytes s, INT64_t [:,::1] out_onehot):
+094:     cdef char* cstr = s
  if (unlikely(__pyx_v_s == Py_None)) {
    PyErr_SetString(PyExc_TypeError, "expected bytes, NoneType found");
    __PYX_ERR(0, 94, __pyx_L1_error)
  }
  __pyx_t_1 = __Pyx_PyBytes_AsWritableString(__pyx_v_s); if (unlikely((!__pyx_t_1) && PyErr_Occurred())) __PYX_ERR(0, 94, __pyx_L1_error)
  __pyx_v_cstr = __pyx_t_1;
+095:     cdef int i, length = len(s)
  if (unlikely(__pyx_v_s == Py_None)) {
    PyErr_SetString(PyExc_TypeError, "object of type 'NoneType' has no len()");
    __PYX_ERR(0, 95, __pyx_L1_error)
  }
  __pyx_t_2 = __Pyx_PyBytes_GET_SIZE(__pyx_v_s); if (unlikely(__pyx_t_2 == ((Py_ssize_t)-1))) __PYX_ERR(0, 95, __pyx_L1_error)
  __pyx_v_length = __pyx_t_2;
 096: 
+097:     for i in range(length):
  __pyx_t_3 = __pyx_v_length;
  __pyx_t_4 = __pyx_t_3;
  for (__pyx_t_5 = 0; __pyx_t_5 < __pyx_t_4; __pyx_t_5+=1) {
    __pyx_v_i = __pyx_t_5;
+098:         if cstr[i] == b'A':
    __pyx_t_6 = ((__pyx_v_cstr[__pyx_v_i]) == 'A');
    if (__pyx_t_6) {
/* … */
      goto __pyx_L5;
    }
+099:             out_onehot[i, 0] = 1
      __pyx_t_7 = __pyx_v_i;
      __pyx_t_8 = 0;
      *((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=1 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ (__pyx_v_out_onehot.data + __pyx_t_7 * __pyx_v_out_onehot.strides[0]) )) + __pyx_t_8)) )) = 1;
+100:         elif cstr[i] == b'C':
    __pyx_t_6 = ((__pyx_v_cstr[__pyx_v_i]) == 'C');
    if (__pyx_t_6) {
/* … */
      goto __pyx_L5;
    }
+101:             out_onehot[i, 1] = 1
      __pyx_t_8 = __pyx_v_i;
      __pyx_t_7 = 1;
      *((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=1 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ (__pyx_v_out_onehot.data + __pyx_t_8 * __pyx_v_out_onehot.strides[0]) )) + __pyx_t_7)) )) = 1;
+102:         elif cstr[i] == b'G':
    __pyx_t_6 = ((__pyx_v_cstr[__pyx_v_i]) == 'G');
    if (__pyx_t_6) {
/* … */
      goto __pyx_L5;
    }
+103:             out_onehot[i, 2] = 1
      __pyx_t_7 = __pyx_v_i;
      __pyx_t_8 = 2;
      *((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=1 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ (__pyx_v_out_onehot.data + __pyx_t_7 * __pyx_v_out_onehot.strides[0]) )) + __pyx_t_8)) )) = 1;
+104:         elif cstr[i] == b'T':
    __pyx_t_6 = ((__pyx_v_cstr[__pyx_v_i]) == 'T');
    if (__pyx_t_6) {
/* … */
      goto __pyx_L5;
    }
+105:             out_onehot[i, 3] = 1
      __pyx_t_8 = __pyx_v_i;
      __pyx_t_7 = 3;
      *((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=1 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ (__pyx_v_out_onehot.data + __pyx_t_8 * __pyx_v_out_onehot.strides[0]) )) + __pyx_t_7)) )) = 1;
 106:         else:
+107:             out_onehot[i, 4] = 1
    /*else*/ {
      __pyx_t_7 = __pyx_v_i;
      __pyx_t_8 = 4;
      *((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=1 */ ((char *) (((__pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t *) ( /* dim=0 */ (__pyx_v_out_onehot.data + __pyx_t_7 * __pyx_v_out_onehot.strides[0]) )) + __pyx_t_8)) )) = 1;
    }
    __pyx_L5:;
  }
 108: 
+109:     return out_onehot
  __Pyx_XDECREF(__pyx_r);
  __pyx_t_9 = __pyx_memoryview_fromslice(__pyx_v_out_onehot, 2, (PyObject *(*)(char *)) __pyx_memview_get_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t, (int (*)(char *, PyObject *)) __pyx_memview_set_nn___pyx_t_11chromatinhd_7loaders_17fragments_helpers_INT64_t, 0);; if (unlikely(!__pyx_t_9)) __PYX_ERR(0, 109, __pyx_L1_error)
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_r = __pyx_t_9;
  __pyx_t_9 = 0;
  goto __pyx_L0;