Generated by Cython 0.19 on Thu Jul 11 09:13:37 2013

Raw output: histogram.c

 1: #!/usr/bin/env python
  /* "histogram.pyx":1
 * #!/usr/bin/env python             # <<<<<<<<<<<<<<
 * # -*- coding: utf-8 -*-
 * #
 */
  __pyx_t_4 = PyDict_New(); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_4));
  if (PyDict_SetItem(__pyx_d, __pyx_n_s____test__, ((PyObject *)__pyx_t_4)) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(((PyObject *)__pyx_t_4)); __pyx_t_4 = 0;
 2: # -*- coding: utf-8 -*-
 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: __author__ = "Jerome Kieffer"
  /* "histogram.pyx":27
 * #
 * 
 * __author__ = "Jerome Kieffer"             # <<<<<<<<<<<<<<
 * __date__ = "20130711"
 * 
 */
  if (PyDict_SetItem(__pyx_d, __pyx_n_s____author__, ((PyObject *)__pyx_kp_s_48)) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
 28: __date__ = "20130711"
  /* "histogram.pyx":28
 * 
 * __author__ = "Jerome Kieffer"
 * __date__ = "20130711"             # <<<<<<<<<<<<<<
 * 
 * import cython
 */
  if (PyDict_SetItem(__pyx_d, __pyx_n_s____date__, ((PyObject *)__pyx_kp_s__20130711)) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
 29: 
 30: import cython
 31: from cython.parallel cimport prange
 32: import numpy
  /* "histogram.pyx":32
 * import cython
 * from cython.parallel cimport prange
 * import numpy             # <<<<<<<<<<<<<<
 * cimport numpy
 * import sys
 */
  __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 (PyDict_SetItem(__pyx_d, __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: cimport numpy
 34: import sys
  /* "histogram.pyx":34
 * import numpy
 * cimport numpy
 * import sys             # <<<<<<<<<<<<<<
 * 
 * from libc.stdlib cimport free, calloc,malloc
 */
  __pyx_t_1 = __Pyx_Import(((PyObject *)__pyx_n_s__sys), 0, -1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  if (PyDict_SetItem(__pyx_d, __pyx_n_s__sys, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
 35: 
 36: from libc.stdlib cimport free, calloc,malloc
 37: from libc.math cimport floor,fabs
 38: from openmp cimport omp_set_num_threads, omp_get_max_threads, omp_get_thread_num
 39: EPS32 = (1.0 + numpy.finfo(numpy.float32).eps)
  /* "histogram.pyx":39
 * from libc.math cimport floor,fabs
 * from openmp cimport omp_set_num_threads, omp_get_max_threads, omp_get_thread_num
 * EPS32 = (1.0 + numpy.finfo(numpy.float32).eps)             # <<<<<<<<<<<<<<
 * 
 * @cython.cdivision(True)
 */
  __pyx_t_1 = PyFloat_FromDouble(1.0); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_2 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __pyx_t_3 = __Pyx_PyObject_GetAttrStr(__pyx_t_2, __pyx_n_s__finfo); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __pyx_t_2 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __pyx_t_4 = __Pyx_PyObject_GetAttrStr(__pyx_t_2, __pyx_n_s__float32); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_4);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __pyx_t_2 = PyTuple_New(1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  PyTuple_SET_ITEM(__pyx_t_2, 0, __pyx_t_4);
  __Pyx_GIVEREF(__pyx_t_4);
  __pyx_t_4 = 0;
  __pyx_t_4 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_t_2), NULL); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_4);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0;
  __pyx_t_2 = __Pyx_PyObject_GetAttrStr(__pyx_t_4, __pyx_n_s__eps); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
  __pyx_t_4 = PyNumber_Add(__pyx_t_1, __pyx_t_2); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_4);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  if (PyDict_SetItem(__pyx_d, __pyx_n_s__EPS32, __pyx_t_4) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
 40: 
 41: @cython.cdivision(True)
 42: @cython.boundscheck(False)
 43: @cython.wraparound(False)
 44: def histogram(numpy.ndarray pos not None, \
/* "histogram.pyx":44
 * @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(CYTHON_UNUSED PyObject *__pyx_self, PyArrayObject *__pyx_v_pos, PyArrayObject *__pyx_v_weights, long __pyx_v_bins, PyObject *__pyx_v_bin_range, PyObject *__pyx_v_pixelSize_in_Pos, PyObject *__pyx_v_nthread, double __pyx_v_dummy) {
  CYTHON_UNUSED long __pyx_v_size;
  __Pyx_memviewslice __pyx_v_cpos = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_cdata = { 0, 0, { 0 }, { 0 }, { 0 } };
  PyArrayObject *__pyx_v_outData = 0;
  PyArrayObject *__pyx_v_outCount = 0;
  PyArrayObject *__pyx_v_outMerge = 0;
  PyArrayObject *__pyx_v_outPos = 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;
  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;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_outCount;
  __Pyx_Buffer __pyx_pybuffer_outCount;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_outData;
  __Pyx_Buffer __pyx_pybuffer_outData;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_outMerge;
  __Pyx_Buffer __pyx_pybuffer_outMerge;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_outPos;
  __Pyx_Buffer __pyx_pybuffer_outPos;
  PyObject *__pyx_r = NULL;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("histogram", 0);
  __pyx_pybuffer_outData.pybuffer.buf = NULL;
  __pyx_pybuffer_outData.refcount = 0;
  __pyx_pybuffernd_outData.data = NULL;
  __pyx_pybuffernd_outData.rcbuffer = &__pyx_pybuffer_outData;
  __pyx_pybuffer_outCount.pybuffer.buf = NULL;
  __pyx_pybuffer_outCount.refcount = 0;
  __pyx_pybuffernd_outCount.data = NULL;
  __pyx_pybuffernd_outCount.rcbuffer = &__pyx_pybuffer_outCount;
  __pyx_pybuffer_outMerge.pybuffer.buf = NULL;
  __pyx_pybuffer_outMerge.refcount = 0;
  __pyx_pybuffernd_outMerge.data = NULL;
  __pyx_pybuffernd_outMerge.rcbuffer = &__pyx_pybuffer_outMerge;
  __pyx_pybuffer_outPos.pybuffer.buf = NULL;
  __pyx_pybuffer_outPos.refcount = 0;
  __pyx_pybuffernd_outPos.data = NULL;
  __pyx_pybuffernd_outPos.rcbuffer = &__pyx_pybuffer_outPos;

  /* "histogram.pyx":44
 * @cython.boundscheck(False)
 * @cython.wraparound(False)
 * def histogram(numpy.ndarray pos not None, \             # <<<<<<<<<<<<<<
 *               numpy.ndarray weights not None, \
 *               long bins=100,
 */
  __pyx_k_tuple_49 = PyTuple_Pack(38, ((PyObject *)__pyx_n_s__pos), ((PyObject *)__pyx_n_s__weights), ((PyObject *)__pyx_n_s__bins), ((PyObject *)__pyx_n_s__bin_range), ((PyObject *)__pyx_n_s__pixelSize_in_Pos), ((PyObject *)__pyx_n_s__nthread), ((PyObject *)__pyx_n_s__dummy), ((PyObject *)__pyx_n_s__size), ((PyObject *)__pyx_n_s__cpos), ((PyObject *)__pyx_n_s__cdata), ((PyObject *)__pyx_n_s__outData), ((PyObject *)__pyx_n_s__outCount), ((PyObject *)__pyx_n_s__outMerge), ((PyObject *)__pyx_n_s__outPos), ((PyObject *)__pyx_n_s__bin_edge_min), ((PyObject *)__pyx_n_s__bin_edge_max), ((PyObject *)__pyx_n_s__bin_width), ((PyObject *)__pyx_n_s__inv_bin_width), ((PyObject *)__pyx_n_s__a), ((PyObject *)__pyx_n_s__d), ((PyObject *)__pyx_n_s__fbin), ((PyObject *)__pyx_n_s__ffbin), ((PyObject *)__pyx_n_s__dInt), ((PyObject *)__pyx_n_s__dIntR), ((PyObject *)__pyx_n_s__dIntL), ((PyObject *)__pyx_n_s__dtmp), ((PyObject *)__pyx_n_s__dbin), ((PyObject *)__pyx_n_s__inv_dbin2), ((PyObject *)__pyx_n_s__tmp_count), ((PyObject *)__pyx_n_s__tmp_data), ((PyObject *)__pyx_n_s__epsilon), ((PyObject *)__pyx_n_s__bin), ((PyObject *)__pyx_n_s__i), ((PyObject *)__pyx_n_s__idx), ((PyObject *)__pyx_n_s__t), ((PyObject *)__pyx_n_s__dest), ((PyObject *)__pyx_n_s__bigCount), ((PyObject *)__pyx_n_s__bigData)); if (unlikely(!__pyx_k_tuple_49)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 44; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_k_tuple_49);
  __Pyx_GIVEREF(((PyObject *)__pyx_k_tuple_49));

  /* "histogram.pyx":44
 * @cython.boundscheck(False)
 * @cython.wraparound(False)
 * def histogram(numpy.ndarray pos not None, \             # <<<<<<<<<<<<<<
 *               numpy.ndarray weights not None, \
 *               long bins=100,
 */
  __pyx_t_4 = PyCFunction_NewEx(&__pyx_mdef_9histogram_1histogram, NULL, __pyx_n_s__histogram); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 44; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_4);
  if (PyDict_SetItem(__pyx_d, __pyx_n_s__histogram, __pyx_t_4) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 44; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
  __pyx_k_codeobj_50 = (PyObject*)__Pyx_PyCode_New(7, 0, 38, 0, 0, __pyx_empty_bytes, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_k_tuple_49, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_kp_s_51, __pyx_n_s__histogram, 44, __pyx_empty_bytes); if (unlikely(!__pyx_k_codeobj_50)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 44; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
 45:               numpy.ndarray weights not None, \
 46:               long bins=100,
 47:               bin_range=None,
    /* "histogram.pyx":47
 *               numpy.ndarray weights not None, \
 *               long bins=100,
 *               bin_range=None,             # <<<<<<<<<<<<<<
 *               pixelSize_in_Pos=None,
 *               nthread=None,
 */
    values[3] = ((PyObject *)Py_None);
 48:               pixelSize_in_Pos=None,
    /* "histogram.pyx":48
 *               long bins=100,
 *               bin_range=None,
 *               pixelSize_in_Pos=None,             # <<<<<<<<<<<<<<
 *               nthread=None,
 *               double dummy=0.0):
 */
    values[4] = ((PyObject *)Py_None);
 49:               nthread=None,
    /* "histogram.pyx":49
 *               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;
      const Py_ssize_t pos_args = PyTuple_GET_SIZE(__pyx_args);
      switch (pos_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 (pos_args) {
        case  0:
        if (likely((values[0] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__pos)) != 0)) kw_args--;
        else goto __pyx_L5_argtuple_error;
        case  1:
        if (likely((values[1] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__weights)) != 0)) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("histogram", 0, 2, 7, 1); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 44; __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, pos_args, "histogram") < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 44; __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 = 46; __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 = 50; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
    } else {
 50:               double dummy=0.0):
      /* "histogram.pyx":50
 *               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 = 44; __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:;
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_pos), __pyx_ptype_5numpy_ndarray, 0, "pos", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 44; __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 = 45; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_r = __pyx_pf_9histogram_histogram(__pyx_self, __pyx_v_pos, __pyx_v_weights, __pyx_v_bins, __pyx_v_bin_range, __pyx_v_pixelSize_in_Pos, __pyx_v_nthread, __pyx_v_dummy);
  goto __pyx_L0;
  __pyx_L1_error:;
  __pyx_r = NULL;
  __pyx_L0:;
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}
 51:     """
 52:     Calculates histogram of pos weighted by weights
 53: 
 54:     @param pos: 2Theta array
 55:     @param weights: array with intensities
 56:     @param bins: number of output bins
 57:     @param pixelSize_in_Pos: size of a pixels in 2theta
 58:     @param nthread: maximum number of thread to use. By default: maximum available.
 59:         One can also limit this with OMP_NUM_THREADS environment variable
 60: 
 61:     @return 2theta, I, weighted histogram, raw histogram
 62:     """
 63: 
 64:     assert pos.size == weights.size
  /* "histogram.pyx":64
 *     """
 * 
 *     assert pos.size == weights.size             # <<<<<<<<<<<<<<
 *     assert  bins > 1
 *     cdef long  size = pos.size
 */
  #ifndef CYTHON_WITHOUT_ASSERTIONS
  __pyx_t_1 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos), __pyx_n_s__size); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 64; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_2 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_weights), __pyx_n_s__size); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 64; __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); __Pyx_XGOTREF(__pyx_t_3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 64; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __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 = 64; __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 = 64; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }
  #endif
 65:     assert  bins > 1
  /* "histogram.pyx":65
 * 
 *     assert pos.size == weights.size
 *     assert  bins > 1             # <<<<<<<<<<<<<<
 *     cdef long  size = pos.size
 *     cdef double[:] cpos = numpy.ascontiguousarray(pos.ravel(),dtype=numpy.float64)
 */
  #ifndef CYTHON_WITHOUT_ASSERTIONS
  if (unlikely(!(__pyx_v_bins > 1))) {
    PyErr_SetNone(PyExc_AssertionError);
    {__pyx_filename = __pyx_f[0]; __pyx_lineno = 65; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }
  #endif
 66:     cdef long  size = pos.size
  /* "histogram.pyx":66
 *     assert pos.size == weights.size
 *     assert  bins > 1
 *     cdef long  size = pos.size             # <<<<<<<<<<<<<<
 *     cdef double[:] cpos = numpy.ascontiguousarray(pos.ravel(),dtype=numpy.float64)
 *     cdef double[:] cdata = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)
 */
  __pyx_t_3 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos), __pyx_n_s__size); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 66; __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 = 66; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_v_size = __pyx_t_5;
 67:     cdef double[:] cpos = numpy.ascontiguousarray(pos.ravel(),dtype=numpy.float64)
  /* "histogram.pyx":67
 *     assert  bins > 1
 *     cdef long  size = pos.size
 *     cdef double[:] cpos = numpy.ascontiguousarray(pos.ravel(),dtype=numpy.float64)             # <<<<<<<<<<<<<<
 *     cdef double[:] cdata = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")
 */
  __pyx_t_3 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_2 = __Pyx_PyObject_GetAttrStr(__pyx_t_3, __pyx_n_s__ascontiguousarray); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __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_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos), __pyx_n_s__ravel); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __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 = 67; __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 = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__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 = PyDict_New(); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_1));
  __pyx_t_6 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_6);
  __pyx_t_7 = __Pyx_PyObject_GetAttrStr(__pyx_t_6, __pyx_n_s__float64); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
  if (PyDict_SetItem(__pyx_t_1, ((PyObject *)__pyx_n_s__dtype), __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
  __pyx_t_7 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_t_3), ((PyObject *)__pyx_t_1)); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
  __pyx_t_8 = __Pyx_PyObject_to_MemoryviewSlice_ds_double(__pyx_t_7);
  if (unlikely(!__pyx_t_8.memview)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
  __pyx_v_cpos = __pyx_t_8;
  __pyx_t_8.memview = NULL;
  __pyx_t_8.data = NULL;
 68:     cdef double[:] cdata = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)
  /* "histogram.pyx":68
 *     cdef long  size = pos.size
 *     cdef double[:] cpos = numpy.ascontiguousarray(pos.ravel(),dtype=numpy.float64)
 *     cdef double[:] cdata = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64")
 */
  __pyx_t_7 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  __pyx_t_1 = __Pyx_PyObject_GetAttrStr(__pyx_t_7, __pyx_n_s__ascontiguousarray); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
  __pyx_t_7 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_weights), __pyx_n_s__ravel); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  __pyx_t_3 = PyObject_Call(__pyx_t_7, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
  __pyx_t_7 = PyTuple_New(1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  PyTuple_SET_ITEM(__pyx_t_7, 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 = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_3));
  __pyx_t_2 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __pyx_t_6 = __Pyx_PyObject_GetAttrStr(__pyx_t_2, __pyx_n_s__float64); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_6);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  if (PyDict_SetItem(__pyx_t_3, ((PyObject *)__pyx_n_s__dtype), __pyx_t_6) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
  __pyx_t_6 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_t_7), ((PyObject *)__pyx_t_3)); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_6);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_7)); __pyx_t_7 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0;
  __pyx_t_9 = __Pyx_PyObject_to_MemoryviewSlice_ds_double(__pyx_t_6);
  if (unlikely(!__pyx_t_9.memview)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
  __pyx_v_cdata = __pyx_t_9;
  __pyx_t_9.memview = NULL;
  __pyx_t_9.data = NULL;
 69:     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")
  /* "histogram.pyx":69
 *     cdef double[:] cpos = numpy.ascontiguousarray(pos.ravel(),dtype=numpy.float64)
 *     cdef double[:] cdata = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outMerge = numpy.zeros(bins, dtype="float64")
 */
  __pyx_t_6 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_6);
  __pyx_t_3 = __Pyx_PyObject_GetAttrStr(__pyx_t_6, __pyx_n_s__zeros); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
  __pyx_t_6 = PyInt_FromLong(__pyx_v_bins); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_6);
  __pyx_t_7 = PyTuple_New(1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  PyTuple_SET_ITEM(__pyx_t_7, 0, __pyx_t_6);
  __Pyx_GIVEREF(__pyx_t_6);
  __pyx_t_6 = 0;
  __pyx_t_6 = PyDict_New(); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_6));
  if (PyDict_SetItem(__pyx_t_6, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64)) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_1 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_t_7), ((PyObject *)__pyx_t_6)); 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);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_7)); __pyx_t_7 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_6)); __pyx_t_6 = 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 = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_10 = ((PyArrayObject *)__pyx_t_1);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_outData.rcbuffer->pybuffer, (PyObject*)__pyx_t_10, &__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) {
      __pyx_v_outData = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 69; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_pybuffernd_outData.diminfo[0].strides = __pyx_pybuffernd_outData.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_outData.diminfo[0].shape = __pyx_pybuffernd_outData.rcbuffer->pybuffer.shape[0];
    }
  }
  __pyx_t_10 = 0;
  __pyx_v_outData = ((PyArrayObject *)__pyx_t_1);
  __pyx_t_1 = 0;
 70:     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64")
  /* "histogram.pyx":70
 *     cdef double[:] cdata = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64")             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outMerge = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outPos = numpy.zeros(bins, dtype="float64")
 */
  __pyx_t_1 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_6 = __Pyx_PyObject_GetAttrStr(__pyx_t_1, __pyx_n_s__zeros); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_6);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_1 = PyInt_FromLong(__pyx_v_bins); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_7 = PyTuple_New(1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  PyTuple_SET_ITEM(__pyx_t_7, 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 = 70; __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 = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_3 = PyObject_Call(__pyx_t_6, ((PyObject *)__pyx_t_7), ((PyObject *)__pyx_t_1)); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_7)); __pyx_t_7 = 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 = 70; __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_pybuffernd_outCount.rcbuffer->pybuffer, (PyObject*)__pyx_t_11, &__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) {
      __pyx_v_outCount = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_pybuffernd_outCount.diminfo[0].strides = __pyx_pybuffernd_outCount.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_outCount.diminfo[0].shape = __pyx_pybuffernd_outCount.rcbuffer->pybuffer.shape[0];
    }
  }
  __pyx_t_11 = 0;
  __pyx_v_outCount = ((PyArrayObject *)__pyx_t_3);
  __pyx_t_3 = 0;
 71:     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outMerge = numpy.zeros(bins, dtype="float64")
  /* "histogram.pyx":71
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outData = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outMerge = numpy.zeros(bins, dtype="float64")             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outPos = numpy.zeros(bins, dtype="float64")
 *     cdef double bin_edge_min, bin_edge_max
 */
  __pyx_t_3 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 71; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_1 = __Pyx_PyObject_GetAttrStr(__pyx_t_3, __pyx_n_s__zeros); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 71; __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_bins); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 71; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_7 = PyTuple_New(1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 71; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  PyTuple_SET_ITEM(__pyx_t_7, 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 = 71; __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 = 71; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_6 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_t_7), ((PyObject *)__pyx_t_3)); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 71; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_6);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_7)); __pyx_t_7 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0;
  if (!(likely(((__pyx_t_6) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_6, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 71; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_12 = ((PyArrayObject *)__pyx_t_6);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_outMerge.rcbuffer->pybuffer, (PyObject*)__pyx_t_12, &__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) {
      __pyx_v_outMerge = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 71; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_pybuffernd_outMerge.diminfo[0].strides = __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_outMerge.diminfo[0].shape = __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.shape[0];
    }
  }
  __pyx_t_12 = 0;
  __pyx_v_outMerge = ((PyArrayObject *)__pyx_t_6);
  __pyx_t_6 = 0;
 72:     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outPos = numpy.zeros(bins, dtype="float64")
  /* "histogram.pyx":72
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outCount = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outMerge = numpy.zeros(bins, dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outPos = numpy.zeros(bins, dtype="float64")             # <<<<<<<<<<<<<<
 *     cdef double bin_edge_min, bin_edge_max
 *     if bin_range is not None:
 */
  __pyx_t_6 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_6);
  __pyx_t_3 = __Pyx_PyObject_GetAttrStr(__pyx_t_6, __pyx_n_s__zeros); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
  __pyx_t_6 = PyInt_FromLong(__pyx_v_bins); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_6);
  __pyx_t_7 = PyTuple_New(1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  PyTuple_SET_ITEM(__pyx_t_7, 0, __pyx_t_6);
  __Pyx_GIVEREF(__pyx_t_6);
  __pyx_t_6 = 0;
  __pyx_t_6 = PyDict_New(); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_6));
  if (PyDict_SetItem(__pyx_t_6, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64)) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_1 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_t_7), ((PyObject *)__pyx_t_6)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_7)); __pyx_t_7 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_6)); __pyx_t_6 = 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 = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_13 = ((PyArrayObject *)__pyx_t_1);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_outPos.rcbuffer->pybuffer, (PyObject*)__pyx_t_13, &__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) {
      __pyx_v_outPos = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_pybuffernd_outPos.rcbuffer->pybuffer.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_pybuffernd_outPos.diminfo[0].strides = __pyx_pybuffernd_outPos.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_outPos.diminfo[0].shape = __pyx_pybuffernd_outPos.rcbuffer->pybuffer.shape[0];
    }
  }
  __pyx_t_13 = 0;
  __pyx_v_outPos = ((PyArrayObject *)__pyx_t_1);
  __pyx_t_1 = 0;
 73:     cdef double bin_edge_min, bin_edge_max
 74:     if bin_range is not None:
  /* "histogram.pyx":74
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] outPos = numpy.zeros(bins, dtype="float64")
 *     cdef double bin_edge_min, bin_edge_max
 *     if bin_range is not None:             # <<<<<<<<<<<<<<
 *         bin_edge_min = bin_range[0]
 *         bin_edge_max = bin_range[1] * EPS32
 */
  __pyx_t_4 = (__pyx_v_bin_range != Py_None);
  if (__pyx_t_4) {
 75:         bin_edge_min = bin_range[0]
    /* "histogram.pyx":75
 *     cdef double bin_edge_min, bin_edge_max
 *     if bin_range is not None:
 *         bin_edge_min = bin_range[0]             # <<<<<<<<<<<<<<
 *         bin_edge_max = bin_range[1] * EPS32
 *     else:
 */
    __pyx_t_1 = __Pyx_GetItemInt(__pyx_v_bin_range, 0, sizeof(long), PyInt_FromLong, 0, 0, 0); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 75; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __pyx_t_14 = __pyx_PyFloat_AsDouble(__pyx_t_1); if (unlikely((__pyx_t_14 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 75; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
    __pyx_v_bin_edge_min = __pyx_t_14;
 76:         bin_edge_max = bin_range[1] * EPS32
    /* "histogram.pyx":76
 *     if bin_range is not None:
 *         bin_edge_min = bin_range[0]
 *         bin_edge_max = bin_range[1] * EPS32             # <<<<<<<<<<<<<<
 *     else:
 *         bin_edge_min = pos.min()
 */
    __pyx_t_1 = __Pyx_GetItemInt(__pyx_v_bin_range, 1, sizeof(long), PyInt_FromLong, 0, 0, 0); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 76; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __pyx_t_6 = __Pyx_GetModuleGlobalName(__pyx_n_s__EPS32); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 76; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_6);
    __pyx_t_7 = PyNumber_Multiply(__pyx_t_1, __pyx_t_6); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 76; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_7);
    __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
    __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
    __pyx_t_14 = __pyx_PyFloat_AsDouble(__pyx_t_7); if (unlikely((__pyx_t_14 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 76; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
    __pyx_v_bin_edge_max = __pyx_t_14;
    goto __pyx_L3;
  }
  /*else*/ {
 77:     else:
 78:         bin_edge_min = pos.min()
    /* "histogram.pyx":78
 *         bin_edge_max = bin_range[1] * EPS32
 *     else:
 *         bin_edge_min = pos.min()             # <<<<<<<<<<<<<<
 *         bin_edge_max = pos.max() * EPS32
 *     cdef double bin_width = (bin_edge_max - bin_edge_min) / (< double > (bins))
 */
    __pyx_t_7 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos), __pyx_n_s__min); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_7);
    __pyx_t_6 = PyObject_Call(__pyx_t_7, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_6);
    __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
    __pyx_t_14 = __pyx_PyFloat_AsDouble(__pyx_t_6); if (unlikely((__pyx_t_14 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
    __pyx_v_bin_edge_min = __pyx_t_14;
 79:         bin_edge_max = pos.max() * EPS32
    /* "histogram.pyx":79
 *     else:
 *         bin_edge_min = pos.min()
 *         bin_edge_max = pos.max() * EPS32             # <<<<<<<<<<<<<<
 *     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_6 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos), __pyx_n_s__max); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_6);
    __pyx_t_7 = PyObject_Call(__pyx_t_6, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_7);
    __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
    __pyx_t_6 = __Pyx_GetModuleGlobalName(__pyx_n_s__EPS32); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_6);
    __pyx_t_1 = PyNumber_Multiply(__pyx_t_7, __pyx_t_6); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
    __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
    __pyx_t_14 = __pyx_PyFloat_AsDouble(__pyx_t_1); if (unlikely((__pyx_t_14 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
    __pyx_v_bin_edge_max = __pyx_t_14;
  }
  __pyx_L3:;
 80:     cdef double bin_width = (bin_edge_max - bin_edge_min) / (< double > (bins))
  /* "histogram.pyx":80
 *         bin_edge_min = pos.min()
 *         bin_edge_max = pos.max() * EPS32
 *     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));
 81:     cdef double inv_bin_width = (< double > (bins)) / (bin_edge_max - bin_edge_min)
  /* "histogram.pyx":81
 *         bin_edge_max = pos.max() * EPS32
 *     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));
 82:     cdef double a = 0.0
  /* "histogram.pyx":82
 *     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;
 83:     cdef double d = 0.0
  /* "histogram.pyx":83
 *     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;
 84:     cdef double fbin = 0.0
  /* "histogram.pyx":84
 *     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;
 85:     cdef double ffbin = 0.0
  /* "histogram.pyx":85
 *     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;
 86:     cdef double dInt = 0.0
  /* "histogram.pyx":86
 *     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;
 87:     cdef double dIntR = 0.0
  /* "histogram.pyx":87
 *     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;
 88:     cdef double dIntL = 0.0
  /* "histogram.pyx":88
 *     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;
 89:     cdef double dtmp = 0.0
  /* "histogram.pyx":89
 *     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;
 90:     cdef double dbin, inv_dbin2 = 0.0
  /* "histogram.pyx":90
 *     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;
 91:     cdef double tmp_count, tmp_data = 0.0
  /* "histogram.pyx":91
 *     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;
 92:     cdef double epsilon = 1e-10
  /* "histogram.pyx":92
 *     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;
 93: 
 94:     cdef long   bin = 0
  /* "histogram.pyx":94
 *     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;
 95:     cdef long   i, idx, t, dest = 0
  /* "histogram.pyx":95
 * 
 *     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;
 96:     if nthread is not None:
  /* "histogram.pyx":96
 *     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) {
 97:         if isinstance(nthread, int) and (nthread > 0):
    /* "histogram.pyx":97
 *     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)
 *     #multi-threaded version of this algorithm fails on MacOSX
 */
    __pyx_t_4 = PyInt_Check(__pyx_v_nthread); 
    if (__pyx_t_4) {
      __pyx_t_1 = PyObject_RichCompare(__pyx_v_nthread, __pyx_int_0, Py_GT); __Pyx_XGOTREF(__pyx_t_1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __pyx_t_15 = __Pyx_PyObject_IsTrue(__pyx_t_1); if (unlikely(__pyx_t_15 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
      __pyx_t_16 = __pyx_t_15;
    } else {
      __pyx_t_16 = __pyx_t_4;
    }
    if (__pyx_t_16) {
 98:             omp_set_num_threads(< int > nthread)
      /* "histogram.pyx":98
 *     if nthread is not None:
 *         if isinstance(nthread, int) and (nthread > 0):
 *             omp_set_num_threads(< int > nthread)             # <<<<<<<<<<<<<<
 *     #multi-threaded version of this algorithm fails on MacOSX
 *     if sys.platform == "darwin":
 */
      __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 = 98; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      omp_set_num_threads(((int)__pyx_t_17));
      goto __pyx_L5;
    }
    __pyx_L5:;
    goto __pyx_L4;
  }
  __pyx_L4:;
 99:     #multi-threaded version of this algorithm fails on MacOSX
 100:     if sys.platform == "darwin":
  /* "histogram.pyx":100
 *             omp_set_num_threads(< int > nthread)
 *     #multi-threaded version of this algorithm fails on MacOSX
 *     if sys.platform == "darwin":             # <<<<<<<<<<<<<<
 *          omp_set_num_threads(< int > 1)
 *     cdef double * bigCount = < double *> calloc(bins * omp_get_max_threads(), sizeof(double))
 */
  __pyx_t_1 = __Pyx_GetModuleGlobalName(__pyx_n_s__sys); 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_t_6 = __Pyx_PyObject_GetAttrStr(__pyx_t_1, __pyx_n_s__platform); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_6);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_1 = PyObject_RichCompare(__pyx_t_6, ((PyObject *)__pyx_n_s__darwin), Py_EQ); __Pyx_XGOTREF(__pyx_t_1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
  __pyx_t_16 = __Pyx_PyObject_IsTrue(__pyx_t_1); if (unlikely(__pyx_t_16 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 100; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  if (__pyx_t_16) {
 101:          omp_set_num_threads(< int > 1)
    /* "histogram.pyx":101
 *     #multi-threaded version of this algorithm fails on MacOSX
 *     if sys.platform == "darwin":
 *          omp_set_num_threads(< int > 1)             # <<<<<<<<<<<<<<
 *     cdef double * bigCount = < double *> calloc(bins * omp_get_max_threads(), sizeof(double))
 *     cdef double * bigData = < double *> calloc(bins * omp_get_max_threads(), sizeof(double))
 */
    omp_set_num_threads(((int)1));
    goto __pyx_L6;
  }
  __pyx_L6:;
 102:     cdef double * bigCount = < double *> calloc(bins * omp_get_max_threads(), sizeof(double))
  /* "histogram.pyx":102
 *     if sys.platform == "darwin":
 *          omp_set_num_threads(< int > 1)
 *     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))));
 103:     cdef double * bigData = < double *> calloc(bins * omp_get_max_threads(), sizeof(double))
  /* "histogram.pyx":103
 *          omp_set_num_threads(< int > 1)
 *     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))));
 104:     if pixelSize_in_Pos is None:
  /* "histogram.pyx":104
 *     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) {
 105:         dbin = 0.5
    /* "histogram.pyx":105
 *     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;
 106:         inv_dbin2 = 4.0
    /* "histogram.pyx":106
 *     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_L7;
  }
 107:     elif isinstance(pixelSize_in_Pos, (int, float)):
  /* "histogram.pyx":107
 *         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_1 = __pyx_v_pixelSize_in_Pos;
  __pyx_t_16 = PyInt_Check(__pyx_t_1); 
  if (!__pyx_t_16) {
    __pyx_t_4 = PyFloat_Check(__pyx_t_1); 
    __pyx_t_15 = __pyx_t_4;
  } else {
    __pyx_t_15 = __pyx_t_16;
  }
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  if (__pyx_t_15) {
 108:         dbin = 0.5 * (< double > pixelSize_in_Pos) * inv_bin_width
    /* "histogram.pyx":108
 *         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 = 108; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __pyx_v_dbin = ((0.5 * ((double)__pyx_t_14)) * __pyx_v_inv_bin_width);
 109:         if dbin > 0.0:
    /* "histogram.pyx":109
 *     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) {
 110:             inv_dbin2 = 1 / dbin / dbin
      /* "histogram.pyx":110
 *         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_L8;
    }
    /*else*/ {
 111:         else:
 112:             inv_dbin2 = 0.0
      /* "histogram.pyx":112
 *             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_L8:;
    goto __pyx_L7;
  }
 113:     elif isinstance(pixelSize_in_Pos, numpy.ndarray):
  /* "histogram.pyx":113
 *         else:
 *             inv_dbin2 = 0.0
 *     elif isinstance(pixelSize_in_Pos, numpy.ndarray):             # <<<<<<<<<<<<<<
 *         pass #TODO
 * 
 */
  __pyx_t_15 = __Pyx_TypeCheck(__pyx_v_pixelSize_in_Pos, ((PyObject*)__pyx_ptype_5numpy_ndarray)); 
  if (__pyx_t_15) {
    goto __pyx_L7;
  }
  __pyx_L7:;
 114:         pass #TODO
 115: 
 116:     if numpy.isnan(dbin) or numpy.isnan(inv_dbin2):
  /* "histogram.pyx":116
 *         pass #TODO
 * 
 *     if numpy.isnan(dbin) or numpy.isnan(inv_dbin2):             # <<<<<<<<<<<<<<
 *         dbin = 0.0
 *         inv_dbin2 = 0.0
 */
  __pyx_t_1 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_6 = __Pyx_PyObject_GetAttrStr(__pyx_t_1, __pyx_n_s__isnan); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_6);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_1 = PyFloat_FromDouble(__pyx_v_dbin); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_7 = PyTuple_New(1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  PyTuple_SET_ITEM(__pyx_t_7, 0, __pyx_t_1);
  __Pyx_GIVEREF(__pyx_t_1);
  __pyx_t_1 = 0;
  __pyx_t_1 = PyObject_Call(__pyx_t_6, ((PyObject *)__pyx_t_7), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_7)); __pyx_t_7 = 0;
  __pyx_t_15 = __Pyx_PyObject_IsTrue(__pyx_t_1); if (unlikely(__pyx_t_15 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  if (!__pyx_t_15) {
    __pyx_t_1 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __pyx_t_7 = __Pyx_PyObject_GetAttrStr(__pyx_t_1, __pyx_n_s__isnan); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_7);
    __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
    __pyx_t_1 = PyFloat_FromDouble(__pyx_v_inv_dbin2); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __pyx_t_6 = PyTuple_New(1); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_6);
    PyTuple_SET_ITEM(__pyx_t_6, 0, __pyx_t_1);
    __Pyx_GIVEREF(__pyx_t_1);
    __pyx_t_1 = 0;
    __pyx_t_1 = PyObject_Call(__pyx_t_7, ((PyObject *)__pyx_t_6), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
    __Pyx_DECREF(((PyObject *)__pyx_t_6)); __pyx_t_6 = 0;
    __pyx_t_16 = __Pyx_PyObject_IsTrue(__pyx_t_1); if (unlikely(__pyx_t_16 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 116; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
    __pyx_t_4 = __pyx_t_16;
  } else {
    __pyx_t_4 = __pyx_t_15;
  }
  if (__pyx_t_4) {
 117:         dbin = 0.0
    /* "histogram.pyx":117
 * 
 *     if numpy.isnan(dbin) or numpy.isnan(inv_dbin2):
 *         dbin = 0.0             # <<<<<<<<<<<<<<
 *         inv_dbin2 = 0.0
 * 
 */
    __pyx_v_dbin = 0.0;
 118:         inv_dbin2 = 0.0
    /* "histogram.pyx":118
 *     if numpy.isnan(dbin) or numpy.isnan(inv_dbin2):
 *         dbin = 0.0
 *         inv_dbin2 = 0.0             # <<<<<<<<<<<<<<
 * 
 *     with nogil:
 */
    __pyx_v_inv_dbin2 = 0.0;
    goto __pyx_L9;
  }
  __pyx_L9:;
 119: 
 120:     with nogil:
  /* "histogram.pyx":120
 *         inv_dbin2 = 0.0
 * 
 *     with nogil:             # <<<<<<<<<<<<<<
 *         for i in prange(size):
 *             d = cdata[i]
 */
  {
      #ifdef WITH_THREAD
      PyThreadState *_save;
      Py_UNBLOCK_THREADS
      #endif
      /*try:*/ {

      /* "histogram.pyx":120
 *         inv_dbin2 = 0.0
 * 
 *     with nogil:             # <<<<<<<<<<<<<<
 *         for i in prange(size):
 *             d = cdata[i]
 */
      /*finally:*/ {
        #ifdef WITH_THREAD
        Py_BLOCK_THREADS
        #endif
      }
  }
 121:         for i in prange(size):
        /* "histogram.pyx":121
 * 
 *     with nogil:
 *         for i in prange(size):             # <<<<<<<<<<<<<<
 *             d = cdata[i]
 *             a = cpos[i]
 */
        __pyx_t_5 = __pyx_v_size;
        if (1 == 0) abort();
        {
            #if ((defined(__APPLE__) || defined(__OSX__)) && (defined(__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)))))
                #undef likely
                #undef unlikely
                #define likely(x)   (x)
                #define unlikely(x) (x)
            #endif
            __pyx_t_19 = (__pyx_t_5 - 0) / 1;
            if (__pyx_t_19 > 0)
            {
                #ifdef _OPENMP
                #pragma omp parallel
                #endif /* _OPENMP */
                {
                    #ifdef _OPENMP
                    #pragma omp for lastprivate(__pyx_v_a) lastprivate(__pyx_v_dtmp) lastprivate(__pyx_v_ffbin) lastprivate(__pyx_v_fbin) lastprivate(__pyx_v_bin) lastprivate(__pyx_v_dIntR) lastprivate(__pyx_v_d) lastprivate(__pyx_v_dInt) firstprivate(__pyx_v_i) lastprivate(__pyx_v_i) lastprivate(__pyx_v_dest) lastprivate(__pyx_v_dIntL)
                    #endif /* _OPENMP */
                    for (__pyx_t_18 = 0; __pyx_t_18 < __pyx_t_19; __pyx_t_18++){
                        {
                            __pyx_v_i = 0 + 1 * __pyx_t_18;
                            /* Initialize private variables to invalid values */
                            __pyx_v_a = ((double)__PYX_NAN());
                            __pyx_v_dtmp = ((double)__PYX_NAN());
                            __pyx_v_ffbin = ((double)__PYX_NAN());
                            __pyx_v_fbin = ((double)__PYX_NAN());
                            __pyx_v_bin = ((long)0xbad0bad0);
                            __pyx_v_dIntR = ((double)__PYX_NAN());
                            __pyx_v_d = ((double)__PYX_NAN());
                            __pyx_v_dInt = ((double)__PYX_NAN());
                            __pyx_v_dest = ((long)0xbad0bad0);
                            __pyx_v_dIntL = ((double)__PYX_NAN());
 122:             d = cdata[i]
                            /* "histogram.pyx":122
 *     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_20 = __pyx_v_i;
                            __pyx_v_d = (*((double *) ( /* dim=0 */ (__pyx_v_cdata.data + __pyx_t_20 * __pyx_v_cdata.strides[0]) )));
 123:             a = cpos[i]
                            /* "histogram.pyx":123
 *         for i in prange(size):
 *             d = cdata[i]
 *             a = cpos[i]             # <<<<<<<<<<<<<<
 *             if (a < bin_edge_min) or (a > bin_edge_max):
 *                 continue
 */
                            __pyx_t_21 = __pyx_v_i;
                            __pyx_v_a = (*((double *) ( /* dim=0 */ (__pyx_v_cpos.data + __pyx_t_21 * __pyx_v_cpos.strides[0]) )));
 124:             if (a < bin_edge_min) or (a > bin_edge_max):
                            /* "histogram.pyx":124
 *             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_4 = (__pyx_v_a < __pyx_v_bin_edge_min);
                            if (!__pyx_t_4) {
                              __pyx_t_15 = (__pyx_v_a > __pyx_v_bin_edge_max);
                              __pyx_t_16 = __pyx_t_15;
                            } else {
                              __pyx_t_16 = __pyx_t_4;
                            }
                            if (__pyx_t_16) {
 125:                 continue
                              /* "histogram.pyx":125
 *             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_L13_continue;
                              goto __pyx_L17;
                            }
                            __pyx_L17:;
 126:             fbin = (a - bin_edge_min) * inv_bin_width
                            /* "histogram.pyx":126
 *             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);
 127:             ffbin = floor(fbin)
                            /* "histogram.pyx":127
 *                 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);
 128:             bin = < long > ffbin
                            /* "histogram.pyx":128
 *             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);
 129:             dest = omp_get_thread_num() * bins + bin
                            /* "histogram.pyx":129
 *             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);
 130:             dInt = 1.0
                            /* "histogram.pyx":130
 *             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;
 131:             if  bin > 0 :
                            /* "histogram.pyx":131
 *             dest = omp_get_thread_num() * bins + bin
 *             dInt = 1.0
 *             if  bin > 0 :             # <<<<<<<<<<<<<<
 *                 dtmp = ffbin - (fbin - dbin)
 *                 if dtmp > 0:
 */
                            __pyx_t_16 = (__pyx_v_bin > 0);
                            if (__pyx_t_16) {
 132:                 dtmp = ffbin - (fbin - dbin)
                              /* "histogram.pyx":132
 *             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));
 133:                 if dtmp > 0:
                              /* "histogram.pyx":133
 *             if  bin > 0 :
 *                 dtmp = ffbin - (fbin - dbin)
 *                 if dtmp > 0:             # <<<<<<<<<<<<<<
 *                     dIntL = 0.5 * dtmp * dtmp * inv_dbin2
 *                     dInt = dInt - dIntL
 */
                              __pyx_t_16 = (__pyx_v_dtmp > 0.0);
                              if (__pyx_t_16) {
 134:                     dIntL = 0.5 * dtmp * dtmp * inv_dbin2
                                /* "histogram.pyx":134
 *                 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);
 135:                     dInt = dInt - dIntL
                                /* "histogram.pyx":135
 *                 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);
 136:                     bigCount[dest - 1] += dIntL
                                /* "histogram.pyx":136
 *                     dIntL = 0.5 * dtmp * dtmp * inv_dbin2
 *                     dInt = dInt - dIntL
 *                     bigCount[dest - 1] += dIntL             # <<<<<<<<<<<<<<
 *                     bigData[dest - 1] += d * dIntL
 * 
 */
                                __pyx_t_22 = (__pyx_v_dest - 1);
                                (__pyx_v_bigCount[__pyx_t_22]) = ((__pyx_v_bigCount[__pyx_t_22]) + __pyx_v_dIntL);
 137:                     bigData[dest - 1] += d * dIntL
                                /* "histogram.pyx":137
 *                     dInt = dInt - dIntL
 *                     bigCount[dest - 1] += dIntL
 *                     bigData[dest - 1] += d * dIntL             # <<<<<<<<<<<<<<
 * 
 *             if bin < bins - 1 :
 */
                                __pyx_t_22 = (__pyx_v_dest - 1);
                                (__pyx_v_bigData[__pyx_t_22]) = ((__pyx_v_bigData[__pyx_t_22]) + (__pyx_v_d * __pyx_v_dIntL));
                                goto __pyx_L19;
                              }
                              __pyx_L19:;
                              goto __pyx_L18;
                            }
                            __pyx_L18:;
 138: 
 139:             if bin < bins - 1 :
                            /* "histogram.pyx":139
 *                     bigData[dest - 1] += d * dIntL
 * 
 *             if bin < bins - 1 :             # <<<<<<<<<<<<<<
 *                 dtmp = fbin + dbin - ffbin - 1
 *                 if dtmp > 0 :
 */
                            __pyx_t_16 = (__pyx_v_bin < (__pyx_v_bins - 1));
                            if (__pyx_t_16) {
 140:                 dtmp = fbin + dbin - ffbin - 1
                              /* "histogram.pyx":140
 * 
 *             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);
 141:                 if dtmp > 0 :
                              /* "histogram.pyx":141
 *             if bin < bins - 1 :
 *                 dtmp = fbin + dbin - ffbin - 1
 *                 if dtmp > 0 :             # <<<<<<<<<<<<<<
 *                     dIntR = 0.5 * dtmp * dtmp * inv_dbin2
 *                     dInt = dInt - dIntR
 */
                              __pyx_t_16 = (__pyx_v_dtmp > 0.0);
                              if (__pyx_t_16) {
 142:                     dIntR = 0.5 * dtmp * dtmp * inv_dbin2
                                /* "histogram.pyx":142
 *                 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);
 143:                     dInt = dInt - dIntR
                                /* "histogram.pyx":143
 *                 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);
 144:                     bigCount[dest + 1] += dIntR
                                /* "histogram.pyx":144
 *                     dIntR = 0.5 * dtmp * dtmp * inv_dbin2
 *                     dInt = dInt - dIntR
 *                     bigCount[dest + 1] += dIntR             # <<<<<<<<<<<<<<
 *                     bigData[dest + 1] += d * dIntR
 *             bigCount[dest] += dInt
 */
                                __pyx_t_22 = (__pyx_v_dest + 1);
                                (__pyx_v_bigCount[__pyx_t_22]) = ((__pyx_v_bigCount[__pyx_t_22]) + __pyx_v_dIntR);
 145:                     bigData[dest + 1] += d * dIntR
                                /* "histogram.pyx":145
 *                     dInt = dInt - dIntR
 *                     bigCount[dest + 1] += dIntR
 *                     bigData[dest + 1] += d * dIntR             # <<<<<<<<<<<<<<
 *             bigCount[dest] += dInt
 *             bigData[dest] += d * dInt
 */
                                __pyx_t_22 = (__pyx_v_dest + 1);
                                (__pyx_v_bigData[__pyx_t_22]) = ((__pyx_v_bigData[__pyx_t_22]) + (__pyx_v_d * __pyx_v_dIntR));
                                goto __pyx_L21;
                              }
                              __pyx_L21:;
                              goto __pyx_L20;
                            }
                            __pyx_L20:;
 146:             bigCount[dest] += dInt
                            /* "histogram.pyx":146
 *                     bigCount[dest + 1] += dIntR
 *                     bigData[dest + 1] += d * dIntR
 *             bigCount[dest] += dInt             # <<<<<<<<<<<<<<
 *             bigData[dest] += d * dInt
 * 
 */
                            __pyx_t_22 = __pyx_v_dest;
                            (__pyx_v_bigCount[__pyx_t_22]) = ((__pyx_v_bigCount[__pyx_t_22]) + __pyx_v_dInt);
 147:             bigData[dest] += d * dInt
                            /* "histogram.pyx":147
 *                     bigData[dest + 1] += d * dIntR
 *             bigCount[dest] += dInt
 *             bigData[dest] += d * dInt             # <<<<<<<<<<<<<<
 * 
 *         for idx in prange(bins):
 */
                            __pyx_t_22 = __pyx_v_dest;
                            (__pyx_v_bigData[__pyx_t_22]) = ((__pyx_v_bigData[__pyx_t_22]) + (__pyx_v_d * __pyx_v_dInt));
                            goto __pyx_L23;
                            __pyx_L13_continue:;
                            goto __pyx_L23;
                            __pyx_L23:;
                        }
                    }
                }
            }
        }
        #if ((defined(__APPLE__) || defined(__OSX__)) && (defined(__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)))))
            #undef likely
            #undef unlikely
            #define likely(x)   __builtin_expect(!!(x), 1)
            #define unlikely(x) __builtin_expect(!!(x), 0)
        #endif
 148: 
 149:         for idx in prange(bins):
        /* "histogram.pyx":149
 *             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_19 = __pyx_v_bins;
        if (1 == 0) abort();
        {
            #if ((defined(__APPLE__) || defined(__OSX__)) && (defined(__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)))))
                #undef likely
                #undef unlikely
                #define likely(x)   (x)
                #define unlikely(x) (x)
            #endif
            __pyx_t_5 = (__pyx_t_19 - 0) / 1;
            if (__pyx_t_5 > 0)
            {
                #ifdef _OPENMP
                #pragma omp parallel
                #endif /* _OPENMP */
                {
                    #ifdef _OPENMP
                    #pragma omp for reduction(+:__pyx_v_tmp_count) reduction(+:__pyx_v_tmp_data) firstprivate(__pyx_v_idx) lastprivate(__pyx_v_idx) lastprivate(__pyx_v_t) lastprivate(__pyx_v_dest)
                    #endif /* _OPENMP */
                    for (__pyx_t_18 = 0; __pyx_t_18 < __pyx_t_5; __pyx_t_18++){
                        {
                            __pyx_v_idx = 0 + 1 * __pyx_t_18;
                            /* Initialize private variables to invalid values */
                            __pyx_v_t = ((long)0xbad0bad0);
                            __pyx_v_dest = ((long)0xbad0bad0);
 150:             outPos[idx] = bin_edge_min + (0.5 +< double > idx) * bin_width
                            /* "histogram.pyx":150
 * 
 *         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_22 = __pyx_v_idx;
                            *__Pyx_BufPtrStrided1d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outPos.rcbuffer->pybuffer.buf, __pyx_t_22, __pyx_pybuffernd_outPos.diminfo[0].strides) = (__pyx_v_bin_edge_min + ((0.5 + ((double)__pyx_v_idx)) * __pyx_v_bin_width));
 151:             tmp_count = 0.0
                            /* "histogram.pyx":151
 *         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;
 152:             tmp_data = 0.0
                            /* "histogram.pyx":152
 *             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;
 153:             for t in range(omp_get_max_threads()):
                            /* "histogram.pyx":153
 *             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_23 = 0; __pyx_t_23 < __pyx_t_17; __pyx_t_23+=1) {
                              __pyx_v_t = __pyx_t_23;
 154:                 dest = t * bins + idx
                              /* "histogram.pyx":154
 *             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);
 155:                 tmp_count += bigCount[dest]
                              /* "histogram.pyx":155
 *             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]));
 156:                 tmp_data += bigData[dest]
                              /* "histogram.pyx":156
 *                 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]));
                            }
 157:             outCount[idx] += tmp_count
                            /* "histogram.pyx":157
 *                 tmp_count += bigCount[dest]
 *                 tmp_data += bigData[dest]
 *             outCount[idx] += tmp_count             # <<<<<<<<<<<<<<
 *             outData[idx] += tmp_data
 *             if outCount[idx] > epsilon:
 */
                            __pyx_t_23 = __pyx_v_idx;
                            *__Pyx_BufPtrStrided1d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_23, __pyx_pybuffernd_outCount.diminfo[0].strides) += __pyx_v_tmp_count;
 158:             outData[idx] += tmp_data
                            /* "histogram.pyx":158
 *                 tmp_data += bigData[dest]
 *             outCount[idx] += tmp_count
 *             outData[idx] += tmp_data             # <<<<<<<<<<<<<<
 *             if outCount[idx] > epsilon:
 *                 outMerge[idx] += tmp_data / tmp_count
 */
                            __pyx_t_24 = __pyx_v_idx;
                            *__Pyx_BufPtrStrided1d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_24, __pyx_pybuffernd_outData.diminfo[0].strides) += __pyx_v_tmp_data;
 159:             if outCount[idx] > epsilon:
                            /* "histogram.pyx":159
 *             outCount[idx] += tmp_count
 *             outData[idx] += tmp_data
 *             if outCount[idx] > epsilon:             # <<<<<<<<<<<<<<
 *                 outMerge[idx] += tmp_data / tmp_count
 *             else:
 */
                            __pyx_t_25 = __pyx_v_idx;
                            __pyx_t_16 = ((*__Pyx_BufPtrStrided1d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_25, __pyx_pybuffernd_outCount.diminfo[0].strides)) > __pyx_v_epsilon);
                            if (__pyx_t_16) {
 160:                 outMerge[idx] += tmp_data / tmp_count
                              /* "histogram.pyx":160
 *             outData[idx] += tmp_data
 *             if outCount[idx] > epsilon:
 *                 outMerge[idx] += tmp_data / tmp_count             # <<<<<<<<<<<<<<
 *             else:
 *                 outMerge[idx] += dummy
 */
                              __pyx_t_26 = __pyx_v_idx;
                              *__Pyx_BufPtrStrided1d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.buf, __pyx_t_26, __pyx_pybuffernd_outMerge.diminfo[0].strides) += (__pyx_v_tmp_data / __pyx_v_tmp_count);
                              goto __pyx_L30;
                            }
                            /*else*/ {
 161:             else:
 162:                 outMerge[idx] += dummy
                              /* "histogram.pyx":162
 *                 outMerge[idx] += tmp_data / tmp_count
 *             else:
 *                 outMerge[idx] += dummy             # <<<<<<<<<<<<<<
 * 
 *     free(bigCount)
 */
                              __pyx_t_27 = __pyx_v_idx;
                              *__Pyx_BufPtrStrided1d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.buf, __pyx_t_27, __pyx_pybuffernd_outMerge.diminfo[0].strides) += __pyx_v_dummy;
                            }
                            __pyx_L30:;
                        }
                    }
                }
            }
        }
        #if ((defined(__APPLE__) || defined(__OSX__)) && (defined(__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)))))
            #undef likely
            #undef unlikely
            #define likely(x)   __builtin_expect(!!(x), 1)
            #define unlikely(x) __builtin_expect(!!(x), 0)
        #endif
      }
 163: 
 164:     free(bigCount)
  /* "histogram.pyx":164
 *                 outMerge[idx] += dummy
 * 
 *     free(bigCount)             # <<<<<<<<<<<<<<
 *     free(bigData)
 *     return  outPos, outMerge, outData, outCount
 */
  free(__pyx_v_bigCount);
 165:     free(bigData)
  /* "histogram.pyx":165
 * 
 *     free(bigCount)
 *     free(bigData)             # <<<<<<<<<<<<<<
 *     return  outPos, outMerge, outData, outCount
 * 
 */
  free(__pyx_v_bigData);
 166:     return  outPos, outMerge, outData, outCount
  /* "histogram.pyx":166
 *     free(bigCount)
 *     free(bigData)
 *     return  outPos, outMerge, outData, outCount             # <<<<<<<<<<<<<<
 * 
 * 
 */
  __Pyx_XDECREF(__pyx_r);
  __pyx_t_1 = PyTuple_New(4); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 166; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_INCREF(((PyObject *)__pyx_v_outPos));
  PyTuple_SET_ITEM(__pyx_t_1, 0, ((PyObject *)__pyx_v_outPos));
  __Pyx_GIVEREF(((PyObject *)__pyx_v_outPos));
  __Pyx_INCREF(((PyObject *)__pyx_v_outMerge));
  PyTuple_SET_ITEM(__pyx_t_1, 1, ((PyObject *)__pyx_v_outMerge));
  __Pyx_GIVEREF(((PyObject *)__pyx_v_outMerge));
  __Pyx_INCREF(((PyObject *)__pyx_v_outData));
  PyTuple_SET_ITEM(__pyx_t_1, 2, ((PyObject *)__pyx_v_outData));
  __Pyx_GIVEREF(((PyObject *)__pyx_v_outData));
  __Pyx_INCREF(((PyObject *)__pyx_v_outCount));
  PyTuple_SET_ITEM(__pyx_t_1, 3, ((PyObject *)__pyx_v_outCount));
  __Pyx_GIVEREF(((PyObject *)__pyx_v_outCount));
  __pyx_r = ((PyObject *)__pyx_t_1);
  __pyx_t_1 = 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_6);
  __Pyx_XDECREF(__pyx_t_7);
  __PYX_XDEC_MEMVIEW(&__pyx_t_8, 1);
  __PYX_XDEC_MEMVIEW(&__pyx_t_9, 1);
  { PyObject *__pyx_type, *__pyx_value, *__pyx_tb;
    __Pyx_ErrFetch(&__pyx_type, &__pyx_value, &__pyx_tb);
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outCount.rcbuffer->pybuffer);
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outData.rcbuffer->pybuffer);
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outMerge.rcbuffer->pybuffer);
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outPos.rcbuffer->pybuffer);
  __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_pybuffernd_outCount.rcbuffer->pybuffer);
  __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outData.rcbuffer->pybuffer);
  __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outMerge.rcbuffer->pybuffer);
  __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outPos.rcbuffer->pybuffer);
  __pyx_L2:;
  __PYX_XDEC_MEMVIEW(&__pyx_v_cpos, 1);
  __PYX_XDEC_MEMVIEW(&__pyx_v_cdata, 1);
  __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_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

/* Python wrapper */
static PyObject *__pyx_pw_9histogram_3histogram2d(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/
static char __pyx_doc_9histogram_2histogram2d[] = "\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_3histogram2d = {__Pyx_NAMESTR("histogram2d"), (PyCFunction)__pyx_pw_9histogram_3histogram2d, METH_VARARGS|METH_KEYWORDS, __Pyx_DOCSTR(__pyx_doc_9histogram_2histogram2d)};
static PyObject *__pyx_pw_9histogram_3histogram2d(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;
  PyObject *__pyx_r = 0;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("histogram2d (wrapper)", 0);
  {
    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};
    PyObject* values[7] = {0,0,0,0,0,0,0};
    values[4] = __pyx_k_1;
 167: 
 168: 
 169: @cython.cdivision(True)
 170: @cython.boundscheck(False)
 171: @cython.wraparound(False)
 172: def histogram2d(numpy.ndarray pos0 not None,
/* "histogram.pyx":172
 * @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_2histogram2d(CYTHON_UNUSED PyObject *__pyx_self, PyArrayObject *__pyx_v_pos0, PyArrayObject *__pyx_v_pos1, PyObject *__pyx_v_bins, PyArrayObject *__pyx_v_weights, PyObject *__pyx_v_split, PyObject *__pyx_v_nthread, 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;
  __Pyx_memviewslice __pyx_v_cpos0 = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_cpos1 = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_memviewslice __pyx_v_data = { 0, 0, { 0 }, { 0 }, { 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_epsilon;
  double __pyx_v_area;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_edges0;
  __Pyx_Buffer __pyx_pybuffer_edges0;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_edges1;
  __Pyx_Buffer __pyx_pybuffer_edges1;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_outCount;
  __Pyx_Buffer __pyx_pybuffer_outCount;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_outData;
  __Pyx_Buffer __pyx_pybuffer_outData;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_outMerge;
  __Pyx_Buffer __pyx_pybuffer_outMerge;
  PyObject *__pyx_r = NULL;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("histogram2d", 0);
  __pyx_pybuffer_outData.pybuffer.buf = NULL;
  __pyx_pybuffer_outData.refcount = 0;
  __pyx_pybuffernd_outData.data = NULL;
  __pyx_pybuffernd_outData.rcbuffer = &__pyx_pybuffer_outData;
  __pyx_pybuffer_outCount.pybuffer.buf = NULL;
  __pyx_pybuffer_outCount.refcount = 0;
  __pyx_pybuffernd_outCount.data = NULL;
  __pyx_pybuffernd_outCount.rcbuffer = &__pyx_pybuffer_outCount;
  __pyx_pybuffer_outMerge.pybuffer.buf = NULL;
  __pyx_pybuffer_outMerge.refcount = 0;
  __pyx_pybuffernd_outMerge.data = NULL;
  __pyx_pybuffernd_outMerge.rcbuffer = &__pyx_pybuffer_outMerge;
  __pyx_pybuffer_edges0.pybuffer.buf = NULL;
  __pyx_pybuffer_edges0.refcount = 0;
  __pyx_pybuffernd_edges0.data = NULL;
  __pyx_pybuffernd_edges0.rcbuffer = &__pyx_pybuffer_edges0;
  __pyx_pybuffer_edges1.pybuffer.buf = NULL;
  __pyx_pybuffer_edges1.refcount = 0;
  __pyx_pybuffernd_edges1.data = NULL;
  __pyx_pybuffernd_edges1.rcbuffer = &__pyx_pybuffer_edges1;

  /* "histogram.pyx":172
 * @cython.boundscheck(False)
 * @cython.wraparound(False)
 * def histogram2d(numpy.ndarray pos0 not None,             # <<<<<<<<<<<<<<
 *                 numpy.ndarray pos1 not None,
 *                 bins not None,
 */
  __pyx_k_tuple_52 = PyTuple_Pack(44, ((PyObject *)__pyx_n_s__pos0), ((PyObject *)__pyx_n_s__pos1), ((PyObject *)__pyx_n_s__bins), ((PyObject *)__pyx_n_s__weights), ((PyObject *)__pyx_n_s__split), ((PyObject *)__pyx_n_s__nthread), ((PyObject *)__pyx_n_s__dummy), ((PyObject *)__pyx_n_s__bin0), ((PyObject *)__pyx_n_s__bin1), ((PyObject *)__pyx_n_s__i), ((PyObject *)__pyx_n_s__j), ((PyObject *)__pyx_n_s__b0), ((PyObject *)__pyx_n_s__b1), ((PyObject *)__pyx_n_s__size), ((PyObject *)__pyx_n_s__csplit), ((PyObject *)__pyx_n_s__cpos0), ((PyObject *)__pyx_n_s__cpos1), ((PyObject *)__pyx_n_s__data), ((PyObject *)__pyx_n_s__outData), ((PyObject *)__pyx_n_s__outCount), ((PyObject *)__pyx_n_s__outMerge), ((PyObject *)__pyx_n_s__edges0), ((PyObject *)__pyx_n_s__edges1), ((PyObject *)__pyx_n_s__min0), ((PyObject *)__pyx_n_s__max0), ((PyObject *)__pyx_n_s__min1), ((PyObject *)__pyx_n_s__max1), ((PyObject *)__pyx_n_s__idp0), ((PyObject *)__pyx_n_s__idp1), ((PyObject *)__pyx_n_s__dbin0), ((PyObject *)__pyx_n_s__dbin1), ((PyObject *)__pyx_n_s__fbin0), ((PyObject *)__pyx_n_s__fbin1), ((PyObject *)__pyx_n_s__p0), ((PyObject *)__pyx_n_s__p1), ((PyObject *)__pyx_n_s__d), ((PyObject *)__pyx_n_s__rest), ((PyObject *)__pyx_n_s__delta0l), ((PyObject *)__pyx_n_s__delta0r), ((PyObject *)__pyx_n_s__delta1l), ((PyObject *)__pyx_n_s__delta1r), ((PyObject *)__pyx_n_s__aera), ((PyObject *)__pyx_n_s__epsilon), ((PyObject *)__pyx_n_s__area)); if (unlikely(!__pyx_k_tuple_52)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 172; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_k_tuple_52);
  __Pyx_GIVEREF(((PyObject *)__pyx_k_tuple_52));

  /* "histogram.pyx":172
 * @cython.boundscheck(False)
 * @cython.wraparound(False)
 * def histogram2d(numpy.ndarray pos0 not None,             # <<<<<<<<<<<<<<
 *                 numpy.ndarray pos1 not None,
 *                 bins not None,
 */
  __pyx_t_4 = PyCFunction_NewEx(&__pyx_mdef_9histogram_3histogram2d, NULL, __pyx_n_s__histogram); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 172; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_4);
  if (PyDict_SetItem(__pyx_d, __pyx_n_s__histogram2d, __pyx_t_4) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 172; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
  __pyx_k_codeobj_53 = (PyObject*)__Pyx_PyCode_New(7, 0, 44, 0, 0, __pyx_empty_bytes, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_k_tuple_52, __pyx_empty_tuple, __pyx_empty_tuple, __pyx_kp_s_51, __pyx_n_s__histogram2d, 172, __pyx_empty_bytes); if (unlikely(!__pyx_k_codeobj_53)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 172; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
 173:                 numpy.ndarray pos1 not None,
 174:                 bins not None,
 175:                 numpy.ndarray weights not None,
 176:                 split=True,
  /* "histogram.pyx":176
 *                 bins not None,
 *                 numpy.ndarray weights not None,
 *                 split=True,             # <<<<<<<<<<<<<<
 *                 nthread=None,
 *                 double dummy=0.0):
 */
  __pyx_t_4 = __Pyx_PyBool_FromLong(1); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 176; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_4);
  __pyx_k_1 = __pyx_t_4;
  __Pyx_GIVEREF(__pyx_t_4);
  __pyx_t_4 = 0;
 177:                 nthread=None,
    /* "histogram.pyx":177
 *                 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;
      const Py_ssize_t pos_args = PyTuple_GET_SIZE(__pyx_args);
      switch (pos_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 (pos_args) {
        case  0:
        if (likely((values[0] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__pos0)) != 0)) kw_args--;
        else goto __pyx_L5_argtuple_error;
        case  1:
        if (likely((values[1] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__pos1)) != 0)) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("histogram2d", 0, 4, 7, 1); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 172; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        }
        case  2:
        if (likely((values[2] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__bins)) != 0)) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("histogram2d", 0, 4, 7, 2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 172; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        }
        case  3:
        if (likely((values[3] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__weights)) != 0)) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("histogram2d", 0, 4, 7, 3); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 172; __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, pos_args, "histogram2d") < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 172; __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 = 178; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
    } else {
 178:                 double dummy=0.0):
      /* "histogram.pyx":178
 *                 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 = 172; __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:;
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_pos0), __pyx_ptype_5numpy_ndarray, 0, "pos0", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 172; __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 = 173; __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 = 174; __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 = 175; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_r = __pyx_pf_9histogram_2histogram2d(__pyx_self, __pyx_v_pos0, __pyx_v_pos1, __pyx_v_bins, __pyx_v_weights, __pyx_v_split, __pyx_v_nthread, __pyx_v_dummy);
  goto __pyx_L0;
  __pyx_L1_error:;
  __pyx_r = NULL;
  __pyx_L0:;
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}
 179:     """
 180:     Calculate 2D histogram of pos0,pos1 weighted by weights
 181: 
 182:     @param pos0: 2Theta array
 183:     @param pos1: Chi array
 184:     @param weights: array with intensities
 185:     @param bins: number of output bins int or 2-tuple of int
 186:     @param nthread: maximum number of thread to use. By default: maximum available.
 187:     One can also limit this with OMP_NUM_THREADS environment variable
 188: 
 189: 
 190:     @return  I, edges0, edges1, weighted histogram(2D), unweighted histogram (2D)
 191:     """
 192:     assert pos0.size == pos1.size
  /* "histogram.pyx":192
 *     @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 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos0), __pyx_n_s__size); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 192; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_2 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos1), __pyx_n_s__size); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 192; __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); __Pyx_XGOTREF(__pyx_t_3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 192; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __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 = 192; __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 = 192; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }
  #endif
 193: #    if weights is not No:
 194:     assert pos0.size == weights.size
  /* "histogram.pyx":194
 *     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 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos0), __pyx_n_s__size); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_2 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_weights), __pyx_n_s__size); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __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); __Pyx_XGOTREF(__pyx_t_1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 194; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __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 = 194; __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 = 194; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }
  #endif
 195:     cdef long  bin0, bin1, i, j, b0, b1
 196:     cdef long  size = pos0.size
  /* "histogram.pyx":196
 *     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 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos0), __pyx_n_s__size); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 196; __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 = 196; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_v_size = __pyx_t_5;
 197:     try:
  /* "histogram.pyx":197
 *     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:*/ {
 198:         bin0, bin1 = tuple(bins)
      /* "histogram.pyx":198
 *     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 = 198; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
      __Pyx_GOTREF(__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 = 198; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
      __Pyx_GOTREF(__pyx_t_2);
      __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
      if (1) {
        PyObject* sequence = __pyx_t_2;
        #if CYTHON_COMPILING_IN_CPYTHON
        Py_ssize_t size = Py_SIZE(sequence);
        #else
        Py_ssize_t size = PySequence_Size(sequence);
        #endif
        if (unlikely(size != 2)) {
          if (size > 2) __Pyx_RaiseTooManyValuesError(2);
          else if (size >= 0) __Pyx_RaiseNeedMoreValuesError(size);
          {__pyx_filename = __pyx_f[0]; __pyx_lineno = 198; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        }
        #if CYTHON_COMPILING_IN_CPYTHON
        __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);
        #else
        __pyx_t_1 = PySequence_ITEM(sequence, 0); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 198; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_3 = PySequence_ITEM(sequence, 1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 198; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        __Pyx_GOTREF(__pyx_t_3);
        #endif
        __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
      } else if (1) {
        __Pyx_RaiseNoneNotIterableError(); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 198; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
      } else
      {
        Py_ssize_t index = -1;
        __pyx_t_9 = PyObject_GetIter(__pyx_t_2); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 198; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        __Pyx_GOTREF(__pyx_t_9);
        __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
        __pyx_t_10 = Py_TYPE(__pyx_t_9)->tp_iternext;
        index = 0; __pyx_t_1 = __pyx_t_10(__pyx_t_9); if (unlikely(!__pyx_t_1)) goto __pyx_L11_unpacking_failed;
        __Pyx_GOTREF(__pyx_t_1);
        index = 1; __pyx_t_3 = __pyx_t_10(__pyx_t_9); if (unlikely(!__pyx_t_3)) goto __pyx_L11_unpacking_failed;
        __Pyx_GOTREF(__pyx_t_3);
        if (__Pyx_IternextUnpackEndCheck(__pyx_t_10(__pyx_t_9), 2) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 198; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        __pyx_t_10 = NULL;
        __Pyx_DECREF(__pyx_t_9); __pyx_t_9 = 0;
        goto __pyx_L12_unpacking_done;
        __pyx_L11_unpacking_failed:;
        __Pyx_DECREF(__pyx_t_9); __pyx_t_9 = 0;
        __pyx_t_10 = NULL;
        if (__Pyx_IterFinish() == 0) __Pyx_RaiseNeedMoreValuesError(index);
        {__pyx_filename = __pyx_f[0]; __pyx_lineno = 198; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        __pyx_L12_unpacking_done:;
      }
      __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 = 198; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
      __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
      __pyx_t_11 = __Pyx_PyInt_AsLong(__pyx_t_3); if (unlikely((__pyx_t_11 == (long)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 198; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_v_bin0 = __pyx_t_5;
      __pyx_v_bin1 = __pyx_t_11;
    }
    __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_L10_try_end;
    __pyx_L3_error:;
    __Pyx_XDECREF(__pyx_t_9); __pyx_t_9 = 0;
    __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;
 199:     except:
    /* "histogram.pyx":199
 *     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 = 199; __pyx_clineno = __LINE__; goto __pyx_L5_except_error;}
      __Pyx_GOTREF(__pyx_t_2);
      __Pyx_GOTREF(__pyx_t_3);
      __Pyx_GOTREF(__pyx_t_1);
 200:         bin0 = bin1 = < long > bins
      /* "histogram.pyx":200
 *         bin0, bin1 = tuple(bins)
 *     except:
 *         bin0 = bin1 = < long > bins             # <<<<<<<<<<<<<<
 *     if bin0 <= 0:
 *         bin0 = 1
 */
      __pyx_t_11 = __Pyx_PyInt_AsLong(__pyx_v_bins); if (unlikely((__pyx_t_11 == (long)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 200; __pyx_clineno = __LINE__; goto __pyx_L5_except_error;}
      __pyx_v_bin0 = ((long)__pyx_t_11);
      __pyx_v_bin1 = ((long)__pyx_t_11);
      __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_L4_exception_handled;
    }
    __pyx_L5_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_L4_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_L10_try_end:;
  }
 201:     if bin0 <= 0:
  /* "histogram.pyx":201
 *     except:
 *         bin0 = bin1 = < long > bins
 *     if bin0 <= 0:             # <<<<<<<<<<<<<<
 *         bin0 = 1
 *     if bin1 <= 0:
 */
  __pyx_t_4 = (__pyx_v_bin0 <= 0);
  if (__pyx_t_4) {
 202:         bin0 = 1
    /* "histogram.pyx":202
 *         bin0 = bin1 = < long > bins
 *     if bin0 <= 0:
 *         bin0 = 1             # <<<<<<<<<<<<<<
 *     if bin1 <= 0:
 *         bin1 = 1
 */
    __pyx_v_bin0 = 1;
    goto __pyx_L15;
  }
  __pyx_L15:;
 203:     if bin1 <= 0:
  /* "histogram.pyx":203
 *     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) {
 204:         bin1 = 1
    /* "histogram.pyx":204
 *         bin0 = 1
 *     if bin1 <= 0:
 *         bin1 = 1             # <<<<<<<<<<<<<<
 *     cdef int csplit = split
 *     cdef double[:] cpos0 = numpy.ascontiguousarray(pos0.ravel(),dtype=numpy.float64)
 */
    __pyx_v_bin1 = 1;
    goto __pyx_L16;
  }
  __pyx_L16:;
 205:     cdef int csplit = split
  /* "histogram.pyx":205
 *     if bin1 <= 0:
 *         bin1 = 1
 *     cdef int csplit = split             # <<<<<<<<<<<<<<
 *     cdef double[:] cpos0 = numpy.ascontiguousarray(pos0.ravel(),dtype=numpy.float64)
 *     cdef double[:] cpos1 = numpy.ascontiguousarray(pos1.ravel(),dtype=numpy.float64)
 */
  __pyx_t_12 = __Pyx_PyInt_AsInt(__pyx_v_split); if (unlikely((__pyx_t_12 == (int)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 205; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_v_csplit = __pyx_t_12;
 206:     cdef double[:] cpos0 = numpy.ascontiguousarray(pos0.ravel(),dtype=numpy.float64)
  /* "histogram.pyx":206
 *         bin1 = 1
 *     cdef int csplit = split
 *     cdef double[:] cpos0 = numpy.ascontiguousarray(pos0.ravel(),dtype=numpy.float64)             # <<<<<<<<<<<<<<
 *     cdef double[:] cpos1 = numpy.ascontiguousarray(pos1.ravel(),dtype=numpy.float64)
 *     cdef double[:] data = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)
 */
  __pyx_t_1 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 206; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_3 = __Pyx_PyObject_GetAttrStr(__pyx_t_1, __pyx_n_s__ascontiguousarray); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 206; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_1 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos0), __pyx_n_s__ravel); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 206; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_2 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 206; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_1 = PyTuple_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 206; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__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 = 206; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_2));
  __pyx_t_9 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 206; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_t_13 = __Pyx_PyObject_GetAttrStr(__pyx_t_9, __pyx_n_s__float64); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 206; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  __Pyx_DECREF(__pyx_t_9); __pyx_t_9 = 0;
  if (PyDict_SetItem(__pyx_t_2, ((PyObject *)__pyx_n_s__dtype), __pyx_t_13) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 206; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0;
  __pyx_t_13 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_t_1), ((PyObject *)__pyx_t_2)); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 206; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0;
  __pyx_t_14 = __Pyx_PyObject_to_MemoryviewSlice_ds_double(__pyx_t_13);
  if (unlikely(!__pyx_t_14.memview)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 206; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0;
  __pyx_v_cpos0 = __pyx_t_14;
  __pyx_t_14.memview = NULL;
  __pyx_t_14.data = NULL;
 207:     cdef double[:] cpos1 = numpy.ascontiguousarray(pos1.ravel(),dtype=numpy.float64)
  /* "histogram.pyx":207
 *     cdef int csplit = split
 *     cdef double[:] cpos0 = numpy.ascontiguousarray(pos0.ravel(),dtype=numpy.float64)
 *     cdef double[:] cpos1 = numpy.ascontiguousarray(pos1.ravel(),dtype=numpy.float64)             # <<<<<<<<<<<<<<
 *     cdef double[:] data = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")
 */
  __pyx_t_13 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 207; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  __pyx_t_2 = __Pyx_PyObject_GetAttrStr(__pyx_t_13, __pyx_n_s__ascontiguousarray); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 207; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0;
  __pyx_t_13 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos1), __pyx_n_s__ravel); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 207; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  __pyx_t_1 = PyObject_Call(__pyx_t_13, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 207; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0;
  __pyx_t_13 = PyTuple_New(1); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 207; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  PyTuple_SET_ITEM(__pyx_t_13, 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 = 207; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_1));
  __pyx_t_3 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 207; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_9 = __Pyx_PyObject_GetAttrStr(__pyx_t_3, __pyx_n_s__float64); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 207; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  if (PyDict_SetItem(__pyx_t_1, ((PyObject *)__pyx_n_s__dtype), __pyx_t_9) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 207; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_9); __pyx_t_9 = 0;
  __pyx_t_9 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_t_13), ((PyObject *)__pyx_t_1)); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 207; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_13)); __pyx_t_13 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
  __pyx_t_15 = __Pyx_PyObject_to_MemoryviewSlice_ds_double(__pyx_t_9);
  if (unlikely(!__pyx_t_15.memview)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 207; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_9); __pyx_t_9 = 0;
  __pyx_v_cpos1 = __pyx_t_15;
  __pyx_t_15.memview = NULL;
  __pyx_t_15.data = NULL;
 208:     cdef double[:] data = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)
  /* "histogram.pyx":208
 *     cdef double[:] cpos0 = numpy.ascontiguousarray(pos0.ravel(),dtype=numpy.float64)
 *     cdef double[:] cpos1 = numpy.ascontiguousarray(pos1.ravel(),dtype=numpy.float64)
 *     cdef double[:] data = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64")
 */
  __pyx_t_9 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 208; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_t_1 = __Pyx_PyObject_GetAttrStr(__pyx_t_9, __pyx_n_s__ascontiguousarray); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 208; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_9); __pyx_t_9 = 0;
  __pyx_t_9 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_weights), __pyx_n_s__ravel); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 208; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_t_13 = PyObject_Call(__pyx_t_9, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 208; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  __Pyx_DECREF(__pyx_t_9); __pyx_t_9 = 0;
  __pyx_t_9 = PyTuple_New(1); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 208; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  PyTuple_SET_ITEM(__pyx_t_9, 0, __pyx_t_13);
  __Pyx_GIVEREF(__pyx_t_13);
  __pyx_t_13 = 0;
  __pyx_t_13 = PyDict_New(); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 208; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_13));
  __pyx_t_2 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 208; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __pyx_t_3 = __Pyx_PyObject_GetAttrStr(__pyx_t_2, __pyx_n_s__float64); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 208; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  if (PyDict_SetItem(__pyx_t_13, ((PyObject *)__pyx_n_s__dtype), __pyx_t_3) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 208; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_t_3 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_t_9), ((PyObject *)__pyx_t_13)); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 208; __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_9)); __pyx_t_9 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_13)); __pyx_t_13 = 0;
  __pyx_t_16 = __Pyx_PyObject_to_MemoryviewSlice_ds_double(__pyx_t_3);
  if (unlikely(!__pyx_t_16.memview)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 208; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_v_data = __pyx_t_16;
  __pyx_t_16.memview = NULL;
  __pyx_t_16.data = NULL;
 209:     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")
  /* "histogram.pyx":209
 *     cdef double[:] cpos1 = numpy.ascontiguousarray(pos1.ravel(),dtype=numpy.float64)
 *     cdef double[:] data = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outMerge = numpy.zeros((bin0, bin1), dtype="float64")
 */
  __pyx_t_3 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 209; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_13 = __Pyx_PyObject_GetAttrStr(__pyx_t_3, __pyx_n_s__zeros); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 209; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  __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 = 209; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_9 = PyInt_FromLong(__pyx_v_bin1); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 209; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_t_1 = PyTuple_New(2); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 209; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__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_9);
  __Pyx_GIVEREF(__pyx_t_9);
  __pyx_t_3 = 0;
  __pyx_t_9 = 0;
  __pyx_t_9 = PyTuple_New(1); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 209; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  PyTuple_SET_ITEM(__pyx_t_9, 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 = 209; __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 = 209; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_3 = PyObject_Call(__pyx_t_13, ((PyObject *)__pyx_t_9), ((PyObject *)__pyx_t_1)); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 209; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_9)); __pyx_t_9 = 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 = 209; __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_pybuffernd_outData.rcbuffer->pybuffer, (PyObject*)__pyx_t_17, &__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 2, 0, __pyx_stack) == -1)) {
      __pyx_v_outData = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 209; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_pybuffernd_outData.diminfo[0].strides = __pyx_pybuffernd_outData.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_outData.diminfo[0].shape = __pyx_pybuffernd_outData.rcbuffer->pybuffer.shape[0]; __pyx_pybuffernd_outData.diminfo[1].strides = __pyx_pybuffernd_outData.rcbuffer->pybuffer.strides[1]; __pyx_pybuffernd_outData.diminfo[1].shape = __pyx_pybuffernd_outData.rcbuffer->pybuffer.shape[1];
    }
  }
  __pyx_t_17 = 0;
  __pyx_v_outData = ((PyArrayObject *)__pyx_t_3);
  __pyx_t_3 = 0;
 210:     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64")
  /* "histogram.pyx":210
 *     cdef double[:] data = numpy.ascontiguousarray(weights.ravel(),dtype=numpy.float64)
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64")             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outMerge = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] edges0, edges1
 */
  __pyx_t_3 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 210; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_1 = __Pyx_PyObject_GetAttrStr(__pyx_t_3, __pyx_n_s__zeros); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 210; __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 = 210; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_9 = PyInt_FromLong(__pyx_v_bin1); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 210; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_t_13 = PyTuple_New(2); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 210; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  PyTuple_SET_ITEM(__pyx_t_13, 0, __pyx_t_3);
  __Pyx_GIVEREF(__pyx_t_3);
  PyTuple_SET_ITEM(__pyx_t_13, 1, __pyx_t_9);
  __Pyx_GIVEREF(__pyx_t_9);
  __pyx_t_3 = 0;
  __pyx_t_9 = 0;
  __pyx_t_9 = PyTuple_New(1); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 210; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  PyTuple_SET_ITEM(__pyx_t_9, 0, ((PyObject *)__pyx_t_13));
  __Pyx_GIVEREF(((PyObject *)__pyx_t_13));
  __pyx_t_13 = 0;
  __pyx_t_13 = PyDict_New(); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 210; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_13));
  if (PyDict_SetItem(__pyx_t_13, ((PyObject *)__pyx_n_s__dtype), ((PyObject *)__pyx_n_s__float64)) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 210; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_3 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_t_9), ((PyObject *)__pyx_t_13)); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 210; __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_9)); __pyx_t_9 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_13)); __pyx_t_13 = 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 = 210; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_18 = ((PyArrayObject *)__pyx_t_3);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_outCount.rcbuffer->pybuffer, (PyObject*)__pyx_t_18, &__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 2, 0, __pyx_stack) == -1)) {
      __pyx_v_outCount = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 210; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_pybuffernd_outCount.diminfo[0].strides = __pyx_pybuffernd_outCount.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_outCount.diminfo[0].shape = __pyx_pybuffernd_outCount.rcbuffer->pybuffer.shape[0]; __pyx_pybuffernd_outCount.diminfo[1].strides = __pyx_pybuffernd_outCount.rcbuffer->pybuffer.strides[1]; __pyx_pybuffernd_outCount.diminfo[1].shape = __pyx_pybuffernd_outCount.rcbuffer->pybuffer.shape[1];
    }
  }
  __pyx_t_18 = 0;
  __pyx_v_outCount = ((PyArrayObject *)__pyx_t_3);
  __pyx_t_3 = 0;
 211:     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outMerge = numpy.zeros((bin0, bin1), dtype="float64")
  /* "histogram.pyx":211
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outData = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outCount = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outMerge = numpy.zeros((bin0, bin1), dtype="float64")             # <<<<<<<<<<<<<<
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] edges0, edges1
 *     cdef double min0 = pos0.min()
 */
  __pyx_t_3 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 211; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_13 = __Pyx_PyObject_GetAttrStr(__pyx_t_3, __pyx_n_s__zeros); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 211; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  __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 = 211; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_9 = PyInt_FromLong(__pyx_v_bin1); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 211; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_t_1 = PyTuple_New(2); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 211; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__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_9);
  __Pyx_GIVEREF(__pyx_t_9);
  __pyx_t_3 = 0;
  __pyx_t_9 = 0;
  __pyx_t_9 = PyTuple_New(1); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 211; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  PyTuple_SET_ITEM(__pyx_t_9, 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 = 211; __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 = 211; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_3 = PyObject_Call(__pyx_t_13, ((PyObject *)__pyx_t_9), ((PyObject *)__pyx_t_1)); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 211; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_9)); __pyx_t_9 = 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 = 211; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_19 = ((PyArrayObject *)__pyx_t_3);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_outMerge.rcbuffer->pybuffer, (PyObject*)__pyx_t_19, &__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 2, 0, __pyx_stack) == -1)) {
      __pyx_v_outMerge = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.buf = NULL;
      {__pyx_filename = __pyx_f[0]; __pyx_lineno = 211; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    } else {__pyx_pybuffernd_outMerge.diminfo[0].strides = __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_outMerge.diminfo[0].shape = __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.shape[0]; __pyx_pybuffernd_outMerge.diminfo[1].strides = __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.strides[1]; __pyx_pybuffernd_outMerge.diminfo[1].shape = __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.shape[1];
    }
  }
  __pyx_t_19 = 0;
  __pyx_v_outMerge = ((PyArrayObject *)__pyx_t_3);
  __pyx_t_3 = 0;
 212:     cdef numpy.ndarray[numpy.float64_t, ndim = 1] edges0, edges1
 213:     cdef double min0 = pos0.min()
  /* "histogram.pyx":213
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 2] outMerge = numpy.zeros((bin0, bin1), dtype="float64")
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] edges0, edges1
 *     cdef double min0 = pos0.min()             # <<<<<<<<<<<<<<
 *     cdef double max0 = pos0.max() * EPS32
 *     cdef double min1 = pos1.min()
 */
  __pyx_t_3 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos0), __pyx_n_s__min); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 213; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __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 = 213; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 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 = 213; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_v_min0 = __pyx_t_20;
 214:     cdef double max0 = pos0.max() * EPS32
  /* "histogram.pyx":214
 *     cdef numpy.ndarray[numpy.float64_t, ndim = 1] edges0, edges1
 *     cdef double min0 = pos0.min()
 *     cdef double max0 = pos0.max() * EPS32             # <<<<<<<<<<<<<<
 *     cdef double min1 = pos1.min()
 *     cdef double max1 = pos1.max() * EPS32
 */
  __pyx_t_1 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos0), __pyx_n_s__max); 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_3 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_empty_tuple), NULL); 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_t_1 = __Pyx_GetModuleGlobalName(__pyx_n_s__EPS32); 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_9 = PyNumber_Multiply(__pyx_t_3, __pyx_t_1); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 214; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_20 = __pyx_PyFloat_AsDouble(__pyx_t_9); if (unlikely((__pyx_t_20 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 214; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_9); __pyx_t_9 = 0;
  __pyx_v_max0 = __pyx_t_20;
 215:     cdef double min1 = pos1.min()
  /* "histogram.pyx":215
 *     cdef double min0 = pos0.min()
 *     cdef double max0 = pos0.max() * EPS32
 *     cdef double min1 = pos1.min()             # <<<<<<<<<<<<<<
 *     cdef double max1 = pos1.max() * EPS32
 *     cdef double idp0 = (< double > bin0) / (max0 - min0)
 */
  __pyx_t_9 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos1), __pyx_n_s__min); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 215; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_t_1 = PyObject_Call(__pyx_t_9, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 215; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_9); __pyx_t_9 = 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 = 215; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_v_min1 = __pyx_t_20;
 216:     cdef double max1 = pos1.max() * EPS32
  /* "histogram.pyx":216
 *     cdef double max0 = pos0.max() * EPS32
 *     cdef double min1 = pos1.min()
 *     cdef double max1 = pos1.max() * EPS32             # <<<<<<<<<<<<<<
 *     cdef double idp0 = (< double > bin0) / (max0 - min0)
 *     cdef double idp1 = (< double > bin1) / (max1 - min1)
 */
  __pyx_t_1 = __Pyx_PyObject_GetAttrStr(((PyObject *)__pyx_v_pos1), __pyx_n_s__max); 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_t_9 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 216; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_1 = __Pyx_GetModuleGlobalName(__pyx_n_s__EPS32); 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_t_3 = PyNumber_Multiply(__pyx_t_9, __pyx_t_1); 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_DECREF(__pyx_t_9); __pyx_t_9 = 0;
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_20 = __pyx_PyFloat_AsDouble(__pyx_t_3); if (unlikely((__pyx_t_20 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 216; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_v_max1 = __pyx_t_20;
 217:     cdef double idp0 = (< double > bin0) / (max0 - min0)
  /* "histogram.pyx":217
 *     cdef double min1 = pos1.min()
 *     cdef double max1 = pos1.max() * EPS32
 *     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));
 218:     cdef double idp1 = (< double > bin1) / (max1 - min1)
  /* "histogram.pyx":218
 *     cdef double max1 = pos1.max() * EPS32
 *     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));
 219:     cdef double dbin0 = 0.5, dbin1 = 0.5
  /* "histogram.pyx":219
 *     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
 *     cdef double epsilon = 1e-10
 */
  __pyx_v_dbin0 = 0.5;
  __pyx_v_dbin1 = 0.5;
 220:     cdef double fbin0, fbin1, p0, p1, d, rest, delta0l, delta0r, delta1l, delta1r, aera
 221:     cdef double epsilon = 1e-10
  /* "histogram.pyx":221
 *     cdef double dbin0 = 0.5, dbin1 = 0.5
 *     cdef double fbin0, fbin1, p0, p1, d, rest, delta0l, delta0r, delta1l, delta1r, aera
 *     cdef double epsilon = 1e-10             # <<<<<<<<<<<<<<
 * 
 *     edges0 = numpy.linspace(min0+(0.5/idp0),max0-(0.5/idp0),bin0)
 */
  __pyx_v_epsilon = 1e-10;
 222: 
 223:     edges0 = numpy.linspace(min0+(0.5/idp0),max0-(0.5/idp0),bin0)
  /* "histogram.pyx":223
 *     cdef double epsilon = 1e-10
 * 
 *     edges0 = numpy.linspace(min0+(0.5/idp0),max0-(0.5/idp0),bin0)             # <<<<<<<<<<<<<<
 *     edges1 = numpy.linspace(min1+(0.5/idp1),max1-(0.5/idp1),bin1)
 *     if nthread is not None:
 */
  __pyx_t_3 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 223; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_1 = __Pyx_PyObject_GetAttrStr(__pyx_t_3, __pyx_n_s__linspace); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 223; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_t_3 = PyFloat_FromDouble((__pyx_v_min0 + (0.5 / __pyx_v_idp0))); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 223; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __pyx_t_9 = PyFloat_FromDouble((__pyx_v_max0 - (0.5 / __pyx_v_idp0))); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 223; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_t_13 = PyInt_FromLong(__pyx_v_bin0); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 223; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  __pyx_t_2 = PyTuple_New(3); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 223; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  PyTuple_SET_ITEM(__pyx_t_2, 0, __pyx_t_3);
  __Pyx_GIVEREF(__pyx_t_3);
  PyTuple_SET_ITEM(__pyx_t_2, 1, __pyx_t_9);
  __Pyx_GIVEREF(__pyx_t_9);
  PyTuple_SET_ITEM(__pyx_t_2, 2, __pyx_t_13);
  __Pyx_GIVEREF(__pyx_t_13);
  __pyx_t_3 = 0;
  __pyx_t_9 = 0;
  __pyx_t_13 = 0;
  __pyx_t_13 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_t_2), NULL); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 223; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0;
  if (!(likely(((__pyx_t_13) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_13, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 223; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_21 = ((PyArrayObject *)__pyx_t_13);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_edges0.rcbuffer->pybuffer);
    __pyx_t_12 = __Pyx_GetBufferAndValidate(&__pyx_pybuffernd_edges0.rcbuffer->pybuffer, (PyObject*)__pyx_t_21, &__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack);
    if (unlikely(__pyx_t_12 < 0)) {
      PyErr_Fetch(&__pyx_t_8, &__pyx_t_7, &__pyx_t_6);
      if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_edges0.rcbuffer->pybuffer, (PyObject*)__pyx_v_edges0, &__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {
        Py_XDECREF(__pyx_t_8); Py_XDECREF(__pyx_t_7); Py_XDECREF(__pyx_t_6);
        __Pyx_RaiseBufferFallbackError();
      } else {
        PyErr_Restore(__pyx_t_8, __pyx_t_7, __pyx_t_6);
      }
    }
    __pyx_pybuffernd_edges0.diminfo[0].strides = __pyx_pybuffernd_edges0.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_edges0.diminfo[0].shape = __pyx_pybuffernd_edges0.rcbuffer->pybuffer.shape[0];
    if (unlikely(__pyx_t_12 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 223; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }
  __pyx_t_21 = 0;
  __pyx_v_edges0 = ((PyArrayObject *)__pyx_t_13);
  __pyx_t_13 = 0;
 224:     edges1 = numpy.linspace(min1+(0.5/idp1),max1-(0.5/idp1),bin1)
  /* "histogram.pyx":224
 * 
 *     edges0 = numpy.linspace(min0+(0.5/idp0),max0-(0.5/idp0),bin0)
 *     edges1 = numpy.linspace(min1+(0.5/idp1),max1-(0.5/idp1),bin1)             # <<<<<<<<<<<<<<
 *     if nthread is not None:
 *         if isinstance(nthread, int) and (nthread > 0):
 */
  __pyx_t_13 = __Pyx_GetModuleGlobalName(__pyx_n_s__numpy); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 224; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  __pyx_t_2 = __Pyx_PyObject_GetAttrStr(__pyx_t_13, __pyx_n_s__linspace); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 224; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0;
  __pyx_t_13 = PyFloat_FromDouble((__pyx_v_min1 + (0.5 / __pyx_v_idp1))); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 224; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_13);
  __pyx_t_1 = PyFloat_FromDouble((__pyx_v_max1 - (0.5 / __pyx_v_idp1))); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 224; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_9 = PyInt_FromLong(__pyx_v_bin1); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 224; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __pyx_t_3 = PyTuple_New(3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 224; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_t_13);
  __Pyx_GIVEREF(__pyx_t_13);
  PyTuple_SET_ITEM(__pyx_t_3, 1, __pyx_t_1);
  __Pyx_GIVEREF(__pyx_t_1);
  PyTuple_SET_ITEM(__pyx_t_3, 2, __pyx_t_9);
  __Pyx_GIVEREF(__pyx_t_9);
  __pyx_t_13 = 0;
  __pyx_t_1 = 0;
  __pyx_t_9 = 0;
  __pyx_t_9 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_t_3), NULL); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 224; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0;
  if (!(likely(((__pyx_t_9) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_9, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 224; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_21 = ((PyArrayObject *)__pyx_t_9);
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_edges1.rcbuffer->pybuffer);
    __pyx_t_12 = __Pyx_GetBufferAndValidate(&__pyx_pybuffernd_edges1.rcbuffer->pybuffer, (PyObject*)__pyx_t_21, &__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack);
    if (unlikely(__pyx_t_12 < 0)) {
      PyErr_Fetch(&__pyx_t_6, &__pyx_t_7, &__pyx_t_8);
      if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_edges1.rcbuffer->pybuffer, (PyObject*)__pyx_v_edges1, &__Pyx_TypeInfo_nn___pyx_t_5numpy_float64_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {
        Py_XDECREF(__pyx_t_6); Py_XDECREF(__pyx_t_7); Py_XDECREF(__pyx_t_8);
        __Pyx_RaiseBufferFallbackError();
      } else {
        PyErr_Restore(__pyx_t_6, __pyx_t_7, __pyx_t_8);
      }
    }
    __pyx_pybuffernd_edges1.diminfo[0].strides = __pyx_pybuffernd_edges1.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_edges1.diminfo[0].shape = __pyx_pybuffernd_edges1.rcbuffer->pybuffer.shape[0];
    if (unlikely(__pyx_t_12 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 224; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }
  __pyx_t_21 = 0;
  __pyx_v_edges1 = ((PyArrayObject *)__pyx_t_9);
  __pyx_t_9 = 0;
 225:     if nthread is not None:
  /* "histogram.pyx":225
 *     edges0 = numpy.linspace(min0+(0.5/idp0),max0-(0.5/idp0),bin0)
 *     edges1 = numpy.linspace(min1+(0.5/idp1),max1-(0.5/idp1),bin1)
 *     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) {
 226:         if isinstance(nthread, int) and (nthread > 0):
    /* "histogram.pyx":226
 *     edges1 = numpy.linspace(min1+(0.5/idp1),max1-(0.5/idp1),bin1)
 *     if nthread is not None:
 *         if isinstance(nthread, int) and (nthread > 0):             # <<<<<<<<<<<<<<
 *             omp_set_num_threads(< int > nthread)
 *     with nogil:
 */
    __pyx_t_4 = PyInt_Check(__pyx_v_nthread); 
    if (__pyx_t_4) {
      __pyx_t_9 = PyObject_RichCompare(__pyx_v_nthread, __pyx_int_0, Py_GT); __Pyx_XGOTREF(__pyx_t_9); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 226; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __pyx_t_22 = __Pyx_PyObject_IsTrue(__pyx_t_9); if (unlikely(__pyx_t_22 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 226; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_9); __pyx_t_9 = 0;
      __pyx_t_23 = __pyx_t_22;
    } else {
      __pyx_t_23 = __pyx_t_4;
    }
    if (__pyx_t_23) {
 227:             omp_set_num_threads(< int > nthread)
      /* "histogram.pyx":227
 *     if nthread is not None:
 *         if isinstance(nthread, int) and (nthread > 0):
 *             omp_set_num_threads(< int > nthread)             # <<<<<<<<<<<<<<
 *     with nogil:
 *         for i in range(size):
 */
      __pyx_t_12 = __Pyx_PyInt_AsInt(__pyx_v_nthread); if (unlikely((__pyx_t_12 == (int)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 227; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      omp_set_num_threads(((int)__pyx_t_12));
      goto __pyx_L18;
    }
    __pyx_L18:;
    goto __pyx_L17;
  }
  __pyx_L17:;
 228:     with nogil:
  /* "histogram.pyx":228
 *         if isinstance(nthread, int) and (nthread > 0):
 *             omp_set_num_threads(< int > nthread)
 *     with nogil:             # <<<<<<<<<<<<<<
 *         for i in range(size):
 *             p0 = cpos0[i]
 */
  {
      #ifdef WITH_THREAD
      PyThreadState *_save;
      Py_UNBLOCK_THREADS
      #endif
      /*try:*/ {

      /* "histogram.pyx":228
 *         if isinstance(nthread, int) and (nthread > 0):
 *             omp_set_num_threads(< int > nthread)
 *     with nogil:             # <<<<<<<<<<<<<<
 *         for i in range(size):
 *             p0 = cpos0[i]
 */
      /*finally:*/ {
        #ifdef WITH_THREAD
        Py_BLOCK_THREADS
        #endif
      }
  }
 229:         for i in range(size):
        /* "histogram.pyx":229
 *             omp_set_num_threads(< int > nthread)
 *     with nogil:
 *         for i in range(size):             # <<<<<<<<<<<<<<
 *             p0 = cpos0[i]
 *             p1 = cpos1[i]
 */
        __pyx_t_11 = __pyx_v_size;
        for (__pyx_t_5 = 0; __pyx_t_5 < __pyx_t_11; __pyx_t_5+=1) {
          __pyx_v_i = __pyx_t_5;
 230:             p0 = cpos0[i]
          /* "histogram.pyx":230
 *     with nogil:
 *         for i in range(size):
 *             p0 = cpos0[i]             # <<<<<<<<<<<<<<
 *             p1 = cpos1[i]
 *             d = data[i]
 */
          __pyx_t_24 = __pyx_v_i;
          __pyx_v_p0 = (*((double *) ( /* dim=0 */ (__pyx_v_cpos0.data + __pyx_t_24 * __pyx_v_cpos0.strides[0]) )));
 231:             p1 = cpos1[i]
          /* "histogram.pyx":231
 *         for i in range(size):
 *             p0 = cpos0[i]
 *             p1 = cpos1[i]             # <<<<<<<<<<<<<<
 *             d = data[i]
 *             fbin0 = (p0 - min0) * idp0
 */
          __pyx_t_25 = __pyx_v_i;
          __pyx_v_p1 = (*((double *) ( /* dim=0 */ (__pyx_v_cpos1.data + __pyx_t_25 * __pyx_v_cpos1.strides[0]) )));
 232:             d = data[i]
          /* "histogram.pyx":232
 *             p0 = cpos0[i]
 *             p1 = cpos1[i]
 *             d = data[i]             # <<<<<<<<<<<<<<
 *             fbin0 = (p0 - min0) * idp0
 *             fbin1 = (p1 - min1) * idp1
 */
          __pyx_t_26 = __pyx_v_i;
          __pyx_v_d = (*((double *) ( /* dim=0 */ (__pyx_v_data.data + __pyx_t_26 * __pyx_v_data.strides[0]) )));
 233:             fbin0 = (p0 - min0) * idp0
          /* "histogram.pyx":233
 *             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);
 234:             fbin1 = (p1 - min1) * idp1
          /* "histogram.pyx":234
 *             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);
 235:             b0 = < long > floor(fbin0)
          /* "histogram.pyx":235
 *             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));
 236:             b1 = < long > floor(fbin1)
          /* "histogram.pyx":236
 *             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));
 237:             if b0 == bin0:
          /* "histogram.pyx":237
 *             b0 = < long > floor(fbin0)
 *             b1 = < long > floor(fbin1)
 *             if b0 == bin0:             # <<<<<<<<<<<<<<
 *                 b0 = bin0 - 1
 *                 fbin0 = (< double > bin0) - 0.5
 */
          __pyx_t_23 = (__pyx_v_b0 == __pyx_v_bin0);
          if (__pyx_t_23) {
 238:                 b0 = bin0 - 1
            /* "histogram.pyx":238
 *             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);
 239:                 fbin0 = (< double > bin0) - 0.5
            /* "histogram.pyx":239
 *             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_L24;
          }
 240:             elif b0 == 0:
          /* "histogram.pyx":240
 *                 b0 = bin0 - 1
 *                 fbin0 = (< double > bin0) - 0.5
 *             elif b0 == 0:             # <<<<<<<<<<<<<<
 *                 fbin0 = 0.5
 *             if b1 == bin1:
 */
          __pyx_t_23 = (__pyx_v_b0 == 0);
          if (__pyx_t_23) {
 241:                 fbin0 = 0.5
            /* "histogram.pyx":241
 *                 fbin0 = (< double > bin0) - 0.5
 *             elif b0 == 0:
 *                 fbin0 = 0.5             # <<<<<<<<<<<<<<
 *             if b1 == bin1:
 *                 b1 = bin1 - 1
 */
            __pyx_v_fbin0 = 0.5;
            goto __pyx_L24;
          }
          __pyx_L24:;
 242:             if b1 == bin1:
          /* "histogram.pyx":242
 *             elif b0 == 0:
 *                 fbin0 = 0.5
 *             if b1 == bin1:             # <<<<<<<<<<<<<<
 *                 b1 = bin1 - 1
 *                 fbin1 = (< double > bin1) - 0.5
 */
          __pyx_t_23 = (__pyx_v_b1 == __pyx_v_bin1);
          if (__pyx_t_23) {
 243:                 b1 = bin1 - 1
            /* "histogram.pyx":243
 *                 fbin0 = 0.5
 *             if b1 == bin1:
 *                 b1 = bin1 - 1             # <<<<<<<<<<<<<<
 *                 fbin1 = (< double > bin1) - 0.5
 *             elif b1 == 0:
 */
            __pyx_v_b1 = (__pyx_v_bin1 - 1);
 244:                 fbin1 = (< double > bin1) - 0.5
            /* "histogram.pyx":244
 *             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_L25;
          }
 245:             elif b1 == 0:
          /* "histogram.pyx":245
 *                 b1 = bin1 - 1
 *                 fbin1 = (< double > bin1) - 0.5
 *             elif b1 == 0:             # <<<<<<<<<<<<<<
 *                 fbin1 = 0.5
 * 
 */
          __pyx_t_23 = (__pyx_v_b1 == 0);
          if (__pyx_t_23) {
 246:                 fbin1 = 0.5
            /* "histogram.pyx":246
 *                 fbin1 = (< double > bin1) - 0.5
 *             elif b1 == 0:
 *                 fbin1 = 0.5             # <<<<<<<<<<<<<<
 * 
 *             delta0l = fbin0 -< double > b0 - dbin0
 */
            __pyx_v_fbin1 = 0.5;
            goto __pyx_L25;
          }
          __pyx_L25:;
 247: 
 248:             delta0l = fbin0 -< double > b0 - dbin0
          /* "histogram.pyx":248
 *                 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);
 249:             delta0r = fbin0 -< double > b0 - 1 + dbin0
          /* "histogram.pyx":249
 * 
 *             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);
 250:             delta1l = fbin1 -< double > b1 - dbin1
          /* "histogram.pyx":250
 *             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);
 251:             delta1r = fbin1 -< double > b1 - 1 + dbin1
          /* "histogram.pyx":251
 *             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);
 252:             rest = 1.0
          /* "histogram.pyx":252
 *             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;
 253:             if csplit == 1:
          /* "histogram.pyx":253
 *             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_23 = (__pyx_v_csplit == 1);
          if (__pyx_t_23) {
 254:                 if delta0l < 0 and b0 > 0:
            /* "histogram.pyx":254
 *             rest = 1.0
 *             if csplit == 1:
 *                 if delta0l < 0 and b0 > 0:             # <<<<<<<<<<<<<<
 *                     if delta1l < 0 and b1 > 0:
 *                         area = delta0l * delta1l
 */
            __pyx_t_23 = (__pyx_v_delta0l < 0.0);
            if (__pyx_t_23) {
              __pyx_t_4 = (__pyx_v_b0 > 0);
              __pyx_t_22 = __pyx_t_4;
            } else {
              __pyx_t_22 = __pyx_t_23;
            }
            if (__pyx_t_22) {
 255:                     if delta1l < 0 and b1 > 0:
              /* "histogram.pyx":255
 *             if csplit == 1:
 *                 if delta0l < 0 and b0 > 0:
 *                     if delta1l < 0 and b1 > 0:             # <<<<<<<<<<<<<<
 *                         area = delta0l * delta1l
 *                         rest -= area
 */
              __pyx_t_22 = (__pyx_v_delta1l < 0.0);
              if (__pyx_t_22) {
                __pyx_t_23 = (__pyx_v_b1 > 0);
                __pyx_t_4 = __pyx_t_23;
              } else {
                __pyx_t_4 = __pyx_t_22;
              }
              if (__pyx_t_4) {
 256:                         area = delta0l * delta1l
                /* "histogram.pyx":256
 *                 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);
 257:                         rest -= area
                /* "histogram.pyx":257
 *                     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);
 258:                         outCount[b0 - 1, b1 - 1] += area
                /* "histogram.pyx":258
 *                         area = delta0l * delta1l
 *                         rest -= area
 *                         outCount[b0 - 1, b1 - 1] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 - 1, b1 - 1] += area * d
 * 
 */
                __pyx_t_27 = (__pyx_v_b0 - 1);
                __pyx_t_28 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_27, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_28, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 259:                         outData[b0 - 1, b1 - 1] += area * d
                /* "histogram.pyx":259
 *                         rest -= area
 *                         outCount[b0 - 1, b1 - 1] += area
 *                         outData[b0 - 1, b1 - 1] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (-delta0l) * (1 + delta1l)
 */
                __pyx_t_29 = (__pyx_v_b0 - 1);
                __pyx_t_30 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_29, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_30, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
 260: 
 261:                         area = (-delta0l) * (1 + delta1l)
                /* "histogram.pyx":261
 *                         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));
 262:                         rest -= area
                /* "histogram.pyx":262
 * 
 *                         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);
 263:                         outCount[b0 - 1, b1 ] += area
                /* "histogram.pyx":263
 *                         area = (-delta0l) * (1 + delta1l)
 *                         rest -= area
 *                         outCount[b0 - 1, b1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 - 1, b1 ] += area * d
 * 
 */
                __pyx_t_31 = (__pyx_v_b0 - 1);
                __pyx_t_32 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_31, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_32, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 264:                         outData[b0 - 1, b1 ] += area * d
                /* "histogram.pyx":264
 *                         rest -= area
 *                         outCount[b0 - 1, b1 ] += area
 *                         outData[b0 - 1, b1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (1 + delta0l) * (-delta1l)
 */
                __pyx_t_33 = (__pyx_v_b0 - 1);
                __pyx_t_34 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_33, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_34, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
 265: 
 266:                         area = (1 + delta0l) * (-delta1l)
                /* "histogram.pyx":266
 *                         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));
 267:                         rest -= area
                /* "histogram.pyx":267
 * 
 *                         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);
 268:                         outCount[b0 , b1 - 1 ] += area
                /* "histogram.pyx":268
 *                         area = (1 + delta0l) * (-delta1l)
 *                         rest -= area
 *                         outCount[b0 , b1 - 1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 , b1 - 1 ] += area * d
 * 
 */
                __pyx_t_35 = __pyx_v_b0;
                __pyx_t_36 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_35, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_36, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 269:                         outData[b0 , b1 - 1 ] += area * d
                /* "histogram.pyx":269
 *                         rest -= area
 *                         outCount[b0 , b1 - 1 ] += area
 *                         outData[b0 , b1 - 1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                     elif delta1r > 0 and b1 < bin1 - 1:
 */
                __pyx_t_37 = __pyx_v_b0;
                __pyx_t_38 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_37, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_38, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
                goto __pyx_L28;
              }
 270: 
 271:                     elif delta1r > 0 and b1 < bin1 - 1:
              /* "histogram.pyx":271
 *                         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_22 = (__pyx_v_b1 < (__pyx_v_bin1 - 1));
                __pyx_t_23 = __pyx_t_22;
              } else {
                __pyx_t_23 = __pyx_t_4;
              }
              if (__pyx_t_23) {
 272:                         area = -delta0l * delta1r
                /* "histogram.pyx":272
 * 
 *                     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);
 273:                         rest -= area
                /* "histogram.pyx":273
 *                     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);
 274:                         outCount[b0 - 1, b1 + 1] += area
                /* "histogram.pyx":274
 *                         area = -delta0l * delta1r
 *                         rest -= area
 *                         outCount[b0 - 1, b1 + 1] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 - 1, b1 + 1] += area * d
 * 
 */
                __pyx_t_39 = (__pyx_v_b0 - 1);
                __pyx_t_40 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_39, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_40, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 275:                         outData[b0 - 1, b1 + 1] += area * d
                /* "histogram.pyx":275
 *                         rest -= area
 *                         outCount[b0 - 1, b1 + 1] += area
 *                         outData[b0 - 1, b1 + 1] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (-delta0l) * (1 - delta1r)
 */
                __pyx_t_41 = (__pyx_v_b0 - 1);
                __pyx_t_42 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_41, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_42, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
 276: 
 277:                         area = (-delta0l) * (1 - delta1r)
                /* "histogram.pyx":277
 *                         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));
 278:                         rest -= area
                /* "histogram.pyx":278
 * 
 *                         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);
 279:                         outCount[b0 - 1, b1 ] += area
                /* "histogram.pyx":279
 *                         area = (-delta0l) * (1 - delta1r)
 *                         rest -= area
 *                         outCount[b0 - 1, b1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 - 1, b1 ] += area * d
 * 
 */
                __pyx_t_43 = (__pyx_v_b0 - 1);
                __pyx_t_44 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_43, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_44, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 280:                         outData[b0 - 1, b1 ] += area * d
                /* "histogram.pyx":280
 *                         rest -= area
 *                         outCount[b0 - 1, b1 ] += area
 *                         outData[b0 - 1, b1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (1 + delta0l) * (delta1r)
 */
                __pyx_t_45 = (__pyx_v_b0 - 1);
                __pyx_t_46 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_45, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_46, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
 281: 
 282:                         area = (1 + delta0l) * (delta1r)
                /* "histogram.pyx":282
 *                         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);
 283:                         rest -= area
                /* "histogram.pyx":283
 * 
 *                         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);
 284:                         outCount[b0 , b1 + 1 ] += area
                /* "histogram.pyx":284
 *                         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_47 = __pyx_v_b0;
                __pyx_t_48 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_47, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_48, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 285:                         outData[b0 , b1 + 1 ] += area * d
                /* "histogram.pyx":285
 *                         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_49 = __pyx_v_b0;
                __pyx_t_50 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_49, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_50, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
                goto __pyx_L28;
              }
              __pyx_L28:;
              goto __pyx_L27;
            }
 286:                 elif delta0r > 0 and b0 < bin0 - 1:
            /* "histogram.pyx":286
 *                         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_23 = (__pyx_v_delta0r > 0.0);
            if (__pyx_t_23) {
              __pyx_t_4 = (__pyx_v_b0 < (__pyx_v_bin0 - 1));
              __pyx_t_22 = __pyx_t_4;
            } else {
              __pyx_t_22 = __pyx_t_23;
            }
            if (__pyx_t_22) {
 287:                     if delta1l < 0 and b1 > 0:
              /* "histogram.pyx":287
 *                         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_22 = (__pyx_v_delta1l < 0.0);
              if (__pyx_t_22) {
                __pyx_t_23 = (__pyx_v_b1 > 0);
                __pyx_t_4 = __pyx_t_23;
              } else {
                __pyx_t_4 = __pyx_t_22;
              }
              if (__pyx_t_4) {
 288:                         area = -delta0r * delta1l
                /* "histogram.pyx":288
 *                 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);
 289:                         rest -= area
                /* "histogram.pyx":289
 *                     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);
 290:                         outCount[b0 + 1, b1 - 1] += area
                /* "histogram.pyx":290
 *                         area = -delta0r * delta1l
 *                         rest -= area
 *                         outCount[b0 + 1, b1 - 1] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 + 1, b1 - 1] += area * d
 * 
 */
                __pyx_t_51 = (__pyx_v_b0 + 1);
                __pyx_t_52 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_51, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_52, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 291:                         outData[b0 + 1, b1 - 1] += area * d
                /* "histogram.pyx":291
 *                         rest -= area
 *                         outCount[b0 + 1, b1 - 1] += area
 *                         outData[b0 + 1, b1 - 1] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (delta0r) * (1 + delta1l)
 */
                __pyx_t_53 = (__pyx_v_b0 + 1);
                __pyx_t_54 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_53, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_54, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
 292: 
 293:                         area = (delta0r) * (1 + delta1l)
                /* "histogram.pyx":293
 *                         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));
 294:                         rest -= area
                /* "histogram.pyx":294
 * 
 *                         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);
 295:                         outCount[b0 + 1, b1 ] += area
                /* "histogram.pyx":295
 *                         area = (delta0r) * (1 + delta1l)
 *                         rest -= area
 *                         outCount[b0 + 1, b1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 + 1, b1 ] += area * d
 * 
 */
                __pyx_t_55 = (__pyx_v_b0 + 1);
                __pyx_t_56 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_55, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_56, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 296:                         outData[b0 + 1, b1 ] += area * d
                /* "histogram.pyx":296
 *                         rest -= area
 *                         outCount[b0 + 1, b1 ] += area
 *                         outData[b0 + 1, b1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (1 - delta0r) * (-delta1l)
 */
                __pyx_t_57 = (__pyx_v_b0 + 1);
                __pyx_t_58 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_57, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_58, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
 297: 
 298:                         area = (1 - delta0r) * (-delta1l)
                /* "histogram.pyx":298
 *                         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));
 299:                         rest -= area
                /* "histogram.pyx":299
 * 
 *                         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);
 300:                         outCount[b0 , b1 - 1 ] += area
                /* "histogram.pyx":300
 *                         area = (1 - delta0r) * (-delta1l)
 *                         rest -= area
 *                         outCount[b0 , b1 - 1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 , b1 - 1 ] += area * d
 * 
 */
                __pyx_t_59 = __pyx_v_b0;
                __pyx_t_60 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_59, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_60, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 301:                         outData[b0 , b1 - 1 ] += area * d
                /* "histogram.pyx":301
 *                         rest -= area
 *                         outCount[b0 , b1 - 1 ] += area
 *                         outData[b0 , b1 - 1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                     elif delta1r > 0 and b1 < bin1 - 1:
 */
                __pyx_t_61 = __pyx_v_b0;
                __pyx_t_62 = (__pyx_v_b1 - 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_61, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_62, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
                goto __pyx_L29;
              }
 302: 
 303:                     elif delta1r > 0 and b1 < bin1 - 1:
              /* "histogram.pyx":303
 *                         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_22 = (__pyx_v_b1 < (__pyx_v_bin1 - 1));
                __pyx_t_23 = __pyx_t_22;
              } else {
                __pyx_t_23 = __pyx_t_4;
              }
              if (__pyx_t_23) {
 304:                         area = delta0r * delta1r
                /* "histogram.pyx":304
 * 
 *                     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);
 305:                         rest -= area
                /* "histogram.pyx":305
 *                     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);
 306:                         outCount[b0 + 1, b1 + 1] += area
                /* "histogram.pyx":306
 *                         area = delta0r * delta1r
 *                         rest -= area
 *                         outCount[b0 + 1, b1 + 1] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 + 1, b1 + 1] += area * d
 * 
 */
                __pyx_t_63 = (__pyx_v_b0 + 1);
                __pyx_t_64 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_63, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_64, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 307:                         outData[b0 + 1, b1 + 1] += area * d
                /* "histogram.pyx":307
 *                         rest -= area
 *                         outCount[b0 + 1, b1 + 1] += area
 *                         outData[b0 + 1, b1 + 1] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (delta0r) * (1 - delta1r)
 */
                __pyx_t_65 = (__pyx_v_b0 + 1);
                __pyx_t_66 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_65, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_66, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
 308: 
 309:                         area = (delta0r) * (1 - delta1r)
                /* "histogram.pyx":309
 *                         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));
 310:                         rest -= area
                /* "histogram.pyx":310
 * 
 *                         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);
 311:                         outCount[b0 + 1, b1 ] += area
                /* "histogram.pyx":311
 *                         area = (delta0r) * (1 - delta1r)
 *                         rest -= area
 *                         outCount[b0 + 1, b1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 + 1, b1 ] += area * d
 * 
 */
                __pyx_t_67 = (__pyx_v_b0 + 1);
                __pyx_t_68 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_67, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_68, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 312:                         outData[b0 + 1, b1 ] += area * d
                /* "histogram.pyx":312
 *                         rest -= area
 *                         outCount[b0 + 1, b1 ] += area
 *                         outData[b0 + 1, b1 ] += area * d             # <<<<<<<<<<<<<<
 * 
 *                         area = (1 - delta0r) * (delta1r)
 */
                __pyx_t_69 = (__pyx_v_b0 + 1);
                __pyx_t_70 = __pyx_v_b1;
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_69, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_70, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
 313: 
 314:                         area = (1 - delta0r) * (delta1r)
                /* "histogram.pyx":314
 *                         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);
 315:                         rest -= area
                /* "histogram.pyx":315
 * 
 *                         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);
 316:                         outCount[b0 , b1 + 1 ] += area
                /* "histogram.pyx":316
 *                         area = (1 - delta0r) * (delta1r)
 *                         rest -= area
 *                         outCount[b0 , b1 + 1 ] += area             # <<<<<<<<<<<<<<
 *                         outData[b0 , b1 + 1 ] += area * d
 *             outCount[b0, b1] += rest
 */
                __pyx_t_71 = __pyx_v_b0;
                __pyx_t_72 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_71, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_72, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_area;
 317:                         outData[b0 , b1 + 1 ] += area * d
                /* "histogram.pyx":317
 *                         rest -= area
 *                         outCount[b0 , b1 + 1 ] += area
 *                         outData[b0 , b1 + 1 ] += area * d             # <<<<<<<<<<<<<<
 *             outCount[b0, b1] += rest
 *             outData[b0, b1] += d * rest
 */
                __pyx_t_73 = __pyx_v_b0;
                __pyx_t_74 = (__pyx_v_b1 + 1);
                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_73, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_74, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_area * __pyx_v_d);
                goto __pyx_L29;
              }
              __pyx_L29:;
              goto __pyx_L27;
            }
            __pyx_L27:;
            goto __pyx_L26;
          }
          __pyx_L26:;
 318:             outCount[b0, b1] += rest
          /* "histogram.pyx":318
 *                         outCount[b0 , b1 + 1 ] += area
 *                         outData[b0 , b1 + 1 ] += area * d
 *             outCount[b0, b1] += rest             # <<<<<<<<<<<<<<
 *             outData[b0, b1] += d * rest
 * 
 */
          __pyx_t_75 = __pyx_v_b0;
          __pyx_t_76 = __pyx_v_b1;
          *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_75, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_76, __pyx_pybuffernd_outCount.diminfo[1].strides) += __pyx_v_rest;
 319:             outData[b0, b1] += d * rest
          /* "histogram.pyx":319
 *                         outData[b0 , b1 + 1 ] += area * d
 *             outCount[b0, b1] += rest
 *             outData[b0, b1] += d * rest             # <<<<<<<<<<<<<<
 * 
 *         for i in prange(bin0):
 */
          __pyx_t_77 = __pyx_v_b0;
          __pyx_t_78 = __pyx_v_b1;
          *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_77, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_78, __pyx_pybuffernd_outData.diminfo[1].strides) += (__pyx_v_d * __pyx_v_rest);
        }
 320: 
 321:         for i in prange(bin0):
        /* "histogram.pyx":321
 *             outData[b0, b1] += d * rest
 * 
 *         for i in prange(bin0):             # <<<<<<<<<<<<<<
 *             for j in range(bin1):
 *                 if outCount[i, j] > epsilon:
 */
        __pyx_t_11 = __pyx_v_bin0;
        if (1 == 0) abort();
        {
            #if ((defined(__APPLE__) || defined(__OSX__)) && (defined(__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)))))
                #undef likely
                #undef unlikely
                #define likely(x)   (x)
                #define unlikely(x) (x)
            #endif
            __pyx_t_79 = (__pyx_t_11 - 0) / 1;
            if (__pyx_t_79 > 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_79; __pyx_t_5++){
                        {
                            __pyx_v_i = 0 + 1 * __pyx_t_5;
                            /* Initialize private variables to invalid values */
                            __pyx_v_j = ((long)0xbad0bad0);
 322:             for j in range(bin1):
                            /* "histogram.pyx":322
 * 
 *         for i in prange(bin0):
 *             for j in range(bin1):             # <<<<<<<<<<<<<<
 *                 if outCount[i, j] > epsilon:
 *                     outMerge[i, j] += outData[i, j] / outCount[i, j]
 */
                            __pyx_t_80 = __pyx_v_bin1;
                            for (__pyx_t_81 = 0; __pyx_t_81 < __pyx_t_80; __pyx_t_81+=1) {
                              __pyx_v_j = __pyx_t_81;
 323:                 if outCount[i, j] > epsilon:
                              /* "histogram.pyx":323
 *         for i in prange(bin0):
 *             for j in range(bin1):
 *                 if outCount[i, j] > epsilon:             # <<<<<<<<<<<<<<
 *                     outMerge[i, j] += outData[i, j] / outCount[i, j]
 *                 else:
 */
                              __pyx_t_82 = __pyx_v_i;
                              __pyx_t_83 = __pyx_v_j;
                              __pyx_t_23 = ((*__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_82, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_83, __pyx_pybuffernd_outCount.diminfo[1].strides)) > __pyx_v_epsilon);
                              if (__pyx_t_23) {
 324:                     outMerge[i, j] += outData[i, j] / outCount[i, j]
                                /* "histogram.pyx":324
 *             for j in range(bin1):
 *                 if outCount[i, j] > epsilon:
 *                     outMerge[i, j] += outData[i, j] / outCount[i, j]             # <<<<<<<<<<<<<<
 *                 else:
 *                     outMerge[i, j] += dummy
 */
                                __pyx_t_84 = __pyx_v_i;
                                __pyx_t_85 = __pyx_v_j;
                                __pyx_t_86 = __pyx_v_i;
                                __pyx_t_87 = __pyx_v_j;
                                __pyx_t_88 = __pyx_v_i;
                                __pyx_t_89 = __pyx_v_j;
                                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.buf, __pyx_t_88, __pyx_pybuffernd_outMerge.diminfo[0].strides, __pyx_t_89, __pyx_pybuffernd_outMerge.diminfo[1].strides) += ((*__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outData.rcbuffer->pybuffer.buf, __pyx_t_84, __pyx_pybuffernd_outData.diminfo[0].strides, __pyx_t_85, __pyx_pybuffernd_outData.diminfo[1].strides)) / (*__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outCount.rcbuffer->pybuffer.buf, __pyx_t_86, __pyx_pybuffernd_outCount.diminfo[0].strides, __pyx_t_87, __pyx_pybuffernd_outCount.diminfo[1].strides)));
                                goto __pyx_L36;
                              }
                              /*else*/ {
 325:                 else:
 326:                     outMerge[i, j] += dummy
                                /* "histogram.pyx":326
 *                     outMerge[i, j] += outData[i, j] / outCount[i, j]
 *                 else:
 *                     outMerge[i, j] += dummy             # <<<<<<<<<<<<<<
 * 
 * 
 */
                                __pyx_t_90 = __pyx_v_i;
                                __pyx_t_91 = __pyx_v_j;
                                *__Pyx_BufPtrStrided2d(__pyx_t_5numpy_float64_t *, __pyx_pybuffernd_outMerge.rcbuffer->pybuffer.buf, __pyx_t_90, __pyx_pybuffernd_outMerge.diminfo[0].strides, __pyx_t_91, __pyx_pybuffernd_outMerge.diminfo[1].strides) += __pyx_v_dummy;
                              }
                              __pyx_L36:;
                            }
                        }
                    }
                }
            }
        }
        #if ((defined(__APPLE__) || defined(__OSX__)) && (defined(__GNUC__) && (__GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)))))
            #undef likely
            #undef unlikely
            #define likely(x)   __builtin_expect(!!(x), 1)
            #define unlikely(x) __builtin_expect(!!(x), 0)
        #endif
      }
 327: 
 328: 
 329:     return outMerge, edges0, edges1, outData, outCount
  /* "histogram.pyx":329
 * 
 * 
 *     return outMerge, edges0, edges1, outData, outCount             # <<<<<<<<<<<<<<
 */
  __Pyx_XDECREF(__pyx_r);
  __pyx_t_9 = PyTuple_New(5); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 329; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_9);
  __Pyx_INCREF(((PyObject *)__pyx_v_outMerge));
  PyTuple_SET_ITEM(__pyx_t_9, 0, ((PyObject *)__pyx_v_outMerge));
  __Pyx_GIVEREF(((PyObject *)__pyx_v_outMerge));
  __Pyx_INCREF(((PyObject *)__pyx_v_edges0));
  PyTuple_SET_ITEM(__pyx_t_9, 1, ((PyObject *)__pyx_v_edges0));
  __Pyx_GIVEREF(((PyObject *)__pyx_v_edges0));
  __Pyx_INCREF(((PyObject *)__pyx_v_edges1));
  PyTuple_SET_ITEM(__pyx_t_9, 2, ((PyObject *)__pyx_v_edges1));
  __Pyx_GIVEREF(((PyObject *)__pyx_v_edges1));
  __Pyx_INCREF(((PyObject *)__pyx_v_outData));
  PyTuple_SET_ITEM(__pyx_t_9, 3, ((PyObject *)__pyx_v_outData));
  __Pyx_GIVEREF(((PyObject *)__pyx_v_outData));
  __Pyx_INCREF(((PyObject *)__pyx_v_outCount));
  PyTuple_SET_ITEM(__pyx_t_9, 4, ((PyObject *)__pyx_v_outCount));
  __Pyx_GIVEREF(((PyObject *)__pyx_v_outCount));
  __pyx_r = ((PyObject *)__pyx_t_9);
  __pyx_t_9 = 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_9);
  __Pyx_XDECREF(__pyx_t_13);
  __PYX_XDEC_MEMVIEW(&__pyx_t_14, 1);
  __PYX_XDEC_MEMVIEW(&__pyx_t_15, 1);
  __PYX_XDEC_MEMVIEW(&__pyx_t_16, 1);
  { PyObject *__pyx_type, *__pyx_value, *__pyx_tb;
    __Pyx_ErrFetch(&__pyx_type, &__pyx_value, &__pyx_tb);
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_edges0.rcbuffer->pybuffer);
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_edges1.rcbuffer->pybuffer);
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outCount.rcbuffer->pybuffer);
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outData.rcbuffer->pybuffer);
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outMerge.rcbuffer->pybuffer);
  __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_pybuffernd_edges0.rcbuffer->pybuffer);
  __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_edges1.rcbuffer->pybuffer);
  __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outCount.rcbuffer->pybuffer);
  __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outData.rcbuffer->pybuffer);
  __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_outMerge.rcbuffer->pybuffer);
  __pyx_L2:;
  __PYX_XDEC_MEMVIEW(&__pyx_v_cpos0, 1);
  __PYX_XDEC_MEMVIEW(&__pyx_v_cpos1, 1);
  __PYX_XDEC_MEMVIEW(&__pyx_v_data, 1);
  __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;
}

/* Python wrapper */
static CYTHON_UNUSED int __pyx_pw_5numpy_7ndarray_1__getbuffer__(PyObject *__pyx_v_self, Py_buffer *__pyx_v_info, int __pyx_v_flags); /*proto*/
static CYTHON_UNUSED int __pyx_pw_5numpy_7ndarray_1__getbuffer__(PyObject *__pyx_v_self, Py_buffer *__pyx_v_info, int __pyx_v_flags) {
  int __pyx_r;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("__getbuffer__ (wrapper)", 0);
  __pyx_r = __pyx_pf_5numpy_7ndarray___getbuffer__(((PyArrayObject *)__pyx_v_self), ((Py_buffer *)__pyx_v_info), ((int)__pyx_v_flags));
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}