Generated by Cython 0.15.1+ on Sat Nov 5 14:14:56 2011
Raw output: histogram.c
1: #!/usr/bin/env python
/* "histogram.pyx":1 * #!/usr/bin/env python #<<<<<<<<<<<<<<* # -*- coding: utf8 -*- * # */ __pyx_t_1 = PyDict_New(); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_1)); if (PyObject_SetAttr(__pyx_m, __pyx_n_s____test__, ((PyObject *)__pyx_t_1))<0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
2: # -*- coding: utf8 -*-
3: #
4: # Project: Azimuthal integration
5: # https://forge.epn-campus.eu/projects/azimuthal
6: #
7: # File: "$Id$"
8: #
9: # Copyright (C) European Synchrotron Radiation Facility, Grenoble, France
10: #
11: # Principal author: Jérôme Kieffer (Jerome.Kieffer@ESRF.eu)
12: #
13: # This program is free software: you can redistribute it and/or modify
14: # it under the terms of the GNU General Public License as published by
15: # the Free Software Foundation, either version 3 of the License, or
16: # (at your option) any later version.
17: #
18: # This program is distributed in the hope that it will be useful,
19: # but WITHOUT ANY WARRANTY; without even the implied warranty of
20: # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
21: # GNU General Public License for more details.
22: #
23: # You should have received a copy of the GNU General Public License
24: # along with this program. If not, see <http://www.gnu.org/licenses/>.
25: #
26:
27:
28:
29: import cython
30: from cython.parallel cimport prange
31: cimport numpy
32: import numpy
/* "histogram.pyx":32 * from cython.parallel cimport prange * cimport numpy * import numpy #<<<<<<<<<<<<<<* cdef extern from "omp.h": * ctypedef struct omp_lock_t: */ __pyx_t_1 = __Pyx_Import(((PyObject *)__pyx_n_s__numpy), 0, -1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); if (PyObject_SetAttr(__pyx_m, __pyx_n_s__numpy, __pyx_t_1)<0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
33: cdef extern from "omp.h":
34: ctypedef struct omp_lock_t:
35: pass
36:
37: extern void omp_set_num_threads(int) nogil
38: extern int omp_get_num_threads() nogil
39: extern int omp_get_max_threads() nogil
40: extern int omp_get_thread_num() nogil
41: extern int omp_get_num_procs() nogil
42:
43: extern int omp_in_parallel() nogil
44: extern void omp_init_lock(omp_lock_t *) nogil
45: extern void omp_destroy_lock(omp_lock_t *) nogil
46: extern void omp_set_lock(omp_lock_t *) nogil
47: extern void omp_unset_lock(omp_lock_t *) nogil
48: extern int omp_test_lock(omp_lock_t *) nogil
49:
50:
51: cdef extern from "stdlib.h":
52: void free(void * ptr)nogil
53: void * calloc(size_t nmemb, size_t size)nogil
54: void * malloc(size_t size)nogil
55: cdef extern from "math.h":
56: double floor(double)nogil
57: double fabs(double)nogil
58: int isnan(double)
59:
60: # "ctypedef" assigns a corresponding compile-time type to DTYPE_t. For
61: # every type in the numpy module there's a corresponding compile-time
62: # type with a _t-suffix.
63: ctypedef numpy.int64_t DTYPE_int64_t
/* "histogram.pyx":63 * # every type in the numpy module there's a corresponding compile-time * # type with a _t-suffix. * ctypedef numpy.int64_t DTYPE_int64_t #<<<<<<<<<<<<<<* ctypedef numpy.float64_t DTYPE_float64_t * */ typedef __pyx_t_5numpy_int64_t __pyx_t_9histogram_DTYPE_int64_t;
64: ctypedef numpy.float64_t DTYPE_float64_t
65:
66: @cython.cdivision(True)
67: @cython.boundscheck(False)
68: @cython.wraparound(False)
69: def histogram(numpy.ndarray pos not None, \
/* "histogram.pyx":69 * @cython.boundscheck(False) * @cython.wraparound(False) * def histogram(numpy.ndarray pos not None, \ #<<<<<<<<<<<<<<* numpy.ndarray weights not None, \ * long bins=100, */ static PyObject *__pyx_pf_9histogram_histogram(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/ static char __pyx_doc_9histogram_histogram[] = "\n Calculates histogram of pos weighted by weights\n \n @param pos: 2Theta array\n @param weights: array with intensities\n @param bins: number of output bins\n @param pixelSize_in_Pos: size of a pixels in 2theta\n @param nthread: maximum number of thread to use. By default: maximum available. \n One can also limit this with OMP_NUM_THREADS environment variable\n \n @return 2theta, I, weighted histogram, raw histogram\n "; static PyMethodDef __pyx_mdef_9histogram_histogram = {__Pyx_NAMESTR("histogram"), (PyCFunction)__pyx_pf_9histogram_histogram, METH_VARARGS|METH_KEYWORDS, __Pyx_DOCSTR(__pyx_doc_9histogram_histogram)}; static PyObject *__pyx_pf_9histogram_histogram(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) { PyArrayObject *__pyx_v_pos = 0; PyArrayObject *__pyx_v_weights = 0; long __pyx_v_bins; PyObject *__pyx_v_bin_range = 0; PyObject *__pyx_v_pixelSize_in_Pos = 0; PyObject *__pyx_v_nthread = 0; double __pyx_v_dummy; long __pyx_v_size; PyArrayObject *__pyx_v_cpos = 0; PyArrayObject *__pyx_v_cdata = 0; PyArrayObject *__pyx_v_outData = 0; PyArrayObject *__pyx_v_outCount = 0; PyArrayObject *__pyx_v_outMerge = 0; PyArrayObject *__pyx_v_outPos = 0; CYTHON_UNUSED PyArrayObject *__pyx_v_temp = 0; double __pyx_v_bin_edge_min; double __pyx_v_bin_edge_max; double __pyx_v_bin_width; double __pyx_v_inv_bin_width; double __pyx_v_a; double __pyx_v_d; double __pyx_v_fbin; double __pyx_v_ffbin; double __pyx_v_dInt; double __pyx_v_dIntR; double __pyx_v_dIntL; CYTHON_UNUSED double __pyx_v_dTmp; double __pyx_v_dbin; double __pyx_v_inv_dbin2; double __pyx_v_tmp_count; double __pyx_v_tmp_data; double __pyx_v_epsilon; long __pyx_v_bin; long __pyx_v_i; long __pyx_v_idx; long __pyx_v_t; long __pyx_v_dest; double *__pyx_v_bigCount; double *__pyx_v_bigData; double __pyx_v_dtmp; Py_buffer __pyx_bstruct_outPos; Py_ssize_t __pyx_bstride_0_outPos = 0; Py_ssize_t __pyx_bshape_0_outPos = 0; Py_buffer __pyx_bstruct_outMerge; Py_ssize_t __pyx_bstride_0_outMerge = 0; Py_ssize_t __pyx_bshape_0_outMerge = 0; Py_buffer __pyx_bstruct_outCount; Py_ssize_t __pyx_bstride_0_outCount = 0; Py_ssize_t __pyx_bshape_0_outCount = 0; Py_buffer __pyx_bstruct_cdata; Py_ssize_t __pyx_bstride_0_cdata = 0; Py_ssize_t __pyx_bshape_0_cdata = 0; Py_buffer __pyx_bstruct_cpos; Py_ssize_t __pyx_bstride_0_cpos = 0; Py_ssize_t __pyx_bshape_0_cpos = 0; Py_buffer __pyx_bstruct_temp; Py_ssize_t __pyx_bstride_0_temp = 0; Py_ssize_t __pyx_bshape_0_temp = 0; Py_buffer __pyx_bstruct_outData; Py_ssize_t __pyx_bstride_0_outData = 0; Py_ssize_t __pyx_bshape_0_outData = 0; PyObject *__pyx_r = NULL; static PyObject **__pyx_pyargnames[] = {&__pyx_n_s__pos,&__pyx_n_s__weights,&__pyx_n_s__bins,&__pyx_n_s__bin_range,&__pyx_n_s__pixelSize_in_Pos,&__pyx_n_s__nthread,&__pyx_n_s__dummy,0}; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("histogram"); __pyx_self = __pyx_self; { PyObject* values[7] = {0,0,0,0,0,0,0}; /* "histogram.pyx":69 * @cython.boundscheck(False) * @cython.wraparound(False) * def histogram(numpy.ndarray pos not None, \ #<<<<<<<<<<<<<<* numpy.ndarray weights not None, \ * long bins=100, */ __pyx_k_tuple_19 = PyTuple_New(40); if (unlikely(!__pyx_k_tuple_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_k_tuple_19)); __Pyx_INCREF(((PyObject *)__pyx_n_s__pos)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 0, ((PyObject *)__pyx_n_s__pos)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__pos)); __Pyx_INCREF(((PyObject *)__pyx_n_s__weights)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 1, ((PyObject *)__pyx_n_s__weights)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__weights)); __Pyx_INCREF(((PyObject *)__pyx_n_s__bins)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 2, ((PyObject *)__pyx_n_s__bins)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__bins)); __Pyx_INCREF(((PyObject *)__pyx_n_s__bin_range)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 3, ((PyObject *)__pyx_n_s__bin_range)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__bin_range)); __Pyx_INCREF(((PyObject *)__pyx_n_s__pixelSize_in_Pos)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 4, ((PyObject *)__pyx_n_s__pixelSize_in_Pos)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__pixelSize_in_Pos)); __Pyx_INCREF(((PyObject *)__pyx_n_s__nthread)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 5, ((PyObject *)__pyx_n_s__nthread)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__nthread)); __Pyx_INCREF(((PyObject *)__pyx_n_s__dummy)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 6, ((PyObject *)__pyx_n_s__dummy)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__dummy)); __Pyx_INCREF(((PyObject *)__pyx_n_s__size)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 7, ((PyObject *)__pyx_n_s__size)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__size)); __Pyx_INCREF(((PyObject *)__pyx_n_s__cpos)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 8, ((PyObject *)__pyx_n_s__cpos)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__cpos)); __Pyx_INCREF(((PyObject *)__pyx_n_s__cdata)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 9, ((PyObject *)__pyx_n_s__cdata)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__cdata)); __Pyx_INCREF(((PyObject *)__pyx_n_s__outData)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 10, ((PyObject *)__pyx_n_s__outData)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__outData)); __Pyx_INCREF(((PyObject *)__pyx_n_s__outCount)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 11, ((PyObject *)__pyx_n_s__outCount)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__outCount)); __Pyx_INCREF(((PyObject *)__pyx_n_s__outMerge)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 12, ((PyObject *)__pyx_n_s__outMerge)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__outMerge)); __Pyx_INCREF(((PyObject *)__pyx_n_s__outPos)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 13, ((PyObject *)__pyx_n_s__outPos)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__outPos)); __Pyx_INCREF(((PyObject *)__pyx_n_s__temp)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 14, ((PyObject *)__pyx_n_s__temp)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__temp)); __Pyx_INCREF(((PyObject *)__pyx_n_s__bin_edge_min)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 15, ((PyObject *)__pyx_n_s__bin_edge_min)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__bin_edge_min)); __Pyx_INCREF(((PyObject *)__pyx_n_s__bin_edge_max)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 16, ((PyObject *)__pyx_n_s__bin_edge_max)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__bin_edge_max)); __Pyx_INCREF(((PyObject *)__pyx_n_s__bin_width)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 17, ((PyObject *)__pyx_n_s__bin_width)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__bin_width)); __Pyx_INCREF(((PyObject *)__pyx_n_s__inv_bin_width)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 18, ((PyObject *)__pyx_n_s__inv_bin_width)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__inv_bin_width)); __Pyx_INCREF(((PyObject *)__pyx_n_s__a)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 19, ((PyObject *)__pyx_n_s__a)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__a)); __Pyx_INCREF(((PyObject *)__pyx_n_s__d)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 20, ((PyObject *)__pyx_n_s__d)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__d)); __Pyx_INCREF(((PyObject *)__pyx_n_s__fbin)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 21, ((PyObject *)__pyx_n_s__fbin)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__fbin)); __Pyx_INCREF(((PyObject *)__pyx_n_s__ffbin)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 22, ((PyObject *)__pyx_n_s__ffbin)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__ffbin)); __Pyx_INCREF(((PyObject *)__pyx_n_s__dInt)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 23, ((PyObject *)__pyx_n_s__dInt)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__dInt)); __Pyx_INCREF(((PyObject *)__pyx_n_s__dIntR)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 24, ((PyObject *)__pyx_n_s__dIntR)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__dIntR)); __Pyx_INCREF(((PyObject *)__pyx_n_s__dIntL)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 25, ((PyObject *)__pyx_n_s__dIntL)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__dIntL)); __Pyx_INCREF(((PyObject *)__pyx_n_s__dTmp)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 26, ((PyObject *)__pyx_n_s__dTmp)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__dTmp)); __Pyx_INCREF(((PyObject *)__pyx_n_s__dbin)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 27, ((PyObject *)__pyx_n_s__dbin)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__dbin)); __Pyx_INCREF(((PyObject *)__pyx_n_s__inv_dbin2)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 28, ((PyObject *)__pyx_n_s__inv_dbin2)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__inv_dbin2)); __Pyx_INCREF(((PyObject *)__pyx_n_s__tmp_count)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 29, ((PyObject *)__pyx_n_s__tmp_count)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__tmp_count)); __Pyx_INCREF(((PyObject *)__pyx_n_s__tmp_data)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 30, ((PyObject *)__pyx_n_s__tmp_data)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__tmp_data)); __Pyx_INCREF(((PyObject *)__pyx_n_s__epsilon)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 31, ((PyObject *)__pyx_n_s__epsilon)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__epsilon)); __Pyx_INCREF(((PyObject *)__pyx_n_s__bin)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 32, ((PyObject *)__pyx_n_s__bin)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__bin)); __Pyx_INCREF(((PyObject *)__pyx_n_s__i)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 33, ((PyObject *)__pyx_n_s__i)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__i)); __Pyx_INCREF(((PyObject *)__pyx_n_s__idx)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 34, ((PyObject *)__pyx_n_s__idx)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__idx)); __Pyx_INCREF(((PyObject *)__pyx_n_s__t)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 35, ((PyObject *)__pyx_n_s__t)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__t)); __Pyx_INCREF(((PyObject *)__pyx_n_s__dest)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 36, ((PyObject *)__pyx_n_s__dest)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__dest)); __Pyx_INCREF(((PyObject *)__pyx_n_s__bigCount)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 37, ((PyObject *)__pyx_n_s__bigCount)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__bigCount)); __Pyx_INCREF(((PyObject *)__pyx_n_s__bigData)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 38, ((PyObject *)__pyx_n_s__bigData)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__bigData)); __Pyx_INCREF(((PyObject *)__pyx_n_s__dtmp)); PyTuple_SET_ITEM(__pyx_k_tuple_19, 39, ((PyObject *)__pyx_n_s__dtmp)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__dtmp)); __Pyx_GIVEREF(((PyObject *)__pyx_k_tuple_19)); /* "histogram.pyx":69 * @cython.boundscheck(False) * @cython.wraparound(False) * def histogram(numpy.ndarray pos not None, \ #<<<<<<<<<<<<<<* numpy.ndarray weights not None, \ * long bins=100, */ __pyx_t_1 = PyCFunction_NewEx(&__pyx_mdef_9histogram_histogram, NULL, __pyx_n_s__histogram); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); if (PyObject_SetAttr(__pyx_m, __pyx_n_s__histogram, __pyx_t_1)<0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_k_codeobj_20 = (PyObject*)__Pyx_PyCode_New(7, 0, 40, 0, 0, __pyx_empty_bytes, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_k_tuple_19, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_kp_s_21, __pyx_n_s__histogram, 69, __pyx_empty_bytes); if (unlikely(!__pyx_k_codeobj_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
70: numpy.ndarray weights not None, \
71: long bins=100,
72: bin_range=None,
    /* "histogram.pyx":72
 *               numpy.ndarray weights not None, \
 *               long bins=100,
 *               bin_range=None,             # <<<<<<<<<<<<<<
 *               pixelSize_in_Pos=None,
 *               nthread=None,
 */
    values[3] = ((PyObject *)Py_None);
73: pixelSize_in_Pos=None,
    /* "histogram.pyx":73
 *               long bins=100,
 *               bin_range=None,
 *               pixelSize_in_Pos=None,             # <<<<<<<<<<<<<<
 *               nthread=None,
 *               double dummy=0.0):
 */
    values[4] = ((PyObject *)Py_None);
74: nthread=None,
    /* "histogram.pyx":74
 *               bin_range=None,
 *               pixelSize_in_Pos=None,
 *               nthread=None,             # <<<<<<<<<<<<<<
 *               double dummy=0.0):
 *     """
 */
    values[5] = ((PyObject *)Py_None);
    if (unlikely(__pyx_kwds)) {
      Py_ssize_t kw_args;
      switch (PyTuple_GET_SIZE(__pyx_args)) {
        case  7: values[6] = PyTuple_GET_ITEM(__pyx_args, 6);
        case  6: values[5] = PyTuple_GET_ITEM(__pyx_args, 5);
        case  5: values[4] = PyTuple_GET_ITEM(__pyx_args, 4);
        case  4: values[3] = PyTuple_GET_ITEM(__pyx_args, 3);
        case  3: values[2] = PyTuple_GET_ITEM(__pyx_args, 2);
        case  2: values[1] = PyTuple_GET_ITEM(__pyx_args, 1);
        case  1: values[0] = PyTuple_GET_ITEM(__pyx_args, 0);
        case  0: break;
        default: goto __pyx_L5_argtuple_error;
      }
      kw_args = PyDict_Size(__pyx_kwds);
      switch (PyTuple_GET_SIZE(__pyx_args)) {
        case  0:
        values[0] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__pos);
        if (likely(values[0])) kw_args--;
        else goto __pyx_L5_argtuple_error;
        case  1:
        values[1] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__weights);
        if (likely(values[1])) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("histogram", 0, 2, 7, 1); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        }
        case  2:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s__bins);
          if (value) { values[2] = value; kw_args--; }
        }
        case  3:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s__bin_range);
          if (value) { values[3] = value; kw_args--; }
        }
        case  4:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s__pixelSize_in_Pos);
          if (value) { values[4] = value; kw_args--; }
        }
        case  5:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s__nthread);
          if (value) { values[5] = value; kw_args--; }
        }
        case  6:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s__dummy);
          if (value) { values[6] = value; kw_args--; }
        }
      }
      if (unlikely(kw_args > 0)) {
        if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_pyargnames, 0, values, PyTuple_GET_SIZE(__pyx_args), "histogram") < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
      }
    } else {
      switch (PyTuple_GET_SIZE(__pyx_args)) {
        case  7: values[6] = PyTuple_GET_ITEM(__pyx_args, 6);
        case  6: values[5] = PyTuple_GET_ITEM(__pyx_args, 5);
        case  5: values[4] = PyTuple_GET_ITEM(__pyx_args, 4);
        case  4: values[3] = PyTuple_GET_ITEM(__pyx_args, 3);
        case  3: values[2] = PyTuple_GET_ITEM(__pyx_args, 2);
        case  2: values[1] = PyTuple_GET_ITEM(__pyx_args, 1);
        values[0] = PyTuple_GET_ITEM(__pyx_args, 0);
        break;
        default: goto __pyx_L5_argtuple_error;
      }
    }
    __pyx_v_pos = ((PyArrayObject *)values[0]);
    __pyx_v_weights = ((PyArrayObject *)values[1]);
    if (values[2]) {
      __pyx_v_bins = __Pyx_PyInt_AsLong(values[2]); if (unlikely((__pyx_v_bins == (long)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 71; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
    } else {
      __pyx_v_bins = ((long)100);
    }
    __pyx_v_bin_range = values[3];
    __pyx_v_pixelSize_in_Pos = values[4];
    __pyx_v_nthread = values[5];
    if (values[6]) {
      __pyx_v_dummy = __pyx_PyFloat_AsDouble(values[6]); if (unlikely((__pyx_v_dummy == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 75; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
    } else {
75: double dummy=0.0):
      /* "histogram.pyx":75
 *               pixelSize_in_Pos=None,
 *               nthread=None,
 *               double dummy=0.0):             # <<<<<<<<<<<<<<
 *     """
 *     Calculates histogram of pos weighted by weights
 */
      __pyx_v_dummy = ((double)0.0);
    }
  }
  goto __pyx_L4_argument_unpacking_done;
  __pyx_L5_argtuple_error:;
  __Pyx_RaiseArgtupleInvalid("histogram", 0, 2, 7, PyTuple_GET_SIZE(__pyx_args)); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
  __pyx_L3_error:;
  __Pyx_AddTraceback("histogram.histogram", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __Pyx_RefNannyFinishContext();
  return NULL;
  __pyx_L4_argument_unpacking_done:;
  __pyx_bstruct_cpos.buf = NULL;
  __pyx_bstruct_cdata.buf = NULL;
  __pyx_bstruct_outData.buf = NULL;
  __pyx_bstruct_outCount.buf = NULL;
  __pyx_bstruct_outMerge.buf = NULL;
  __pyx_bstruct_outPos.buf = NULL;
  __pyx_bstruct_temp.buf = NULL;
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_pos), __pyx_ptype_5numpy_ndarray, 0, "pos", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_weights), __pyx_ptype_5numpy_ndarray, 0, "weights", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
76: """
77: Calculates histogram of pos weighted by weights
78:
79: @param pos: 2Theta array
80: @param weights: array with intensities
81: @param bins: number of output bins
82: @param pixelSize_in_Pos: size of a pixels in 2theta
83: @param nthread: maximum number of thread to use. By default: maximum available.
84: One can also limit this with OMP_NUM_THREADS environment variable
85:
86: @return 2theta, I, weighted histogram, raw histogram
87: """
88:
89: assert pos.size == weights.size
/* "histogram.pyx":89 * """ * * assert pos.size == weights.size #<<<<<<<<<<<<<<* assert bins > 1 * cdef long size = pos.size */ #ifndef CYTHON_WITHOUT_ASSERTIONS __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_pos), __pyx_n_s__size); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 89; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_2 = PyObject_GetAttr(((PyObject *)__pyx_v_weights), __pyx_n_s__size); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 89; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __pyx_t_3 = PyObject_RichCompare(__pyx_t_1, __pyx_t_2, Py_EQ); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 89; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __pyx_t_4 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_4<0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 89; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; if (unlikely(!__pyx_t_4)) { PyErr_SetNone(PyExc_AssertionError); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 89; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } #endif
90: assert bins > 1
/* "histogram.pyx":90 * * assert pos.size == weights.size * assert bins > 1 #<<<<<<<<<<<<<<* cdef long size = pos.size * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos = pos.astype("float64").ravel() */ #ifndef CYTHON_WITHOUT_ASSERTIONS if (unlikely(!(__pyx_v_bins > 1))) { PyErr_SetNone(PyExc_AssertionError); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 90; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } #endif
91: cdef long size = pos.size
/* "histogram.pyx":91 * assert pos.size == weights.size * assert bins > 1 * cdef long size = pos.size #<<<<<<<<<<<<<<* cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos = pos.astype("float64").ravel() * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cdata = weights.astype("float64").ravel() */ __pyx_t_3 = PyObject_GetAttr(((PyObject *)__pyx_v_pos), __pyx_n_s__size); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 91; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_5 = __Pyx_PyInt_AsLong(__pyx_t_3); if (unlikely((__pyx_t_5 == (long)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 91; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_v_size = __pyx_t_5;
 92:     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos = pos.astype("float64").ravel()
  /* "histogram.pyx":92
 *     assert  bins > 1
 *     cdef long  size = pos.size
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos = pos.astype("float64").ravel()             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cdata = weights.astype("float64").ravel()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")
 */
  __pyx_t_3 = PyObject_GetAttr(((PyObject *)__pyx_v_pos), __pyx_n_s__astype); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 92; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_2 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_k_tuple_1), NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 92; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_t_3 = PyObject_GetAttr(__pyx_t_2, __pyx_n_s__ravel); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 92; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __pyx_t_2 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 92; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  if (!(likely(((__pyx_t_2) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_2, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 92; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_6 = ((PyArrayObject *)__pyx_t_2);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_cpos, (PyObject*)__pyx_t_6, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {
      __pyx_v_cpos = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_cpos.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 92; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_bstride_0_cpos = __pyx_bstruct_cpos.strides[0];
      __pyx_bshape_0_cpos = __pyx_bstruct_cpos.shape[0];
    }
  }
  __pyx_t_6 = 0;
  __pyx_v_cpos = ((PyArrayObject *)__pyx_t_2);
  __pyx_t_2 = 0;
  /* "histogram.pyx":92
 *     assert  bins > 1
 *     cdef long  size = pos.size
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos = pos.astype("float64").ravel()             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cdata = weights.astype("float64").ravel()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")
 */
  __pyx_k_tuple_1 = PyTuple_New(1); if (unlikely(!__pyx_k_tuple_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 92; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_k_tuple_1));
  __Pyx_INCREF(((PyObject *)__pyx_n_s__float64));
  PyTuple_SET_ITEM(__pyx_k_tuple_1, 0, ((PyObject *)__pyx_n_s__float64));
  __Pyx_GIVEREF(((PyObject *)__pyx_n_s__float64));
  __Pyx_GIVEREF(((PyObject *)__pyx_k_tuple_1));
 93:     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cdata = weights.astype("float64").ravel()
  /* "histogram.pyx":93
 *     cdef long  size = pos.size
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos = pos.astype("float64").ravel()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cdata = weights.astype("float64").ravel()             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64")
 */
  __pyx_t_2 = PyObject_GetAttr(((PyObject *)__pyx_v_weights), __pyx_n_s__astype); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 93; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __pyx_t_3 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_k_tuple_2), NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 93; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __pyx_t_2 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__ravel); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 93; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_t_3 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 93; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  if (!(likely(((__pyx_t_3) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_3, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 93; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_7 = ((PyArrayObject *)__pyx_t_3);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_cdata, (PyObject*)__pyx_t_7, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {
      __pyx_v_cdata = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_cdata.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 93; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_bstride_0_cdata = __pyx_bstruct_cdata.strides[0];
      __pyx_bshape_0_cdata = __pyx_bstruct_cdata.shape[0];
    }
  }
  __pyx_t_7 = 0;
  __pyx_v_cdata = ((PyArrayObject *)__pyx_t_3);
  __pyx_t_3 = 0;
  /* "histogram.pyx":93
 *     cdef long  size = pos.size
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos = pos.astype("float64").ravel()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cdata = weights.astype("float64").ravel()             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64")
 */
  __pyx_k_tuple_2 = PyTuple_New(1); if (unlikely(!__pyx_k_tuple_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 93; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_k_tuple_2));
  __Pyx_INCREF(((PyObject *)__pyx_n_s__float64));
  PyTuple_SET_ITEM(__pyx_k_tuple_2, 0, ((PyObject *)__pyx_n_s__float64));
  __Pyx_GIVEREF(((PyObject *)__pyx_n_s__float64));
  __Pyx_GIVEREF(((PyObject *)__pyx_k_tuple_2));
94: cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")
  /* "histogram.pyx":94
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos = pos.astype("float64").ravel()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cdata = weights.astype("float64").ravel()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outMerge = numpy.zeros(bins, dtype="float64")
 */
  __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 94; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_2 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__zeros); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 94; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_t_3 = PyInt_FromLong(__pyx_v_bins); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 94; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_1 = PyTuple_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 94; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_1));
  PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_t_3);
  __Pyx_GIVEREF(__pyx_t_3);
  __pyx_t_3 = 0;
  __pyx_t_3 = PyDict_New(); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 94; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_3));
  if (PyDict_SetItem(__pyx_t_3, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64)) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 94; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_8 = PyEval_CallObjectWithKeywords(__pyx_t_2, ((PyObject *)__pyx_t_1), ((PyObject *)__pyx_t_3)); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 94; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_8);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0;
  if (!(likely(((__pyx_t_8) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_8, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 94; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_9 = ((PyArrayObject *)__pyx_t_8);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_outData, (PyObject*)__pyx_t_9, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) {
      __pyx_v_outData = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_outData.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 94; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_bstride_0_outData = __pyx_bstruct_outData.strides[0];
      __pyx_bshape_0_outData = __pyx_bstruct_outData.shape[0];
    }
  }
  __pyx_t_9 = 0;
  __pyx_v_outData = ((PyArrayObject *)__pyx_t_8);
  __pyx_t_8 = 0;
95: cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64")
  /* "histogram.pyx":95
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cdata = weights.astype("float64").ravel()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64")             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outMerge = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outPos = numpy.zeros(bins, dtype="float64")
 */
  __pyx_t_8 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 95; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_8);
  __pyx_t_3 = PyObject_GetAttr(__pyx_t_8, __pyx_n_s__zeros); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 95; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0;
  __pyx_t_8 = PyInt_FromLong(__pyx_v_bins); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 95; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_8);
  __pyx_t_1 = PyTuple_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 95; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_1));
  PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_t_8);
  __Pyx_GIVEREF(__pyx_t_8);
  __pyx_t_8 = 0;
  __pyx_t_8 = PyDict_New(); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 95; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_8));
  if (PyDict_SetItem(__pyx_t_8, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64)) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 95; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_2 = PyEval_CallObjectWithKeywords(__pyx_t_3, ((PyObject *)__pyx_t_1), ((PyObject *)__pyx_t_8)); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 95; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_8)); __pyx_t_8 = 0;
  if (!(likely(((__pyx_t_2) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_2, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 95; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_10 = ((PyArrayObject *)__pyx_t_2);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_outCount, (PyObject*)__pyx_t_10, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) {
      __pyx_v_outCount = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_outCount.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 95; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_bstride_0_outCount = __pyx_bstruct_outCount.strides[0];
      __pyx_bshape_0_outCount = __pyx_bstruct_outCount.shape[0];
    }
  }
  __pyx_t_10 = 0;
  __pyx_v_outCount = ((PyArrayObject *)__pyx_t_2);
  __pyx_t_2 = 0;
96: cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outMerge = numpy.zeros(bins, dtype="float64")
/* "histogram.pyx":96 * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outMerge = numpy.zeros(bins, dtype="float64") #<<<<<<<<<<<<<<* cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outPos = numpy.zeros(bins, dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] temp = numpy.zeros(size - 1, dtype="float64") */ __pyx_t_2 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 96; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __pyx_t_8 = PyObject_GetAttr(__pyx_t_2, __pyx_n_s__zeros); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 96; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __pyx_t_2 = PyInt_FromLong(__pyx_v_bins); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 96; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __pyx_t_1 = PyTuple_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 96; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_1)); PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_t_2); __Pyx_GIVEREF(__pyx_t_2); __pyx_t_2 = 0; __pyx_t_2 = PyDict_New(); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 96; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_2)); if (PyDict_SetItem(__pyx_t_2, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64))<0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 96; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_3 = PyEval_CallObjectWithKeywords(__pyx_t_8, ((PyObject *)__pyx_t_1), ((PyObject *)__pyx_t_2)); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 96; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0; if (!(likely(((__pyx_t_3) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_3, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 96; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_11 = ((PyArrayObject *)__pyx_t_3); { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_outMerge, (PyObject*)__pyx_t_11, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) { __pyx_v_outMerge = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_outMerge.buf = NULL; {__pyx_filename = __pyx_f[0]; __pyx_lineno = 96; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } else {__pyx_bstride_0_outMerge = __pyx_bstruct_outMerge.strides[0]; __pyx_bshape_0_outMerge = __pyx_bstruct_outMerge.shape[0]; } } __pyx_t_11 = 0; __pyx_v_outMerge = ((PyArrayObject *)__pyx_t_3); __pyx_t_3 = 0;
97: cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outPos = numpy.zeros(bins, dtype="float64")
/* "histogram.pyx":97 * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outMerge = numpy.zeros(bins, dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outPos = numpy.zeros(bins, dtype="float64") #<<<<<<<<<<<<<<* cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] temp = numpy.zeros(size - 1, dtype="float64") * cdef double bin_edge_min = pos.min() */ __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_2 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__zeros); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = PyInt_FromLong(__pyx_v_bins); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_1 = PyTuple_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_1)); PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_t_3); __Pyx_GIVEREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = PyDict_New(); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_3)); if (PyDict_SetItem(__pyx_t_3, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64))<0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_8 = PyEval_CallObjectWithKeywords(__pyx_t_2, ((PyObject *)__pyx_t_1), ((PyObject *)__pyx_t_3)); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0; if (!(likely(((__pyx_t_8) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_8, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_12 = ((PyArrayObject *)__pyx_t_8); { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_outPos, (PyObject*)__pyx_t_12, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) { __pyx_v_outPos = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_outPos.buf = NULL; {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } else {__pyx_bstride_0_outPos = __pyx_bstruct_outPos.strides[0]; __pyx_bshape_0_outPos = __pyx_bstruct_outPos.shape[0]; } } __pyx_t_12 = 0; __pyx_v_outPos = ((PyArrayObject *)__pyx_t_8); __pyx_t_8 = 0;
98: cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] temp = numpy.zeros(size - 1, dtype="float64")
/* "histogram.pyx":98 * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outMerge = numpy.zeros(bins, dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outPos = numpy.zeros(bins, dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] temp = numpy.zeros(size - 1, dtype="float64") #<<<<<<<<<<<<<<* cdef double bin_edge_min = pos.min() * cdef double bin_edge_max = pos.max() * (1 + numpy.finfo(numpy.double).eps) */ __pyx_t_8 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 98; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __pyx_t_3 = PyObject_GetAttr(__pyx_t_8, __pyx_n_s__zeros); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 98; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = PyInt_FromLong((__pyx_v_size - 1)); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 98; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __pyx_t_1 = PyTuple_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 98; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_1)); PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_t_8); __Pyx_GIVEREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = PyDict_New(); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 98; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_8)); if (PyDict_SetItem(__pyx_t_8, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64))<0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 98; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_2 = PyEval_CallObjectWithKeywords(__pyx_t_3, ((PyObject *)__pyx_t_1), ((PyObject *)__pyx_t_8)); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 98; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_8)); __pyx_t_8 = 0; if (!(likely(((__pyx_t_2) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_2, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 98; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_13 = ((PyArrayObject *)__pyx_t_2); { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_temp, (PyObject*)__pyx_t_13, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) { __pyx_v_temp = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_temp.buf = NULL; {__pyx_filename = __pyx_f[0]; __pyx_lineno = 98; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } else {__pyx_bstride_0_temp = __pyx_bstruct_temp.strides[0]; __pyx_bshape_0_temp = __pyx_bstruct_temp.shape[0]; } } __pyx_t_13 = 0; __pyx_v_temp = ((PyArrayObject *)__pyx_t_2); __pyx_t_2 = 0;
99: cdef double bin_edge_min = pos.min()
/* "histogram.pyx":99 * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] outPos = numpy.zeros(bins, dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] temp = numpy.zeros(size - 1, dtype="float64") * cdef double bin_edge_min = pos.min() #<<<<<<<<<<<<<<* cdef double bin_edge_max = pos.max() * (1 + numpy.finfo(numpy.double).eps) * if bin_range is not None: */ __pyx_t_2 = PyObject_GetAttr(((PyObject *)__pyx_v_pos), __pyx_n_s__min); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 99; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __pyx_t_8 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 99; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __pyx_t_14 = __pyx_PyFloat_AsDouble(__pyx_t_8); if (unlikely((__pyx_t_14 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 99; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_v_bin_edge_min = __pyx_t_14;
100: cdef double bin_edge_max = pos.max() * (1 + numpy.finfo(numpy.double).eps)
/* "histogram.pyx":100 * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] temp = numpy.zeros(size - 1, dtype="float64") * cdef double bin_edge_min = pos.min() * cdef double bin_edge_max = pos.max() * (1 + numpy.finfo(numpy.double).eps) #<<<<<<<<<<<<<<* if bin_range is not None: * bin_edge_min = bin_range[0] */ __pyx_t_8 = PyObject_GetAttr(((PyObject *)__pyx_v_pos), __pyx_n_s__max); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __pyx_t_2 = PyObject_Call(__pyx_t_8, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __pyx_t_1 = PyObject_GetAttr(__pyx_t_8, __pyx_n_s__finfo); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __pyx_t_3 = PyObject_GetAttr(__pyx_t_8, __pyx_n_s__double); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = PyTuple_New(1); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_8)); PyTuple_SET_ITEM(__pyx_t_8, 0, __pyx_t_3); __Pyx_GIVEREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_t_8), NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_8)); __pyx_t_8 = 0; __pyx_t_8 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__eps); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = PyNumber_Add(__pyx_int_1, __pyx_t_8); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = PyNumber_Multiply(__pyx_t_2, __pyx_t_3); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_14 = __pyx_PyFloat_AsDouble(__pyx_t_8); if (unlikely((__pyx_t_14 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_v_bin_edge_max = __pyx_t_14;
101: if bin_range is not None:
/* "histogram.pyx":101 * cdef double bin_edge_min = pos.min() * cdef double bin_edge_max = pos.max() * (1 + numpy.finfo(numpy.double).eps) * if bin_range is not None: #<<<<<<<<<<<<<<* bin_edge_min = bin_range[0] * bin_edge_max = bin_range[1] * (1 + numpy.finfo(numpy.double).eps) */ __pyx_t_4 = (__pyx_v_bin_range != Py_None); if (__pyx_t_4) {
102: bin_edge_min = bin_range[0]
    /* "histogram.pyx":102
 *     cdef double bin_edge_max = pos.max() * (1 + numpy.finfo(numpy.double).eps)
 *     if bin_range is not None:
 *         bin_edge_min = bin_range[0]             # <<<<<<<<<<<<<<
 *         bin_edge_max = bin_range[1] * (1 + numpy.finfo(numpy.double).eps)
 *     cdef double bin_width = (bin_edge_max - bin_edge_min) / (< double > (bins))
 */
    __pyx_t_8 = __Pyx_GetItemInt(__pyx_v_bin_range, 0, sizeof(long), PyInt_FromLong); if (!__pyx_t_8) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 102; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_8);
    __pyx_t_14 = __pyx_PyFloat_AsDouble(__pyx_t_8); if (unlikely((__pyx_t_14 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 102; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0;
    __pyx_v_bin_edge_min = __pyx_t_14;
103: bin_edge_max = bin_range[1] * (1 + numpy.finfo(numpy.double).eps)
    /* "histogram.pyx":103
 *     if bin_range is not None:
 *         bin_edge_min = bin_range[0]
 *         bin_edge_max = bin_range[1] * (1 + numpy.finfo(numpy.double).eps)             # <<<<<<<<<<<<<<
 *     cdef double bin_width = (bin_edge_max - bin_edge_min) / (< double > (bins))
 *     cdef double inv_bin_width = (< double > (bins)) / (bin_edge_max - bin_edge_min)
 */
    __pyx_t_8 = __Pyx_GetItemInt(__pyx_v_bin_range, 1, sizeof(long), PyInt_FromLong); if (!__pyx_t_8) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 103; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_8);
    __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 103; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_3);
    __pyx_t_2 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__finfo); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 103; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_2);
    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
    __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 103; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_3);
    __pyx_t_1 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__double); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 103; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
    __pyx_t_3 = PyTuple_New(1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 103; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(((PyObject *)__pyx_t_3));
    PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_t_1);
    __Pyx_GIVEREF(__pyx_t_1);
    __pyx_t_1 = 0;
    __pyx_t_1 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_t_3), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 103; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
    __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0;
    __pyx_t_3 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__eps); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 103; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_3);
    __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
    __pyx_t_1 = PyNumber_Add(__pyx_int_1, __pyx_t_3); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 103; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
    __pyx_t_3 = PyNumber_Multiply(__pyx_t_8, __pyx_t_1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 103; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_3);
    __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0;
    __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
    __pyx_t_14 = __pyx_PyFloat_AsDouble(__pyx_t_3); if (unlikely((__pyx_t_14 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 103; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
    __pyx_v_bin_edge_max = __pyx_t_14;
    goto __pyx_L6;
  }
  __pyx_L6:;
104: cdef double bin_width = (bin_edge_max - bin_edge_min) / (< double > (bins))
/* "histogram.pyx":104 * bin_edge_min = bin_range[0] * bin_edge_max = bin_range[1] * (1 + numpy.finfo(numpy.double).eps) * cdef double bin_width = (bin_edge_max - bin_edge_min) / (<double > (bins)) #<<<<<<<<<<<<<<* cdef double inv_bin_width = (<double > (bins)) / (bin_edge_max - bin_edge_min) * cdef double a = 0.0 */ __pyx_v_bin_width = ((__pyx_v_bin_edge_max - __pyx_v_bin_edge_min) / ((double)__pyx_v_bins));
105: cdef double inv_bin_width = (< double > (bins)) / (bin_edge_max - bin_edge_min)
/* "histogram.pyx":105 * bin_edge_max = bin_range[1] * (1 + numpy.finfo(numpy.double).eps) * cdef double bin_width = (bin_edge_max - bin_edge_min) / (<double > (bins)) * cdef double inv_bin_width = (<double > (bins)) / (bin_edge_max - bin_edge_min) #<<<<<<<<<<<<<<* cdef double a = 0.0 * cdef double d = 0.0 */ __pyx_v_inv_bin_width = (((double)__pyx_v_bins) / (__pyx_v_bin_edge_max - __pyx_v_bin_edge_min));
106: cdef double a = 0.0
/* "histogram.pyx":106 * cdef double bin_width = (bin_edge_max - bin_edge_min) / (<double > (bins)) * cdef double inv_bin_width = (<double > (bins)) / (bin_edge_max - bin_edge_min) * cdef double a = 0.0 #<<<<<<<<<<<<<<* cdef double d = 0.0 * cdef double fbin = 0.0 */ __pyx_v_a = 0.0;
107: cdef double d = 0.0
/* "histogram.pyx":107 * cdef double inv_bin_width = (<double > (bins)) / (bin_edge_max - bin_edge_min) * cdef double a = 0.0 * cdef double d = 0.0 #<<<<<<<<<<<<<<* cdef double fbin = 0.0 * cdef double ffbin = 0.0 */ __pyx_v_d = 0.0;
108: cdef double fbin = 0.0
/* "histogram.pyx":108 * cdef double a = 0.0 * cdef double d = 0.0 * cdef double fbin = 0.0 #<<<<<<<<<<<<<<* cdef double ffbin = 0.0 * cdef double dInt = 0.0 */ __pyx_v_fbin = 0.0;
109: cdef double ffbin = 0.0
/* "histogram.pyx":109 * cdef double d = 0.0 * cdef double fbin = 0.0 * cdef double ffbin = 0.0 #<<<<<<<<<<<<<<* cdef double dInt = 0.0 * cdef double dIntR = 0.0 */ __pyx_v_ffbin = 0.0;
110: cdef double dInt = 0.0
/* "histogram.pyx":110 * cdef double fbin = 0.0 * cdef double ffbin = 0.0 * cdef double dInt = 0.0 #<<<<<<<<<<<<<<* cdef double dIntR = 0.0 * cdef double dIntL = 0.0 */ __pyx_v_dInt = 0.0;
111: cdef double dIntR = 0.0
/* "histogram.pyx":111 * cdef double ffbin = 0.0 * cdef double dInt = 0.0 * cdef double dIntR = 0.0 #<<<<<<<<<<<<<<* cdef double dIntL = 0.0 * cdef double dTmp = 0.0 */ __pyx_v_dIntR = 0.0;
112: cdef double dIntL = 0.0
/* "histogram.pyx":112 * cdef double dInt = 0.0 * cdef double dIntR = 0.0 * cdef double dIntL = 0.0 #<<<<<<<<<<<<<<* cdef double dTmp = 0.0 * cdef double dbin, inv_dbin2 = 0.0 */ __pyx_v_dIntL = 0.0;
113: cdef double dTmp = 0.0
/* "histogram.pyx":113 * cdef double dIntR = 0.0 * cdef double dIntL = 0.0 * cdef double dTmp = 0.0 #<<<<<<<<<<<<<<* cdef double dbin, inv_dbin2 = 0.0 * cdef double tmp_count, tmp_data = 0.0 */ __pyx_v_dTmp = 0.0;
114: cdef double dbin, inv_dbin2 = 0.0
/* "histogram.pyx":114 * cdef double dIntL = 0.0 * cdef double dTmp = 0.0 * cdef double dbin, inv_dbin2 = 0.0 #<<<<<<<<<<<<<<* cdef double tmp_count, tmp_data = 0.0 * cdef double epsilon = 1e-10 */ __pyx_v_inv_dbin2 = 0.0;
115: cdef double tmp_count, tmp_data = 0.0
/* "histogram.pyx":115 * cdef double dTmp = 0.0 * cdef double dbin, inv_dbin2 = 0.0 * cdef double tmp_count, tmp_data = 0.0 #<<<<<<<<<<<<<<* cdef double epsilon = 1e-10 * */ __pyx_v_tmp_data = 0.0;
116: cdef double epsilon = 1e-10
/* "histogram.pyx":116 * cdef double dbin, inv_dbin2 = 0.0 * cdef double tmp_count, tmp_data = 0.0 * cdef double epsilon = 1e-10 #<<<<<<<<<<<<<<* * cdef long bin = 0 */ __pyx_v_epsilon = 1e-10;
117:
118: cdef long bin = 0
/* "histogram.pyx":118 * cdef double epsilon = 1e-10 * * cdef long bin = 0 #<<<<<<<<<<<<<<* cdef long i, idx, t, dest = 0 * if nthread is not None: */ __pyx_v_bin = 0;
119: cdef long i, idx, t, dest = 0
/* "histogram.pyx":119 * * cdef long bin = 0 * cdef long i, idx, t, dest = 0 #<<<<<<<<<<<<<<* if nthread is not None: * if isinstance(nthread, int) and (nthread > 0): */ __pyx_v_dest = 0;
120: if nthread is not None:
/* "histogram.pyx":120 * cdef long bin = 0 * cdef long i, idx, t, dest = 0 * if nthread is not None: #<<<<<<<<<<<<<<* if isinstance(nthread, int) and (nthread > 0): * omp_set_num_threads(<int > nthread) */ __pyx_t_4 = (__pyx_v_nthread != Py_None); if (__pyx_t_4) {
121: if isinstance(nthread, int) and (nthread > 0):
    /* "histogram.pyx":121
 *     cdef long   i, idx, t, dest = 0
 *     if nthread is not None:
 *         if isinstance(nthread, int) and (nthread > 0):             # <<<<<<<<<<<<<<
 *             omp_set_num_threads(< int > nthread)
 * 
 */
    __pyx_t_3 = ((PyObject *)((PyObject*)(&PyInt_Type)));
    __Pyx_INCREF(__pyx_t_3);
    __pyx_t_4 = __Pyx_TypeCheck(__pyx_v_nthread, __pyx_t_3); 
    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
    if (__pyx_t_4) {
      __pyx_t_3 = PyObject_RichCompare(__pyx_v_nthread, __pyx_int_0, Py_GT); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 121; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_15 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_15 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 121; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_t_16 = __pyx_t_15;
    } else {
      __pyx_t_16 = __pyx_t_4;
    }
    if (__pyx_t_16) {
122: omp_set_num_threads(< int > nthread)
      /* "histogram.pyx":122
 *     if nthread is not None:
 *         if isinstance(nthread, int) and (nthread > 0):
 *             omp_set_num_threads(< int > nthread)             # <<<<<<<<<<<<<<
 * 
 *     cdef double * bigCount = < double *> calloc(bins * omp_get_max_threads(), sizeof(double))
 */
      __pyx_t_17 = __Pyx_PyInt_AsInt(__pyx_v_nthread); if (unlikely((__pyx_t_17 == (int)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 122; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      omp_set_num_threads(((int)__pyx_t_17));
      goto __pyx_L8;
    }
    __pyx_L8:;
    goto __pyx_L7;
  }
  __pyx_L7:;
123:
124: cdef double * bigCount = < double *> calloc(bins * omp_get_max_threads(), sizeof(double))
/* "histogram.pyx":124 * omp_set_num_threads(<int > nthread) * * cdef double * bigCount =<double *> calloc(bins * omp_get_max_threads(), sizeof(double)) #<<<<<<<<<<<<<<* cdef double * bigData =<double *> calloc(bins * omp_get_max_threads(), sizeof(double)) * if pixelSize_in_Pos is None: */ __pyx_v_bigCount = ((double *)calloc((__pyx_v_bins * omp_get_max_threads()), (sizeof(double))));
125: cdef double * bigData = < double *> calloc(bins * omp_get_max_threads(), sizeof(double))
/* "histogram.pyx":125 * * cdef double * bigCount =<double *> calloc(bins * omp_get_max_threads(), sizeof(double)) * cdef double * bigData =<double *> calloc(bins * omp_get_max_threads(), sizeof(double)) #<<<<<<<<<<<<<<* if pixelSize_in_Pos is None: * dbin = 0.5 */ __pyx_v_bigData = ((double *)calloc((__pyx_v_bins * omp_get_max_threads()), (sizeof(double))));
126: if pixelSize_in_Pos is None:
/* "histogram.pyx":126 * cdef double * bigCount =<double *> calloc(bins * omp_get_max_threads(), sizeof(double)) * cdef double * bigData =<double *> calloc(bins * omp_get_max_threads(), sizeof(double)) * if pixelSize_in_Pos is None: #<<<<<<<<<<<<<<* dbin = 0.5 * inv_dbin2 = 4.0 */ __pyx_t_16 = (__pyx_v_pixelSize_in_Pos == Py_None); if (__pyx_t_16) {
127: dbin = 0.5
    /* "histogram.pyx":127
 *     cdef double * bigData = < double *> calloc(bins * omp_get_max_threads(), sizeof(double))
 *     if pixelSize_in_Pos is None:
 *         dbin = 0.5             # <<<<<<<<<<<<<<
 *         inv_dbin2 = 4.0
 *     elif isinstance(pixelSize_in_Pos, (int, float)):
 */
    __pyx_v_dbin = 0.5;
128: inv_dbin2 = 4.0
    /* "histogram.pyx":128
 *     if pixelSize_in_Pos is None:
 *         dbin = 0.5
 *         inv_dbin2 = 4.0             # <<<<<<<<<<<<<<
 *     elif isinstance(pixelSize_in_Pos, (int, float)):
 *         dbin = 0.5 * (< double > pixelSize_in_Pos) * inv_bin_width
 */
    __pyx_v_inv_dbin2 = 4.0;
    goto __pyx_L9;
  }
129: elif isinstance(pixelSize_in_Pos, (int, float)):
/* "histogram.pyx":129 * dbin = 0.5 * inv_dbin2 = 4.0 * elif isinstance(pixelSize_in_Pos, (int, float)): #<<<<<<<<<<<<<<* dbin = 0.5 * (<double > pixelSize_in_Pos) * inv_bin_width * if dbin > 0.0: */ __Pyx_INCREF(__pyx_v_pixelSize_in_Pos); __pyx_t_3 = __pyx_v_pixelSize_in_Pos; __pyx_t_16 = PyInt_Check(__pyx_t_3); if (!__pyx_t_16) { __pyx_t_4 = PyFloat_Check(__pyx_t_3); __pyx_t_15 = __pyx_t_4; } else { __pyx_t_15 = __pyx_t_16; } __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; if (__pyx_t_15) {
130: dbin = 0.5 * (< double > pixelSize_in_Pos) * inv_bin_width
    /* "histogram.pyx":130
 *         inv_dbin2 = 4.0
 *     elif isinstance(pixelSize_in_Pos, (int, float)):
 *         dbin = 0.5 * (< double > pixelSize_in_Pos) * inv_bin_width             # <<<<<<<<<<<<<<
 *         if dbin > 0.0:
 *             inv_dbin2 = 1 / dbin / dbin
 */
    __pyx_t_14 = __pyx_PyFloat_AsDouble(__pyx_v_pixelSize_in_Pos); if (unlikely((__pyx_t_14 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 130; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __pyx_v_dbin = ((0.5 * ((double)__pyx_t_14)) * __pyx_v_inv_bin_width);
131: if dbin > 0.0:
    /* "histogram.pyx":131
 *     elif isinstance(pixelSize_in_Pos, (int, float)):
 *         dbin = 0.5 * (< double > pixelSize_in_Pos) * inv_bin_width
 *         if dbin > 0.0:             # <<<<<<<<<<<<<<
 *             inv_dbin2 = 1 / dbin / dbin
 *         else:
 */
    __pyx_t_15 = (__pyx_v_dbin > 0.0);
    if (__pyx_t_15) {
132: inv_dbin2 = 1 / dbin / dbin
      /* "histogram.pyx":132
 *         dbin = 0.5 * (< double > pixelSize_in_Pos) * inv_bin_width
 *         if dbin > 0.0:
 *             inv_dbin2 = 1 / dbin / dbin             # <<<<<<<<<<<<<<
 *         else:
 *             inv_dbin2 = 0.0
 */
      __pyx_v_inv_dbin2 = ((1.0 / __pyx_v_dbin) / __pyx_v_dbin);
      goto __pyx_L10;
    }
    /*else*/ {
133: else:
134: inv_dbin2 = 0.0
      /* "histogram.pyx":134
 *             inv_dbin2 = 1 / dbin / dbin
 *         else:
 *             inv_dbin2 = 0.0             # <<<<<<<<<<<<<<
 *     elif isinstance(pixelSize_in_Pos, numpy.ndarray):
 *         pass #TODO
 */
      __pyx_v_inv_dbin2 = 0.0;
    }
    __pyx_L10:;
    goto __pyx_L9;
  }
135: elif isinstance(pixelSize_in_Pos, numpy.ndarray):
/* "histogram.pyx":135 * else: * inv_dbin2 = 0.0 * elif isinstance(pixelSize_in_Pos, numpy.ndarray): #<<<<<<<<<<<<<<* pass #TODO * */ __pyx_t_3 = ((PyObject *)((PyObject*)__pyx_ptype_5numpy_ndarray)); __Pyx_INCREF(__pyx_t_3); __pyx_t_15 = __Pyx_TypeCheck(__pyx_v_pixelSize_in_Pos, __pyx_t_3); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; if (__pyx_t_15) { goto __pyx_L9; } __pyx_L9:;
136: pass #TODO
137:
138: if isnan(dbin) or isnan(inv_dbin2):
/* "histogram.pyx":138 * pass #TODO * * if isnan(dbin) or isnan(inv_dbin2): #<<<<<<<<<<<<<<* dbin = 0.0 * inv_dbin2 = 0.0 */ __pyx_t_17 = isnan(__pyx_v_dbin); if (!__pyx_t_17) { __pyx_t_18 = isnan(__pyx_v_inv_dbin2); __pyx_t_15 = __pyx_t_18; } else { __pyx_t_15 = __pyx_t_17; } if (__pyx_t_15) {
139: dbin = 0.0
    /* "histogram.pyx":139
 * 
 *     if isnan(dbin) or isnan(inv_dbin2):
 *         dbin = 0.0             # <<<<<<<<<<<<<<
 *         inv_dbin2 = 0.0
 * 
 */
    __pyx_v_dbin = 0.0;
140: inv_dbin2 = 0.0
    /* "histogram.pyx":140
 *     if isnan(dbin) or isnan(inv_dbin2):
 *         dbin = 0.0
 *         inv_dbin2 = 0.0             # <<<<<<<<<<<<<<
 * 
 *     with nogil:
 */
    __pyx_v_inv_dbin2 = 0.0;
    goto __pyx_L11;
  }
  __pyx_L11:;
141:
142: with nogil:
/* "histogram.pyx":142 * inv_dbin2 = 0.0 * * with nogil: #<<<<<<<<<<<<<<* for i in prange(size): * d = cdata[i] */ { #ifdef WITH_THREAD PyThreadState *_save = NULL; #endif Py_UNBLOCK_THREADS /*try:*/ { /* "histogram.pyx":142 * inv_dbin2 = 0.0 * * with nogil: #<<<<<<<<<<<<<<* for i in prange(size): * d = cdata[i] */ /*finally:*/ { Py_BLOCK_THREADS } }
143: for i in prange(size):
        /* "histogram.pyx":143
 * 
 *     with nogil:
 *         for i in prange(size):             # <<<<<<<<<<<<<<
 *             d = cdata[i]
 *             a = cpos[i]
 */
        __pyx_t_5 = __pyx_v_size;
        if (1 == 0) abort();
        {
            __pyx_t_20 = (__pyx_t_5 - 0) / 1;
            if (__pyx_t_20 > 0)
            {
                __pyx_v_i = 0;
                #ifdef _OPENMP
                #pragma omp parallel
                #endif /* _OPENMP */
                {
                    #ifdef _OPENMP
                    #pragma omp for lastprivate(__pyx_v_dest) lastprivate(__pyx_v_dIntL) firstprivate(__pyx_v_i) lastprivate(__pyx_v_i) lastprivate(__pyx_v_dIntR) lastprivate(__pyx_v_bin) lastprivate(__pyx_v_dtmp) lastprivate(__pyx_v_a) lastprivate(__pyx_v_ffbin) lastprivate(__pyx_v_d) lastprivate(__pyx_v_fbin) lastprivate(__pyx_v_dInt)
                    #endif /* _OPENMP */
                    for (__pyx_t_19 = 0; __pyx_t_19 < __pyx_t_20; __pyx_t_19++){
                        {
                            __pyx_v_i = 0 + 1 * __pyx_t_19;
                            /* Initialize private variables to invalid values */
                            __pyx_v_dest = ((long)0xbad0bad0);
                            __pyx_v_dIntL = ((double)__PYX_NAN);
                            __pyx_v_dIntR = ((double)__PYX_NAN);
                            __pyx_v_bin = ((long)0xbad0bad0);
                            __pyx_v_dtmp = ((double)__PYX_NAN);
                            __pyx_v_a = ((double)__PYX_NAN);
                            __pyx_v_ffbin = ((double)__PYX_NAN);
                            __pyx_v_d = ((double)__PYX_NAN);
                            __pyx_v_fbin = ((double)__PYX_NAN);
                            __pyx_v_dInt = ((double)__PYX_NAN);
144: d = cdata[i]
                            /* "histogram.pyx":144
 *     with nogil:
 *         for i in prange(size):
 *             d = cdata[i]             # <<<<<<<<<<<<<<
 *             a = cpos[i]
 *             if (a < bin_edge_min) or (a > bin_edge_max):
 */
                            __pyx_t_21 = __pyx_v_i;
                            __pyx_v_d = (*__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_cdata.buf, __pyx_t_21, __pyx_bstride_0_cdata));
145: a = cpos[i]
                            /* "histogram.pyx":145
 *         for i in prange(size):
 *             d = cdata[i]
 *             a = cpos[i]             # <<<<<<<<<<<<<<
 *             if (a < bin_edge_min) or (a > bin_edge_max):
 *                 continue
 */
                            __pyx_t_22 = __pyx_v_i;
                            __pyx_v_a = (*__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_cpos.buf, __pyx_t_22, __pyx_bstride_0_cpos));
146: if (a < bin_edge_min) or (a > bin_edge_max):
                            /* "histogram.pyx":146
 *             d = cdata[i]
 *             a = cpos[i]
 *             if (a < bin_edge_min) or (a > bin_edge_max):             # <<<<<<<<<<<<<<
 *                 continue
 *             fbin = (a - bin_edge_min) * inv_bin_width
 */
                            __pyx_t_15 = (__pyx_v_a < __pyx_v_bin_edge_min);
                            if (!__pyx_t_15) {
                              __pyx_t_16 = (__pyx_v_a > __pyx_v_bin_edge_max);
                              __pyx_t_4 = __pyx_t_16;
                            } else {
                              __pyx_t_4 = __pyx_t_15;
                            }
                            if (__pyx_t_4) {
147: continue
                              /* "histogram.pyx":147
 *             a = cpos[i]
 *             if (a < bin_edge_min) or (a > bin_edge_max):
 *                 continue             # <<<<<<<<<<<<<<
 *             fbin = (a - bin_edge_min) * inv_bin_width
 *             ffbin = floor(fbin)
 */
                              goto __pyx_L15_continue;
                              goto __pyx_L19;
                            }
                            __pyx_L19:;
148: fbin = (a - bin_edge_min) * inv_bin_width
                            /* "histogram.pyx":148
 *             if (a < bin_edge_min) or (a > bin_edge_max):
 *                 continue
 *             fbin = (a - bin_edge_min) * inv_bin_width             # <<<<<<<<<<<<<<
 *             ffbin = floor(fbin)
 *             bin = < long > ffbin
 */
                            __pyx_v_fbin = ((__pyx_v_a - __pyx_v_bin_edge_min) * __pyx_v_inv_bin_width);
149: ffbin = floor(fbin)
                            /* "histogram.pyx":149
 *                 continue
 *             fbin = (a - bin_edge_min) * inv_bin_width
 *             ffbin = floor(fbin)             # <<<<<<<<<<<<<<
 *             bin = < long > ffbin
 *             dest = omp_get_thread_num() * bins + bin
 */
                            __pyx_v_ffbin = floor(__pyx_v_fbin);
150: bin = < long > ffbin
                            /* "histogram.pyx":150
 *             fbin = (a - bin_edge_min) * inv_bin_width
 *             ffbin = floor(fbin)
 *             bin = < long > ffbin             # <<<<<<<<<<<<<<
 *             dest = omp_get_thread_num() * bins + bin
 *             dInt = 1.0
 */
                            __pyx_v_bin = ((long)__pyx_v_ffbin);
151: dest = omp_get_thread_num() * bins + bin
                            /* "histogram.pyx":151
 *             ffbin = floor(fbin)
 *             bin = < long > ffbin
 *             dest = omp_get_thread_num() * bins + bin             # <<<<<<<<<<<<<<
 *             dInt = 1.0
 *             if  bin > 0 :
 */
                            __pyx_v_dest = ((omp_get_thread_num() * __pyx_v_bins) + __pyx_v_bin);
152: dInt = 1.0
                            /* "histogram.pyx":152
 *             bin = < long > ffbin
 *             dest = omp_get_thread_num() * bins + bin
 *             dInt = 1.0             # <<<<<<<<<<<<<<
 *             if  bin > 0 :
 *                 dtmp = ffbin - (fbin - dbin)
 */
                            __pyx_v_dInt = 1.0;
153: if bin > 0 :
                            /* "histogram.pyx":153
 *             dest = omp_get_thread_num() * bins + bin
 *             dInt = 1.0
 *             if  bin > 0 :             # <<<<<<<<<<<<<<
 *                 dtmp = ffbin - (fbin - dbin)
 *                 if dtmp > 0:
 */
                            __pyx_t_4 = (__pyx_v_bin > 0);
                            if (__pyx_t_4) {
154: dtmp = ffbin - (fbin - dbin)
                              /* "histogram.pyx":154
 *             dInt = 1.0
 *             if  bin > 0 :
 *                 dtmp = ffbin - (fbin - dbin)             # <<<<<<<<<<<<<<
 *                 if dtmp > 0:
 *                     dIntL = 0.5 * dtmp * dtmp * inv_dbin2
 */
                              __pyx_v_dtmp = (__pyx_v_ffbin - (__pyx_v_fbin - __pyx_v_dbin));
155: if dtmp > 0:
                              /* "histogram.pyx":155
 *             if  bin > 0 :
 *                 dtmp = ffbin - (fbin - dbin)
 *                 if dtmp > 0:             # <<<<<<<<<<<<<<
 *                     dIntL = 0.5 * dtmp * dtmp * inv_dbin2
 *                     dInt = dInt - dIntL
 */
                              __pyx_t_4 = (__pyx_v_dtmp > 0.0);
                              if (__pyx_t_4) {
156: dIntL = 0.5 * dtmp * dtmp * inv_dbin2
                                /* "histogram.pyx":156
 *                 dtmp = ffbin - (fbin - dbin)
 *                 if dtmp > 0:
 *                     dIntL = 0.5 * dtmp * dtmp * inv_dbin2             # <<<<<<<<<<<<<<
 *                     dInt = dInt - dIntL
 *                     bigCount[dest - 1] += dIntL
 */
                                __pyx_v_dIntL = (((0.5 * __pyx_v_dtmp) * __pyx_v_dtmp) * __pyx_v_inv_dbin2);
157: dInt = dInt - dIntL
                                /* "histogram.pyx":157
 *                 if dtmp > 0:
 *                     dIntL = 0.5 * dtmp * dtmp * inv_dbin2
 *                     dInt = dInt - dIntL             # <<<<<<<<<<<<<<
 *                     bigCount[dest - 1] += dIntL
 *                     bigData[dest - 1] += d * dIntL
 */
                                __pyx_v_dInt = (__pyx_v_dInt - __pyx_v_dIntL);
158: bigCount[dest - 1] += dIntL
                                /* "histogram.pyx":158
 *                     dIntL = 0.5 * dtmp * dtmp * inv_dbin2
 *                     dInt = dInt - dIntL
 *                     bigCount[dest - 1] += dIntL             # <<<<<<<<<<<<<<
 *                     bigData[dest - 1] += d * dIntL
 * 
 */
                                __pyx_t_23 = (__pyx_v_dest - 1);
                                (__pyx_v_bigCount[__pyx_t_23]) = ((__pyx_v_bigCount[__pyx_t_23]) + __pyx_v_dIntL);
159: bigData[dest - 1] += d * dIntL
                                /* "histogram.pyx":159
 *                     dInt = dInt - dIntL
 *                     bigCount[dest - 1] += dIntL
 *                     bigData[dest - 1] += d * dIntL             # <<<<<<<<<<<<<<
 * 
 *             if bin < bins - 1 :
 */
                                __pyx_t_23 = (__pyx_v_dest - 1);
                                (__pyx_v_bigData[__pyx_t_23]) = ((__pyx_v_bigData[__pyx_t_23]) + (__pyx_v_d * __pyx_v_dIntL));
                                goto __pyx_L21;
                              }
                              __pyx_L21:;
                              goto __pyx_L20;
                            }
                            __pyx_L20:;
160:
161: if bin < bins - 1 :
                            /* "histogram.pyx":161
 *                     bigData[dest - 1] += d * dIntL
 * 
 *             if bin < bins - 1 :             # <<<<<<<<<<<<<<
 *                 dtmp = fbin + dbin - ffbin - 1
 *                 if dtmp > 0 :
 */
                            __pyx_t_4 = (__pyx_v_bin < (__pyx_v_bins - 1));
                            if (__pyx_t_4) {
162: dtmp = fbin + dbin - ffbin - 1
                              /* "histogram.pyx":162
 * 
 *             if bin < bins - 1 :
 *                 dtmp = fbin + dbin - ffbin - 1             # <<<<<<<<<<<<<<
 *                 if dtmp > 0 :
 *                     dIntR = 0.5 * dtmp * dtmp * inv_dbin2
 */
                              __pyx_v_dtmp = (((__pyx_v_fbin + __pyx_v_dbin) - __pyx_v_ffbin) - 1.0);
163: if dtmp > 0 :
                              /* "histogram.pyx":163
 *             if bin < bins - 1 :
 *                 dtmp = fbin + dbin - ffbin - 1
 *                 if dtmp > 0 :             # <<<<<<<<<<<<<<
 *                     dIntR = 0.5 * dtmp * dtmp * inv_dbin2
 *                     dInt = dInt - dIntR
 */
                              __pyx_t_4 = (__pyx_v_dtmp > 0.0);
                              if (__pyx_t_4) {
164: dIntR = 0.5 * dtmp * dtmp * inv_dbin2
                                /* "histogram.pyx":164
 *                 dtmp = fbin + dbin - ffbin - 1
 *                 if dtmp > 0 :
 *                     dIntR = 0.5 * dtmp * dtmp * inv_dbin2             # <<<<<<<<<<<<<<
 *                     dInt = dInt - dIntR
 *                     bigCount[dest + 1] += dIntR
 */
                                __pyx_v_dIntR = (((0.5 * __pyx_v_dtmp) * __pyx_v_dtmp) * __pyx_v_inv_dbin2);
165: dInt = dInt - dIntR
                                /* "histogram.pyx":165
 *                 if dtmp > 0 :
 *                     dIntR = 0.5 * dtmp * dtmp * inv_dbin2
 *                     dInt = dInt - dIntR             # <<<<<<<<<<<<<<
 *                     bigCount[dest + 1] += dIntR
 *                     bigData[dest + 1] += d * dIntR
 */
                                __pyx_v_dInt = (__pyx_v_dInt - __pyx_v_dIntR);
166: bigCount[dest + 1] += dIntR
                                /* "histogram.pyx":166
 *                     dIntR = 0.5 * dtmp * dtmp * inv_dbin2
 *                     dInt = dInt - dIntR
 *                     bigCount[dest + 1] += dIntR             # <<<<<<<<<<<<<<
 *                     bigData[dest + 1] += d * dIntR
 *             bigCount[dest] += dInt
 */
                                __pyx_t_23 = (__pyx_v_dest + 1);
                                (__pyx_v_bigCount[__pyx_t_23]) = ((__pyx_v_bigCount[__pyx_t_23]) + __pyx_v_dIntR);
167: bigData[dest + 1] += d * dIntR
                                /* "histogram.pyx":167
 *                     dInt = dInt - dIntR
 *                     bigCount[dest + 1] += dIntR
 *                     bigData[dest + 1] += d * dIntR             # <<<<<<<<<<<<<<
 *             bigCount[dest] += dInt
 *             bigData[dest] += d * dInt
 */
                                __pyx_t_23 = (__pyx_v_dest + 1);
                                (__pyx_v_bigData[__pyx_t_23]) = ((__pyx_v_bigData[__pyx_t_23]) + (__pyx_v_d * __pyx_v_dIntR));
                                goto __pyx_L23;
                              }
                              __pyx_L23:;
                              goto __pyx_L22;
                            }
                            __pyx_L22:;
168: bigCount[dest] += dInt
                            /* "histogram.pyx":168
 *                     bigCount[dest + 1] += dIntR
 *                     bigData[dest + 1] += d * dIntR
 *             bigCount[dest] += dInt             # <<<<<<<<<<<<<<
 *             bigData[dest] += d * dInt
 * 
 */
                            __pyx_t_23 = __pyx_v_dest;
                            (__pyx_v_bigCount[__pyx_t_23]) = ((__pyx_v_bigCount[__pyx_t_23]) + __pyx_v_dInt);
169: bigData[dest] += d * dInt
                            /* "histogram.pyx":169
 *                     bigData[dest + 1] += d * dIntR
 *             bigCount[dest] += dInt
 *             bigData[dest] += d * dInt             # <<<<<<<<<<<<<<
 * 
 *         for idx in prange(bins):
 */
                            __pyx_t_23 = __pyx_v_dest;
                            (__pyx_v_bigData[__pyx_t_23]) = ((__pyx_v_bigData[__pyx_t_23]) + (__pyx_v_d * __pyx_v_dInt));
                            goto __pyx_L25;
                            __pyx_L15_continue:;
                            goto __pyx_L25;
                            __pyx_L25:;
                        }
                    }
                }
            }
        }
170:
171: for idx in prange(bins):
        /* "histogram.pyx":171
 *             bigData[dest] += d * dInt
 * 
 *         for idx in prange(bins):             # <<<<<<<<<<<<<<
 *             outPos[idx] = bin_edge_min + (0.5 +< double > idx) * bin_width
 *             tmp_count = 0.0
 */
        __pyx_t_20 = __pyx_v_bins;
        if (1 == 0) abort();
        {
            __pyx_t_5 = (__pyx_t_20 - 0) / 1;
            if (__pyx_t_5 > 0)
            {
                __pyx_v_idx = 0;
                #ifdef _OPENMP
                #pragma omp parallel
                #endif /* _OPENMP */
                {
                    #ifdef _OPENMP
                    #pragma omp for reduction(+:__pyx_v_tmp_count) reduction(+:__pyx_v_tmp_data) lastprivate(__pyx_v_t) firstprivate(__pyx_v_idx) lastprivate(__pyx_v_idx) lastprivate(__pyx_v_dest)
                    #endif /* _OPENMP */
                    for (__pyx_t_19 = 0; __pyx_t_19 < __pyx_t_5; __pyx_t_19++){
                        {
                            __pyx_v_idx = 0 + 1 * __pyx_t_19;
                            /* Initialize private variables to invalid values */
                            __pyx_v_t = ((long)0xbad0bad0);
                            __pyx_v_dest = ((long)0xbad0bad0);
172: outPos[idx] = bin_edge_min + (0.5 +< double > idx) * bin_width
                            /* "histogram.pyx":172
 * 
 *         for idx in prange(bins):
 *             outPos[idx] = bin_edge_min + (0.5 +< double > idx) * bin_width             # <<<<<<<<<<<<<<
 *             tmp_count = 0.0
 *             tmp_data = 0.0
 */
                            __pyx_t_23 = __pyx_v_idx;
                            *__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outPos.buf, __pyx_t_23, __pyx_bstride_0_outPos) = (__pyx_v_bin_edge_min + ((0.5 + ((double)__pyx_v_idx)) * __pyx_v_bin_width));
173: tmp_count = 0.0
                            /* "histogram.pyx":173
 *         for idx in prange(bins):
 *             outPos[idx] = bin_edge_min + (0.5 +< double > idx) * bin_width
 *             tmp_count = 0.0             # <<<<<<<<<<<<<<
 *             tmp_data = 0.0
 *             for t in range(omp_get_max_threads()):
 */
                            __pyx_v_tmp_count = 0.0;
174: tmp_data = 0.0
                            /* "histogram.pyx":174
 *             outPos[idx] = bin_edge_min + (0.5 +< double > idx) * bin_width
 *             tmp_count = 0.0
 *             tmp_data = 0.0             # <<<<<<<<<<<<<<
 *             for t in range(omp_get_max_threads()):
 *                 dest = t * bins + idx
 */
                            __pyx_v_tmp_data = 0.0;
175: for t in range(omp_get_max_threads()):
                            /* "histogram.pyx":175
 *             tmp_count = 0.0
 *             tmp_data = 0.0
 *             for t in range(omp_get_max_threads()):             # <<<<<<<<<<<<<<
 *                 dest = t * bins + idx
 *                 tmp_count += bigCount[dest]
 */
                            __pyx_t_17 = omp_get_max_threads();
                            for (__pyx_t_24 = 0; __pyx_t_24 < __pyx_t_17; __pyx_t_24+=1) {
                              __pyx_v_t = __pyx_t_24;
176: dest = t * bins + idx
                              /* "histogram.pyx":176
 *             tmp_data = 0.0
 *             for t in range(omp_get_max_threads()):
 *                 dest = t * bins + idx             # <<<<<<<<<<<<<<
 *                 tmp_count += bigCount[dest]
 *                 tmp_data += bigData[dest]
 */
                              __pyx_v_dest = ((__pyx_v_t * __pyx_v_bins) + __pyx_v_idx);
177: tmp_count += bigCount[dest]
                              /* "histogram.pyx":177
 *             for t in range(omp_get_max_threads()):
 *                 dest = t * bins + idx
 *                 tmp_count += bigCount[dest]             # <<<<<<<<<<<<<<
 *                 tmp_data += bigData[dest]
 *             outCount[idx] += tmp_count
 */
                              __pyx_v_tmp_count = (__pyx_v_tmp_count + (__pyx_v_bigCount[__pyx_v_dest]));
178: tmp_data += bigData[dest]
                              /* "histogram.pyx":178
 *                 dest = t * bins + idx
 *                 tmp_count += bigCount[dest]
 *                 tmp_data += bigData[dest]             # <<<<<<<<<<<<<<
 *             outCount[idx] += tmp_count
 *             outData[idx] += tmp_data
 */
                              __pyx_v_tmp_data = (__pyx_v_tmp_data + (__pyx_v_bigData[__pyx_v_dest]));
                            }
179: outCount[idx] += tmp_count
                            /* "histogram.pyx":179
 *                 tmp_count += bigCount[dest]
 *                 tmp_data += bigData[dest]
 *             outCount[idx] += tmp_count             # <<<<<<<<<<<<<<
 *             outData[idx] += tmp_data
 *             if outCount[idx] > epsilon:
 */
                            __pyx_t_24 = __pyx_v_idx;
                            *__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_24, __pyx_bstride_0_outCount) += __pyx_v_tmp_count;
180: outData[idx] += tmp_data
                            /* "histogram.pyx":180
 *                 tmp_data += bigData[dest]
 *             outCount[idx] += tmp_count
 *             outData[idx] += tmp_data             # <<<<<<<<<<<<<<
 *             if outCount[idx] > epsilon:
 *                 outMerge[idx] += tmp_data / tmp_count
 */
                            __pyx_t_25 = __pyx_v_idx;
                            *__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_25, __pyx_bstride_0_outData) += __pyx_v_tmp_data;
181: if outCount[idx] > epsilon:
                            /* "histogram.pyx":181
 *             outCount[idx] += tmp_count
 *             outData[idx] += tmp_data
 *             if outCount[idx] > epsilon:             # <<<<<<<<<<<<<<
 *                 outMerge[idx] += tmp_data / tmp_count
 *             else:
 */
                            __pyx_t_26 = __pyx_v_idx;
                            __pyx_t_4 = ((*__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_26, __pyx_bstride_0_outCount)) > __pyx_v_epsilon);
                            if (__pyx_t_4) {
182: outMerge[idx] += tmp_data / tmp_count
                              /* "histogram.pyx":182
 *             outData[idx] += tmp_data
 *             if outCount[idx] > epsilon:
 *                 outMerge[idx] += tmp_data / tmp_count             # <<<<<<<<<<<<<<
 *             else:
 *                 outMerge[idx] += dummy
 */
                              __pyx_t_27 = __pyx_v_idx;
                              *__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outMerge.buf, __pyx_t_27, __pyx_bstride_0_outMerge) += (__pyx_v_tmp_data / __pyx_v_tmp_count);
                              goto __pyx_L32;
                            }
                            /*else*/ {
183: else:
184: outMerge[idx] += dummy
                              /* "histogram.pyx":184
 *                 outMerge[idx] += tmp_data / tmp_count
 *             else:
 *                 outMerge[idx] += dummy             # <<<<<<<<<<<<<<
 * 
 *     free(bigCount)
 */
                              __pyx_t_28 = __pyx_v_idx;
                              *__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outMerge.buf, __pyx_t_28, __pyx_bstride_0_outMerge) += __pyx_v_dummy;
                            }
                            __pyx_L32:;
                        }
                    }
                }
            }
        }
      }
185:
186: free(bigCount)
/* "histogram.pyx":186 * outMerge[idx] += dummy * * free(bigCount) #<<<<<<<<<<<<<<* free(bigData) * return outPos, outMerge, outData, outCount */ free(__pyx_v_bigCount);
187: free(bigData)
/* "histogram.pyx":187 * * free(bigCount) * free(bigData) #<<<<<<<<<<<<<<* return outPos, outMerge, outData, outCount * */ free(__pyx_v_bigData);
188: return outPos, outMerge, outData, outCount
/* "histogram.pyx":188 * free(bigCount) * free(bigData) * return outPos, outMerge, outData, outCount #<<<<<<<<<<<<<<* * */ __Pyx_XDECREF(__pyx_r); __pyx_t_3 = PyTuple_New(4); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 188; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_3)); __Pyx_INCREF(((PyObject *)__pyx_v_outPos)); PyTuple_SET_ITEM(__pyx_t_3, 0, ((PyObject *)__pyx_v_outPos)); __Pyx_GIVEREF(((PyObject *)__pyx_v_outPos)); __Pyx_INCREF(((PyObject *)__pyx_v_outMerge)); PyTuple_SET_ITEM(__pyx_t_3, 1, ((PyObject *)__pyx_v_outMerge)); __Pyx_GIVEREF(((PyObject *)__pyx_v_outMerge)); __Pyx_INCREF(((PyObject *)__pyx_v_outData)); PyTuple_SET_ITEM(__pyx_t_3, 2, ((PyObject *)__pyx_v_outData)); __Pyx_GIVEREF(((PyObject *)__pyx_v_outData)); __Pyx_INCREF(((PyObject *)__pyx_v_outCount)); PyTuple_SET_ITEM(__pyx_t_3, 3, ((PyObject *)__pyx_v_outCount)); __Pyx_GIVEREF(((PyObject *)__pyx_v_outCount)); __pyx_r = ((PyObject *)__pyx_t_3); __pyx_t_3 = 0; goto __pyx_L0; __pyx_r = Py_None; __Pyx_INCREF(Py_None); goto __pyx_L0; __pyx_L1_error:; __Pyx_XDECREF(__pyx_t_1); __Pyx_XDECREF(__pyx_t_2); __Pyx_XDECREF(__pyx_t_3); __Pyx_XDECREF(__pyx_t_8); { PyObject *__pyx_type, *__pyx_value, *__pyx_tb; __Pyx_ErrFetch(&__pyx_type, &__pyx_value, &__pyx_tb); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outPos); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outMerge); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outCount); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_cdata); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_cpos); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_temp); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outData); __Pyx_ErrRestore(__pyx_type, __pyx_value, __pyx_tb);} __Pyx_AddTraceback("histogram.histogram", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = NULL; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outPos); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outMerge); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outCount); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_cdata); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_cpos); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_temp); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outData); __pyx_L2:; __Pyx_XDECREF((PyObject *)__pyx_v_cpos); __Pyx_XDECREF((PyObject *)__pyx_v_cdata); __Pyx_XDECREF((PyObject *)__pyx_v_outData); __Pyx_XDECREF((PyObject *)__pyx_v_outCount); __Pyx_XDECREF((PyObject *)__pyx_v_outMerge); __Pyx_XDECREF((PyObject *)__pyx_v_outPos); __Pyx_XDECREF((PyObject *)__pyx_v_temp); __Pyx_XGIVEREF(__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }
189:
190:
191: @cython.cdivision(True)
192: @cython.boundscheck(False)
193: @cython.wraparound(False)
194: def histogram2d(numpy.ndarray pos0 not None,
/* "histogram.pyx":194 * @cython.boundscheck(False) * @cython.wraparound(False) * def histogram2d(numpy.ndarray pos0 not None, #<<<<<<<<<<<<<<* numpy.ndarray pos1 not None, * bins not None, */ static PyObject *__pyx_pf_9histogram_1histogram2d(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/ static char __pyx_doc_9histogram_1histogram2d[] = "\n Calculate 2D histogram of pos0,pos1 weighted by weights\n\n @param pos0: 2Theta array\n @param pos1: Chi array\n @param weights: array with intensities\n @param bins: number of output bins int or 2-tuple of int\n @param nthread: maximum number of thread to use. By default: maximum available. \n One can also limit this with OMP_NUM_THREADS environment variable\n\n \n @return I, edges0, edges1, weighted histogram(2D), unweighted histogram (2D)\n "; static PyMethodDef __pyx_mdef_9histogram_1histogram2d = {__Pyx_NAMESTR("histogram2d"), (PyCFunction)__pyx_pf_9histogram_1histogram2d, METH_VARARGS|METH_KEYWORDS, __Pyx_DOCSTR(__pyx_doc_9histogram_1histogram2d)}; static PyObject *__pyx_pf_9histogram_1histogram2d(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) { PyArrayObject *__pyx_v_pos0 = 0; PyArrayObject *__pyx_v_pos1 = 0; PyObject *__pyx_v_bins = 0; PyArrayObject *__pyx_v_weights = 0; PyObject *__pyx_v_split = 0; PyObject *__pyx_v_nthread = 0; double __pyx_v_dummy; long __pyx_v_bin0; long __pyx_v_bin1; long __pyx_v_i; long __pyx_v_j; long __pyx_v_b0; long __pyx_v_b1; long __pyx_v_size; int __pyx_v_csplit; PyArrayObject *__pyx_v_cpos0 = 0; PyArrayObject *__pyx_v_cpos1 = 0; PyArrayObject *__pyx_v_data = 0; PyArrayObject *__pyx_v_outData = 0; PyArrayObject *__pyx_v_outCount = 0; PyArrayObject *__pyx_v_outMerge = 0; PyArrayObject *__pyx_v_edges0 = 0; PyArrayObject *__pyx_v_edges1 = 0; double __pyx_v_min0; double __pyx_v_max0; double __pyx_v_min1; double __pyx_v_max1; double __pyx_v_idp0; double __pyx_v_idp1; double __pyx_v_dbin0; double __pyx_v_dbin1; double __pyx_v_fbin0; double __pyx_v_fbin1; double __pyx_v_p0; double __pyx_v_p1; double __pyx_v_d; double __pyx_v_rest; double __pyx_v_delta0l; double __pyx_v_delta0r; double __pyx_v_delta1l; double __pyx_v_delta1r; double __pyx_v_area; Py_buffer __pyx_bstruct_edges0; Py_ssize_t __pyx_bstride_0_edges0 = 0; Py_ssize_t __pyx_bshape_0_edges0 = 0; Py_buffer __pyx_bstruct_edges1; Py_ssize_t __pyx_bstride_0_edges1 = 0; Py_ssize_t __pyx_bshape_0_edges1 = 0; Py_buffer __pyx_bstruct_outMerge; Py_ssize_t __pyx_bstride_0_outMerge = 0; Py_ssize_t __pyx_bstride_1_outMerge = 0; Py_ssize_t __pyx_bshape_0_outMerge = 0; Py_ssize_t __pyx_bshape_1_outMerge = 0; Py_buffer __pyx_bstruct_outCount; Py_ssize_t __pyx_bstride_0_outCount = 0; Py_ssize_t __pyx_bstride_1_outCount = 0; Py_ssize_t __pyx_bshape_0_outCount = 0; Py_ssize_t __pyx_bshape_1_outCount = 0; Py_buffer __pyx_bstruct_data; Py_ssize_t __pyx_bstride_0_data = 0; Py_ssize_t __pyx_bshape_0_data = 0; Py_buffer __pyx_bstruct_cpos1; Py_ssize_t __pyx_bstride_0_cpos1 = 0; Py_ssize_t __pyx_bshape_0_cpos1 = 0; Py_buffer __pyx_bstruct_cpos0; Py_ssize_t __pyx_bstride_0_cpos0 = 0; Py_ssize_t __pyx_bshape_0_cpos0 = 0; Py_buffer __pyx_bstruct_outData; Py_ssize_t __pyx_bstride_0_outData = 0; Py_ssize_t __pyx_bstride_1_outData = 0; Py_ssize_t __pyx_bshape_0_outData = 0; Py_ssize_t __pyx_bshape_1_outData = 0; PyObject *__pyx_r = NULL; static PyObject **__pyx_pyargnames[] = {&__pyx_n_s__pos0,&__pyx_n_s__pos1,&__pyx_n_s__bins,&__pyx_n_s__weights,&__pyx_n_s__split,&__pyx_n_s__nthread,&__pyx_n_s__dummy,0}; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("histogram2d"); __pyx_self = __pyx_self; { PyObject* values[7] = {0,0,0,0,0,0,0}; values[4] = __pyx_k_3; /* "histogram.pyx":194 * @cython.boundscheck(False) * @cython.wraparound(False) * def histogram2d(numpy.ndarray pos0 not None, #<<<<<<<<<<<<<<* numpy.ndarray pos1 not None, * bins not None, */ __pyx_k_tuple_22 = PyTuple_New(43); if (unlikely(!__pyx_k_tuple_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_k_tuple_22)); __Pyx_INCREF(((PyObject *)__pyx_n_s__pos0)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 0, ((PyObject *)__pyx_n_s__pos0)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__pos0)); __Pyx_INCREF(((PyObject *)__pyx_n_s__pos1)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 1, ((PyObject *)__pyx_n_s__pos1)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__pos1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__bins)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 2, ((PyObject *)__pyx_n_s__bins)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__bins)); __Pyx_INCREF(((PyObject *)__pyx_n_s__weights)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 3, ((PyObject *)__pyx_n_s__weights)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__weights)); __Pyx_INCREF(((PyObject *)__pyx_n_s__split)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 4, ((PyObject *)__pyx_n_s__split)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__split)); __Pyx_INCREF(((PyObject *)__pyx_n_s__nthread)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 5, ((PyObject *)__pyx_n_s__nthread)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__nthread)); __Pyx_INCREF(((PyObject *)__pyx_n_s__dummy)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 6, ((PyObject *)__pyx_n_s__dummy)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__dummy)); __Pyx_INCREF(((PyObject *)__pyx_n_s__bin0)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 7, ((PyObject *)__pyx_n_s__bin0)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__bin0)); __Pyx_INCREF(((PyObject *)__pyx_n_s__bin1)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 8, ((PyObject *)__pyx_n_s__bin1)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__bin1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__i)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 9, ((PyObject *)__pyx_n_s__i)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__i)); __Pyx_INCREF(((PyObject *)__pyx_n_s__j)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 10, ((PyObject *)__pyx_n_s__j)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__j)); __Pyx_INCREF(((PyObject *)__pyx_n_s__b0)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 11, ((PyObject *)__pyx_n_s__b0)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__b0)); __Pyx_INCREF(((PyObject *)__pyx_n_s__b1)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 12, ((PyObject *)__pyx_n_s__b1)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__b1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__size)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 13, ((PyObject *)__pyx_n_s__size)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__size)); __Pyx_INCREF(((PyObject *)__pyx_n_s__csplit)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 14, ((PyObject *)__pyx_n_s__csplit)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__csplit)); __Pyx_INCREF(((PyObject *)__pyx_n_s__cpos0)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 15, ((PyObject *)__pyx_n_s__cpos0)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__cpos0)); __Pyx_INCREF(((PyObject *)__pyx_n_s__cpos1)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 16, ((PyObject *)__pyx_n_s__cpos1)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__cpos1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__data)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 17, ((PyObject *)__pyx_n_s__data)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__data)); __Pyx_INCREF(((PyObject *)__pyx_n_s__outData)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 18, ((PyObject *)__pyx_n_s__outData)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__outData)); __Pyx_INCREF(((PyObject *)__pyx_n_s__outCount)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 19, ((PyObject *)__pyx_n_s__outCount)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__outCount)); __Pyx_INCREF(((PyObject *)__pyx_n_s__outMerge)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 20, ((PyObject *)__pyx_n_s__outMerge)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__outMerge)); __Pyx_INCREF(((PyObject *)__pyx_n_s__edges0)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 21, ((PyObject *)__pyx_n_s__edges0)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__edges0)); __Pyx_INCREF(((PyObject *)__pyx_n_s__edges1)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 22, ((PyObject *)__pyx_n_s__edges1)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__edges1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__min0)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 23, ((PyObject *)__pyx_n_s__min0)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__min0)); __Pyx_INCREF(((PyObject *)__pyx_n_s__max0)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 24, ((PyObject *)__pyx_n_s__max0)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__max0)); __Pyx_INCREF(((PyObject *)__pyx_n_s__min1)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 25, ((PyObject *)__pyx_n_s__min1)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__min1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__max1)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 26, ((PyObject *)__pyx_n_s__max1)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__max1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__idp0)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 27, ((PyObject *)__pyx_n_s__idp0)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__idp0)); __Pyx_INCREF(((PyObject *)__pyx_n_s__idp1)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 28, ((PyObject *)__pyx_n_s__idp1)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__idp1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__dbin0)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 29, ((PyObject *)__pyx_n_s__dbin0)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__dbin0)); __Pyx_INCREF(((PyObject *)__pyx_n_s__dbin1)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 30, ((PyObject *)__pyx_n_s__dbin1)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__dbin1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__fbin0)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 31, ((PyObject *)__pyx_n_s__fbin0)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__fbin0)); __Pyx_INCREF(((PyObject *)__pyx_n_s__fbin1)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 32, ((PyObject *)__pyx_n_s__fbin1)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__fbin1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__p0)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 33, ((PyObject *)__pyx_n_s__p0)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__p0)); __Pyx_INCREF(((PyObject *)__pyx_n_s__p1)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 34, ((PyObject *)__pyx_n_s__p1)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__p1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__d)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 35, ((PyObject *)__pyx_n_s__d)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__d)); __Pyx_INCREF(((PyObject *)__pyx_n_s__rest)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 36, ((PyObject *)__pyx_n_s__rest)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__rest)); __Pyx_INCREF(((PyObject *)__pyx_n_s__delta0l)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 37, ((PyObject *)__pyx_n_s__delta0l)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__delta0l)); __Pyx_INCREF(((PyObject *)__pyx_n_s__delta0r)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 38, ((PyObject *)__pyx_n_s__delta0r)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__delta0r)); __Pyx_INCREF(((PyObject *)__pyx_n_s__delta1l)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 39, ((PyObject *)__pyx_n_s__delta1l)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__delta1l)); __Pyx_INCREF(((PyObject *)__pyx_n_s__delta1r)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 40, ((PyObject *)__pyx_n_s__delta1r)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__delta1r)); __Pyx_INCREF(((PyObject *)__pyx_n_s__aera)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 41, ((PyObject *)__pyx_n_s__aera)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__aera)); __Pyx_INCREF(((PyObject *)__pyx_n_s__area)); PyTuple_SET_ITEM(__pyx_k_tuple_22, 42, ((PyObject *)__pyx_n_s__area)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__area)); __Pyx_GIVEREF(((PyObject *)__pyx_k_tuple_22)); /* "histogram.pyx":194 * @cython.boundscheck(False) * @cython.wraparound(False) * def histogram2d(numpy.ndarray pos0 not None, #<<<<<<<<<<<<<<* numpy.ndarray pos1 not None, * bins not None, */ __pyx_t_1 = PyCFunction_NewEx(&__pyx_mdef_9histogram_1histogram2d, NULL, __pyx_n_s__histogram); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); if (PyObject_SetAttr(__pyx_m, __pyx_n_s__histogram2d, __pyx_t_1)<0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
195: numpy.ndarray pos1 not None,
196: bins not None,
197: numpy.ndarray weights not None,
198: split=True,
/* "histogram.pyx":198 * bins not None, * numpy.ndarray weights not None, * split=True, #<<<<<<<<<<<<<<* nthread=None, * double dummy=0.0): */ __pyx_t_1 = __Pyx_PyBool_FromLong(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 198; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_k_3 = __pyx_t_1; __Pyx_GIVEREF(__pyx_t_1); __pyx_t_1 = 0;
199: nthread=None,
    /* "histogram.pyx":199
 *                 numpy.ndarray weights not None,
 *                 split=True,
 *                 nthread=None,             # <<<<<<<<<<<<<<
 *                 double dummy=0.0):
 *     """
 */
    values[5] = ((PyObject *)Py_None);
    if (unlikely(__pyx_kwds)) {
      Py_ssize_t kw_args;
      switch (PyTuple_GET_SIZE(__pyx_args)) {
        case  7: values[6] = PyTuple_GET_ITEM(__pyx_args, 6);
        case  6: values[5] = PyTuple_GET_ITEM(__pyx_args, 5);
        case  5: values[4] = PyTuple_GET_ITEM(__pyx_args, 4);
        case  4: values[3] = PyTuple_GET_ITEM(__pyx_args, 3);
        case  3: values[2] = PyTuple_GET_ITEM(__pyx_args, 2);
        case  2: values[1] = PyTuple_GET_ITEM(__pyx_args, 1);
        case  1: values[0] = PyTuple_GET_ITEM(__pyx_args, 0);
        case  0: break;
        default: goto __pyx_L5_argtuple_error;
      }
      kw_args = PyDict_Size(__pyx_kwds);
      switch (PyTuple_GET_SIZE(__pyx_args)) {
        case  0:
        values[0] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__pos0);
        if (likely(values[0])) kw_args--;
        else goto __pyx_L5_argtuple_error;
        case  1:
        values[1] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__pos1);
        if (likely(values[1])) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("histogram2d", 0, 4, 7, 1); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        }
        case  2:
        values[2] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__bins);
        if (likely(values[2])) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("histogram2d", 0, 4, 7, 2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        }
        case  3:
        values[3] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__weights);
        if (likely(values[3])) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("histogram2d", 0, 4, 7, 3); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        }
        case  4:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s__split);
          if (value) { values[4] = value; kw_args--; }
        }
        case  5:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s__nthread);
          if (value) { values[5] = value; kw_args--; }
        }
        case  6:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s__dummy);
          if (value) { values[6] = value; kw_args--; }
        }
      }
      if (unlikely(kw_args > 0)) {
        if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_pyargnames, 0, values, PyTuple_GET_SIZE(__pyx_args), "histogram2d") < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
      }
    } else {
      switch (PyTuple_GET_SIZE(__pyx_args)) {
        case  7: values[6] = PyTuple_GET_ITEM(__pyx_args, 6);
        case  6: values[5] = PyTuple_GET_ITEM(__pyx_args, 5);
        case  5: values[4] = PyTuple_GET_ITEM(__pyx_args, 4);
        case  4: values[3] = PyTuple_GET_ITEM(__pyx_args, 3);
        values[2] = PyTuple_GET_ITEM(__pyx_args, 2);
        values[1] = PyTuple_GET_ITEM(__pyx_args, 1);
        values[0] = PyTuple_GET_ITEM(__pyx_args, 0);
        break;
        default: goto __pyx_L5_argtuple_error;
      }
    }
    __pyx_v_pos0 = ((PyArrayObject *)values[0]);
    __pyx_v_pos1 = ((PyArrayObject *)values[1]);
    __pyx_v_bins = values[2];
    __pyx_v_weights = ((PyArrayObject *)values[3]);
    __pyx_v_split = values[4];
    __pyx_v_nthread = values[5];
    if (values[6]) {
      __pyx_v_dummy = __pyx_PyFloat_AsDouble(values[6]); if (unlikely((__pyx_v_dummy == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 200; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
    } else {
200: double dummy=0.0):
      /* "histogram.pyx":200
 *                 split=True,
 *                 nthread=None,
 *                 double dummy=0.0):             # <<<<<<<<<<<<<<
 *     """
 *     Calculate 2D histogram of pos0,pos1 weighted by weights
 */
      __pyx_v_dummy = ((double)0.0);
    }
  }
  goto __pyx_L4_argument_unpacking_done;
  __pyx_L5_argtuple_error:;
  __Pyx_RaiseArgtupleInvalid("histogram2d", 0, 4, 7, PyTuple_GET_SIZE(__pyx_args)); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
  __pyx_L3_error:;
  __Pyx_AddTraceback("histogram.histogram2d", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __Pyx_RefNannyFinishContext();
  return NULL;
  __pyx_L4_argument_unpacking_done:;
  __pyx_bstruct_cpos0.buf = NULL;
  __pyx_bstruct_cpos1.buf = NULL;
  __pyx_bstruct_data.buf = NULL;
  __pyx_bstruct_outData.buf = NULL;
  __pyx_bstruct_outCount.buf = NULL;
  __pyx_bstruct_outMerge.buf = NULL;
  __pyx_bstruct_edges0.buf = NULL;
  __pyx_bstruct_edges1.buf = NULL;
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_pos0), __pyx_ptype_5numpy_ndarray, 0, "pos0", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_pos1), __pyx_ptype_5numpy_ndarray, 0, "pos1", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 195; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(((PyObject *)__pyx_v_bins) == Py_None)) {
    PyErr_Format(PyExc_TypeError, "Argument 'bins' must not be None"); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 196; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_weights), __pyx_ptype_5numpy_ndarray, 0, "weights", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 197; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
201: """
202: Calculate 2D histogram of pos0,pos1 weighted by weights
203:
204: @param pos0: 2Theta array
205: @param pos1: Chi array
206: @param weights: array with intensities
207: @param bins: number of output bins int or 2-tuple of int
208: @param nthread: maximum number of thread to use. By default: maximum available.
209: One can also limit this with OMP_NUM_THREADS environment variable
210:
211:
212: @return I, edges0, edges1, weighted histogram(2D), unweighted histogram (2D)
213: """
214: assert pos0.size == pos1.size
/* "histogram.pyx":214 * @return I, edges0, edges1, weighted histogram(2D), unweighted histogram (2D) * """ * assert pos0.size == pos1.size #<<<<<<<<<<<<<<* # if weights is not No: * assert pos0.size == weights.size */ #ifndef CYTHON_WITHOUT_ASSERTIONS __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_pos0), __pyx_n_s__size); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 214; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_2 = PyObject_GetAttr(((PyObject *)__pyx_v_pos1), __pyx_n_s__size); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 214; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __pyx_t_3 = PyObject_RichCompare(__pyx_t_1, __pyx_t_2, Py_EQ); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 214; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __pyx_t_4 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_4<0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 214; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; if (unlikely(!__pyx_t_4)) { PyErr_SetNone(PyExc_AssertionError); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 214; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } #endif
215: # if weights is not No:
216: assert pos0.size == weights.size
/* "histogram.pyx":216 * assert pos0.size == pos1.size * # if weights is not No: * assert pos0.size == weights.size #<<<<<<<<<<<<<<* cdef long bin0, bin1, i, j, b0, b1 * cdef long size = pos0.size */ #ifndef CYTHON_WITHOUT_ASSERTIONS __pyx_t_3 = PyObject_GetAttr(((PyObject *)__pyx_v_pos0), __pyx_n_s__size); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 216; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_2 = PyObject_GetAttr(((PyObject *)__pyx_v_weights), __pyx_n_s__size); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 216; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __pyx_t_1 = PyObject_RichCompare(__pyx_t_3, __pyx_t_2, Py_EQ); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 216; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __pyx_t_4 = __Pyx_PyObject_IsTrue(__pyx_t_1); if (unlikely(__pyx_t_4<0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 216; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; if (unlikely(!__pyx_t_4)) { PyErr_SetNone(PyExc_AssertionError); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 216; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } #endif
217: cdef long bin0, bin1, i, j, b0, b1
218: cdef long size = pos0.size
/* "histogram.pyx":218 * assert pos0.size == weights.size * cdef long bin0, bin1, i, j, b0, b1 * cdef long size = pos0.size #<<<<<<<<<<<<<<* try: * bin0, bin1 = tuple(bins) */ __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_pos0), __pyx_n_s__size); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 218; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_5 = __Pyx_PyInt_AsLong(__pyx_t_1); if (unlikely((__pyx_t_5 == (long)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 218; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_v_size = __pyx_t_5;
219: try:
/* "histogram.pyx":219 * cdef long bin0, bin1, i, j, b0, b1 * cdef long size = pos0.size * try: #<<<<<<<<<<<<<<* bin0, bin1 = tuple(bins) * except: */ { __Pyx_ExceptionSave(&__pyx_t_6, &__pyx_t_7, &__pyx_t_8); __Pyx_XGOTREF(__pyx_t_6); __Pyx_XGOTREF(__pyx_t_7); __Pyx_XGOTREF(__pyx_t_8); /*try:*/ {
220: bin0, bin1 = tuple(bins)
      /* "histogram.pyx":220
 *     cdef long  size = pos0.size
 *     try:
 *         bin0, bin1 = tuple(bins)             # <<<<<<<<<<<<<<
 *     except:
 *         bin0 = bin1 = < long > bins
 */
      __pyx_t_1 = PyTuple_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 220; __pyx_clineno = __LINE__; goto __pyx_L6_error;}
      __Pyx_GOTREF(((PyObject *)__pyx_t_1));
      __Pyx_INCREF(__pyx_v_bins);
      PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_v_bins);
      __Pyx_GIVEREF(__pyx_v_bins);
      __pyx_t_2 = PyObject_Call(((PyObject *)((PyObject*)(&PyTuple_Type))), ((PyObject *)__pyx_t_1), NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 220; __pyx_clineno = __LINE__; goto __pyx_L6_error;}
      __Pyx_GOTREF(__pyx_t_2);
      __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
      if (likely(PyTuple_CheckExact(__pyx_t_2))) {
        PyObject* sequence = __pyx_t_2;
        if (unlikely(PyTuple_GET_SIZE(sequence) != 2)) {
          if (PyTuple_GET_SIZE(sequence) > 2) __Pyx_RaiseTooManyValuesError(2);
          else __Pyx_RaiseNeedMoreValuesError(PyTuple_GET_SIZE(sequence));
          {__pyx_filename = __pyx_f[0]; __pyx_lineno = 220; __pyx_clineno = __LINE__; goto __pyx_L6_error;}
        }
        __pyx_t_1 = PyTuple_GET_ITEM(sequence, 0); 
        __pyx_t_3 = PyTuple_GET_ITEM(sequence, 1); 
        __Pyx_INCREF(__pyx_t_1);
        __Pyx_INCREF(__pyx_t_3);
        __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
      } else {
        __Pyx_UnpackTupleError(__pyx_t_2, 2);
        {__pyx_filename = __pyx_f[0]; __pyx_lineno = 220; __pyx_clineno = __LINE__; goto __pyx_L6_error;}
      }
      __pyx_t_5 = __Pyx_PyInt_AsLong(__pyx_t_1); if (unlikely((__pyx_t_5 == (long)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 220; __pyx_clineno = __LINE__; goto __pyx_L6_error;}
      __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
      __pyx_t_9 = __Pyx_PyInt_AsLong(__pyx_t_3); if (unlikely((__pyx_t_9 == (long)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 220; __pyx_clineno = __LINE__; goto __pyx_L6_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_v_bin0 = __pyx_t_5;
      __pyx_v_bin1 = __pyx_t_9;
    }
    __Pyx_XDECREF(__pyx_t_6); __pyx_t_6 = 0;
    __Pyx_XDECREF(__pyx_t_7); __pyx_t_7 = 0;
    __Pyx_XDECREF(__pyx_t_8); __pyx_t_8 = 0;
    goto __pyx_L13_try_end;
    __pyx_L6_error:;
    __Pyx_XDECREF(__pyx_t_1); __pyx_t_1 = 0;
    __Pyx_XDECREF(__pyx_t_3); __pyx_t_3 = 0;
    __Pyx_XDECREF(__pyx_t_2); __pyx_t_2 = 0;
221: except:
    /* "histogram.pyx":221
 *     try:
 *         bin0, bin1 = tuple(bins)
 *     except:             # <<<<<<<<<<<<<<
 *         bin0 = bin1 = < long > bins
 *     if bin0 <= 0:
 */
    /*except:*/ {
      __Pyx_AddTraceback("histogram.histogram2d", __pyx_clineno, __pyx_lineno, __pyx_filename);
      if (__Pyx_GetException(&__pyx_t_2, &__pyx_t_3, &__pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 221; __pyx_clineno = __LINE__; goto __pyx_L8_except_error;}
      __Pyx_GOTREF(__pyx_t_2);
      __Pyx_GOTREF(__pyx_t_3);
      __Pyx_GOTREF(__pyx_t_1);
222: bin0 = bin1 = < long > bins
      /* "histogram.pyx":222
 *         bin0, bin1 = tuple(bins)
 *     except:
 *         bin0 = bin1 = < long > bins             # <<<<<<<<<<<<<<
 *     if bin0 <= 0:
 *         bin0 = 1
 */
      __pyx_t_9 = __Pyx_PyInt_AsLong(__pyx_v_bins); if (unlikely((__pyx_t_9 == (long)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 222; __pyx_clineno = __LINE__; goto __pyx_L8_except_error;}
      __pyx_v_bin0 = ((long)__pyx_t_9);
      __pyx_v_bin1 = ((long)__pyx_t_9);
      __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
      goto __pyx_L7_exception_handled;
    }
    __pyx_L8_except_error:;
    __Pyx_XGIVEREF(__pyx_t_6);
    __Pyx_XGIVEREF(__pyx_t_7);
    __Pyx_XGIVEREF(__pyx_t_8);
    __Pyx_ExceptionReset(__pyx_t_6, __pyx_t_7, __pyx_t_8);
    goto __pyx_L1_error;
    __pyx_L7_exception_handled:;
    __Pyx_XGIVEREF(__pyx_t_6);
    __Pyx_XGIVEREF(__pyx_t_7);
    __Pyx_XGIVEREF(__pyx_t_8);
    __Pyx_ExceptionReset(__pyx_t_6, __pyx_t_7, __pyx_t_8);
    __pyx_L13_try_end:;
  }
223: if bin0 <= 0:
/* "histogram.pyx":223 * except: * bin0 = bin1 =<long > bins * if bin0<= 0: #<<<<<<<<<<<<<<* bin0 = 1 * if bin1<= 0: */ __pyx_t_4 = (__pyx_v_bin0<= 0); if (__pyx_t_4) {
224: bin0 = 1
    /* "histogram.pyx":224
 *         bin0 = bin1 = < long > bins
 *     if bin0 <= 0:
 *         bin0 = 1             # <<<<<<<<<<<<<<
 *     if bin1 <= 0:
 *         bin1 = 1
 */
    __pyx_v_bin0 = 1;
    goto __pyx_L16;
  }
  __pyx_L16:;
225: if bin1 <= 0:
/* "histogram.pyx":225 * if bin0<= 0: * bin0 = 1 * if bin1<= 0: #<<<<<<<<<<<<<<* bin1 = 1 * cdef int csplit = split */ __pyx_t_4 = (__pyx_v_bin1<= 0); if (__pyx_t_4) {
226: bin1 = 1
    /* "histogram.pyx":226
 *         bin0 = 1
 *     if bin1 <= 0:
 *         bin1 = 1             # <<<<<<<<<<<<<<
 *     cdef int csplit = split
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos0 = pos0.astype("float64").flatten()
 */
    __pyx_v_bin1 = 1;
    goto __pyx_L17;
  }
  __pyx_L17:;
227: cdef int csplit = split
/* "histogram.pyx":227 * if bin1<= 0: * bin1 = 1 * cdef int csplit = split #<<<<<<<<<<<<<<* cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos0 = pos0.astype("float64").flatten() * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos1 = pos1.astype("float64").flatten() */ __pyx_t_10 = __Pyx_PyInt_AsInt(__pyx_v_split); if (unlikely((__pyx_t_10 == (int)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 227; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_v_csplit = __pyx_t_10;
 228:     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos0 = pos0.astype("float64").flatten()
  /* "histogram.pyx":228
 *         bin1 = 1
 *     cdef int csplit = split
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos0 = pos0.astype("float64").flatten()             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos1 = pos1.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] data = weights.astype("float64").flatten()
 */
  __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_pos0), __pyx_n_s__astype); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 228; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_3 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_k_tuple_4), NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 228; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_1 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__flatten); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 228; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_t_3 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 228; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  if (!(likely(((__pyx_t_3) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_3, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 228; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_11 = ((PyArrayObject *)__pyx_t_3);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_cpos0, (PyObject*)__pyx_t_11, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {
      __pyx_v_cpos0 = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_cpos0.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 228; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_bstride_0_cpos0 = __pyx_bstruct_cpos0.strides[0];
      __pyx_bshape_0_cpos0 = __pyx_bstruct_cpos0.shape[0];
    }
  }
  __pyx_t_11 = 0;
  __pyx_v_cpos0 = ((PyArrayObject *)__pyx_t_3);
  __pyx_t_3 = 0;
  /* "histogram.pyx":228
 *         bin1 = 1
 *     cdef int csplit = split
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos0 = pos0.astype("float64").flatten()             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos1 = pos1.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] data = weights.astype("float64").flatten()
 */
  __pyx_k_tuple_4 = PyTuple_New(1); if (unlikely(!__pyx_k_tuple_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 228; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_k_tuple_4));
  __Pyx_INCREF(((PyObject *)__pyx_n_s__float64));
  PyTuple_SET_ITEM(__pyx_k_tuple_4, 0, ((PyObject *)__pyx_n_s__float64));
  __Pyx_GIVEREF(((PyObject *)__pyx_n_s__float64));
  __Pyx_GIVEREF(((PyObject *)__pyx_k_tuple_4));
 229:     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos1 = pos1.astype("float64").flatten()
  /* "histogram.pyx":229
 *     cdef int csplit = split
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos0 = pos0.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos1 = pos1.astype("float64").flatten()             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] data = weights.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")
 */
  __pyx_t_3 = PyObject_GetAttr(((PyObject *)__pyx_v_pos1), __pyx_n_s__astype); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 229; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_1 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_k_tuple_5), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 229; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_t_3 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__flatten); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 229; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_1 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 229; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  if (!(likely(((__pyx_t_1) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_1, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 229; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_12 = ((PyArrayObject *)__pyx_t_1);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_cpos1, (PyObject*)__pyx_t_12, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {
      __pyx_v_cpos1 = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_cpos1.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 229; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_bstride_0_cpos1 = __pyx_bstruct_cpos1.strides[0];
      __pyx_bshape_0_cpos1 = __pyx_bstruct_cpos1.shape[0];
    }
  }
  __pyx_t_12 = 0;
  __pyx_v_cpos1 = ((PyArrayObject *)__pyx_t_1);
  __pyx_t_1 = 0;
  /* "histogram.pyx":229
 *     cdef int csplit = split
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos0 = pos0.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos1 = pos1.astype("float64").flatten()             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] data = weights.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")
 */
  __pyx_k_tuple_5 = PyTuple_New(1); if (unlikely(!__pyx_k_tuple_5)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 229; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_k_tuple_5));
  __Pyx_INCREF(((PyObject *)__pyx_n_s__float64));
  PyTuple_SET_ITEM(__pyx_k_tuple_5, 0, ((PyObject *)__pyx_n_s__float64));
  __Pyx_GIVEREF(((PyObject *)__pyx_n_s__float64));
  __Pyx_GIVEREF(((PyObject *)__pyx_k_tuple_5));
 230:     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] data = weights.astype("float64").flatten()
  /* "histogram.pyx":230
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos0 = pos0.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos1 = pos1.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] data = weights.astype("float64").flatten()             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64")
 */
  __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_weights), __pyx_n_s__astype); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 230; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_3 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_k_tuple_6), NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 230; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_1 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__flatten); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 230; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_t_3 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 230; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  if (!(likely(((__pyx_t_3) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_3, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 230; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_13 = ((PyArrayObject *)__pyx_t_3);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_data, (PyObject*)__pyx_t_13, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {
      __pyx_v_data = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_data.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 230; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_bstride_0_data = __pyx_bstruct_data.strides[0];
      __pyx_bshape_0_data = __pyx_bstruct_data.shape[0];
    }
  }
  __pyx_t_13 = 0;
  __pyx_v_data = ((PyArrayObject *)__pyx_t_3);
  __pyx_t_3 = 0;
  /* "histogram.pyx":230
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos0 = pos0.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos1 = pos1.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] data = weights.astype("float64").flatten()             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64")
 */
  __pyx_k_tuple_6 = PyTuple_New(1); if (unlikely(!__pyx_k_tuple_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 230; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_k_tuple_6));
  __Pyx_INCREF(((PyObject *)__pyx_n_s__float64));
  PyTuple_SET_ITEM(__pyx_k_tuple_6, 0, ((PyObject *)__pyx_n_s__float64));
  __Pyx_GIVEREF(((PyObject *)__pyx_n_s__float64));
  __Pyx_GIVEREF(((PyObject *)__pyx_k_tuple_6));
231: cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")
  /* "histogram.pyx":231
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] cpos1 = pos1.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] data = weights.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outMerge = numpy.zeros((bin0, bin1), dtype="float64")
 */
  __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 231; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_1 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__zeros); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 231; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_t_3 = PyInt_FromLong(__pyx_v_bin0); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 231; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_2 = PyInt_FromLong(__pyx_v_bin1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 231; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __pyx_t_14 = PyTuple_New(2); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 231; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_14));
  PyTuple_SET_ITEM(__pyx_t_14, 0, __pyx_t_3);
  __Pyx_GIVEREF(__pyx_t_3);
  PyTuple_SET_ITEM(__pyx_t_14, 1, __pyx_t_2);
  __Pyx_GIVEREF(__pyx_t_2);
  __pyx_t_3 = 0;
  __pyx_t_2 = 0;
  __pyx_t_2 = PyTuple_New(1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 231; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_2));
  PyTuple_SET_ITEM(__pyx_t_2, 0, ((PyObject *)__pyx_t_14));
  __Pyx_GIVEREF(((PyObject *)__pyx_t_14));
  __pyx_t_14 = 0;
  __pyx_t_14 = PyDict_New(); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 231; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_14));
  if (PyDict_SetItem(__pyx_t_14, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64)) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 231; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_3 = PyEval_CallObjectWithKeywords(__pyx_t_1, ((PyObject *)__pyx_t_2), ((PyObject *)__pyx_t_14)); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 231; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_14)); __pyx_t_14 = 0;
  if (!(likely(((__pyx_t_3) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_3, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 231; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_15 = ((PyArrayObject *)__pyx_t_3);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_outData, (PyObject*)__pyx_t_15, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 2, 0, __pyx_stack) == -1)) {
      __pyx_v_outData = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_outData.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 231; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_bstride_0_outData = __pyx_bstruct_outData.strides[0]; __pyx_bstride_1_outData = __pyx_bstruct_outData.strides[1];
      __pyx_bshape_0_outData = __pyx_bstruct_outData.shape[0]; __pyx_bshape_1_outData = __pyx_bstruct_outData.shape[1];
    }
  }
  __pyx_t_15 = 0;
  __pyx_v_outData = ((PyArrayObject *)__pyx_t_3);
  __pyx_t_3 = 0;
232: cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64")
  /* "histogram.pyx":232
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] data = weights.astype("float64").flatten()
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64")             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outMerge = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges0 = numpy.zeros(bin0, dtype="float64")
 */
  __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 232; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_14 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__zeros); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 232; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_14);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_t_3 = PyInt_FromLong(__pyx_v_bin0); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 232; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_2 = PyInt_FromLong(__pyx_v_bin1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 232; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __pyx_t_1 = PyTuple_New(2); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 232; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_1));
  PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_t_3);
  __Pyx_GIVEREF(__pyx_t_3);
  PyTuple_SET_ITEM(__pyx_t_1, 1, __pyx_t_2);
  __Pyx_GIVEREF(__pyx_t_2);
  __pyx_t_3 = 0;
  __pyx_t_2 = 0;
  __pyx_t_2 = PyTuple_New(1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 232; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_2));
  PyTuple_SET_ITEM(__pyx_t_2, 0, ((PyObject *)__pyx_t_1));
  __Pyx_GIVEREF(((PyObject *)__pyx_t_1));
  __pyx_t_1 = 0;
  __pyx_t_1 = PyDict_New(); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 232; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_1));
  if (PyDict_SetItem(__pyx_t_1, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64)) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 232; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_3 = PyEval_CallObjectWithKeywords(__pyx_t_14, ((PyObject *)__pyx_t_2), ((PyObject *)__pyx_t_1)); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 232; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_14); __pyx_t_14 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
  if (!(likely(((__pyx_t_3) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_3, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 232; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_16 = ((PyArrayObject *)__pyx_t_3);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_outCount, (PyObject*)__pyx_t_16, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 2, 0, __pyx_stack) == -1)) {
      __pyx_v_outCount = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_outCount.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 232; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_bstride_0_outCount = __pyx_bstruct_outCount.strides[0]; __pyx_bstride_1_outCount = __pyx_bstruct_outCount.strides[1];
      __pyx_bshape_0_outCount = __pyx_bstruct_outCount.shape[0]; __pyx_bshape_1_outCount = __pyx_bstruct_outCount.shape[1];
    }
  }
  __pyx_t_16 = 0;
  __pyx_v_outCount = ((PyArrayObject *)__pyx_t_3);
  __pyx_t_3 = 0;
233: cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outMerge = numpy.zeros((bin0, bin1), dtype="float64")
/* "histogram.pyx":233 * cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outMerge = numpy.zeros((bin0, bin1), dtype="float64") #<<<<<<<<<<<<<<* cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges0 = numpy.zeros(bin0, dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges1 = numpy.zeros(bin1, dtype="float64") */ __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 233; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_1 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__zeros); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 233; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = PyInt_FromLong(__pyx_v_bin0); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 233; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_2 = PyInt_FromLong(__pyx_v_bin1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 233; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __pyx_t_14 = PyTuple_New(2); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 233; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_14)); PyTuple_SET_ITEM(__pyx_t_14, 0, __pyx_t_3); __Pyx_GIVEREF(__pyx_t_3); PyTuple_SET_ITEM(__pyx_t_14, 1, __pyx_t_2); __Pyx_GIVEREF(__pyx_t_2); __pyx_t_3 = 0; __pyx_t_2 = 0; __pyx_t_2 = PyTuple_New(1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 233; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_2)); PyTuple_SET_ITEM(__pyx_t_2, 0, ((PyObject *)__pyx_t_14)); __Pyx_GIVEREF(((PyObject *)__pyx_t_14)); __pyx_t_14 = 0; __pyx_t_14 = PyDict_New(); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 233; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_14)); if (PyDict_SetItem(__pyx_t_14, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64))<0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 233; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_3 = PyEval_CallObjectWithKeywords(__pyx_t_1, ((PyObject *)__pyx_t_2), ((PyObject *)__pyx_t_14)); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 233; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_14)); __pyx_t_14 = 0; if (!(likely(((__pyx_t_3) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_3, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 233; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_17 = ((PyArrayObject *)__pyx_t_3); { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_outMerge, (PyObject*)__pyx_t_17, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 2, 0, __pyx_stack) == -1)) { __pyx_v_outMerge = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_outMerge.buf = NULL; {__pyx_filename = __pyx_f[0]; __pyx_lineno = 233; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } else {__pyx_bstride_0_outMerge = __pyx_bstruct_outMerge.strides[0]; __pyx_bstride_1_outMerge = __pyx_bstruct_outMerge.strides[1]; __pyx_bshape_0_outMerge = __pyx_bstruct_outMerge.shape[0]; __pyx_bshape_1_outMerge = __pyx_bstruct_outMerge.shape[1]; } } __pyx_t_17 = 0; __pyx_v_outMerge = ((PyArrayObject *)__pyx_t_3); __pyx_t_3 = 0;
234: cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges0 = numpy.zeros(bin0, dtype="float64")
/* "histogram.pyx":234 * cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outMerge = numpy.zeros((bin0, bin1), dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges0 = numpy.zeros(bin0, dtype="float64") #<<<<<<<<<<<<<<* cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges1 = numpy.zeros(bin1, dtype="float64") * cdef double min0 = pos0.min() */ __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 234; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_14 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__zeros); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 234; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_14); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = PyInt_FromLong(__pyx_v_bin0); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 234; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_2 = PyTuple_New(1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 234; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_2)); PyTuple_SET_ITEM(__pyx_t_2, 0, __pyx_t_3); __Pyx_GIVEREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = PyDict_New(); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 234; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_3)); if (PyDict_SetItem(__pyx_t_3, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64))<0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 234; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_1 = PyEval_CallObjectWithKeywords(__pyx_t_14, ((PyObject *)__pyx_t_2), ((PyObject *)__pyx_t_3)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 234; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_14); __pyx_t_14 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0; if (!(likely(((__pyx_t_1) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_1, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 234; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_18 = ((PyArrayObject *)__pyx_t_1); { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_edges0, (PyObject*)__pyx_t_18, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) { __pyx_v_edges0 = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_edges0.buf = NULL; {__pyx_filename = __pyx_f[0]; __pyx_lineno = 234; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } else {__pyx_bstride_0_edges0 = __pyx_bstruct_edges0.strides[0]; __pyx_bshape_0_edges0 = __pyx_bstruct_edges0.shape[0]; } } __pyx_t_18 = 0; __pyx_v_edges0 = ((PyArrayObject *)__pyx_t_1); __pyx_t_1 = 0;
235: cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges1 = numpy.zeros(bin1, dtype="float64")
/* "histogram.pyx":235 * cdef numpy.ndarray[DTYPE_float64_t, ndim = 2] outMerge = numpy.zeros((bin0, bin1), dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges0 = numpy.zeros(bin0, dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges1 = numpy.zeros(bin1, dtype="float64") #<<<<<<<<<<<<<<* cdef double min0 = pos0.min() * cdef double max0 = pos0.max() */ __pyx_t_1 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 235; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_3 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__zeros); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 235; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = PyInt_FromLong(__pyx_v_bin1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 235; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_2 = PyTuple_New(1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 235; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_2)); PyTuple_SET_ITEM(__pyx_t_2, 0, __pyx_t_1); __Pyx_GIVEREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = PyDict_New(); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 235; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_1)); if (PyDict_SetItem(__pyx_t_1, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64))<0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 235; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_14 = PyEval_CallObjectWithKeywords(__pyx_t_3, ((PyObject *)__pyx_t_2), ((PyObject *)__pyx_t_1)); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 235; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_14); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0; if (!(likely(((__pyx_t_14) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_14, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 235; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_19 = ((PyArrayObject *)__pyx_t_14); { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_edges1, (PyObject*)__pyx_t_19, &__Pyx_TypeInfo_nn___pyx_t_9histogram_DTYPE_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) { __pyx_v_edges1 = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_edges1.buf = NULL; {__pyx_filename = __pyx_f[0]; __pyx_lineno = 235; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } else {__pyx_bstride_0_edges1 = __pyx_bstruct_edges1.strides[0]; __pyx_bshape_0_edges1 = __pyx_bstruct_edges1.shape[0]; } } __pyx_t_19 = 0; __pyx_v_edges1 = ((PyArrayObject *)__pyx_t_14); __pyx_t_14 = 0;
236: cdef double min0 = pos0.min()
/* "histogram.pyx":236 * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges0 = numpy.zeros(bin0, dtype="float64") * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges1 = numpy.zeros(bin1, dtype="float64") * cdef double min0 = pos0.min() #<<<<<<<<<<<<<<* cdef double max0 = pos0.max() * cdef double min1 = pos1.min() */ __pyx_t_14 = PyObject_GetAttr(((PyObject *)__pyx_v_pos0), __pyx_n_s__min); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 236; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_14); __pyx_t_1 = PyObject_Call(__pyx_t_14, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 236; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_14); __pyx_t_14 = 0; __pyx_t_20 = __pyx_PyFloat_AsDouble(__pyx_t_1); if (unlikely((__pyx_t_20 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 236; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_v_min0 = __pyx_t_20;
237: cdef double max0 = pos0.max()
/* "histogram.pyx":237 * cdef numpy.ndarray[DTYPE_float64_t, ndim = 1] edges1 = numpy.zeros(bin1, dtype="float64") * cdef double min0 = pos0.min() * cdef double max0 = pos0.max() #<<<<<<<<<<<<<<* cdef double min1 = pos1.min() * cdef double max1 = pos1.max() */ __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_pos0), __pyx_n_s__max); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 237; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_14 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 237; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_14); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_20 = __pyx_PyFloat_AsDouble(__pyx_t_14); if (unlikely((__pyx_t_20 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 237; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_14); __pyx_t_14 = 0; __pyx_v_max0 = __pyx_t_20;
238: cdef double min1 = pos1.min()
/* "histogram.pyx":238 * cdef double min0 = pos0.min() * cdef double max0 = pos0.max() * cdef double min1 = pos1.min() #<<<<<<<<<<<<<<* cdef double max1 = pos1.max() * cdef double idp0 = (<double > bin0) / (max0 - min0) */ __pyx_t_14 = PyObject_GetAttr(((PyObject *)__pyx_v_pos1), __pyx_n_s__min); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 238; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_14); __pyx_t_1 = PyObject_Call(__pyx_t_14, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 238; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_14); __pyx_t_14 = 0; __pyx_t_20 = __pyx_PyFloat_AsDouble(__pyx_t_1); if (unlikely((__pyx_t_20 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 238; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_v_min1 = __pyx_t_20;
239: cdef double max1 = pos1.max()
/* "histogram.pyx":239 * cdef double max0 = pos0.max() * cdef double min1 = pos1.min() * cdef double max1 = pos1.max() #<<<<<<<<<<<<<<* cdef double idp0 = (<double > bin0) / (max0 - min0) * cdef double idp1 = (<double > bin1) / (max1 - min1) */ __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_pos1), __pyx_n_s__max); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 239; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_14 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 239; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_14); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_20 = __pyx_PyFloat_AsDouble(__pyx_t_14); if (unlikely((__pyx_t_20 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 239; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_14); __pyx_t_14 = 0; __pyx_v_max1 = __pyx_t_20;
240: cdef double idp0 = (< double > bin0) / (max0 - min0)
/* "histogram.pyx":240 * cdef double min1 = pos1.min() * cdef double max1 = pos1.max() * cdef double idp0 = (<double > bin0) / (max0 - min0) #<<<<<<<<<<<<<<* cdef double idp1 = (<double > bin1) / (max1 - min1) * cdef double dbin0 = 0.5, dbin1 = 0.5 */ __pyx_v_idp0 = (((double)__pyx_v_bin0) / (__pyx_v_max0 - __pyx_v_min0));
241: cdef double idp1 = (< double > bin1) / (max1 - min1)
/* "histogram.pyx":241 * cdef double max1 = pos1.max() * cdef double idp0 = (<double > bin0) / (max0 - min0) * cdef double idp1 = (<double > bin1) / (max1 - min1) #<<<<<<<<<<<<<<* cdef double dbin0 = 0.5, dbin1 = 0.5 * cdef double fbin0, fbin1, p0, p1, d, rest, delta0l, delta0r, delta1l, delta1r, aera */ __pyx_v_idp1 = (((double)__pyx_v_bin1) / (__pyx_v_max1 - __pyx_v_min1));
242: cdef double dbin0 = 0.5, dbin1 = 0.5
/* "histogram.pyx":242 * cdef double idp0 = (<double > bin0) / (max0 - min0) * cdef double idp1 = (<double > bin1) / (max1 - min1) * cdef double dbin0 = 0.5, dbin1 = 0.5 #<<<<<<<<<<<<<<* cdef double fbin0, fbin1, p0, p1, d, rest, delta0l, delta0r, delta1l, delta1r, aera * if nthread is not None: */ __pyx_v_dbin0 = 0.5; __pyx_v_dbin1 = 0.5;
243: cdef double fbin0, fbin1, p0, p1, d, rest, delta0l, delta0r, delta1l, delta1r, aera
244: if nthread is not None:
/* "histogram.pyx":244 * cdef double dbin0 = 0.5, dbin1 = 0.5 * cdef double fbin0, fbin1, p0, p1, d, rest, delta0l, delta0r, delta1l, delta1r, aera * if nthread is not None: #<<<<<<<<<<<<<<* if isinstance(nthread, int) and (nthread > 0): * omp_set_num_threads(<int > nthread) */ __pyx_t_4 = (__pyx_v_nthread != Py_None); if (__pyx_t_4) {
245: if isinstance(nthread, int) and (nthread > 0):
    /* "histogram.pyx":245
 *     cdef double fbin0, fbin1, p0, p1, d, rest, delta0l, delta0r, delta1l, delta1r, aera
 *     if nthread is not None:
 *         if isinstance(nthread, int) and (nthread > 0):             # <<<<<<<<<<<<<<
 *             omp_set_num_threads(< int > nthread)
 *     with nogil:
 */
    __pyx_t_14 = ((PyObject *)((PyObject*)(&PyInt_Type)));
    __Pyx_INCREF(__pyx_t_14);
    __pyx_t_4 = __Pyx_TypeCheck(__pyx_v_nthread, __pyx_t_14); 
    __Pyx_DECREF(__pyx_t_14); __pyx_t_14 = 0;
    if (__pyx_t_4) {
      __pyx_t_14 = PyObject_RichCompare(__pyx_v_nthread, __pyx_int_0, Py_GT); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 245; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_14);
      __pyx_t_21 = __Pyx_PyObject_IsTrue(__pyx_t_14); if (unlikely(__pyx_t_21 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 245; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_14); __pyx_t_14 = 0;
      __pyx_t_22 = __pyx_t_21;
    } else {
      __pyx_t_22 = __pyx_t_4;
    }
    if (__pyx_t_22) {
246: omp_set_num_threads(< int > nthread)
      /* "histogram.pyx":246
 *     if nthread is not None:
 *         if isinstance(nthread, int) and (nthread > 0):
 *             omp_set_num_threads(< int > nthread)             # <<<<<<<<<<<<<<
 *     with nogil:
 *         for i in prange(bin0):
 */
      __pyx_t_10 = __Pyx_PyInt_AsInt(__pyx_v_nthread); if (unlikely((__pyx_t_10 == (int)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 246; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      omp_set_num_threads(((int)__pyx_t_10));
      goto __pyx_L19;
    }
    __pyx_L19:;
    goto __pyx_L18;
  }
  __pyx_L18:;
247: with nogil:
/* "histogram.pyx":247 * if isinstance(nthread, int) and (nthread > 0): * omp_set_num_threads(<int > nthread) * with nogil: #<<<<<<<<<<<<<<* for i in prange(bin0): * edges0[i] = min0 + (0.5 +<double > i) / idp0 */ { #ifdef WITH_THREAD PyThreadState *_save = NULL; #endif Py_UNBLOCK_THREADS /*try:*/ { /* "histogram.pyx":247 * if isinstance(nthread, int) and (nthread > 0): * omp_set_num_threads(<int > nthread) * with nogil: #<<<<<<<<<<<<<<* for i in prange(bin0): * edges0[i] = min0 + (0.5 +<double > i) / idp0 */ /*finally:*/ { Py_BLOCK_THREADS } }
248: for i in prange(bin0):
        /* "histogram.pyx":248
 *             omp_set_num_threads(< int > nthread)
 *     with nogil:
 *         for i in prange(bin0):             # <<<<<<<<<<<<<<
 *             edges0[i] = min0 + (0.5 +< double > i) / idp0
 *         for i in prange(bin1):
 */
        __pyx_t_9 = __pyx_v_bin0;
        if (1 == 0) abort();
        {
            __pyx_t_23 = (__pyx_t_9 - 0) / 1;
            if (__pyx_t_23 > 0)
            {
                __pyx_v_i = 0;
                #ifdef _OPENMP
                #pragma omp parallel
                #endif /* _OPENMP */
                {
                    #ifdef _OPENMP
                    #pragma omp for firstprivate(__pyx_v_i) lastprivate(__pyx_v_i)
                    #endif /* _OPENMP */
                    for (__pyx_t_5 = 0; __pyx_t_5 < __pyx_t_23; __pyx_t_5++){
                        {
                            __pyx_v_i = 0 + 1 * __pyx_t_5;
249: edges0[i] = min0 + (0.5 +< double > i) / idp0
                            /* "histogram.pyx":249
 *     with nogil:
 *         for i in prange(bin0):
 *             edges0[i] = min0 + (0.5 +< double > i) / idp0             # <<<<<<<<<<<<<<
 *         for i in prange(bin1):
 *             edges1[i] = min1 + (0.5 +< double > i) / idp1
 */
                            __pyx_t_24 = __pyx_v_i;
                            *__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_edges0.buf, __pyx_t_24, __pyx_bstride_0_edges0) = (__pyx_v_min0 + ((0.5 + ((double)__pyx_v_i)) / __pyx_v_idp0));
                        }
                    }
                }
            }
        }
250: for i in prange(bin1):
        /* "histogram.pyx":250
 *         for i in prange(bin0):
 *             edges0[i] = min0 + (0.5 +< double > i) / idp0
 *         for i in prange(bin1):             # <<<<<<<<<<<<<<
 *             edges1[i] = min1 + (0.5 +< double > i) / idp1
 *         for i in range(size):
 */
        __pyx_t_23 = __pyx_v_bin1;
        if (1 == 0) abort();
        {
            __pyx_t_9 = (__pyx_t_23 - 0) / 1;
            if (__pyx_t_9 > 0)
            {
                __pyx_v_i = 0;
                #ifdef _OPENMP
                #pragma omp parallel
                #endif /* _OPENMP */
                {
                    #ifdef _OPENMP
                    #pragma omp for firstprivate(__pyx_v_i) lastprivate(__pyx_v_i)
                    #endif /* _OPENMP */
                    for (__pyx_t_5 = 0; __pyx_t_5 < __pyx_t_9; __pyx_t_5++){
                        {
                            __pyx_v_i = 0 + 1 * __pyx_t_5;
251: edges1[i] = min1 + (0.5 +< double > i) / idp1
                            /* "histogram.pyx":251
 *             edges0[i] = min0 + (0.5 +< double > i) / idp0
 *         for i in prange(bin1):
 *             edges1[i] = min1 + (0.5 +< double > i) / idp1             # <<<<<<<<<<<<<<
 *         for i in range(size):
 *             p0 = cpos0[i]
 */
                            __pyx_t_25 = __pyx_v_i;
                            *__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_edges1.buf, __pyx_t_25, __pyx_bstride_0_edges1) = (__pyx_v_min1 + ((0.5 + ((double)__pyx_v_i)) / __pyx_v_idp1));
                        }
                    }
                }
            }
        }
252: for i in range(size):
        /* "histogram.pyx":252
 *         for i in prange(bin1):
 *             edges1[i] = min1 + (0.5 +< double > i) / idp1
 *         for i in range(size):             # <<<<<<<<<<<<<<
 *             p0 = cpos0[i]
 *             p1 = cpos1[i]
 */
        __pyx_t_9 = __pyx_v_size;
        for (__pyx_t_5 = 0; __pyx_t_5 < __pyx_t_9; __pyx_t_5+=1) {
          __pyx_v_i = __pyx_t_5;
253: p0 = cpos0[i]
          /* "histogram.pyx":253
 *             edges1[i] = min1 + (0.5 +< double > i) / idp1
 *         for i in range(size):
 *             p0 = cpos0[i]             # <<<<<<<<<<<<<<
 *             p1 = cpos1[i]
 *             d = data[i]
 */
          __pyx_t_23 = __pyx_v_i;
          __pyx_v_p0 = (*__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_cpos0.buf, __pyx_t_23, __pyx_bstride_0_cpos0));
254: p1 = cpos1[i]
          /* "histogram.pyx":254
 *         for i in range(size):
 *             p0 = cpos0[i]
 *             p1 = cpos1[i]             # <<<<<<<<<<<<<<
 *             d = data[i]
 *             fbin0 = (p0 - min0) * idp0
 */
          __pyx_t_26 = __pyx_v_i;
          __pyx_v_p1 = (*__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_cpos1.buf, __pyx_t_26, __pyx_bstride_0_cpos1));
255: d = data[i]
          /* "histogram.pyx":255
 *             p0 = cpos0[i]
 *             p1 = cpos1[i]
 *             d = data[i]             # <<<<<<<<<<<<<<
 *             fbin0 = (p0 - min0) * idp0
 *             fbin1 = (p1 - min1) * idp1
 */
          __pyx_t_27 = __pyx_v_i;
          __pyx_v_d = (*__Pyx_BufPtrStrided1d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_data.buf, __pyx_t_27, __pyx_bstride_0_data));
256: fbin0 = (p0 - min0) * idp0
          /* "histogram.pyx":256
 *             p1 = cpos1[i]
 *             d = data[i]
 *             fbin0 = (p0 - min0) * idp0             # <<<<<<<<<<<<<<
 *             fbin1 = (p1 - min1) * idp1
 *             b0 = < long > floor(fbin0)
 */
          __pyx_v_fbin0 = ((__pyx_v_p0 - __pyx_v_min0) * __pyx_v_idp0);
257: fbin1 = (p1 - min1) * idp1
          /* "histogram.pyx":257
 *             d = data[i]
 *             fbin0 = (p0 - min0) * idp0
 *             fbin1 = (p1 - min1) * idp1             # <<<<<<<<<<<<<<
 *             b0 = < long > floor(fbin0)
 *             b1 = < long > floor(fbin1)
 */
          __pyx_v_fbin1 = ((__pyx_v_p1 - __pyx_v_min1) * __pyx_v_idp1);
258: b0 = < long > floor(fbin0)
          /* "histogram.pyx":258
 *             fbin0 = (p0 - min0) * idp0
 *             fbin1 = (p1 - min1) * idp1
 *             b0 = < long > floor(fbin0)             # <<<<<<<<<<<<<<
 *             b1 = < long > floor(fbin1)
 *             if b0 == bin0:
 */
          __pyx_v_b0 = ((long)floor(__pyx_v_fbin0));
259: b1 = < long > floor(fbin1)
          /* "histogram.pyx":259
 *             fbin1 = (p1 - min1) * idp1
 *             b0 = < long > floor(fbin0)
 *             b1 = < long > floor(fbin1)             # <<<<<<<<<<<<<<
 *             if b0 == bin0:
 *                 b0 = bin0 - 1
 */
          __pyx_v_b1 = ((long)floor(__pyx_v_fbin1));
260: if b0 == bin0:
          /* "histogram.pyx":260
 *             b0 = < long > floor(fbin0)
 *             b1 = < long > floor(fbin1)
 *             if b0 == bin0:             # <<<<<<<<<<<<<<
 *                 b0 = bin0 - 1
 *                 fbin0 = (< double > bin0) - 0.5
 */
          __pyx_t_22 = (__pyx_v_b0 == __pyx_v_bin0);
          if (__pyx_t_22) {
261: b0 = bin0 - 1
            /* "histogram.pyx":261
 *             b1 = < long > floor(fbin1)
 *             if b0 == bin0:
 *                 b0 = bin0 - 1             # <<<<<<<<<<<<<<
 *                 fbin0 = (< double > bin0) - 0.5
 *             elif b0 == 0:
 */
            __pyx_v_b0 = (__pyx_v_bin0 - 1);
262: fbin0 = (< double > bin0) - 0.5
            /* "histogram.pyx":262
 *             if b0 == bin0:
 *                 b0 = bin0 - 1
 *                 fbin0 = (< double > bin0) - 0.5             # <<<<<<<<<<<<<<
 *             elif b0 == 0:
 *                 fbin0 = 0.5
 */
            __pyx_v_fbin0 = (((double)__pyx_v_bin0) - 0.5);
            goto __pyx_L37;
          }
263: elif b0 == 0:
          /* "histogram.pyx":263
 *                 b0 = bin0 - 1
 *                 fbin0 = (< double > bin0) - 0.5
 *             elif b0 == 0:             # <<<<<<<<<<<<<<
 *                 fbin0 = 0.5
 *             if b1 == bin1:
 */
          __pyx_t_22 = (__pyx_v_b0 == 0);
          if (__pyx_t_22) {
264: fbin0 = 0.5
            /* "histogram.pyx":264
 *                 fbin0 = (< double > bin0) - 0.5
 *             elif b0 == 0:
 *                 fbin0 = 0.5             # <<<<<<<<<<<<<<
 *             if b1 == bin1:
 *                 b1 = bin1 - 1
 */
            __pyx_v_fbin0 = 0.5;
            goto __pyx_L37;
          }
          __pyx_L37:;
265: if b1 == bin1:
          /* "histogram.pyx":265
 *             elif b0 == 0:
 *                 fbin0 = 0.5
 *             if b1 == bin1:             # <<<<<<<<<<<<<<
 *                 b1 = bin1 - 1
 *                 fbin1 = (< double > bin1) - 0.5
 */
          __pyx_t_22 = (__pyx_v_b1 == __pyx_v_bin1);
          if (__pyx_t_22) {
266: b1 = bin1 - 1
            /* "histogram.pyx":266
 *                 fbin0 = 0.5
 *             if b1 == bin1:
 *                 b1 = bin1 - 1             # <<<<<<<<<<<<<<
 *                 fbin1 = (< double > bin1) - 0.5
 *             elif b1 == 0:
 */
            __pyx_v_b1 = (__pyx_v_bin1 - 1);
267: fbin1 = (< double > bin1) - 0.5
            /* "histogram.pyx":267
 *             if b1 == bin1:
 *                 b1 = bin1 - 1
 *                 fbin1 = (< double > bin1) - 0.5             # <<<<<<<<<<<<<<
 *             elif b1 == 0:
 *                 fbin1 = 0.5
 */
            __pyx_v_fbin1 = (((double)__pyx_v_bin1) - 0.5);
            goto __pyx_L38;
          }
268: elif b1 == 0:
          /* "histogram.pyx":268
 *                 b1 = bin1 - 1
 *                 fbin1 = (< double > bin1) - 0.5
 *             elif b1 == 0:             # <<<<<<<<<<<<<<
 *                 fbin1 = 0.5
 * 
 */
          __pyx_t_22 = (__pyx_v_b1 == 0);
          if (__pyx_t_22) {
269: fbin1 = 0.5
            /* "histogram.pyx":269
 *                 fbin1 = (< double > bin1) - 0.5
 *             elif b1 == 0:
 *                 fbin1 = 0.5             # <<<<<<<<<<<<<<
 * 
 *             delta0l = fbin0 -< double > b0 - dbin0
 */
            __pyx_v_fbin1 = 0.5;
            goto __pyx_L38;
          }
          __pyx_L38:;
270:
271: delta0l = fbin0 -< double > b0 - dbin0
          /* "histogram.pyx":271
 *                 fbin1 = 0.5
 * 
 *             delta0l = fbin0 -< double > b0 - dbin0             # <<<<<<<<<<<<<<
 *             delta0r = fbin0 -< double > b0 - 1 + dbin0
 *             delta1l = fbin1 -< double > b1 - dbin1
 */
          __pyx_v_delta0l = ((__pyx_v_fbin0 - ((double)__pyx_v_b0)) - __pyx_v_dbin0);
272: delta0r = fbin0 -< double > b0 - 1 + dbin0
          /* "histogram.pyx":272
 * 
 *             delta0l = fbin0 -< double > b0 - dbin0
 *             delta0r = fbin0 -< double > b0 - 1 + dbin0             # <<<<<<<<<<<<<<
 *             delta1l = fbin1 -< double > b1 - dbin1
 *             delta1r = fbin1 -< double > b1 - 1 + dbin1
 */
          __pyx_v_delta0r = (((__pyx_v_fbin0 - ((double)__pyx_v_b0)) - 1.0) + __pyx_v_dbin0);
273: delta1l = fbin1 -< double > b1 - dbin1
          /* "histogram.pyx":273
 *             delta0l = fbin0 -< double > b0 - dbin0
 *             delta0r = fbin0 -< double > b0 - 1 + dbin0
 *             delta1l = fbin1 -< double > b1 - dbin1             # <<<<<<<<<<<<<<
 *             delta1r = fbin1 -< double > b1 - 1 + dbin1
 *             rest = 1.0
 */
          __pyx_v_delta1l = ((__pyx_v_fbin1 - ((double)__pyx_v_b1)) - __pyx_v_dbin1);
274: delta1r = fbin1 -< double > b1 - 1 + dbin1
          /* "histogram.pyx":274
 *             delta0r = fbin0 -< double > b0 - 1 + dbin0
 *             delta1l = fbin1 -< double > b1 - dbin1
 *             delta1r = fbin1 -< double > b1 - 1 + dbin1             # <<<<<<<<<<<<<<
 *             rest = 1.0
 *             if csplit == 1:
 */
          __pyx_v_delta1r = (((__pyx_v_fbin1 - ((double)__pyx_v_b1)) - 1.0) + __pyx_v_dbin1);
275: rest = 1.0
          /* "histogram.pyx":275
 *             delta1l = fbin1 -< double > b1 - dbin1
 *             delta1r = fbin1 -< double > b1 - 1 + dbin1
 *             rest = 1.0             # <<<<<<<<<<<<<<
 *             if csplit == 1:
 *                 if delta0l < 0 and b0 > 0:
 */
          __pyx_v_rest = 1.0;
276: if csplit == 1:
          /* "histogram.pyx":276
 *             delta1r = fbin1 -< double > b1 - 1 + dbin1
 *             rest = 1.0
 *             if csplit == 1:             # <<<<<<<<<<<<<<
 *                 if delta0l < 0 and b0 > 0:
 *                     if delta1l < 0 and b1 > 0:
 */
          __pyx_t_22 = (__pyx_v_csplit == 1);
          if (__pyx_t_22) {
277: if delta0l < 0 and b0 > 0:
            /* "histogram.pyx":277
 *             rest = 1.0
 *             if csplit == 1:
 *                 if delta0l < 0 and b0 > 0:             # <<<<<<<<<<<<<<
 *                     if delta1l < 0 and b1 > 0:
 *                         area = delta0l * delta1l
 */
            __pyx_t_22 = (__pyx_v_delta0l < 0.0);
            if (__pyx_t_22) {
              __pyx_t_4 = (__pyx_v_b0 > 0);
              __pyx_t_21 = __pyx_t_4;
            } else {
              __pyx_t_21 = __pyx_t_22;
            }
            if (__pyx_t_21) {
278: if delta1l < 0 and b1 > 0:
              /* "histogram.pyx":278
 *             if csplit == 1:
 *                 if delta0l < 0 and b0 > 0:
 *                     if delta1l < 0 and b1 > 0:             # <<<<<<<<<<<<<<
 *                         area = delta0l * delta1l
 *                         rest -= area
 */
              __pyx_t_21 = (__pyx_v_delta1l < 0.0);
              if (__pyx_t_21) {
                __pyx_t_22 = (__pyx_v_b1 > 0);
                __pyx_t_4 = __pyx_t_22;
              } else {
                __pyx_t_4 = __pyx_t_21;
              }
              if (__pyx_t_4) {
279: area = delta0l * delta1l
                /* "histogram.pyx":279
 *                 if delta0l < 0 and b0 > 0:
 *                     if delta1l < 0 and b1 > 0:
 *                         area = delta0l * delta1l             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 - 1, b1 - 1] += area
 */
                __pyx_v_area = (__pyx_v_delta0l * __pyx_v_delta1l);
280: rest -= area
                /* "histogram.pyx":280
 *                     if delta1l < 0 and b1 > 0:
 *                         area = delta0l * delta1l
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 - 1, b1 - 1] += area
 *                         outData[b0 - 1, b1 - 1] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
281: outCount[b0 - 1, b1 - 1] += area
                /* "histogram.pyx":281
 *                         area = delta0l * delta1l
 *                         rest -= area
 *                         outCount[b0 - 1, b1 - 1] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 - 1, b1 - 1] += area * d
 * 
 */
                __pyx_t_28 = (__pyx_v_b0 - 1);
                __pyx_t_29 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_28, __pyx_bstride_0_outCount, __pyx_t_29, __pyx_bstride_1_outCount) += __pyx_v_area;
282: outData[b0 - 1, b1 - 1] += area * d
                /* "histogram.pyx":282
 *                         rest -= area
 *                         outCount[b0 - 1, b1 - 1] += area
 *                         outData[b0 - 1, b1 - 1] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (-delta0l) * (1 + delta1l)
 */
                __pyx_t_30 = (__pyx_v_b0 - 1);
                __pyx_t_31 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_30, __pyx_bstride_0_outData, __pyx_t_31, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
283:
284: area = (-delta0l) * (1 + delta1l)
                /* "histogram.pyx":284
 *                         outData[b0 - 1, b1 - 1] += area * d
 * 
 *                         area = (-delta0l) * (1 + delta1l)             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 - 1, b1 ] += area
 */
                __pyx_v_area = ((-__pyx_v_delta0l) * (1.0 + __pyx_v_delta1l));
285: rest -= area
                /* "histogram.pyx":285
 * 
 *                         area = (-delta0l) * (1 + delta1l)
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 - 1, b1 ] += area
 *                         outData[b0 - 1, b1 ] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
286: outCount[b0 - 1, b1 ] += area
                /* "histogram.pyx":286
 *                         area = (-delta0l) * (1 + delta1l)
 *                         rest -= area
 *                         outCount[b0 - 1, b1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 - 1, b1 ] += area * d
 * 
 */
                __pyx_t_32 = (__pyx_v_b0 - 1);
                __pyx_t_33 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_32, __pyx_bstride_0_outCount, __pyx_t_33, __pyx_bstride_1_outCount) += __pyx_v_area;
287: outData[b0 - 1, b1 ] += area * d
                /* "histogram.pyx":287
 *                         rest -= area
 *                         outCount[b0 - 1, b1 ] += area
 *                         outData[b0 - 1, b1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (1 + delta0l) * (-delta1l)
 */
                __pyx_t_34 = (__pyx_v_b0 - 1);
                __pyx_t_35 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_34, __pyx_bstride_0_outData, __pyx_t_35, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
288:
289: area = (1 + delta0l) * (-delta1l)
                /* "histogram.pyx":289
 *                         outData[b0 - 1, b1 ] += area * d
 * 
 *                         area = (1 + delta0l) * (-delta1l)             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 , b1 - 1 ] += area
 */
                __pyx_v_area = ((1.0 + __pyx_v_delta0l) * (-__pyx_v_delta1l));
290: rest -= area
                /* "histogram.pyx":290
 * 
 *                         area = (1 + delta0l) * (-delta1l)
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 , b1 - 1 ] += area
 *                         outData[b0 , b1 - 1 ] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
291: outCount[b0 , b1 - 1 ] += area
                /* "histogram.pyx":291
 *                         area = (1 + delta0l) * (-delta1l)
 *                         rest -= area
 *                         outCount[b0 , b1 - 1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 , b1 - 1 ] += area * d
 * 
 */
                __pyx_t_36 = __pyx_v_b0;
                __pyx_t_37 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_36, __pyx_bstride_0_outCount, __pyx_t_37, __pyx_bstride_1_outCount) += __pyx_v_area;
292: outData[b0 , b1 - 1 ] += area * d
                /* "histogram.pyx":292
 *                         rest -= area
 *                         outCount[b0 , b1 - 1 ] += area
 *                         outData[b0 , b1 - 1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                     elif delta1r > 0 and b1 < bin1 - 1:
 */
                __pyx_t_38 = __pyx_v_b0;
                __pyx_t_39 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_38, __pyx_bstride_0_outData, __pyx_t_39, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
                goto __pyx_L41;
              }
293:
294: elif delta1r > 0 and b1 < bin1 - 1:
              /* "histogram.pyx":294
 *                         outData[b0 , b1 - 1 ] += area * d
 * 
 *                     elif delta1r > 0 and b1 < bin1 - 1:             # <<<<<<<<<<<<<<
 *                         area = -delta0l * delta1r
 *                         rest -= area
 */
              __pyx_t_4 = (__pyx_v_delta1r > 0.0);
              if (__pyx_t_4) {
                __pyx_t_21 = (__pyx_v_b1 < (__pyx_v_bin1 - 1));
                __pyx_t_22 = __pyx_t_21;
              } else {
                __pyx_t_22 = __pyx_t_4;
              }
              if (__pyx_t_22) {
295: area = -delta0l * delta1r
                /* "histogram.pyx":295
 * 
 *                     elif delta1r > 0 and b1 < bin1 - 1:
 *                         area = -delta0l * delta1r             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 - 1, b1 + 1] += area
 */
                __pyx_v_area = ((-__pyx_v_delta0l) * __pyx_v_delta1r);
296: rest -= area
                /* "histogram.pyx":296
 *                     elif delta1r > 0 and b1 < bin1 - 1:
 *                         area = -delta0l * delta1r
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 - 1, b1 + 1] += area
 *                         outData[b0 - 1, b1 + 1] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
297: outCount[b0 - 1, b1 + 1] += area
                /* "histogram.pyx":297
 *                         area = -delta0l * delta1r
 *                         rest -= area
 *                         outCount[b0 - 1, b1 + 1] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 - 1, b1 + 1] += area * d
 * 
 */
                __pyx_t_40 = (__pyx_v_b0 - 1);
                __pyx_t_41 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_40, __pyx_bstride_0_outCount, __pyx_t_41, __pyx_bstride_1_outCount) += __pyx_v_area;
298: outData[b0 - 1, b1 + 1] += area * d
                /* "histogram.pyx":298
 *                         rest -= area
 *                         outCount[b0 - 1, b1 + 1] += area
 *                         outData[b0 - 1, b1 + 1] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (-delta0l) * (1 - delta1r)
 */
                __pyx_t_42 = (__pyx_v_b0 - 1);
                __pyx_t_43 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_42, __pyx_bstride_0_outData, __pyx_t_43, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
299:
300: area = (-delta0l) * (1 - delta1r)
                /* "histogram.pyx":300
 *                         outData[b0 - 1, b1 + 1] += area * d
 * 
 *                         area = (-delta0l) * (1 - delta1r)             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 - 1, b1 ] += area
 */
                __pyx_v_area = ((-__pyx_v_delta0l) * (1.0 - __pyx_v_delta1r));
301: rest -= area
                /* "histogram.pyx":301
 * 
 *                         area = (-delta0l) * (1 - delta1r)
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 - 1, b1 ] += area
 *                         outData[b0 - 1, b1 ] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
302: outCount[b0 - 1, b1 ] += area
                /* "histogram.pyx":302
 *                         area = (-delta0l) * (1 - delta1r)
 *                         rest -= area
 *                         outCount[b0 - 1, b1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 - 1, b1 ] += area * d
 * 
 */
                __pyx_t_44 = (__pyx_v_b0 - 1);
                __pyx_t_45 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_44, __pyx_bstride_0_outCount, __pyx_t_45, __pyx_bstride_1_outCount) += __pyx_v_area;
303: outData[b0 - 1, b1 ] += area * d
                /* "histogram.pyx":303
 *                         rest -= area
 *                         outCount[b0 - 1, b1 ] += area
 *                         outData[b0 - 1, b1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (1 + delta0l) * (delta1r)
 */
                __pyx_t_46 = (__pyx_v_b0 - 1);
                __pyx_t_47 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_46, __pyx_bstride_0_outData, __pyx_t_47, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
304:
305: area = (1 + delta0l) * (delta1r)
                /* "histogram.pyx":305
 *                         outData[b0 - 1, b1 ] += area * d
 * 
 *                         area = (1 + delta0l) * (delta1r)             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 , b1 + 1 ] += area
 */
                __pyx_v_area = ((1.0 + __pyx_v_delta0l) * __pyx_v_delta1r);
306: rest -= area
                /* "histogram.pyx":306
 * 
 *                         area = (1 + delta0l) * (delta1r)
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 , b1 + 1 ] += area
 *                         outData[b0 , b1 + 1 ] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
307: outCount[b0 , b1 + 1 ] += area
                /* "histogram.pyx":307
 *                         area = (1 + delta0l) * (delta1r)
 *                         rest -= area
 *                         outCount[b0 , b1 + 1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 , b1 + 1 ] += area * d
 *                 elif delta0r > 0 and b0 < bin0 - 1:
 */
                __pyx_t_48 = __pyx_v_b0;
                __pyx_t_49 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_48, __pyx_bstride_0_outCount, __pyx_t_49, __pyx_bstride_1_outCount) += __pyx_v_area;
308: outData[b0 , b1 + 1 ] += area * d
                /* "histogram.pyx":308
 *                         rest -= area
 *                         outCount[b0 , b1 + 1 ] += area
 *                         outData[b0 , b1 + 1 ] += area * d             # <<<<<<<<<<<<<<
 *                 elif delta0r > 0 and b0 < bin0 - 1:
 *                     if delta1l < 0 and b1 > 0:
 */
                __pyx_t_50 = __pyx_v_b0;
                __pyx_t_51 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_50, __pyx_bstride_0_outData, __pyx_t_51, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
                goto __pyx_L41;
              }
              __pyx_L41:;
              goto __pyx_L40;
            }
309: elif delta0r > 0 and b0 < bin0 - 1:
            /* "histogram.pyx":309
 *                         outCount[b0 , b1 + 1 ] += area
 *                         outData[b0 , b1 + 1 ] += area * d
 *                 elif delta0r > 0 and b0 < bin0 - 1:             # <<<<<<<<<<<<<<
 *                     if delta1l < 0 and b1 > 0:
 *                         area = -delta0r * delta1l
 */
            __pyx_t_22 = (__pyx_v_delta0r > 0.0);
            if (__pyx_t_22) {
              __pyx_t_4 = (__pyx_v_b0 < (__pyx_v_bin0 - 1));
              __pyx_t_21 = __pyx_t_4;
            } else {
              __pyx_t_21 = __pyx_t_22;
            }
            if (__pyx_t_21) {
310: if delta1l < 0 and b1 > 0:
              /* "histogram.pyx":310
 *                         outData[b0 , b1 + 1 ] += area * d
 *                 elif delta0r > 0 and b0 < bin0 - 1:
 *                     if delta1l < 0 and b1 > 0:             # <<<<<<<<<<<<<<
 *                         area = -delta0r * delta1l
 *                         rest -= area
 */
              __pyx_t_21 = (__pyx_v_delta1l < 0.0);
              if (__pyx_t_21) {
                __pyx_t_22 = (__pyx_v_b1 > 0);
                __pyx_t_4 = __pyx_t_22;
              } else {
                __pyx_t_4 = __pyx_t_21;
              }
              if (__pyx_t_4) {
311: area = -delta0r * delta1l
                /* "histogram.pyx":311
 *                 elif delta0r > 0 and b0 < bin0 - 1:
 *                     if delta1l < 0 and b1 > 0:
 *                         area = -delta0r * delta1l             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 + 1, b1 - 1] += area
 */
                __pyx_v_area = ((-__pyx_v_delta0r) * __pyx_v_delta1l);
312: rest -= area
                /* "histogram.pyx":312
 *                     if delta1l < 0 and b1 > 0:
 *                         area = -delta0r * delta1l
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 + 1, b1 - 1] += area
 *                         outData[b0 + 1, b1 - 1] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
313: outCount[b0 + 1, b1 - 1] += area
                /* "histogram.pyx":313
 *                         area = -delta0r * delta1l
 *                         rest -= area
 *                         outCount[b0 + 1, b1 - 1] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 + 1, b1 - 1] += area * d
 * 
 */
                __pyx_t_52 = (__pyx_v_b0 + 1);
                __pyx_t_53 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_52, __pyx_bstride_0_outCount, __pyx_t_53, __pyx_bstride_1_outCount) += __pyx_v_area;
314: outData[b0 + 1, b1 - 1] += area * d
                /* "histogram.pyx":314
 *                         rest -= area
 *                         outCount[b0 + 1, b1 - 1] += area
 *                         outData[b0 + 1, b1 - 1] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (delta0r) * (1 + delta1l)
 */
                __pyx_t_54 = (__pyx_v_b0 + 1);
                __pyx_t_55 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_54, __pyx_bstride_0_outData, __pyx_t_55, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
315:
316: area = (delta0r) * (1 + delta1l)
                /* "histogram.pyx":316
 *                         outData[b0 + 1, b1 - 1] += area * d
 * 
 *                         area = (delta0r) * (1 + delta1l)             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 + 1, b1 ] += area
 */
                __pyx_v_area = (__pyx_v_delta0r * (1.0 + __pyx_v_delta1l));
317: rest -= area
                /* "histogram.pyx":317
 * 
 *                         area = (delta0r) * (1 + delta1l)
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 + 1, b1 ] += area
 *                         outData[b0 + 1, b1 ] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
318: outCount[b0 + 1, b1 ] += area
                /* "histogram.pyx":318
 *                         area = (delta0r) * (1 + delta1l)
 *                         rest -= area
 *                         outCount[b0 + 1, b1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 + 1, b1 ] += area * d
 * 
 */
                __pyx_t_56 = (__pyx_v_b0 + 1);
                __pyx_t_57 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_56, __pyx_bstride_0_outCount, __pyx_t_57, __pyx_bstride_1_outCount) += __pyx_v_area;
319: outData[b0 + 1, b1 ] += area * d
                /* "histogram.pyx":319
 *                         rest -= area
 *                         outCount[b0 + 1, b1 ] += area
 *                         outData[b0 + 1, b1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (1 - delta0r) * (-delta1l)
 */
                __pyx_t_58 = (__pyx_v_b0 + 1);
                __pyx_t_59 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_58, __pyx_bstride_0_outData, __pyx_t_59, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
320:
321: area = (1 - delta0r) * (-delta1l)
                /* "histogram.pyx":321
 *                         outData[b0 + 1, b1 ] += area * d
 * 
 *                         area = (1 - delta0r) * (-delta1l)             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 , b1 - 1 ] += area
 */
                __pyx_v_area = ((1.0 - __pyx_v_delta0r) * (-__pyx_v_delta1l));
322: rest -= area
                /* "histogram.pyx":322
 * 
 *                         area = (1 - delta0r) * (-delta1l)
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 , b1 - 1 ] += area
 *                         outData[b0 , b1 - 1 ] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
323: outCount[b0 , b1 - 1 ] += area
                /* "histogram.pyx":323
 *                         area = (1 - delta0r) * (-delta1l)
 *                         rest -= area
 *                         outCount[b0 , b1 - 1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 , b1 - 1 ] += area * d
 * 
 */
                __pyx_t_60 = __pyx_v_b0;
                __pyx_t_61 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_60, __pyx_bstride_0_outCount, __pyx_t_61, __pyx_bstride_1_outCount) += __pyx_v_area;
324: outData[b0 , b1 - 1 ] += area * d
                /* "histogram.pyx":324
 *                         rest -= area
 *                         outCount[b0 , b1 - 1 ] += area
 *                         outData[b0 , b1 - 1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                     elif delta1r > 0 and b1 < bin1 - 1:
 */
                __pyx_t_62 = __pyx_v_b0;
                __pyx_t_63 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_62, __pyx_bstride_0_outData, __pyx_t_63, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
                goto __pyx_L42;
              }
325:
326: elif delta1r > 0 and b1 < bin1 - 1:
              /* "histogram.pyx":326
 *                         outData[b0 , b1 - 1 ] += area * d
 * 
 *                     elif delta1r > 0 and b1 < bin1 - 1:             # <<<<<<<<<<<<<<
 *                         area = delta0r * delta1r
 *                         rest -= area
 */
              __pyx_t_4 = (__pyx_v_delta1r > 0.0);
              if (__pyx_t_4) {
                __pyx_t_21 = (__pyx_v_b1 < (__pyx_v_bin1 - 1));
                __pyx_t_22 = __pyx_t_21;
              } else {
                __pyx_t_22 = __pyx_t_4;
              }
              if (__pyx_t_22) {
327: area = delta0r * delta1r
                /* "histogram.pyx":327
 * 
 *                     elif delta1r > 0 and b1 < bin1 - 1:
 *                         area = delta0r * delta1r             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 + 1, b1 + 1] += area
 */
                __pyx_v_area = (__pyx_v_delta0r * __pyx_v_delta1r);
328: rest -= area
                /* "histogram.pyx":328
 *                     elif delta1r > 0 and b1 < bin1 - 1:
 *                         area = delta0r * delta1r
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 + 1, b1 + 1] += area
 *                         outData[b0 + 1, b1 + 1] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
329: outCount[b0 + 1, b1 + 1] += area
                /* "histogram.pyx":329
 *                         area = delta0r * delta1r
 *                         rest -= area
 *                         outCount[b0 + 1, b1 + 1] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 + 1, b1 + 1] += area * d
 * 
 */
                __pyx_t_64 = (__pyx_v_b0 + 1);
                __pyx_t_65 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_64, __pyx_bstride_0_outCount, __pyx_t_65, __pyx_bstride_1_outCount) += __pyx_v_area;
330: outData[b0 + 1, b1 + 1] += area * d
                /* "histogram.pyx":330
 *                         rest -= area
 *                         outCount[b0 + 1, b1 + 1] += area
 *                         outData[b0 + 1, b1 + 1] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (delta0r) * (1 - delta1r)
 */
                __pyx_t_66 = (__pyx_v_b0 + 1);
                __pyx_t_67 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_66, __pyx_bstride_0_outData, __pyx_t_67, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
331:
332: area = (delta0r) * (1 - delta1r)
                /* "histogram.pyx":332
 *                         outData[b0 + 1, b1 + 1] += area * d
 * 
 *                         area = (delta0r) * (1 - delta1r)             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 + 1, b1 ] += area
 */
                __pyx_v_area = (__pyx_v_delta0r * (1.0 - __pyx_v_delta1r));
333: rest -= area
                /* "histogram.pyx":333
 * 
 *                         area = (delta0r) * (1 - delta1r)
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 + 1, b1 ] += area
 *                         outData[b0 + 1, b1 ] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
334: outCount[b0 + 1, b1 ] += area
                /* "histogram.pyx":334
 *                         area = (delta0r) * (1 - delta1r)
 *                         rest -= area
 *                         outCount[b0 + 1, b1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 + 1, b1 ] += area * d
 * 
 */
                __pyx_t_68 = (__pyx_v_b0 + 1);
                __pyx_t_69 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_68, __pyx_bstride_0_outCount, __pyx_t_69, __pyx_bstride_1_outCount) += __pyx_v_area;
335: outData[b0 + 1, b1 ] += area * d
                /* "histogram.pyx":335
 *                         rest -= area
 *                         outCount[b0 + 1, b1 ] += area
 *                         outData[b0 + 1, b1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (1 - delta0r) * (delta1r)
 */
                __pyx_t_70 = (__pyx_v_b0 + 1);
                __pyx_t_71 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_70, __pyx_bstride_0_outData, __pyx_t_71, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
336:
337: area = (1 - delta0r) * (delta1r)
                /* "histogram.pyx":337
 *                         outData[b0 + 1, b1 ] += area * d
 * 
 *                         area = (1 - delta0r) * (delta1r)             # <<<<<<<<<<<<<<
 *                         rest -= area
 *                         outCount[b0 , b1 + 1 ] += area
 */
                __pyx_v_area = ((1.0 - __pyx_v_delta0r) * __pyx_v_delta1r);
338: rest -= area
                /* "histogram.pyx":338
 * 
 *                         area = (1 - delta0r) * (delta1r)
 *                         rest -= area             # <<<<<<<<<<<<<<
 *                         outCount[b0 , b1 + 1 ] += area
 *                         outData[b0 , b1 + 1 ] += area * d
 */
                __pyx_v_rest = (__pyx_v_rest - __pyx_v_area);
339: outCount[b0 , b1 + 1 ] += area
                /* "histogram.pyx":339
 *                         area = (1 - delta0r) * (delta1r)
 *                         rest -= area
 *                         outCount[b0 , b1 + 1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 , b1 + 1 ] += area * d
 *             outCount[b0, b1] += rest
 */
                __pyx_t_72 = __pyx_v_b0;
                __pyx_t_73 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_72, __pyx_bstride_0_outCount, __pyx_t_73, __pyx_bstride_1_outCount) += __pyx_v_area;
340: outData[b0 , b1 + 1 ] += area * d
                /* "histogram.pyx":340
 *                         rest -= area
 *                         outCount[b0 , b1 + 1 ] += area
 *                         outData[b0 , b1 + 1 ] += area * d             # <<<<<<<<<<<<<<
 *             outCount[b0, b1] += rest
 *             outData[b0, b1] += d * rest
 */
                __pyx_t_74 = __pyx_v_b0;
                __pyx_t_75 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_74, __pyx_bstride_0_outData, __pyx_t_75, __pyx_bstride_1_outData) += (__pyx_v_area * __pyx_v_d);
                goto __pyx_L42;
              }
              __pyx_L42:;
              goto __pyx_L40;
            }
            __pyx_L40:;
            goto __pyx_L39;
          }
          __pyx_L39:;
341: outCount[b0, b1] += rest
          /* "histogram.pyx":341
 *                         outCount[b0 , b1 + 1 ] += area
 *                         outData[b0 , b1 + 1 ] += area * d
 *             outCount[b0, b1] += rest             # <<<<<<<<<<<<<<
 *             outData[b0, b1] += d * rest
 * 
 */
          __pyx_t_76 = __pyx_v_b0;
          __pyx_t_77 = __pyx_v_b1;
          *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_76, __pyx_bstride_0_outCount, __pyx_t_77, __pyx_bstride_1_outCount) += __pyx_v_rest;
342: outData[b0, b1] += d * rest
          /* "histogram.pyx":342
 *                         outData[b0 , b1 + 1 ] += area * d
 *             outCount[b0, b1] += rest
 *             outData[b0, b1] += d * rest             # <<<<<<<<<<<<<<
 * 
 *         for i in prange(bin0):
 */
          __pyx_t_78 = __pyx_v_b0;
          __pyx_t_79 = __pyx_v_b1;
          *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_78, __pyx_bstride_0_outData, __pyx_t_79, __pyx_bstride_1_outData) += (__pyx_v_d * __pyx_v_rest);
        }
343:
344: for i in prange(bin0):
        /* "histogram.pyx":344
 *             outData[b0, b1] += d * rest
 * 
 *         for i in prange(bin0):             # <<<<<<<<<<<<<<
 *             for j in range(bin1):
 *                 if outCount[i, j] > 1e-10:
 */
        __pyx_t_9 = __pyx_v_bin0;
        if (1 == 0) abort();
        {
            __pyx_t_80 = (__pyx_t_9 - 0) / 1;
            if (__pyx_t_80 > 0)
            {
                __pyx_v_i = 0;
                #ifdef _OPENMP
                #pragma omp parallel
                #endif /* _OPENMP */
                {
                    #ifdef _OPENMP
                    #pragma omp for lastprivate(__pyx_v_j) firstprivate(__pyx_v_i) lastprivate(__pyx_v_i)
                    #endif /* _OPENMP */
                    for (__pyx_t_5 = 0; __pyx_t_5 < __pyx_t_80; __pyx_t_5++){
                        {
                            __pyx_v_i = 0 + 1 * __pyx_t_5;
                            /* Initialize private variables to invalid values */
                            __pyx_v_j = ((long)0xbad0bad0);
345: for j in range(bin1):
                            /* "histogram.pyx":345
 * 
 *         for i in prange(bin0):
 *             for j in range(bin1):             # <<<<<<<<<<<<<<
 *                 if outCount[i, j] > 1e-10:
 *                     outMerge[i, j] += outData[i, j] / outCount[i, j]
 */
                            __pyx_t_81 = __pyx_v_bin1;
                            for (__pyx_t_82 = 0; __pyx_t_82 < __pyx_t_81; __pyx_t_82+=1) {
                              __pyx_v_j = __pyx_t_82;
346: if outCount[i, j] > 1e-10:
                              /* "histogram.pyx":346
 *         for i in prange(bin0):
 *             for j in range(bin1):
 *                 if outCount[i, j] > 1e-10:             # <<<<<<<<<<<<<<
 *                     outMerge[i, j] += outData[i, j] / outCount[i, j]
 *                 else:
 */
                              __pyx_t_83 = __pyx_v_i;
                              __pyx_t_84 = __pyx_v_j;
                              __pyx_t_22 = ((*__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_83, __pyx_bstride_0_outCount, __pyx_t_84, __pyx_bstride_1_outCount)) > 1e-10);
                              if (__pyx_t_22) {
347: outMerge[i, j] += outData[i, j] / outCount[i, j]
                                /* "histogram.pyx":347
 *             for j in range(bin1):
 *                 if outCount[i, j] > 1e-10:
 *                     outMerge[i, j] += outData[i, j] / outCount[i, j]             # <<<<<<<<<<<<<<
 *                 else:
 *                     outMerge[i, j] += dummy
 */
                                __pyx_t_85 = __pyx_v_i;
                                __pyx_t_86 = __pyx_v_j;
                                __pyx_t_87 = __pyx_v_i;
                                __pyx_t_88 = __pyx_v_j;
                                __pyx_t_89 = __pyx_v_i;
                                __pyx_t_90 = __pyx_v_j;
                                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outMerge.buf, __pyx_t_89, __pyx_bstride_0_outMerge, __pyx_t_90, __pyx_bstride_1_outMerge) += ((*__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outData.buf, __pyx_t_85, __pyx_bstride_0_outData, __pyx_t_86, __pyx_bstride_1_outData)) / (*__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outCount.buf, __pyx_t_87, __pyx_bstride_0_outCount, __pyx_t_88, __pyx_bstride_1_outCount)));
                                goto __pyx_L49;
                              }
                              /*else*/ {
348: else:
349: outMerge[i, j] += dummy
                                /* "histogram.pyx":349
 *                     outMerge[i, j] += outData[i, j] / outCount[i, j]
 *                 else:
 *                     outMerge[i, j] += dummy             # <<<<<<<<<<<<<<
 * 
 * 
 */
                                __pyx_t_91 = __pyx_v_i;
                                __pyx_t_92 = __pyx_v_j;
                                *__Pyx_BufPtrStrided2d(__pyx_t_9histogram_DTYPE_float64_t *, __pyx_bstruct_outMerge.buf, __pyx_t_91, __pyx_bstride_0_outMerge, __pyx_t_92, __pyx_bstride_1_outMerge) += __pyx_v_dummy;
                              }
                              __pyx_L49:;
                            }
                        }
                    }
                }
            }
        }
      }
350:
351:
352: return outMerge, edges0, edges1, outData, outCount
/* "histogram.pyx":352 * * * return outMerge, edges0, edges1, outData, outCount #<<<<<<<<<<<<<<*/ __Pyx_XDECREF(__pyx_r); __pyx_t_14 = PyTuple_New(5); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 352; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_14)); __Pyx_INCREF(((PyObject *)__pyx_v_outMerge)); PyTuple_SET_ITEM(__pyx_t_14, 0, ((PyObject *)__pyx_v_outMerge)); __Pyx_GIVEREF(((PyObject *)__pyx_v_outMerge)); __Pyx_INCREF(((PyObject *)__pyx_v_edges0)); PyTuple_SET_ITEM(__pyx_t_14, 1, ((PyObject *)__pyx_v_edges0)); __Pyx_GIVEREF(((PyObject *)__pyx_v_edges0)); __Pyx_INCREF(((PyObject *)__pyx_v_edges1)); PyTuple_SET_ITEM(__pyx_t_14, 2, ((PyObject *)__pyx_v_edges1)); __Pyx_GIVEREF(((PyObject *)__pyx_v_edges1)); __Pyx_INCREF(((PyObject *)__pyx_v_outData)); PyTuple_SET_ITEM(__pyx_t_14, 3, ((PyObject *)__pyx_v_outData)); __Pyx_GIVEREF(((PyObject *)__pyx_v_outData)); __Pyx_INCREF(((PyObject *)__pyx_v_outCount)); PyTuple_SET_ITEM(__pyx_t_14, 4, ((PyObject *)__pyx_v_outCount)); __Pyx_GIVEREF(((PyObject *)__pyx_v_outCount)); __pyx_r = ((PyObject *)__pyx_t_14); __pyx_t_14 = 0; goto __pyx_L0; __pyx_r = Py_None; __Pyx_INCREF(Py_None); goto __pyx_L0; __pyx_L1_error:; __Pyx_XDECREF(__pyx_t_1); __Pyx_XDECREF(__pyx_t_2); __Pyx_XDECREF(__pyx_t_3); __Pyx_XDECREF(__pyx_t_14); { PyObject *__pyx_type, *__pyx_value, *__pyx_tb; __Pyx_ErrFetch(&__pyx_type, &__pyx_value, &__pyx_tb); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_edges0); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_edges1); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outMerge); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outCount); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_data); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_cpos1); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_cpos0); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outData); __Pyx_ErrRestore(__pyx_type, __pyx_value, __pyx_tb);} __Pyx_AddTraceback("histogram.histogram2d", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = NULL; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_edges0); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_edges1); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outMerge); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outCount); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_data); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_cpos1); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_cpos0); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_outData); __pyx_L2:; __Pyx_XDECREF((PyObject *)__pyx_v_cpos0); __Pyx_XDECREF((PyObject *)__pyx_v_cpos1); __Pyx_XDECREF((PyObject *)__pyx_v_data); __Pyx_XDECREF((PyObject *)__pyx_v_outData); __Pyx_XDECREF((PyObject *)__pyx_v_outCount); __Pyx_XDECREF((PyObject *)__pyx_v_outMerge); __Pyx_XDECREF((PyObject *)__pyx_v_edges0); __Pyx_XDECREF((PyObject *)__pyx_v_edges1); __Pyx_XGIVEREF(__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }