#

Generated by Cython 0.21

#define PY_SSIZE_T_CLEAN
#ifndef CYTHON_USE_PYLONG_INTERNALS
#ifdef PYLONG_BITS_IN_DIGIT
#define CYTHON_USE_PYLONG_INTERNALS 0
#else
#include "pyconfig.h"
#ifdef PYLONG_BITS_IN_DIGIT
#define CYTHON_USE_PYLONG_INTERNALS 1
#else
#define CYTHON_USE_PYLONG_INTERNALS 0
#endif
#endif
#endif
#include "Python.h"
#ifndef Py_PYTHON_H
    #error Python headers needed to compile C extensions, please install development version of Python.
#elif PY_VERSION_HEX < 0x02060000 || (0x03000000 <= PY_VERSION_HEX && PY_VERSION_HEX < 0x03020000)
    #error Cython requires Python 2.6+ or Python 3.2+.
#else
#define CYTHON_ABI "0_21"
#include <stddef.h>
#ifndef offsetof
#define offsetof(type, member) ( (size_t) & ((type*)0) -> member )
#endif
#if !defined(WIN32) && !defined(MS_WINDOWS)
  #ifndef __stdcall
    #define __stdcall
  #endif
  #ifndef __cdecl
    #define __cdecl
  #endif
  #ifndef __fastcall
    #define __fastcall
  #endif
#endif
#ifndef DL_IMPORT
  #define DL_IMPORT(t) t
#endif
#ifndef DL_EXPORT
  #define DL_EXPORT(t) t
#endif
#ifndef PY_LONG_LONG
  #define PY_LONG_LONG LONG_LONG
#endif
#ifndef Py_HUGE_VAL
  #define Py_HUGE_VAL HUGE_VAL
#endif
#ifdef PYPY_VERSION
#define CYTHON_COMPILING_IN_PYPY 1
#define CYTHON_COMPILING_IN_CPYTHON 0
#else
#define CYTHON_COMPILING_IN_PYPY 0
#define CYTHON_COMPILING_IN_CPYTHON 1
#endif
#if CYTHON_COMPILING_IN_PYPY && PY_VERSION_HEX < 0x02070600
#define Py_OptimizeFlag 0
#endif
#define __PYX_BUILD_PY_SSIZE_T "n"
#define CYTHON_FORMAT_SSIZE_T "z"
#if PY_MAJOR_VERSION < 3
  #define __Pyx_BUILTIN_MODULE_NAME "__builtin__"
  #define __Pyx_PyCode_New(a, k, l, s, f, code, c, n, v, fv, cell, fn, name, fline, lnos) \
          PyCode_New(a+k, l, s, f, code, c, n, v, fv, cell, fn, name, fline, lnos)
  #define __Pyx_DefaultClassType PyClass_Type
#else
  #define __Pyx_BUILTIN_MODULE_NAME "builtins"
  #define __Pyx_PyCode_New(a, k, l, s, f, code, c, n, v, fv, cell, fn, name, fline, lnos) \
          PyCode_New(a, k, l, s, f, code, c, n, v, fv, cell, fn, name, fline, lnos)
  #define __Pyx_DefaultClassType PyType_Type
#endif
#if PY_MAJOR_VERSION >= 3
  #define Py_TPFLAGS_CHECKTYPES 0
  #define Py_TPFLAGS_HAVE_INDEX 0
#endif
#if PY_MAJOR_VERSION >= 3
  #define Py_TPFLAGS_HAVE_NEWBUFFER 0
#endif
#if PY_VERSION_HEX < 0x030400a1 && !defined(Py_TPFLAGS_HAVE_FINALIZE)
  #define Py_TPFLAGS_HAVE_FINALIZE 0
#endif
#if PY_VERSION_HEX > 0x03030000 && defined(PyUnicode_KIND)
  #define CYTHON_PEP393_ENABLED 1
  #define __Pyx_PyUnicode_READY(op)       (likely(PyUnicode_IS_READY(op)) ? \
                                              0 : _PyUnicode_Ready((PyObject *)(op)))
  #define __Pyx_PyUnicode_GET_LENGTH(u)   PyUnicode_GET_LENGTH(u)
  #define __Pyx_PyUnicode_READ_CHAR(u, i) PyUnicode_READ_CHAR(u, i)
  #define __Pyx_PyUnicode_KIND(u)         PyUnicode_KIND(u)
  #define __Pyx_PyUnicode_DATA(u)         PyUnicode_DATA(u)
  #define __Pyx_PyUnicode_READ(k, d, i)   PyUnicode_READ(k, d, i)
#else
  #define CYTHON_PEP393_ENABLED 0
  #define __Pyx_PyUnicode_READY(op)       (0)
  #define __Pyx_PyUnicode_GET_LENGTH(u)   PyUnicode_GET_SIZE(u)
  #define __Pyx_PyUnicode_READ_CHAR(u, i) ((Py_UCS4)(PyUnicode_AS_UNICODE(u)[i]))
  #define __Pyx_PyUnicode_KIND(u)         (sizeof(Py_UNICODE))
  #define __Pyx_PyUnicode_DATA(u)         ((void*)PyUnicode_AS_UNICODE(u))
  #define __Pyx_PyUnicode_READ(k, d, i)   ((void)(k), (Py_UCS4)(((Py_UNICODE*)d)[i]))
#endif
#if CYTHON_COMPILING_IN_PYPY
  #define __Pyx_PyUnicode_Concat(a, b)      PyNumber_Add(a, b)
  #define __Pyx_PyUnicode_ConcatSafe(a, b)  PyNumber_Add(a, b)
#else
  #define __Pyx_PyUnicode_Concat(a, b)      PyUnicode_Concat(a, b)
  #define __Pyx_PyUnicode_ConcatSafe(a, b)  ((unlikely((a) == Py_None) || unlikely((b) == Py_None)) ? \
      PyNumber_Add(a, b) : __Pyx_PyUnicode_Concat(a, b))
#endif
#define __Pyx_PyString_FormatSafe(a, b)   ((unlikely((a) == Py_None)) ? PyNumber_Remainder(a, b) : __Pyx_PyString_Format(a, b))
#define __Pyx_PyUnicode_FormatSafe(a, b)  ((unlikely((a) == Py_None)) ? PyNumber_Remainder(a, b) : PyUnicode_Format(a, b))
#if PY_MAJOR_VERSION >= 3
  #define __Pyx_PyString_Format(a, b)  PyUnicode_Format(a, b)
#else
  #define __Pyx_PyString_Format(a, b)  PyString_Format(a, b)
#endif
#if PY_MAJOR_VERSION >= 3
  #define PyBaseString_Type            PyUnicode_Type
  #define PyStringObject               PyUnicodeObject
  #define PyString_Type                PyUnicode_Type
  #define PyString_Check               PyUnicode_Check
  #define PyString_CheckExact          PyUnicode_CheckExact
#endif
#if PY_MAJOR_VERSION >= 3
  #define __Pyx_PyBaseString_Check(obj) PyUnicode_Check(obj)
  #define __Pyx_PyBaseString_CheckExact(obj) PyUnicode_CheckExact(obj)
#else
  #define __Pyx_PyBaseString_Check(obj) (PyString_Check(obj) || PyUnicode_Check(obj))
  #define __Pyx_PyBaseString_CheckExact(obj) (PyString_CheckExact(obj) || PyUnicode_CheckExact(obj))
#endif
#ifndef PySet_CheckExact
  #define PySet_CheckExact(obj)        (Py_TYPE(obj) == &PySet_Type)
#endif
#define __Pyx_TypeCheck(obj, type) PyObject_TypeCheck(obj, (PyTypeObject *)type)
#if PY_MAJOR_VERSION >= 3
  #define PyIntObject                  PyLongObject
  #define PyInt_Type                   PyLong_Type
  #define PyInt_Check(op)              PyLong_Check(op)
  #define PyInt_CheckExact(op)         PyLong_CheckExact(op)
  #define PyInt_FromString             PyLong_FromString
  #define PyInt_FromUnicode            PyLong_FromUnicode
  #define PyInt_FromLong               PyLong_FromLong
  #define PyInt_FromSize_t             PyLong_FromSize_t
  #define PyInt_FromSsize_t            PyLong_FromSsize_t
  #define PyInt_AsLong                 PyLong_AsLong
  #define PyInt_AS_LONG                PyLong_AS_LONG
  #define PyInt_AsSsize_t              PyLong_AsSsize_t
  #define PyInt_AsUnsignedLongMask     PyLong_AsUnsignedLongMask
  #define PyInt_AsUnsignedLongLongMask PyLong_AsUnsignedLongLongMask
  #define PyNumber_Int                 PyNumber_Long
#endif
#if PY_MAJOR_VERSION >= 3
  #define PyBoolObject                 PyLongObject
#endif
#if PY_VERSION_HEX < 0x030200A4
  typedef long Py_hash_t;
  #define __Pyx_PyInt_FromHash_t PyInt_FromLong
  #define __Pyx_PyInt_AsHash_t   PyInt_AsLong
#else
  #define __Pyx_PyInt_FromHash_t PyInt_FromSsize_t
  #define __Pyx_PyInt_AsHash_t   PyInt_AsSsize_t
#endif
#if PY_MAJOR_VERSION >= 3
  #define PyMethod_New(func, self, klass) ((self) ? PyMethod_New(func, self) : PyInstanceMethod_New(func))
#endif
#ifndef CYTHON_INLINE
  #if defined(__GNUC__)
    #define CYTHON_INLINE __inline__
  #elif defined(_MSC_VER)
    #define CYTHON_INLINE __inline
  #elif defined (__STDC_VERSION__) && __STDC_VERSION__ >= 199901L
    #define CYTHON_INLINE inline
  #else
    #define CYTHON_INLINE
  #endif
#endif
#ifndef CYTHON_RESTRICT
  #if defined(__GNUC__)
    #define CYTHON_RESTRICT __restrict__
  #elif defined(_MSC_VER) && _MSC_VER >= 1400
    #define CYTHON_RESTRICT __restrict
  #elif defined (__STDC_VERSION__) && __STDC_VERSION__ >= 199901L
    #define CYTHON_RESTRICT restrict
  #else
    #define CYTHON_RESTRICT
  #endif
#endif
#ifdef NAN
#define __PYX_NAN() ((float) NAN)
#else
static CYTHON_INLINE float __PYX_NAN() {

#

Initialize NaN. The sign is irrelevant, an exponent with all bits 1 and a nonzero mantissa means NaN. If the first bit in the mantissa is 1, it is a quiet NaN.

  float value;
  memset(&value, 0xFF, sizeof(value));
  return value;
}
#endif
#ifdef __cplusplus
template<typename T>

#

void __Pyx_call_destructor(T* x) {
    x->~T();
}
#endif


#if PY_MAJOR_VERSION >= 3
  #define __Pyx_PyNumber_Divide(x,y)         PyNumber_TrueDivide(x,y)
  #define __Pyx_PyNumber_InPlaceDivide(x,y)  PyNumber_InPlaceTrueDivide(x,y)
#else
  #define __Pyx_PyNumber_Divide(x,y)         PyNumber_TrueDivide(x,y)
  #define __Pyx_PyNumber_InPlaceDivide(x,y)  PyNumber_InPlaceTrueDivide(x,y)
#endif

#ifndef __PYX_EXTERN_C
  #ifdef __cplusplus
    #define __PYX_EXTERN_C extern "C"
  #else
    #define __PYX_EXTERN_C extern
  #endif
#endif

#if defined(WIN32) || defined(MS_WINDOWS)
#define _USE_MATH_DEFINES
#endif
#include <math.h>
#define __PYX_HAVE__src__particles__interpolation__interpolation
#define __PYX_HAVE_API__src__particles__interpolation__interpolation
#include "string.h"
#include "stdio.h"
#include "stdlib.h"
#include "numpy/arrayobject.h"
#include "numpy/ufuncobject.h"
#ifdef _OPENMP
#include <omp.h>
#endif /* _OPENMP */

#ifdef PYREX_WITHOUT_ASSERTIONS
#define CYTHON_WITHOUT_ASSERTIONS
#endif

#ifndef CYTHON_UNUSED
# if defined(__GNUC__)
#   if !(defined(__cplusplus)) || (__GNUC__ > 3 || (__GNUC__ == 3 && __GNUC_MINOR__ >= 4))
#     define CYTHON_UNUSED __attribute__ ((__unused__))
#   else
#     define CYTHON_UNUSED
#   endif
# elif defined(__ICC) || (defined(__INTEL_COMPILER) && !defined(_MSC_VER))
#   define CYTHON_UNUSED __attribute__ ((__unused__))
# else
#   define CYTHON_UNUSED
# endif
#endif
typedef struct {PyObject **p; char *s; const Py_ssize_t n; const char* encoding;
                const char is_unicode; const char is_str; const char intern; } __Pyx_StringTabEntry;

#define __PYX_DEFAULT_STRING_ENCODING_IS_ASCII 0
#define __PYX_DEFAULT_STRING_ENCODING_IS_DEFAULT 0
#define __PYX_DEFAULT_STRING_ENCODING ""
#define __Pyx_PyObject_FromString __Pyx_PyBytes_FromString
#define __Pyx_PyObject_FromStringAndSize __Pyx_PyBytes_FromStringAndSize
#define __Pyx_fits_Py_ssize_t(v, type, is_signed)  (    \
    (sizeof(type) < sizeof(Py_ssize_t))  ||             \
    (sizeof(type) > sizeof(Py_ssize_t) &&               \
          likely(v < (type)PY_SSIZE_T_MAX ||            \
                 v == (type)PY_SSIZE_T_MAX)  &&         \
          (!is_signed || likely(v > (type)PY_SSIZE_T_MIN ||       \
                                v == (type)PY_SSIZE_T_MIN)))  ||  \
    (sizeof(type) == sizeof(Py_ssize_t) &&              \
          (is_signed || likely(v < (type)PY_SSIZE_T_MAX ||        \
                               v == (type)PY_SSIZE_T_MAX)))  )
static CYTHON_INLINE char* __Pyx_PyObject_AsString(PyObject*);
static CYTHON_INLINE char* __Pyx_PyObject_AsStringAndSize(PyObject*, Py_ssize_t* length);
#define __Pyx_PyByteArray_FromString(s) PyByteArray_FromStringAndSize((const char*)s, strlen((const char*)s))
#define __Pyx_PyByteArray_FromStringAndSize(s, l) PyByteArray_FromStringAndSize((const char*)s, l)
#define __Pyx_PyBytes_FromString        PyBytes_FromString
#define __Pyx_PyBytes_FromStringAndSize PyBytes_FromStringAndSize
static CYTHON_INLINE PyObject* __Pyx_PyUnicode_FromString(const char*);
#if PY_MAJOR_VERSION < 3
    #define __Pyx_PyStr_FromString        __Pyx_PyBytes_FromString
    #define __Pyx_PyStr_FromStringAndSize __Pyx_PyBytes_FromStringAndSize
#else
    #define __Pyx_PyStr_FromString        __Pyx_PyUnicode_FromString
    #define __Pyx_PyStr_FromStringAndSize __Pyx_PyUnicode_FromStringAndSize
#endif
#define __Pyx_PyObject_AsSString(s)    ((signed char*) __Pyx_PyObject_AsString(s))
#define __Pyx_PyObject_AsUString(s)    ((unsigned char*) __Pyx_PyObject_AsString(s))
#define __Pyx_PyObject_FromUString(s)  __Pyx_PyObject_FromString((const char*)s)
#define __Pyx_PyBytes_FromUString(s)   __Pyx_PyBytes_FromString((const char*)s)
#define __Pyx_PyByteArray_FromUString(s)   __Pyx_PyByteArray_FromString((const char*)s)
#define __Pyx_PyStr_FromUString(s)     __Pyx_PyStr_FromString((const char*)s)
#define __Pyx_PyUnicode_FromUString(s) __Pyx_PyUnicode_FromString((const char*)s)
#if PY_MAJOR_VERSION < 3
static CYTHON_INLINE size_t __Pyx_Py_UNICODE_strlen(const Py_UNICODE *u)
{
    const Py_UNICODE *u_end = u;
    while (*u_end++) ;
    return (size_t)(u_end - u - 1);
}
#else
#define __Pyx_Py_UNICODE_strlen Py_UNICODE_strlen
#endif
#define __Pyx_PyUnicode_FromUnicode(u)       PyUnicode_FromUnicode(u, __Pyx_Py_UNICODE_strlen(u))
#define __Pyx_PyUnicode_FromUnicodeAndLength PyUnicode_FromUnicode
#define __Pyx_PyUnicode_AsUnicode            PyUnicode_AsUnicode
#define __Pyx_Owned_Py_None(b) (Py_INCREF(Py_None), Py_None)
#define __Pyx_PyBool_FromLong(b) ((b) ? (Py_INCREF(Py_True), Py_True) : (Py_INCREF(Py_False), Py_False))
static CYTHON_INLINE int __Pyx_PyObject_IsTrue(PyObject*);
static CYTHON_INLINE PyObject* __Pyx_PyNumber_Int(PyObject* x);
static CYTHON_INLINE Py_ssize_t __Pyx_PyIndex_AsSsize_t(PyObject*);
static CYTHON_INLINE PyObject * __Pyx_PyInt_FromSize_t(size_t);
#if CYTHON_COMPILING_IN_CPYTHON
#define __pyx_PyFloat_AsDouble(x) (PyFloat_CheckExact(x) ? PyFloat_AS_DOUBLE(x) : PyFloat_AsDouble(x))
#else
#define __pyx_PyFloat_AsDouble(x) PyFloat_AsDouble(x)
#endif
#define __pyx_PyFloat_AsFloat(x) ((float) __pyx_PyFloat_AsDouble(x))
#if PY_MAJOR_VERSION < 3 && __PYX_DEFAULT_STRING_ENCODING_IS_ASCII
static int __Pyx_sys_getdefaultencoding_not_ascii;
static int __Pyx_init_sys_getdefaultencoding_params(void) {
    PyObject* sys;
    PyObject* default_encoding = NULL;
    PyObject* ascii_chars_u = NULL;
    PyObject* ascii_chars_b = NULL;
    const char* default_encoding_c;
    sys = PyImport_ImportModule("sys");
    if (!sys) goto bad;
    default_encoding = PyObject_CallMethod(sys, (char*) (const char*) "getdefaultencoding", NULL);
    Py_DECREF(sys);
    if (!default_encoding) goto bad;
    default_encoding_c = PyBytes_AsString(default_encoding);
    if (!default_encoding_c) goto bad;
    if (strcmp(default_encoding_c, "ascii") == 0) {
        __Pyx_sys_getdefaultencoding_not_ascii = 0;
    } else {
        char ascii_chars[128];

        int c;
        for (c = 0; c < 128; c++) {
            ascii_chars[c] = c;
        }
        __Pyx_sys_getdefaultencoding_not_ascii = 1;
        ascii_chars_u = PyUnicode_DecodeASCII(ascii_chars, 128, NULL);
        if (!ascii_chars_u) goto bad;
        ascii_chars_b = PyUnicode_AsEncodedString(ascii_chars_u, default_encoding_c, NULL);
        if (!ascii_chars_b || !PyBytes_Check(ascii_chars_b) || memcmp(ascii_chars, PyBytes_AS_STRING(ascii_chars_b), 128) != 0) {
            PyErr_Format(
                PyExc_ValueError,
                "This module compiled with c_string_encoding=ascii, but default encoding '%.200s' is not a superset of ascii.",
                default_encoding_c);
            goto bad;
        }
        Py_DECREF(ascii_chars_u);
        Py_DECREF(ascii_chars_b);
    }
    Py_DECREF(default_encoding);
    return 0;
bad:
    Py_XDECREF(default_encoding);
    Py_XDECREF(ascii_chars_u);
    Py_XDECREF(ascii_chars_b);
    return -1;
}
#endif
#if __PYX_DEFAULT_STRING_ENCODING_IS_DEFAULT && PY_MAJOR_VERSION >= 3
#define __Pyx_PyUnicode_FromStringAndSize(c_str, size) PyUnicode_DecodeUTF8(c_str, size, NULL)
#else
#define __Pyx_PyUnicode_FromStringAndSize(c_str, size) PyUnicode_Decode(c_str, size, __PYX_DEFAULT_STRING_ENCODING, NULL)
#if __PYX_DEFAULT_STRING_ENCODING_IS_DEFAULT
static char* __PYX_DEFAULT_STRING_ENCODING;
static int __Pyx_init_sys_getdefaultencoding_params(void) {
    PyObject* sys;
    PyObject* default_encoding = NULL;
    char* default_encoding_c;
    sys = PyImport_ImportModule("sys");
    if (!sys) goto bad;
    default_encoding = PyObject_CallMethod(sys, (char*) (const char*) "getdefaultencoding", NULL);
    Py_DECREF(sys);
    if (!default_encoding) goto bad;
    default_encoding_c = PyBytes_AsString(default_encoding);
    if (!default_encoding_c) goto bad;
    __PYX_DEFAULT_STRING_ENCODING = (char*) malloc(strlen(default_encoding_c));
    if (!__PYX_DEFAULT_STRING_ENCODING) goto bad;
    strcpy(__PYX_DEFAULT_STRING_ENCODING, default_encoding_c);
    Py_DECREF(default_encoding);
    return 0;
bad:
    Py_XDECREF(default_encoding);
    return -1;
}
#endif
#endif

#

Test for GCC > 2.95

#if defined(__GNUC__)     && (__GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)))
  #define likely(x)   __builtin_expect(!!(x), 1)
  #define unlikely(x) __builtin_expect(!!(x), 0)
#else /* !__GNUC__ or GCC < 2.95 */
  #define likely(x)   (x)
  #define unlikely(x) (x)
#endif /* __GNUC__ */

static PyObject *__pyx_m;
static PyObject *__pyx_d;
static PyObject *__pyx_b;
static PyObject *__pyx_empty_tuple;
static PyObject *__pyx_empty_bytes;
static int __pyx_lineno;
static int __pyx_clineno = 0;
static const char * __pyx_cfilenm= __FILE__;
static const char *__pyx_filename;

#if !defined(CYTHON_CCOMPLEX)
  #if defined(__cplusplus)
    #define CYTHON_CCOMPLEX 1
  #elif defined(_Complex_I)
    #define CYTHON_CCOMPLEX 1
  #else
    #define CYTHON_CCOMPLEX 0
  #endif
#endif
#if CYTHON_CCOMPLEX
  #ifdef __cplusplus
    #include <complex>
  #else
    #include <complex.h>
  #endif
#endif
#if CYTHON_CCOMPLEX && !defined(__cplusplus) && defined(__sun__) && defined(__GNUC__)
  #undef _Complex_I
  #define _Complex_I 1.0fj
#endif


static const char *__pyx_f[] = {
  "src/particles/interpolation/interpolation.pyx",
  "__init__.pxd",
  "type.pxd",
};
#define IS_UNSIGNED(type) (((type) -1) > 0)
struct __Pyx_StructField_;
#define __PYX_BUF_FLAGS_PACKED_STRUCT (1 << 0)
typedef struct {
  const char* name;
  struct __Pyx_StructField_* fields;
  size_t size;
  size_t arraysize[8];

  int ndim;
  char typegroup;
  char is_unsigned;

#

  int flags;
} __Pyx_TypeInfo;
typedef struct __Pyx_StructField_ {
  __Pyx_TypeInfo* type;
  const char* name;
  size_t offset;
} __Pyx_StructField;
typedef struct {
  __Pyx_StructField* field;
  size_t parent_offset;
} __Pyx_BufFmt_StackElem;
typedef struct {
  __Pyx_StructField root;
  __Pyx_BufFmt_StackElem* head;
  size_t fmt_offset;
  size_t new_count, enc_count;
  size_t struct_alignment;

#

  int is_complex;
  char enc_type;
  char new_packmode;
  char enc_packmode;
  char is_valid_array;
} __Pyx_BufFmt_Context;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":723

in Cython to enable them only on the right systems.

ctypedef npy_int8 int8_t # <<<<<<<<<<<<<< ctypedef npy_int16 int16_t ctypedef npy_int32 int32_t

typedef npy_int8 __pyx_t_5numpy_int8_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":724

ctypedef npy_int8 int8_t ctypedef npy_int16 int16_t # <<<<<<<<<<<<<< ctypedef npy_int32 int32_t ctypedef npy_int64 int64_t

typedef npy_int16 __pyx_t_5numpy_int16_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":725 ctypedef npy_int8 int8_t ctypedef npy_int16 int16_t ctypedef npy_int32 int32_t # <<<<<<<<<<<<<< ctypedef npy_int64 int64_t

ctypedef npy_int96 int96_t

typedef npy_int32 __pyx_t_5numpy_int32_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":726 ctypedef npy_int16 int16_t ctypedef npy_int32 int32_t ctypedef npy_int64 int64_t # <<<<<<<<<<<<<<

ctypedef npy_int96 int96_t

ctypedef npy_int128 int128_t

typedef npy_int64 __pyx_t_5numpy_int64_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":730

ctypedef npy_int128 int128_t

ctypedef npy_uint8 uint8_t # <<<<<<<<<<<<<< ctypedef npy_uint16 uint16_t ctypedef npy_uint32 uint32_t

typedef npy_uint8 __pyx_t_5numpy_uint8_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":731

ctypedef npy_uint8 uint8_t ctypedef npy_uint16 uint16_t # <<<<<<<<<<<<<< ctypedef npy_uint32 uint32_t ctypedef npy_uint64 uint64_t

typedef npy_uint16 __pyx_t_5numpy_uint16_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":732 ctypedef npy_uint8 uint8_t ctypedef npy_uint16 uint16_t ctypedef npy_uint32 uint32_t # <<<<<<<<<<<<<< ctypedef npy_uint64 uint64_t

ctypedef npy_uint96 uint96_t

typedef npy_uint32 __pyx_t_5numpy_uint32_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":733 ctypedef npy_uint16 uint16_t ctypedef npy_uint32 uint32_t ctypedef npy_uint64 uint64_t # <<<<<<<<<<<<<<

ctypedef npy_uint96 uint96_t

ctypedef npy_uint128 uint128_t

typedef npy_uint64 __pyx_t_5numpy_uint64_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":737

ctypedef npy_uint128 uint128_t

ctypedef npy_float32 float32_t # <<<<<<<<<<<<<< ctypedef npy_float64 float64_t

ctypedef npy_float80 float80_t

typedef npy_float32 __pyx_t_5numpy_float32_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":738

ctypedef npy_float32 float32_t ctypedef npy_float64 float64_t # <<<<<<<<<<<<<<

ctypedef npy_float80 float80_t

ctypedef npy_float128 float128_t

typedef npy_float64 __pyx_t_5numpy_float64_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":747

The int types are mapped a bit surprising --

numpy.int corresponds to 'l' and numpy.long to 'q'

ctypedef npy_long int_t # <<<<<<<<<<<<<< ctypedef npy_longlong long_t ctypedef npy_longlong longlong_t

typedef npy_long __pyx_t_5numpy_int_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":748

numpy.int corresponds to 'l' and numpy.long to 'q'

ctypedef npy_long int_t ctypedef npy_longlong long_t # <<<<<<<<<<<<<< ctypedef npy_longlong longlong_t

typedef npy_longlong __pyx_t_5numpy_long_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":749 ctypedef npy_long int_t ctypedef npy_longlong long_t ctypedef npy_longlong longlong_t # <<<<<<<<<<<<<<

ctypedef npy_ulong uint_t

typedef npy_longlong __pyx_t_5numpy_longlong_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":751 ctypedef npy_longlong longlong_t

ctypedef npy_ulong uint_t # <<<<<<<<<<<<<< ctypedef npy_ulonglong ulong_t ctypedef npy_ulonglong ulonglong_t

typedef npy_ulong __pyx_t_5numpy_uint_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":752

ctypedef npy_ulong uint_t ctypedef npy_ulonglong ulong_t # <<<<<<<<<<<<<< ctypedef npy_ulonglong ulonglong_t

typedef npy_ulonglong __pyx_t_5numpy_ulong_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":753 ctypedef npy_ulong uint_t ctypedef npy_ulonglong ulong_t ctypedef npy_ulonglong ulonglong_t # <<<<<<<<<<<<<<

ctypedef npy_intp intp_t

typedef npy_ulonglong __pyx_t_5numpy_ulonglong_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":755 ctypedef npy_ulonglong ulonglong_t

ctypedef npy_intp intp_t # <<<<<<<<<<<<<< ctypedef npy_uintp uintp_t

typedef npy_intp __pyx_t_5numpy_intp_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":756

ctypedef npy_intp intp_t ctypedef npy_uintp uintp_t # <<<<<<<<<<<<<<

ctypedef npy_double float_t

typedef npy_uintp __pyx_t_5numpy_uintp_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":758 ctypedef npy_uintp uintp_t

ctypedef npy_double float_t # <<<<<<<<<<<<<< ctypedef npy_double double_t ctypedef npy_longdouble longdouble_t

typedef npy_double __pyx_t_5numpy_float_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":759

ctypedef npy_double float_t ctypedef npy_double double_t # <<<<<<<<<<<<<< ctypedef npy_longdouble longdouble_t

typedef npy_double __pyx_t_5numpy_double_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":760 ctypedef npy_double float_t ctypedef npy_double double_t ctypedef npy_longdouble longdouble_t # <<<<<<<<<<<<<<

ctypedef npy_cfloat cfloat_t

typedef npy_longdouble __pyx_t_5numpy_longdouble_t;

#

"src/particles/interpolation/interpolation.pyx":19

every type in the numpy module there's a corresponding compile-time

type with a _t-suffix.

ctypedef numpy.float_t DTYPE_t # <<<<<<<<<<<<<<

cdef int exponential_interpolation = True

typedef __pyx_t_5numpy_float_t __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t;
#if CYTHON_CCOMPLEX
  #ifdef __cplusplus
    typedef ::std::complex< float > __pyx_t_float_complex;
  #else
    typedef float _Complex __pyx_t_float_complex;
  #endif
#else
    typedef struct { float real, imag; } __pyx_t_float_complex;
#endif

#if CYTHON_CCOMPLEX
  #ifdef __cplusplus
    typedef ::std::complex< double > __pyx_t_double_complex;
  #else
    typedef double _Complex __pyx_t_double_complex;
  #endif
#else
    typedef struct { double real, imag; } __pyx_t_double_complex;
#endif

#

--- Type declarations ---

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":762 ctypedef npy_longdouble longdouble_t

ctypedef npy_cfloat cfloat_t # <<<<<<<<<<<<<< ctypedef npy_cdouble cdouble_t ctypedef npy_clongdouble clongdouble_t

typedef npy_cfloat __pyx_t_5numpy_cfloat_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":763

ctypedef npy_cfloat cfloat_t ctypedef npy_cdouble cdouble_t # <<<<<<<<<<<<<< ctypedef npy_clongdouble clongdouble_t

typedef npy_cdouble __pyx_t_5numpy_cdouble_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":764 ctypedef npy_cfloat cfloat_t ctypedef npy_cdouble cdouble_t ctypedef npy_clongdouble clongdouble_t # <<<<<<<<<<<<<<

ctypedef npy_cdouble complex_t

typedef npy_clongdouble __pyx_t_5numpy_clongdouble_t;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":766 ctypedef npy_clongdouble clongdouble_t

ctypedef npy_cdouble complex_t # <<<<<<<<<<<<<<

cdef inline object PyArray_MultiIterNew1(a):

typedef npy_cdouble __pyx_t_5numpy_complex_t;
#ifndef CYTHON_REFNANNY
  #define CYTHON_REFNANNY 0
#endif
#if CYTHON_REFNANNY
  typedef struct {
    void (*INCREF)(void*, PyObject*, int);
    void (*DECREF)(void*, PyObject*, int);
    void (*GOTREF)(void*, PyObject*, int);
    void (*GIVEREF)(void*, PyObject*, int);
    void* (*SetupContext)(const char*, int, const char*);
    void (*FinishContext)(void**);
  } __Pyx_RefNannyAPIStruct;
  static __Pyx_RefNannyAPIStruct *__Pyx_RefNanny = NULL;
  static __Pyx_RefNannyAPIStruct *__Pyx_RefNannyImportAPI(const char *modname);
  #define __Pyx_RefNannyDeclarations void *__pyx_refnanny = NULL;
#ifdef WITH_THREAD
  #define __Pyx_RefNannySetupContext(name, acquire_gil) \
          if (acquire_gil) { \
              PyGILState_STATE __pyx_gilstate_save = PyGILState_Ensure(); \
              __pyx_refnanny = __Pyx_RefNanny->SetupContext((name), __LINE__, __FILE__); \
              PyGILState_Release(__pyx_gilstate_save); \
          } else { \
              __pyx_refnanny = __Pyx_RefNanny->SetupContext((name), __LINE__, __FILE__); \
          }
#else
  #define __Pyx_RefNannySetupContext(name, acquire_gil) \
          __pyx_refnanny = __Pyx_RefNanny->SetupContext((name), __LINE__, __FILE__)
#endif
  #define __Pyx_RefNannyFinishContext() \
          __Pyx_RefNanny->FinishContext(&__pyx_refnanny)
  #define __Pyx_INCREF(r)  __Pyx_RefNanny->INCREF(__pyx_refnanny, (PyObject *)(r), __LINE__)
  #define __Pyx_DECREF(r)  __Pyx_RefNanny->DECREF(__pyx_refnanny, (PyObject *)(r), __LINE__)
  #define __Pyx_GOTREF(r)  __Pyx_RefNanny->GOTREF(__pyx_refnanny, (PyObject *)(r), __LINE__)
  #define __Pyx_GIVEREF(r) __Pyx_RefNanny->GIVEREF(__pyx_refnanny, (PyObject *)(r), __LINE__)
  #define __Pyx_XINCREF(r)  do { if((r) != NULL) {__Pyx_INCREF(r); }} while(0)
  #define __Pyx_XDECREF(r)  do { if((r) != NULL) {__Pyx_DECREF(r); }} while(0)
  #define __Pyx_XGOTREF(r)  do { if((r) != NULL) {__Pyx_GOTREF(r); }} while(0)
  #define __Pyx_XGIVEREF(r) do { if((r) != NULL) {__Pyx_GIVEREF(r);}} while(0)
#else
  #define __Pyx_RefNannyDeclarations
  #define __Pyx_RefNannySetupContext(name, acquire_gil)
  #define __Pyx_RefNannyFinishContext()
  #define __Pyx_INCREF(r) Py_INCREF(r)
  #define __Pyx_DECREF(r) Py_DECREF(r)
  #define __Pyx_GOTREF(r)
  #define __Pyx_GIVEREF(r)
  #define __Pyx_XINCREF(r) Py_XINCREF(r)
  #define __Pyx_XDECREF(r) Py_XDECREF(r)
  #define __Pyx_XGOTREF(r)
  #define __Pyx_XGIVEREF(r)
#endif
#define __Pyx_XDECREF_SET(r, v) do {                            \
        PyObject *tmp = (PyObject *) r;                         \
        r = v; __Pyx_XDECREF(tmp);                              \
    } while (0)
#define __Pyx_DECREF_SET(r, v) do {                             \
        PyObject *tmp = (PyObject *) r;                         \
        r = v; __Pyx_DECREF(tmp);                               \
    } while (0)
#define __Pyx_CLEAR(r)    do { PyObject* tmp = ((PyObject*)(r)); r = NULL; __Pyx_DECREF(tmp);} while(0)
#define __Pyx_XCLEAR(r)   do { if((r) != NULL) {PyObject* tmp = ((PyObject*)(r)); r = NULL; __Pyx_DECREF(tmp);}} while(0)

#if CYTHON_COMPILING_IN_CPYTHON
static CYTHON_INLINE PyObject* __Pyx_PyObject_GetAttrStr(PyObject* obj, PyObject* attr_name) {
    PyTypeObject* tp = Py_TYPE(obj);
    if (likely(tp->tp_getattro))
        return tp->tp_getattro(obj, attr_name);
#if PY_MAJOR_VERSION < 3
    if (likely(tp->tp_getattr))
        return tp->tp_getattr(obj, PyString_AS_STRING(attr_name));
#endif
    return PyObject_GetAttr(obj, attr_name);
}
#else
#define __Pyx_PyObject_GetAttrStr(o,n) PyObject_GetAttr(o,n)
#endif

static PyObject *__Pyx_GetBuiltinName(PyObject *name);

static CYTHON_INLINE PyObject *__Pyx_GetModuleGlobalName(PyObject *name);

#if CYTHON_COMPILING_IN_CPYTHON
static CYTHON_INLINE PyObject* __Pyx_PyObject_Call(PyObject *func, PyObject *arg, PyObject *kw);
#else
#define __Pyx_PyObject_Call(func, arg, kw) PyObject_Call(func, arg, kw)
#endif

#if CYTHON_COMPILING_IN_CPYTHON
static CYTHON_INLINE PyObject* __Pyx_PyObject_CallMethO(PyObject *func, PyObject *arg);
#endif

static CYTHON_INLINE PyObject* __Pyx_PyObject_CallOneArg(PyObject *func, PyObject *arg);

static CYTHON_INLINE void __Pyx_ErrRestore(PyObject *type, PyObject *value, PyObject *tb);
static CYTHON_INLINE void __Pyx_ErrFetch(PyObject **type, PyObject **value, PyObject **tb);

static void __Pyx_WriteUnraisable(const char *name, int clineno,

                                  int lineno, const char *filename,

#

                                  int full_traceback);

static void __Pyx_RaiseArgtupleInvalid(const char* func_name, int exact,
    Py_ssize_t num_min, Py_ssize_t num_max, Py_ssize_t num_found);

static void __Pyx_RaiseDoubleKeywordsError(const char* func_name, PyObject* kw_name);

static int __Pyx_ParseOptionalKeywords(PyObject *kwds, PyObject **argnames[], \
    PyObject *kwds2, PyObject *values[], Py_ssize_t num_pos_args, \
    const char* function_name);

static CYTHON_INLINE int __Pyx_ArgTypeTest(PyObject *obj, PyTypeObject *type, int none_allowed,
    const char *name, int exact);

static CYTHON_INLINE int  __Pyx_GetBufferAndValidate(Py_buffer* buf, PyObject* obj,
    __Pyx_TypeInfo* dtype, int flags, int nd, int cast, __Pyx_BufFmt_StackElem* stack);
static CYTHON_INLINE void __Pyx_SafeReleaseBuffer(Py_buffer* info);

#define __Pyx_BufPtrStrided1d(type, buf, i0, s0) (type)((char*)buf + i0 * s0)
static void __Pyx_Raise(PyObject *type, PyObject *value, PyObject *tb, PyObject *cause);

static CYTHON_INLINE void __Pyx_RaiseTooManyValuesError(Py_ssize_t expected);

static CYTHON_INLINE void __Pyx_RaiseNeedMoreValuesError(Py_ssize_t index);

static CYTHON_INLINE void __Pyx_RaiseNoneNotIterableError(void);

static CYTHON_INLINE int __Pyx_TypeTest(PyObject *obj, PyTypeObject *type);

typedef struct {

#

    int code_line;
    PyCodeObject* code_object;
} __Pyx_CodeObjectCacheEntry;
struct __Pyx_CodeObjectCache {

#

    int count;

#

    int max_count;
    __Pyx_CodeObjectCacheEntry* entries;
};
static struct __Pyx_CodeObjectCache __pyx_code_cache = {0,0,NULL};
static int __pyx_bisect_code_objects(__Pyx_CodeObjectCacheEntry* entries, int count, int code_line);
static PyCodeObject *__pyx_find_code_object(int code_line);
static void __pyx_insert_code_object(int code_line, PyCodeObject* code_object);

static void __Pyx_AddTraceback(const char *funcname, int c_line,

                               int py_line, const char *filename);

typedef struct {
  Py_ssize_t shape, strides, suboffsets;
} __Pyx_Buf_DimInfo;
typedef struct {
  size_t refcount;
  Py_buffer pybuffer;
} __Pyx_Buffer;
typedef struct {
  __Pyx_Buffer *rcbuffer;
  char *data;
  __Pyx_Buf_DimInfo diminfo[8];
} __Pyx_LocalBuf_ND;

#if PY_MAJOR_VERSION < 3
    static int __Pyx_GetBuffer(PyObject *obj, Py_buffer *view, int flags);
    static void __Pyx_ReleaseBuffer(Py_buffer *view);
#else
    #define __Pyx_GetBuffer PyObject_GetBuffer
    #define __Pyx_ReleaseBuffer PyBuffer_Release
#endif


static Py_ssize_t __Pyx_zeros[] = {0, 0, 0, 0, 0, 0, 0, 0};
static Py_ssize_t __Pyx_minusones[] = {-1, -1, -1, -1, -1, -1, -1, -1};

static PyObject *__Pyx_Import(PyObject *name, PyObject *from_list, int level);

static CYTHON_INLINE int __Pyx_PyInt_As_int(PyObject *);

static CYTHON_INLINE unsigned int __Pyx_PyInt_As_unsigned_int(PyObject *);

static CYTHON_INLINE PyObject* __Pyx_PyInt_From_int(int value);

#if CYTHON_CCOMPLEX
  #ifdef __cplusplus
    #define __Pyx_CREAL(z) ((z).real())
    #define __Pyx_CIMAG(z) ((z).imag())
  #else
    #define __Pyx_CREAL(z) (__real__(z))
    #define __Pyx_CIMAG(z) (__imag__(z))
  #endif
#else
    #define __Pyx_CREAL(z) ((z).real)
    #define __Pyx_CIMAG(z) ((z).imag)
#endif
#if (defined(_WIN32) || defined(__clang__)) && defined(__cplusplus) && CYTHON_CCOMPLEX
    #define __Pyx_SET_CREAL(z,x) ((z).real(x))
    #define __Pyx_SET_CIMAG(z,y) ((z).imag(y))
#else
    #define __Pyx_SET_CREAL(z,x) __Pyx_CREAL(z) = (x)
    #define __Pyx_SET_CIMAG(z,y) __Pyx_CIMAG(z) = (y)
#endif

static CYTHON_INLINE __pyx_t_float_complex __pyx_t_float_complex_from_parts(float, float);

#if CYTHON_CCOMPLEX
    #define __Pyx_c_eqf(a, b)   ((a)==(b))
    #define __Pyx_c_sumf(a, b)  ((a)+(b))
    #define __Pyx_c_difff(a, b) ((a)-(b))
    #define __Pyx_c_prodf(a, b) ((a)*(b))
    #define __Pyx_c_quotf(a, b) ((a)/(b))
    #define __Pyx_c_negf(a)     (-(a))
  #ifdef __cplusplus
    #define __Pyx_c_is_zerof(z) ((z)==(float)0)
    #define __Pyx_c_conjf(z)    (::std::conj(z))
    #if 1
        #define __Pyx_c_absf(z)     (::std::abs(z))
        #define __Pyx_c_powf(a, b)  (::std::pow(a, b))
    #endif
  #else
    #define __Pyx_c_is_zerof(z) ((z)==0)
    #define __Pyx_c_conjf(z)    (conjf(z))
    #if 1
        #define __Pyx_c_absf(z)     (cabsf(z))
        #define __Pyx_c_powf(a, b)  (cpowf(a, b))
    #endif
 #endif
#else
    static CYTHON_INLINE int __Pyx_c_eqf(__pyx_t_float_complex, __pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_sumf(__pyx_t_float_complex, __pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_difff(__pyx_t_float_complex, __pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_prodf(__pyx_t_float_complex, __pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_quotf(__pyx_t_float_complex, __pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_negf(__pyx_t_float_complex);
    static CYTHON_INLINE int __Pyx_c_is_zerof(__pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_conjf(__pyx_t_float_complex);
    #if 1
        static CYTHON_INLINE float __Pyx_c_absf(__pyx_t_float_complex);
        static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_powf(__pyx_t_float_complex, __pyx_t_float_complex);
    #endif
#endif

static CYTHON_INLINE __pyx_t_double_complex __pyx_t_double_complex_from_parts(double, double);

#if CYTHON_CCOMPLEX
    #define __Pyx_c_eq(a, b)   ((a)==(b))
    #define __Pyx_c_sum(a, b)  ((a)+(b))
    #define __Pyx_c_diff(a, b) ((a)-(b))
    #define __Pyx_c_prod(a, b) ((a)*(b))
    #define __Pyx_c_quot(a, b) ((a)/(b))
    #define __Pyx_c_neg(a)     (-(a))
  #ifdef __cplusplus
    #define __Pyx_c_is_zero(z) ((z)==(double)0)
    #define __Pyx_c_conj(z)    (::std::conj(z))
    #if 1
        #define __Pyx_c_abs(z)     (::std::abs(z))
        #define __Pyx_c_pow(a, b)  (::std::pow(a, b))
    #endif
  #else
    #define __Pyx_c_is_zero(z) ((z)==0)
    #define __Pyx_c_conj(z)    (conj(z))
    #if 1
        #define __Pyx_c_abs(z)     (cabs(z))
        #define __Pyx_c_pow(a, b)  (cpow(a, b))
    #endif
 #endif
#else
    static CYTHON_INLINE int __Pyx_c_eq(__pyx_t_double_complex, __pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_sum(__pyx_t_double_complex, __pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_diff(__pyx_t_double_complex, __pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_prod(__pyx_t_double_complex, __pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_quot(__pyx_t_double_complex, __pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_neg(__pyx_t_double_complex);
    static CYTHON_INLINE int __Pyx_c_is_zero(__pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_conj(__pyx_t_double_complex);
    #if 1
        static CYTHON_INLINE double __Pyx_c_abs(__pyx_t_double_complex);
        static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_pow(__pyx_t_double_complex, __pyx_t_double_complex);
    #endif
#endif

static CYTHON_INLINE PyObject* __Pyx_PyInt_From_long(long value);

static CYTHON_INLINE long __Pyx_PyInt_As_long(PyObject *);

static int __Pyx_check_binary_version(void);

#if !defined(__Pyx_PyIdentifier_FromString)
#if PY_MAJOR_VERSION < 3
  #define __Pyx_PyIdentifier_FromString(s) PyString_FromString(s)
#else
  #define __Pyx_PyIdentifier_FromString(s) PyUnicode_FromString(s)
#endif
#endif

static PyObject *__Pyx_ImportModule(const char *name);

static PyTypeObject *__Pyx_ImportType(const char *module_name, const char *class_name, size_t size, int strict);

static int __Pyx_InitStrings(__Pyx_StringTabEntry *t);

#

Module declarations from 'cpython.buffer'

#

Module declarations from 'cpython.ref'

#

Module declarations from 'libc.string'

#

Module declarations from 'libc.stdio'

#

Module declarations from 'cpython.object'

#

Module declarations from 'builtin'

#

Module declarations from 'cpython.type'

static PyTypeObject *__pyx_ptype_7cpython_4type_type = 0;

#

Module declarations from 'libc.stdlib'

#

Module declarations from 'numpy'

#

Module declarations from 'numpy'

static PyTypeObject *__pyx_ptype_5numpy_dtype = 0;
static PyTypeObject *__pyx_ptype_5numpy_flatiter = 0;
static PyTypeObject *__pyx_ptype_5numpy_broadcast = 0;
static PyTypeObject *__pyx_ptype_5numpy_ndarray = 0;
static PyTypeObject *__pyx_ptype_5numpy_ufunc = 0;
static CYTHON_INLINE char *__pyx_f_5numpy__util_dtypestring(PyArray_Descr *, char *, char *, int *); /*proto*/

#

Module declarations from 'cython'

#

Module declarations from 'src.particles.interpolation.interpolation'

static int __pyx_v_3src_9particles_13interpolation_13interpolation_exponential_interpolation;
static __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_f_3src_9particles_13interpolation_13interpolation_energy(__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t, __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t); /*proto*/
static __Pyx_TypeInfo __Pyx_TypeInfo_nn___pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t = { "DTYPE_t", NULL, sizeof(__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t), { 0 }, 0, 'R', 0, 0 };
#define __Pyx_MODULE_NAME "src.particles.interpolation.interpolation"

#

Implementation of 'src.particles.interpolation.interpolation'

int __pyx_module_is_main_src__particles__interpolation__interpolation = 0;

#

static PyObject *__pyx_builtin_ValueError;
static PyObject *__pyx_builtin_range;
static PyObject *__pyx_builtin_RuntimeError;
static PyObject *__pyx_pf_3src_9particles_13interpolation_13interpolation_distribution_interpolation(CYTHON_UNUSED PyObject *__pyx_self, PyArrayObject *__pyx_v_grid, PyArrayObject *__pyx_v_distribution, __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_p, __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_m, int __pyx_v_eta, __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_MIN_MOMENTUM, __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_MOMENTUM_STEP); /* proto */
static int __pyx_pf_5numpy_7ndarray___getbuffer__(PyArrayObject *__pyx_v_self, Py_buffer *__pyx_v_info, int __pyx_v_flags); /* proto */
static void __pyx_pf_5numpy_7ndarray_2__releasebuffer__(PyArrayObject *__pyx_v_self, Py_buffer *__pyx_v_info); /* proto */
static char __pyx_k_B[] = "B";
static char __pyx_k_H[] = "H";
static char __pyx_k_I[] = "I";
static char __pyx_k_L[] = "L";
static char __pyx_k_O[] = "O";
static char __pyx_k_Q[] = "Q";
static char __pyx_k_b[] = "b";
static char __pyx_k_d[] = "d";
static char __pyx_k_f[] = "f";
static char __pyx_k_g[] = "g";
static char __pyx_k_h[] = "h";
static char __pyx_k_i[] = "i";
static char __pyx_k_l[] = "l";
static char __pyx_k_m[] = "m";
static char __pyx_k_p[] = "p";
static char __pyx_k_q[] = "q";
static char __pyx_k_Zd[] = "Zd";
static char __pyx_k_Zf[] = "Zf";
static char __pyx_k_Zg[] = "Zg";
static char __pyx_k_E_p[] = "E_p";
static char __pyx_k_eta[] = "eta";
static char __pyx_k_exp[] = "exp";
static char __pyx_k_log[] = "log";
static char __pyx_k_grid[] = "grid";
static char __pyx_k_main[] = "__main__";
static char __pyx_k_math[] = "math";
static char __pyx_k_sqrt[] = "sqrt";
static char __pyx_k_test[] = "__test__";
static char __pyx_k_DTYPE[] = "DTYPE";
static char __pyx_k_E_low[] = "E_low";
static char __pyx_k_float[] = "float_";
static char __pyx_k_g_low[] = "g_low";
static char __pyx_k_i_low[] = "i_low";
static char __pyx_k_index[] = "index";
static char __pyx_k_numpy[] = "numpy";
static char __pyx_k_p_low[] = "p_low";
static char __pyx_k_range[] = "range";
static char __pyx_k_E_high[] = "E_high";
static char __pyx_k_g_high[] = "g_high";
static char __pyx_k_i_high[] = "i_high";
static char __pyx_k_import[] = "__import__";
static char __pyx_k_p_high[] = "p_high";
static char __pyx_k_remnant[] = "remnant";
static char __pyx_k_ValueError[] = "ValueError";
static char __pyx_k_MIN_MOMENTUM[] = "MIN_MOMENTUM";
static char __pyx_k_RuntimeError[] = "RuntimeError";
static char __pyx_k_distribution[] = "distribution";
static char __pyx_k_MOMENTUM_STEP[] = "MOMENTUM_STEP";
static char __pyx_k_distribution_interpolation[] = "distribution_interpolation";
static char __pyx_k_ndarray_is_not_C_contiguous[] = "ndarray is not C contiguous";
static char __pyx_k_Users_b_Repositories_pyBBN_src[] = "/Users/b/Repositories/pyBBN/src/particles/interpolation/interpolation.pyx";
static char __pyx_k_Cython_version_of_the_distribut[] = " Cython version of the distribution interpolation.\n\n    TODO: currently massless only\n    TODO: no significant performance gain ";
static char __pyx_k_unknown_dtype_code_in_numpy_pxd[] = "unknown dtype code in numpy.pxd (%d)";
static char __pyx_k_Format_string_allocated_too_shor[] = "Format string allocated too short, see comment in numpy.pxd";
static char __pyx_k_Non_native_byte_order_not_suppor[] = "Non-native byte order not supported";
static char __pyx_k_ndarray_is_not_Fortran_contiguou[] = "ndarray is not Fortran contiguous";
static char __pyx_k_src_particles_interpolation_inte[] = "src.particles.interpolation.interpolation";
static char __pyx_k_Format_string_allocated_too_shor_2[] = "Format string allocated too short.";
static PyObject *__pyx_n_s_DTYPE;
static PyObject *__pyx_n_s_E_high;
static PyObject *__pyx_n_s_E_low;
static PyObject *__pyx_n_s_E_p;
static PyObject *__pyx_kp_u_Format_string_allocated_too_shor;
static PyObject *__pyx_kp_u_Format_string_allocated_too_shor_2;
static PyObject *__pyx_n_s_MIN_MOMENTUM;
static PyObject *__pyx_n_s_MOMENTUM_STEP;
static PyObject *__pyx_kp_u_Non_native_byte_order_not_suppor;
static PyObject *__pyx_n_s_RuntimeError;
static PyObject *__pyx_kp_s_Users_b_Repositories_pyBBN_src;
static PyObject *__pyx_n_s_ValueError;
static PyObject *__pyx_n_s_distribution;
static PyObject *__pyx_n_s_distribution_interpolation;
static PyObject *__pyx_n_s_eta;
static PyObject *__pyx_n_s_exp;
static PyObject *__pyx_n_s_float;
static PyObject *__pyx_n_s_g;
static PyObject *__pyx_n_s_g_high;
static PyObject *__pyx_n_s_g_low;
static PyObject *__pyx_n_s_grid;
static PyObject *__pyx_n_s_i;
static PyObject *__pyx_n_s_i_high;
static PyObject *__pyx_n_s_i_low;
static PyObject *__pyx_n_s_import;
static PyObject *__pyx_n_s_index;
static PyObject *__pyx_n_s_log;
static PyObject *__pyx_n_s_m;
static PyObject *__pyx_n_s_main;
static PyObject *__pyx_n_s_math;
static PyObject *__pyx_kp_u_ndarray_is_not_C_contiguous;
static PyObject *__pyx_kp_u_ndarray_is_not_Fortran_contiguou;
static PyObject *__pyx_n_s_numpy;
static PyObject *__pyx_n_s_p;
static PyObject *__pyx_n_s_p_high;
static PyObject *__pyx_n_s_p_low;
static PyObject *__pyx_n_s_range;
static PyObject *__pyx_n_s_remnant;
static PyObject *__pyx_n_s_sqrt;
static PyObject *__pyx_n_s_src_particles_interpolation_inte;
static PyObject *__pyx_n_s_test;
static PyObject *__pyx_kp_u_unknown_dtype_code_in_numpy_pxd;
static PyObject *__pyx_float_1_;
static PyObject *__pyx_int_1;
static PyObject *__pyx_tuple_;
static PyObject *__pyx_tuple__2;
static PyObject *__pyx_tuple__3;
static PyObject *__pyx_tuple__4;
static PyObject *__pyx_tuple__5;
static PyObject *__pyx_tuple__6;
static PyObject *__pyx_tuple__7;
static PyObject *__pyx_codeobj__8;

#

"src/particles/interpolation/interpolation.pyx":23 cdef int exponential_interpolation = True

cdef DTYPE_t energy(DTYPE_t p, DTYPE_t m): <<<<<<<<<<<<<<

return numpy.sqrt(p2 + m2)

static __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_f_3src_9particles_13interpolation_13interpolation_energy(__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_p, __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_m) {
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_r;
  __Pyx_RefNannyDeclarations
  PyObject *__pyx_t_1 = NULL;
  PyObject *__pyx_t_2 = NULL;
  PyObject *__pyx_t_3 = NULL;
  PyObject *__pyx_t_4 = NULL;
  PyObject *__pyx_t_5 = NULL;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_t_6;

  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;

#

  int __pyx_clineno = 0;
  __Pyx_RefNannySetupContext("energy", 0);

#

"src/particles/interpolation/interpolation.pyx":24

cdef DTYPE_t energy(DTYPE_t p, DTYPE_t m): return numpy.sqrt(p2 + m2) # <<<<<<<<<<<<<<

  __pyx_t_2 = __Pyx_GetModuleGlobalName(__pyx_n_s_numpy); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __pyx_t_3 = __Pyx_PyObject_GetAttrStr(__pyx_t_2, __pyx_n_s_sqrt); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_3);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __pyx_t_2 = PyFloat_FromDouble((pow(__pyx_v_p, 2.0) + pow(__pyx_v_m, 2.0))); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __pyx_t_4 = NULL;
  if (CYTHON_COMPILING_IN_CPYTHON && unlikely(PyMethod_Check(__pyx_t_3))) {
    __pyx_t_4 = PyMethod_GET_SELF(__pyx_t_3);
    if (likely(__pyx_t_4)) {
      PyObject* function = PyMethod_GET_FUNCTION(__pyx_t_3);
      __Pyx_INCREF(__pyx_t_4);
      __Pyx_INCREF(function);
      __Pyx_DECREF_SET(__pyx_t_3, function);
    }
  }
  if (!__pyx_t_4) {
    __pyx_t_1 = __Pyx_PyObject_CallOneArg(__pyx_t_3, __pyx_t_2); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
    __Pyx_GOTREF(__pyx_t_1);
  } else {
    __pyx_t_5 = PyTuple_New(1+1); if (unlikely(!__pyx_t_5)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_5);
    PyTuple_SET_ITEM(__pyx_t_5, 0, __pyx_t_4); __Pyx_GIVEREF(__pyx_t_4); __pyx_t_4 = NULL;
    PyTuple_SET_ITEM(__pyx_t_5, 0+1, __pyx_t_2);
    __Pyx_GIVEREF(__pyx_t_2);
    __pyx_t_2 = 0;
    __pyx_t_1 = __Pyx_PyObject_Call(__pyx_t_3, __pyx_t_5, NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __Pyx_DECREF(__pyx_t_5); __pyx_t_5 = 0;
  }
  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
  __pyx_t_6 = __pyx_PyFloat_AsDouble(__pyx_t_1); if (unlikely((__pyx_t_6 == (npy_double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_r = __pyx_t_6;
  goto __pyx_L0;

#

"src/particles/interpolation/interpolation.pyx":23 cdef int exponential_interpolation = True

cdef DTYPE_t energy(DTYPE_t p, DTYPE_t m): <<<<<<<<<<<<<<

return numpy.sqrt(p2 + m2)

#

function exit code

  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_1);
  __Pyx_XDECREF(__pyx_t_2);
  __Pyx_XDECREF(__pyx_t_3);
  __Pyx_XDECREF(__pyx_t_4);
  __Pyx_XDECREF(__pyx_t_5);
  __Pyx_WriteUnraisable("src.particles.interpolation.interpolation.energy", __pyx_clineno, __pyx_lineno, __pyx_filename, 0);
  __pyx_r = 0;
  __pyx_L0:;
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

#

"src/particles/interpolation/interpolation.pyx":30 cython .wraparound(False)
cython .cdivision(True)
def distribution_interpolation(numpy.ndarray[DTYPE_t, ndim=1] grid, # <<<<<<<<<<<<<< numpy.ndarray[DTYPE_t, ndim=1] distribution, DTYPE_t p,

#

Python wrapper

static PyObject *__pyx_pw_3src_9particles_13interpolation_13interpolation_1distribution_interpolation(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/
static PyMethodDef __pyx_mdef_3src_9particles_13interpolation_13interpolation_1distribution_interpolation = {"distribution_interpolation", (PyCFunction)__pyx_pw_3src_9particles_13interpolation_13interpolation_1distribution_interpolation, METH_VARARGS|METH_KEYWORDS, 0};
static PyObject *__pyx_pw_3src_9particles_13interpolation_13interpolation_1distribution_interpolation(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) {
  PyArrayObject *__pyx_v_grid = 0;
  PyArrayObject *__pyx_v_distribution = 0;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_p;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_m;

  int __pyx_v_eta;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_MIN_MOMENTUM;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_MOMENTUM_STEP;

#

  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;

#

  int __pyx_clineno = 0;
  PyObject *__pyx_r = 0;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("distribution_interpolation (wrapper)", 0);
  {
    static PyObject **__pyx_pyargnames[] = {&__pyx_n_s_grid,&__pyx_n_s_distribution,&__pyx_n_s_p,&__pyx_n_s_m,&__pyx_n_s_eta,&__pyx_n_s_MIN_MOMENTUM,&__pyx_n_s_MOMENTUM_STEP,0};
    PyObject* values[7] = {0,0,0,0,0,0,0};
    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_grid)) != 0)) kw_args--;
        else goto __pyx_L5_argtuple_error;
        case  1:
        if (likely((values[1] = PyDict_GetItem(__pyx_kwds, __pyx_n_s_distribution)) != 0)) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("distribution_interpolation", 0, 4, 7, 1); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        }
        case  2:
        if (likely((values[2] = PyDict_GetItem(__pyx_kwds, __pyx_n_s_p)) != 0)) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("distribution_interpolation", 0, 4, 7, 2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        }
        case  3:
        if (likely((values[3] = PyDict_GetItem(__pyx_kwds, __pyx_n_s_m)) != 0)) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("distribution_interpolation", 0, 4, 7, 3); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
        }
        case  4:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s_eta);
          if (value) { values[4] = value; kw_args--; }
        }
        case  5:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s_MIN_MOMENTUM);
          if (value) { values[5] = value; kw_args--; }
        }
        case  6:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s_MOMENTUM_STEP);
          if (value) { values[6] = value; kw_args--; }
        }
      }
      if (unlikely(kw_args > 0)) {
        if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_pyargnames, 0, values, pos_args, "distribution_interpolation") < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __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_grid = ((PyArrayObject *)values[0]);
    __pyx_v_distribution = ((PyArrayObject *)values[1]);
    __pyx_v_p = __pyx_PyFloat_AsDouble(values[2]); if (unlikely((__pyx_v_p == (npy_double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
    __pyx_v_m = __pyx_PyFloat_AsDouble(values[3]); if (unlikely((__pyx_v_m == (npy_double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
    if (values[4]) {
      __pyx_v_eta = __Pyx_PyInt_As_int(values[4]); if (unlikely((__pyx_v_eta == (int)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
    } else {
      __pyx_v_eta = ((int)1);
    }
    if (values[5]) {
      __pyx_v_MIN_MOMENTUM = __pyx_PyFloat_AsDouble(values[5]); if (unlikely((__pyx_v_MIN_MOMENTUM == (npy_double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
    } else {
      __pyx_v_MIN_MOMENTUM = ((__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t)0.);
    }
    if (values[6]) {
      __pyx_v_MOMENTUM_STEP = __pyx_PyFloat_AsDouble(values[6]); if (unlikely((__pyx_v_MOMENTUM_STEP == (npy_double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
    } else {
      __pyx_v_MOMENTUM_STEP = ((__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t)0.);
    }
  }
  goto __pyx_L4_argument_unpacking_done;
  __pyx_L5_argtuple_error:;
  __Pyx_RaiseArgtupleInvalid("distribution_interpolation", 0, 4, 7, PyTuple_GET_SIZE(__pyx_args)); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
  __pyx_L3_error:;
  __Pyx_AddTraceback("src.particles.interpolation.interpolation.distribution_interpolation", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __Pyx_RefNannyFinishContext();
  return NULL;
  __pyx_L4_argument_unpacking_done:;
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_grid), __pyx_ptype_5numpy_ndarray, 1, "grid", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_distribution), __pyx_ptype_5numpy_ndarray, 1, "distribution", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_r = __pyx_pf_3src_9particles_13interpolation_13interpolation_distribution_interpolation(__pyx_self, __pyx_v_grid, __pyx_v_distribution, __pyx_v_p, __pyx_v_m, __pyx_v_eta, __pyx_v_MIN_MOMENTUM, __pyx_v_MOMENTUM_STEP);

#

function exit code

  goto __pyx_L0;
  __pyx_L1_error:;
  __pyx_r = NULL;
  __pyx_L0:;
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

static PyObject *__pyx_pf_3src_9particles_13interpolation_13interpolation_distribution_interpolation(CYTHON_UNUSED PyObject *__pyx_self, PyArrayObject *__pyx_v_grid, PyArrayObject *__pyx_v_distribution, __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_p, __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_m, int __pyx_v_eta, __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_MIN_MOMENTUM, __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_MOMENTUM_STEP) {
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_p_low;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_p_high;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_remnant;
  CYTHON_UNUSED unsigned int __pyx_v_i;
  unsigned int __pyx_v_i_low;
  unsigned int __pyx_v_i_high;
  PyObject *__pyx_v_index = NULL;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_E_p;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_E_low;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_E_high;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_g_high;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_g_low;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_v_g;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_distribution;
  __Pyx_Buffer __pyx_pybuffer_distribution;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_grid;
  __Pyx_Buffer __pyx_pybuffer_grid;
  PyObject *__pyx_r = NULL;
  __Pyx_RefNannyDeclarations
  PyObject *__pyx_t_1 = NULL;
  PyObject *__pyx_t_2 = NULL;

#

  int __pyx_t_3;
  Py_ssize_t __pyx_t_4;
  Py_ssize_t __pyx_t_5;
  unsigned int __pyx_t_6;
  unsigned int __pyx_t_7;
  unsigned int __pyx_t_8;
  unsigned int __pyx_t_9;
  PyObject *__pyx_t_10 = NULL;
  PyObject *__pyx_t_11 = NULL;
  PyObject *__pyx_t_12 = NULL;
  __pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t __pyx_t_13;
  unsigned int __pyx_t_14;
  unsigned int __pyx_t_15;

#

  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;

#

  int __pyx_clineno = 0;
  __Pyx_RefNannySetupContext("distribution_interpolation", 0);
  __pyx_pybuffer_grid.pybuffer.buf = NULL;
  __pyx_pybuffer_grid.refcount = 0;
  __pyx_pybuffernd_grid.data = NULL;
  __pyx_pybuffernd_grid.rcbuffer = &__pyx_pybuffer_grid;
  __pyx_pybuffer_distribution.pybuffer.buf = NULL;
  __pyx_pybuffer_distribution.refcount = 0;
  __pyx_pybuffernd_distribution.data = NULL;
  __pyx_pybuffernd_distribution.rcbuffer = &__pyx_pybuffer_distribution;
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_grid.rcbuffer->pybuffer, (PyObject*)__pyx_v_grid, &__Pyx_TypeInfo_nn___pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }
  __pyx_pybuffernd_grid.diminfo[0].strides = __pyx_pybuffernd_grid.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_grid.diminfo[0].shape = __pyx_pybuffernd_grid.rcbuffer->pybuffer.shape[0];
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_distribution.rcbuffer->pybuffer, (PyObject*)__pyx_v_distribution, &__Pyx_TypeInfo_nn___pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }
  __pyx_pybuffernd_distribution.diminfo[0].strides = __pyx_pybuffernd_distribution.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_distribution.diminfo[0].shape = __pyx_pybuffernd_distribution.rcbuffer->pybuffer.shape[0];

#

"src/particles/interpolation/interpolation.pyx":37 DTYPE_t MIN_MOMENTUM=0., DTYPE_t MOMENTUM_STEP=0.):

cdef DTYPE_t p_low = -1, p_high = -1, remnant # <<<<<<<<<<<<<<

cdef unsigned int i = 0, i_low, i_high

  __pyx_v_p_low = -1.0;
  __pyx_v_p_high = -1.0;

#

"src/particles/interpolation/interpolation.pyx":39 cdef DTYPE_t p_low = -1, p_high = -1, remnant

cdef unsigned int i = 0, i_low, i_high # <<<<<<<<<<<<<<

remnant = (p - MIN_MOMENTUM) % MOMENTUM_STEP

  __pyx_v_i = 0;

#

"src/particles/interpolation/interpolation.pyx":41 cdef unsigned int i = 0, i_low, i_high

remnant = (p - MIN_MOMENTUM) % MOMENTUM_STEP # <<<<<<<<<<<<<< index = int((p - MIN_MOMENTUM) / MOMENTUM_STEP - remnant)

  __pyx_v_remnant = fmod((__pyx_v_p - __pyx_v_MIN_MOMENTUM), __pyx_v_MOMENTUM_STEP);

#

"src/particles/interpolation/interpolation.pyx":42

remnant = (p - MIN_MOMENTUM) % MOMENTUM_STEP index = int((p - MIN_MOMENTUM) / MOMENTUM_STEP - remnant) # <<<<<<<<<<<<<<

if remnant == 0:

  __pyx_t_1 = PyFloat_FromDouble((((__pyx_v_p - __pyx_v_MIN_MOMENTUM) / __pyx_v_MOMENTUM_STEP) - __pyx_v_remnant)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 42; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_2 = PyTuple_New(1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 42; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  PyTuple_SET_ITEM(__pyx_t_2, 0, __pyx_t_1);
  __Pyx_GIVEREF(__pyx_t_1);
  __pyx_t_1 = 0;
  __pyx_t_1 = __Pyx_PyObject_Call(((PyObject *)((PyObject*)(&PyInt_Type))), __pyx_t_2, NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 42; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __pyx_v_index = __pyx_t_1;
  __pyx_t_1 = 0;

#

"src/particles/interpolation/interpolation.pyx":44 index = int((p - MIN_MOMENTUM) / MOMENTUM_STEP - remnant)

if remnant == 0: <<<<<<<<<<<<<<

return distribution[index]

  __pyx_t_3 = ((__pyx_v_remnant == 0.0) != 0);
  if (__pyx_t_3) {

#

"src/particles/interpolation/interpolation.pyx":45

if remnant == 0: return distribution[index] # <<<<<<<<<<<<<<

if index >= len(grid) - 1:

    __Pyx_XDECREF(__pyx_r);
    __pyx_t_1 = PyObject_GetItem(((PyObject *)__pyx_v_distribution), __pyx_v_index); if (unlikely(__pyx_t_1 == NULL)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 45; __pyx_clineno = __LINE__; goto __pyx_L1_error;};
    __Pyx_GOTREF(__pyx_t_1);
    __pyx_r = __pyx_t_1;
    __pyx_t_1 = 0;
    goto __pyx_L0;
  }

#

"src/particles/interpolation/interpolation.pyx":47 return distribution[index]

if index >= len(grid) - 1: <<<<<<<<<<<<<<

return distribution[len(grid) - 1]

  __pyx_t_4 = PyObject_Length(((PyObject *)__pyx_v_grid)); if (unlikely(__pyx_t_4 == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_1 = PyInt_FromSsize_t((__pyx_t_4 - 1)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_t_2 = PyObject_RichCompare(__pyx_v_index, __pyx_t_1, Py_GE); __Pyx_XGOTREF(__pyx_t_2); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
  __pyx_t_3 = __Pyx_PyObject_IsTrue(__pyx_t_2); if (unlikely(__pyx_t_3 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  if (__pyx_t_3) {

#

"src/particles/interpolation/interpolation.pyx":48

if index >= len(grid) - 1: return distribution[len(grid) - 1] # <<<<<<<<<<<<<<

Determine the closest grid points

    __Pyx_XDECREF(__pyx_r);
    __pyx_t_4 = PyObject_Length(((PyObject *)__pyx_v_grid)); if (unlikely(__pyx_t_4 == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 48; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __pyx_t_5 = (__pyx_t_4 - 1);
    __pyx_t_2 = PyFloat_FromDouble((*__Pyx_BufPtrStrided1d(__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t *, __pyx_pybuffernd_distribution.rcbuffer->pybuffer.buf, __pyx_t_5, __pyx_pybuffernd_distribution.diminfo[0].strides))); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 48; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_2);
    __pyx_r = __pyx_t_2;
    __pyx_t_2 = 0;
    goto __pyx_L0;
  }

#

"src/particles/interpolation/interpolation.pyx":52

Determine the closest grid points

i_low = index # <<<<<<<<<<<<<< p_low = grid[i_low]

  __pyx_t_6 = __Pyx_PyInt_As_unsigned_int(__pyx_v_index); if (unlikely((__pyx_t_6 == (unsigned int)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 52; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_v_i_low = __pyx_t_6;

#

"src/particles/interpolation/interpolation.pyx":53

i_low = index p_low = grid[i_low] # <<<<<<<<<<<<<<

i_high = index + 1

  __pyx_t_6 = __pyx_v_i_low;
  __pyx_v_p_low = (*__Pyx_BufPtrStrided1d(__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t *, __pyx_pybuffernd_grid.rcbuffer->pybuffer.buf, __pyx_t_6, __pyx_pybuffernd_grid.diminfo[0].strides));

#

"src/particles/interpolation/interpolation.pyx":55 p_low = grid[i_low]

i_high = index + 1 # <<<<<<<<<<<<<< p_high = grid[i_high]

  __pyx_t_2 = PyNumber_Add(__pyx_v_index, __pyx_int_1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 55; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_2);
  __pyx_t_7 = __Pyx_PyInt_As_unsigned_int(__pyx_t_2); if (unlikely((__pyx_t_7 == (unsigned int)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 55; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __pyx_v_i_high = __pyx_t_7;

#

"src/particles/interpolation/interpolation.pyx":56

i_high = index + 1 p_high = grid[i_high] # <<<<<<<<<<<<<<

cdef DTYPE_t E_p, E_low, E_high, g_high, g_low, g

  __pyx_t_7 = __pyx_v_i_high;
  __pyx_v_p_high = (*__Pyx_BufPtrStrided1d(__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t *, __pyx_pybuffernd_grid.rcbuffer->pybuffer.buf, __pyx_t_7, __pyx_pybuffernd_grid.diminfo[0].strides));

#

"src/particles/interpolation/interpolation.pyx":60 cdef DTYPE_t E_p, E_low, E_high, g_high, g_low, g

if exponential_interpolation: <<<<<<<<<<<<<<

""" === Exponential interpolation === """ E_p = energy(p, m)

  __pyx_t_3 = (__pyx_v_3src_9particles_13interpolation_13interpolation_exponential_interpolation != 0);
  if (__pyx_t_3) {

#

"src/particles/interpolation/interpolation.pyx":62 if exponential_interpolation: """ === Exponential interpolation === """ E_p = energy(p, m) # <<<<<<<<<<<<<< E_low = energy(p_low, m) E_high = energy(p_high, m)

    __pyx_v_E_p = __pyx_f_3src_9particles_13interpolation_13interpolation_energy(__pyx_v_p, __pyx_v_m);

#

"src/particles/interpolation/interpolation.pyx":63 """ === Exponential interpolation === """ E_p = energy(p, m) E_low = energy(p_low, m) # <<<<<<<<<<<<<< E_high = energy(p_high, m)

    __pyx_v_E_low = __pyx_f_3src_9particles_13interpolation_13interpolation_energy(__pyx_v_p_low, __pyx_v_m);

#

"src/particles/interpolation/interpolation.pyx":64 E_p = energy(p, m) E_low = energy(p_low, m) E_high = energy(p_high, m) # <<<<<<<<<<<<<<

"""

    __pyx_v_E_high = __pyx_f_3src_9particles_13interpolation_13interpolation_energy(__pyx_v_p_high, __pyx_v_m);

#

"src/particles/interpolation/interpolation.pyx":73 """

g_high = distribution[i_high]             # <<<<<<<<<<<<<<
g_low = distribution[i_low]

    __pyx_t_8 = __pyx_v_i_high;
    __pyx_v_g_high = (*__Pyx_BufPtrStrided1d(__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t *, __pyx_pybuffernd_distribution.rcbuffer->pybuffer.buf, __pyx_t_8, __pyx_pybuffernd_distribution.diminfo[0].strides));

#

"src/particles/interpolation/interpolation.pyx":74

g_high = distribution[i_high]
g_low = distribution[i_low]             # <<<<<<<<<<<<<<

if g_high > 0:

    __pyx_t_9 = __pyx_v_i_low;
    __pyx_v_g_low = (*__Pyx_BufPtrStrided1d(__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t *, __pyx_pybuffernd_distribution.rcbuffer->pybuffer.buf, __pyx_t_9, __pyx_pybuffernd_distribution.diminfo[0].strides));

#

"src/particles/interpolation/interpolation.pyx":76 g_low = distribution[i_low]

if g_high > 0
:   # <<<<<<<<<<<<<<

    g_high = (1. / g_high - 1.)
    if g_high > 0:

    __pyx_t_3 = ((__pyx_v_g_high > 0.0) != 0);
    if (__pyx_t_3) {

#

"src/particles/interpolation/interpolation.pyx":77

if g_high > 0: g_high = (1. / g_high - 1.) # <<<<<<<<<<<<<< if g_high > 0: g_high = math.log(g_high)

      __pyx_v_g_high = ((((__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t)1.) / __pyx_v_g_high) - 1.);

#

"src/particles/interpolation/interpolation.pyx":78 if g_high > 0: g_high = (1. / g_high - 1.) if g_high > 0 : # <<<<<<<<<<<<<<

      g_high = math.log(g_high)

      __pyx_t_3 = ((__pyx_v_g_high > 0.0) != 0);
      if (__pyx_t_3) {

#

"src/particles/interpolation/interpolation.pyx":79 g_high = (1. / g_high - 1.) if g_high > 0: g_high = math.log(g_high) # <<<<<<<<<<<<<<

if g_low > 0:

        __pyx_t_1 = __Pyx_GetModuleGlobalName(__pyx_n_s_math); 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_t_10 = __Pyx_PyObject_GetAttrStr(__pyx_t_1, __pyx_n_s_log); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_10);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyFloat_FromDouble(__pyx_v_g_high); 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_t_11 = NULL;
        if (CYTHON_COMPILING_IN_CPYTHON && unlikely(PyMethod_Check(__pyx_t_10))) {
          __pyx_t_11 = PyMethod_GET_SELF(__pyx_t_10);
          if (likely(__pyx_t_11)) {
            PyObject* function = PyMethod_GET_FUNCTION(__pyx_t_10);
            __Pyx_INCREF(__pyx_t_11);
            __Pyx_INCREF(function);
            __Pyx_DECREF_SET(__pyx_t_10, function);
          }
        }
        if (!__pyx_t_11) {
          __pyx_t_2 = __Pyx_PyObject_CallOneArg(__pyx_t_10, __pyx_t_1); if (unlikely(!__pyx_t_2)) {__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_GOTREF(__pyx_t_2);
        } else {
          __pyx_t_12 = PyTuple_New(1+1); if (unlikely(!__pyx_t_12)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_12);
          PyTuple_SET_ITEM(__pyx_t_12, 0, __pyx_t_11); __Pyx_GIVEREF(__pyx_t_11); __pyx_t_11 = NULL;
          PyTuple_SET_ITEM(__pyx_t_12, 0+1, __pyx_t_1);
          __Pyx_GIVEREF(__pyx_t_1);
          __pyx_t_1 = 0;
          __pyx_t_2 = __Pyx_PyObject_Call(__pyx_t_10, __pyx_t_12, NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_2);
          __Pyx_DECREF(__pyx_t_12); __pyx_t_12 = 0;
        }
        __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
        __pyx_t_13 = __pyx_PyFloat_AsDouble(__pyx_t_2); if (unlikely((__pyx_t_13 == (npy_double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
        __pyx_v_g_high = __pyx_t_13;
        goto __pyx_L7;
      }
      __pyx_L7:;
      goto __pyx_L6;
    }
    __pyx_L6:;

#

"src/particles/interpolation/interpolation.pyx":81 g_high = math.log(g_high)

if g_low > 0
:   # <<<<<<<<<<<<<<

    g_low = (1. / g_low - 1.)
    if g_low > 0:

    __pyx_t_3 = ((__pyx_v_g_low > 0.0) != 0);
    if (__pyx_t_3) {

#

"src/particles/interpolation/interpolation.pyx":82

if g_low > 0: g_low = (1. / g_low - 1.) # <<<<<<<<<<<<<< if g_low > 0: g_low = math.log(g_low)

      __pyx_v_g_low = ((((__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t)1.) / __pyx_v_g_low) - 1.);

#

"src/particles/interpolation/interpolation.pyx":83 if g_low > 0: g_low = (1. / g_low - 1.) if g_low > 0 : # <<<<<<<<<<<<<<

      g_low = math.log(g_low)

      __pyx_t_3 = ((__pyx_v_g_low > 0.0) != 0);
      if (__pyx_t_3) {

#

"src/particles/interpolation/interpolation.pyx":84 g_low = (1. / g_low - 1.) if g_low > 0: g_low = math.log(g_low) # <<<<<<<<<<<<<<

g = ((E_p - E_low) * g_high + (E_high - E_p) * g_low) / (E_high - E_low)

        __pyx_t_10 = __Pyx_GetModuleGlobalName(__pyx_n_s_math); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 84; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_10);
        __pyx_t_12 = __Pyx_PyObject_GetAttrStr(__pyx_t_10, __pyx_n_s_log); if (unlikely(!__pyx_t_12)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 84; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_12);
        __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
        __pyx_t_10 = PyFloat_FromDouble(__pyx_v_g_low); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 84; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_10);
        __pyx_t_1 = NULL;
        if (CYTHON_COMPILING_IN_CPYTHON && unlikely(PyMethod_Check(__pyx_t_12))) {
          __pyx_t_1 = PyMethod_GET_SELF(__pyx_t_12);
          if (likely(__pyx_t_1)) {
            PyObject* function = PyMethod_GET_FUNCTION(__pyx_t_12);
            __Pyx_INCREF(__pyx_t_1);
            __Pyx_INCREF(function);
            __Pyx_DECREF_SET(__pyx_t_12, function);
          }
        }
        if (!__pyx_t_1) {
          __pyx_t_2 = __Pyx_PyObject_CallOneArg(__pyx_t_12, __pyx_t_10); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 84; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __Pyx_GOTREF(__pyx_t_2);
        } else {
          __pyx_t_11 = PyTuple_New(1+1); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 84; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_11);
          PyTuple_SET_ITEM(__pyx_t_11, 0, __pyx_t_1); __Pyx_GIVEREF(__pyx_t_1); __pyx_t_1 = NULL;
          PyTuple_SET_ITEM(__pyx_t_11, 0+1, __pyx_t_10);
          __Pyx_GIVEREF(__pyx_t_10);
          __pyx_t_10 = 0;
          __pyx_t_2 = __Pyx_PyObject_Call(__pyx_t_12, __pyx_t_11, NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 84; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_2);
          __Pyx_DECREF(__pyx_t_11); __pyx_t_11 = 0;
        }
        __Pyx_DECREF(__pyx_t_12); __pyx_t_12 = 0;
        __pyx_t_13 = __pyx_PyFloat_AsDouble(__pyx_t_2); if (unlikely((__pyx_t_13 == (npy_double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 84; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
        __pyx_v_g_low = __pyx_t_13;
        goto __pyx_L9;
      }
      __pyx_L9:;
      goto __pyx_L8;
    }
    __pyx_L8:;

#

"src/particles/interpolation/interpolation.pyx":86 g_low = math.log(g_low)

g = ((E_p - E_low) * g_high + (E_high - E_p) * g_low) / (E_high - E_low)             # <<<<<<<<<<<<<<

return 1. / (numpy.exp(g) + eta)

    __pyx_v_g = ((((__pyx_v_E_p - __pyx_v_E_low) * __pyx_v_g_high) + ((__pyx_v_E_high - __pyx_v_E_p) * __pyx_v_g_low)) / (__pyx_v_E_high - __pyx_v_E_low));

#

"src/particles/interpolation/interpolation.pyx":88 g = ((E_p - E_low) * g_high + (E_high - E_p) * g_low) / (E_high - E_low)

return 1. / (numpy.exp(g) + eta)             # <<<<<<<<<<<<<<

else:

    __Pyx_XDECREF(__pyx_r);
    __pyx_t_12 = __Pyx_GetModuleGlobalName(__pyx_n_s_numpy); if (unlikely(!__pyx_t_12)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_12);
    __pyx_t_11 = __Pyx_PyObject_GetAttrStr(__pyx_t_12, __pyx_n_s_exp); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_11);
    __Pyx_DECREF(__pyx_t_12); __pyx_t_12 = 0;
    __pyx_t_12 = PyFloat_FromDouble(__pyx_v_g); if (unlikely(!__pyx_t_12)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_12);
    __pyx_t_10 = NULL;
    if (CYTHON_COMPILING_IN_CPYTHON && unlikely(PyMethod_Check(__pyx_t_11))) {
      __pyx_t_10 = PyMethod_GET_SELF(__pyx_t_11);
      if (likely(__pyx_t_10)) {
        PyObject* function = PyMethod_GET_FUNCTION(__pyx_t_11);
        __Pyx_INCREF(__pyx_t_10);
        __Pyx_INCREF(function);
        __Pyx_DECREF_SET(__pyx_t_11, function);
      }
    }
    if (!__pyx_t_10) {
      __pyx_t_2 = __Pyx_PyObject_CallOneArg(__pyx_t_11, __pyx_t_12); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_12); __pyx_t_12 = 0;
      __Pyx_GOTREF(__pyx_t_2);
    } else {
      __pyx_t_1 = PyTuple_New(1+1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_1);
      PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_t_10); __Pyx_GIVEREF(__pyx_t_10); __pyx_t_10 = NULL;
      PyTuple_SET_ITEM(__pyx_t_1, 0+1, __pyx_t_12);
      __Pyx_GIVEREF(__pyx_t_12);
      __pyx_t_12 = 0;
      __pyx_t_2 = __Pyx_PyObject_Call(__pyx_t_11, __pyx_t_1, NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_2);
      __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
    }
    __Pyx_DECREF(__pyx_t_11); __pyx_t_11 = 0;
    __pyx_t_11 = __Pyx_PyInt_From_int(__pyx_v_eta); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_11);
    __pyx_t_1 = PyNumber_Add(__pyx_t_2, __pyx_t_11); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
    __Pyx_DECREF(__pyx_t_11); __pyx_t_11 = 0;
    __pyx_t_11 = __Pyx_PyNumber_Divide(__pyx_float_1_, __pyx_t_1); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_11);
    __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
    __pyx_r = __pyx_t_11;
    __pyx_t_11 = 0;
    goto __pyx_L0;
  }

#

else

#

"src/particles/interpolation/interpolation.pyx":92 else: """ === Linear interpolation === """ return ( # <<<<<<<<<<<<<< distribution[i_low] * (p_high - p) + distribution[i_high] * (p - p_low) ) / (p_high - p_low)

    __Pyx_XDECREF(__pyx_r);

#

"src/particles/interpolation/interpolation.pyx":93 """ === Linear interpolation === """ return ( distribution[i_low] * (p_high - p) + distribution[i_high] * (p - p_low) # <<<<<<<<<<<<<< ) / (p_high - p_low)

    __pyx_t_14 = __pyx_v_i_low;
    __pyx_t_15 = __pyx_v_i_high;

#

"src/particles/interpolation/interpolation.pyx":94 return ( distribution[i_low] * (p_high - p) + distribution[i_high] * (p - p_low) ) / (p_high - p_low) # <<<<<<<<<<<<<<

    __pyx_t_11 = PyFloat_FromDouble(((((*__Pyx_BufPtrStrided1d(__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t *, __pyx_pybuffernd_distribution.rcbuffer->pybuffer.buf, __pyx_t_14, __pyx_pybuffernd_distribution.diminfo[0].strides)) * (__pyx_v_p_high - __pyx_v_p)) + ((*__Pyx_BufPtrStrided1d(__pyx_t_3src_9particles_13interpolation_13interpolation_DTYPE_t *, __pyx_pybuffernd_distribution.rcbuffer->pybuffer.buf, __pyx_t_15, __pyx_pybuffernd_distribution.diminfo[0].strides)) * (__pyx_v_p - __pyx_v_p_low))) / (__pyx_v_p_high - __pyx_v_p_low))); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 94; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_11);
    __pyx_r = __pyx_t_11;
    __pyx_t_11 = 0;
    goto __pyx_L0;
  }

#

"src/particles/interpolation/interpolation.pyx":30 cython .wraparound(False)
cython .cdivision(True)
def distribution_interpolation(numpy.ndarray[DTYPE_t, ndim=1] grid, # <<<<<<<<<<<<<< numpy.ndarray[DTYPE_t, ndim=1] distribution, DTYPE_t p,

#

function exit code

  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_1);
  __Pyx_XDECREF(__pyx_t_2);
  __Pyx_XDECREF(__pyx_t_10);
  __Pyx_XDECREF(__pyx_t_11);
  __Pyx_XDECREF(__pyx_t_12);
  { PyObject *__pyx_type, *__pyx_value, *__pyx_tb;
    __Pyx_ErrFetch(&__pyx_type, &__pyx_value, &__pyx_tb);
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_distribution.rcbuffer->pybuffer);
    __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_grid.rcbuffer->pybuffer);
  __Pyx_ErrRestore(__pyx_type, __pyx_value, __pyx_tb);}
  __Pyx_AddTraceback("src.particles.interpolation.interpolation.distribution_interpolation", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = NULL;
  goto __pyx_L2;
  __pyx_L0:;
  __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_distribution.rcbuffer->pybuffer);
  __Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_grid.rcbuffer->pybuffer);
  __pyx_L2:;
  __Pyx_XDECREF(__pyx_v_index);
  __Pyx_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":194

experimental exception made for getbuffer and releasebuffer

-- the details of this may change.

    def __getbuffer__(ndarray self, Py_buffer* info, int flags)
    :   # <<<<<<<<<<<<<<

This implementation of getbuffer is geared towards Cython

requirements, and does not yet fullfill the PEP.

#

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));

#

function exit code

  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

static int __pyx_pf_5numpy_7ndarray___getbuffer__(PyArrayObject *__pyx_v_self, Py_buffer *__pyx_v_info, int __pyx_v_flags) {

#

  int __pyx_v_copy_shape;

#

  int __pyx_v_i;

#

  int __pyx_v_ndim;

#

  int __pyx_v_endian_detector;

#

  int __pyx_v_little_endian;

#

  int __pyx_v_t;
  char *__pyx_v_f;
  PyArray_Descr *__pyx_v_descr = 0;

#

  int __pyx_v_offset;

#

  int __pyx_v_hasfields;

#

  int __pyx_r;
  __Pyx_RefNannyDeclarations

#

  int __pyx_t_1;

#

  int __pyx_t_2;
  PyObject *__pyx_t_3 = NULL;

#

  int __pyx_t_4;

#

  int __pyx_t_5;
  PyObject *__pyx_t_6 = NULL;
  char *__pyx_t_7;

#

  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;

#

  int __pyx_clineno = 0;
  __Pyx_RefNannySetupContext("__getbuffer__", 0);
  if (__pyx_v_info != NULL) {
    __pyx_v_info->obj = Py_None; __Pyx_INCREF(Py_None);
    __Pyx_GIVEREF(__pyx_v_info->obj);
  }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":200

of flags

      if info == NULL: return             # <<<<<<<<<<<<<<

      cdef int copy_shape, i, ndim

  __pyx_t_1 = ((__pyx_v_info == NULL) != 0);
  if (__pyx_t_1) {
    __pyx_r = 0;
    goto __pyx_L0;
  }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":203

      cdef int copy_shape, i, ndim
      cdef int endian_detector = 1             # <<<<<<<<<<<<<<
      cdef bint little_endian = ((<char*>&endian_detector)[0] != 0)

  __pyx_v_endian_detector = 1;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":204 cdef int copy_shape, i, ndim cdef int endian_detector = 1 cdef bint little_endian = ((&endian_detector)[0] != 0) # <<<<<<<<<<<<<<

      ndim = PyArray_NDIM(self)

  __pyx_v_little_endian = ((((char *)(&__pyx_v_endian_detector))[0]) != 0);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":206 cdef bint little_endian = ((&endian_detector)[0] != 0)

      ndim = PyArray_NDIM(self)             # <<<<<<<<<<<<<<

      if sizeof(npy_intp) != sizeof(Py_ssize_t):

  __pyx_v_ndim = PyArray_NDIM(__pyx_v_self);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":208 ndim = PyArray_NDIM(self)

      if sizeof(npy_intp) != sizeof(Py_ssize_t)
      :   # <<<<<<<<<<<<<<

          copy_shape = 1
      else:

  __pyx_t_1 = (((sizeof(npy_intp)) != (sizeof(Py_ssize_t))) != 0);
  if (__pyx_t_1) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":209

    if sizeof(npy_intp) != sizeof(Py_ssize_t):
        copy_shape = 1             # <<<<<<<<<<<<<<
    else:
        copy_shape = 0

    __pyx_v_copy_shape = 1;
    goto __pyx_L4;
  }

#

else

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":211 copy_shape = 1 else: copy_shape = 0 # <<<<<<<<<<<<<<

    if ((flags & pybuf.PyBUF_C_CONTIGUOUS == pybuf.PyBUF_C_CONTIGUOUS)

    __pyx_v_copy_shape = 0;
  }
  __pyx_L4:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":213 copy_shape = 0

      if ((flags & pybuf.PyBUF_C_CONTIGUOUS == pybuf.PyBUF_C_CONTIGUOUS)             # <<<<<<<<<<<<<<
          and not PyArray_CHKFLAGS(self, NPY_C_CONTIGUOUS)):
          raise ValueError(u"ndarray is not C contiguous")

  __pyx_t_2 = (((__pyx_v_flags & PyBUF_C_CONTIGUOUS) == PyBUF_C_CONTIGUOUS) != 0);
  if (__pyx_t_2) {
    goto __pyx_L7_next_and;
  } else {
    __pyx_t_1 = __pyx_t_2;
    goto __pyx_L6_bool_binop_done;
  }
  __pyx_L7_next_and:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":214

      if ((flags & pybuf.PyBUF_C_CONTIGUOUS == pybuf.PyBUF_C_CONTIGUOUS)
          and not PyArray_CHKFLAGS(self, NPY_C_CONTIGUOUS))
          :   # <<<<<<<<<<<<<<

          raise ValueError(u"ndarray is not C contiguous")

  __pyx_t_2 = ((!(PyArray_CHKFLAGS(__pyx_v_self, NPY_C_CONTIGUOUS) != 0)) != 0);
  __pyx_t_1 = __pyx_t_2;
  __pyx_L6_bool_binop_done:;
  if (__pyx_t_1) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":215 if ((flags & pybuf.PyBUF_C_CONTIGUOUS == pybuf.PyBUF_C_CONTIGUOUS) and not PyArray_CHKFLAGS(self, NPY_C_CONTIGUOUS)): raise ValueError(u"ndarray is not C contiguous") # <<<<<<<<<<<<<<

    if ((flags & pybuf.PyBUF_F_CONTIGUOUS == pybuf.PyBUF_F_CONTIGUOUS)

    __pyx_t_3 = __Pyx_PyObject_Call(__pyx_builtin_ValueError, __pyx_tuple_, NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 215; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_3);
    __Pyx_Raise(__pyx_t_3, 0, 0, 0);
    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
    {__pyx_filename = __pyx_f[1]; __pyx_lineno = 215; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":217 raise ValueError(u"ndarray is not C contiguous")

      if ((flags & pybuf.PyBUF_F_CONTIGUOUS == pybuf.PyBUF_F_CONTIGUOUS)             # <<<<<<<<<<<<<<
          and not PyArray_CHKFLAGS(self, NPY_F_CONTIGUOUS)):
          raise ValueError(u"ndarray is not Fortran contiguous")

  __pyx_t_2 = (((__pyx_v_flags & PyBUF_F_CONTIGUOUS) == PyBUF_F_CONTIGUOUS) != 0);
  if (__pyx_t_2) {
    goto __pyx_L10_next_and;
  } else {
    __pyx_t_1 = __pyx_t_2;
    goto __pyx_L9_bool_binop_done;
  }
  __pyx_L10_next_and:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":218

      if ((flags & pybuf.PyBUF_F_CONTIGUOUS == pybuf.PyBUF_F_CONTIGUOUS)
          and not PyArray_CHKFLAGS(self, NPY_F_CONTIGUOUS))
          :   # <<<<<<<<<<<<<<

          raise ValueError(u"ndarray is not Fortran contiguous")

  __pyx_t_2 = ((!(PyArray_CHKFLAGS(__pyx_v_self, NPY_F_CONTIGUOUS) != 0)) != 0);
  __pyx_t_1 = __pyx_t_2;
  __pyx_L9_bool_binop_done:;
  if (__pyx_t_1) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":219 if ((flags & pybuf.PyBUF_F_CONTIGUOUS == pybuf.PyBUF_F_CONTIGUOUS) and not PyArray_CHKFLAGS(self, NPY_F_CONTIGUOUS)): raise ValueError(u"ndarray is not Fortran contiguous") # <<<<<<<<<<<<<<

    info.buf = PyArray_DATA(self)

    __pyx_t_3 = __Pyx_PyObject_Call(__pyx_builtin_ValueError, __pyx_tuple__2, NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 219; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_3);
    __Pyx_Raise(__pyx_t_3, 0, 0, 0);
    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
    {__pyx_filename = __pyx_f[1]; __pyx_lineno = 219; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":221 raise ValueError(u"ndarray is not Fortran contiguous")

      info.buf = PyArray_DATA(self)             # <<<<<<<<<<<<<<
      info.ndim = ndim
      if copy_shape:

  __pyx_v_info->buf = PyArray_DATA(__pyx_v_self);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":222

      info.buf = PyArray_DATA(self)
      info.ndim = ndim             # <<<<<<<<<<<<<<
      if copy_shape:

Allocate new buffer for strides and shape info.

  __pyx_v_info->ndim = __pyx_v_ndim;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":223 info.buf = PyArray_DATA(self) info.ndim = ndim if copy_shape : # <<<<<<<<<<<<<<

Allocate new buffer for strides and shape info.

This is allocated as one block, strides first.

  __pyx_t_1 = (__pyx_v_copy_shape != 0);
  if (__pyx_t_1) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":226

Allocate new buffer for strides and shape info.

This is allocated as one block, strides first.

        info.strides = <Py_ssize_t*>stdlib.malloc(sizeof(Py_ssize_t) * <size_t>ndim * 2)             # <<<<<<<<<<<<<<
        info.shape = info.strides + ndim
        for i in range(ndim):

    __pyx_v_info->strides = ((Py_ssize_t *)malloc((((sizeof(Py_ssize_t)) * ((size_t)__pyx_v_ndim)) * 2)));

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":227

This is allocated as one block, strides first.

        info.strides = <Py_ssize_t*>stdlib.malloc(sizeof(Py_ssize_t) * <size_t>ndim * 2)
        info.shape = info.strides + ndim             # <<<<<<<<<<<<<<
        for i in range(ndim):
            info.strides[i] = PyArray_STRIDES(self)[i]

    __pyx_v_info->shape = (__pyx_v_info->strides + __pyx_v_ndim);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":228 info.strides = stdlib.malloc(sizeof(Py_ssize_t) * ndim * 2) info.shape = info.strides + ndim for i in range(ndim) : # <<<<<<<<<<<<<<

            info.strides[i] = PyArray_STRIDES(self)[i]
            info.shape[i] = PyArray_DIMS(self)[i]

    __pyx_t_4 = __pyx_v_ndim;
    for (__pyx_t_5 = 0; __pyx_t_5 < __pyx_t_4; __pyx_t_5+=1) {
      __pyx_v_i = __pyx_t_5;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":229 info.shape = info.strides + ndim for i in range(ndim): info.strides[i] = PyArray_STRIDES(self)[i] # <<<<<<<<<<<<<< info.shape[i] = PyArray_DIMS(self)[i] else:

      (__pyx_v_info->strides[__pyx_v_i]) = (PyArray_STRIDES(__pyx_v_self)[__pyx_v_i]);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":230 for i in range(ndim): info.strides[i] = PyArray_STRIDES(self)[i] info.shape[i] = PyArray_DIMS(self)[i] # <<<<<<<<<<<<<< else: info.strides = PyArray_STRIDES(self)

      (__pyx_v_info->shape[__pyx_v_i]) = (PyArray_DIMS(__pyx_v_self)[__pyx_v_i]);
    }
    goto __pyx_L11;
  }

#

else

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":232 info.shape[i] = PyArray_DIMS(self)[i] else: info.strides = PyArray_STRIDES(self) # <<<<<<<<<<<<<< info.shape = PyArray_DIMS(self) info.suboffsets = NULL

    __pyx_v_info->strides = ((Py_ssize_t *)PyArray_STRIDES(__pyx_v_self));

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":233 else: info.strides = PyArray_STRIDES(self) info.shape = PyArray_DIMS(self) # <<<<<<<<<<<<<< info.suboffsets = NULL info.itemsize = PyArray_ITEMSIZE(self)

    __pyx_v_info->shape = ((Py_ssize_t *)PyArray_DIMS(__pyx_v_self));
  }
  __pyx_L11:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":234 info.strides = PyArray_STRIDES(self) info.shape = PyArray_DIMS(self) info.suboffsets = NULL # <<<<<<<<<<<<<< info.itemsize = PyArray_ITEMSIZE(self) info.readonly = not PyArray_ISWRITEABLE(self)

  __pyx_v_info->suboffsets = NULL;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":235 info.shape = PyArray_DIMS(self) info.suboffsets = NULL info.itemsize = PyArray_ITEMSIZE(self) # <<<<<<<<<<<<<< info.readonly = not PyArray_ISWRITEABLE(self)

  __pyx_v_info->itemsize = PyArray_ITEMSIZE(__pyx_v_self);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":236 info.suboffsets = NULL info.itemsize = PyArray_ITEMSIZE(self) info.readonly = not PyArray_ISWRITEABLE(self) # <<<<<<<<<<<<<<

      cdef int t

  __pyx_v_info->readonly = (!(PyArray_ISWRITEABLE(__pyx_v_self) != 0));

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":239

      cdef int t
      cdef char* f = NULL             # <<<<<<<<<<<<<<
      cdef dtype descr = self.descr
      cdef list stack

  __pyx_v_f = NULL;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":240 cdef int t cdef char* f = NULL cdef dtype descr = self.descr # <<<<<<<<<<<<<< cdef list stack cdef int offset

  __pyx_t_3 = ((PyObject *)__pyx_v_self->descr);
  __Pyx_INCREF(__pyx_t_3);
  __pyx_v_descr = ((PyArray_Descr *)__pyx_t_3);
  __pyx_t_3 = 0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":244 cdef int offset

      cdef bint hasfields = PyDataType_HASFIELDS(descr)             # <<<<<<<<<<<<<<

      if not hasfields and not copy_shape:

  __pyx_v_hasfields = PyDataType_HASFIELDS(__pyx_v_descr);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":246 cdef bint hasfields = PyDataType_HASFIELDS(descr)

      if not hasfields and not copy_shape
      :   # <<<<<<<<<<<<<<

do not call releasebuffer

          info.obj = None

  __pyx_t_2 = ((!(__pyx_v_hasfields != 0)) != 0);
  if (__pyx_t_2) {
    goto __pyx_L16_next_and;
  } else {
    __pyx_t_1 = __pyx_t_2;
    goto __pyx_L15_bool_binop_done;
  }
  __pyx_L16_next_and:;
  __pyx_t_2 = ((!(__pyx_v_copy_shape != 0)) != 0);
  __pyx_t_1 = __pyx_t_2;
  __pyx_L15_bool_binop_done:;
  if (__pyx_t_1) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":248 if not hasfields and not copy_shape:

do not call releasebuffer

        info.obj = None             # <<<<<<<<<<<<<<
    else:

need to call releasebuffer

    __Pyx_INCREF(Py_None);
    __Pyx_GIVEREF(Py_None);
    __Pyx_GOTREF(__pyx_v_info->obj);
    __Pyx_DECREF(__pyx_v_info->obj);
    __pyx_v_info->obj = Py_None;
    goto __pyx_L14;
  }

#

else

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":251 else:

need to call releasebuffer

        info.obj = self             # <<<<<<<<<<<<<<

    if not hasfields:

    __Pyx_INCREF(((PyObject *)__pyx_v_self));
    __Pyx_GIVEREF(((PyObject *)__pyx_v_self));
    __Pyx_GOTREF(__pyx_v_info->obj);
    __Pyx_DECREF(__pyx_v_info->obj);
    __pyx_v_info->obj = ((PyObject *)__pyx_v_self);
  }
  __pyx_L14:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":253 info.obj = self

      if not hasfields
      :   # <<<<<<<<<<<<<<

          t = descr.type_num
          if ((descr.byteorder == c'>' and little_endian) or

  __pyx_t_1 = ((!(__pyx_v_hasfields != 0)) != 0);
  if (__pyx_t_1) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":254

    if not hasfields:
        t = descr.type_num             # <<<<<<<<<<<<<<
        if ((descr.byteorder == c'>' and little_endian) or
            (descr.byteorder == c'<' and not little_endian)):

    __pyx_t_4 = __pyx_v_descr->type_num;
    __pyx_v_t = __pyx_t_4;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":255 if not hasfields: t = descr.type_num if ((descr.byteorder == c'>' and little_endian) or # <<<<<<<<<<<<<< (descr.byteorder == c'<' and not little_endian)): raise ValueError(u"Non-native byte order not supported")

    __pyx_t_2 = ((__pyx_v_descr->byteorder == '>') != 0);
    if (!__pyx_t_2) {
      goto __pyx_L20_next_or;
    } else {
      goto __pyx_L21_next_and;
    }
    __pyx_L21_next_and:;
    __pyx_t_2 = (__pyx_v_little_endian != 0);
    if (!__pyx_t_2) {
      goto __pyx_L20_next_or;
    } else {
      __pyx_t_1 = __pyx_t_2;
      goto __pyx_L19_bool_binop_done;
    }
    __pyx_L20_next_or:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":256 t = descr.type_num if ((descr.byteorder == c'>' and little_endian) or (descr.byteorder == c'<' and not little_endian)) : # <<<<<<<<<<<<<<

            raise ValueError(u"Non-native byte order not supported")
        if   t == NPY_BYTE
        :   f = "b"

    __pyx_t_2 = ((__pyx_v_descr->byteorder == '<') != 0);
    if (__pyx_t_2) {
      goto __pyx_L22_next_and;
    } else {
      __pyx_t_1 = __pyx_t_2;
      goto __pyx_L19_bool_binop_done;
    }
    __pyx_L22_next_and:;
    __pyx_t_2 = ((!(__pyx_v_little_endian != 0)) != 0);
    __pyx_t_1 = __pyx_t_2;
    __pyx_L19_bool_binop_done:;
    if (__pyx_t_1) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":257 if ((descr.byteorder == c'>' and little_endian) or (descr.byteorder == c'<' and not little_endian)): raise ValueError(u"Non-native byte order not supported") # <<<<<<<<<<<<<< if t == NPY_BYTE : f = "b"

      elif t == NPY_UBYTE
      :   f = "B"

      __pyx_t_3 = __Pyx_PyObject_Call(__pyx_builtin_ValueError, __pyx_tuple__3, NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 257; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __Pyx_Raise(__pyx_t_3, 0, 0, 0);
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      {__pyx_filename = __pyx_f[1]; __pyx_lineno = 257; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":274 elif t == NPY_CDOUBLE : f = "Zd"

        elif t == NPY_CLONGDOUBLE: f = "Zg"
        elif t == NPY_OBJECT
        :   f = "O"             # <<<<<<<<<<<<<<

        else:
            raise ValueError(u"unknown dtype code in numpy.pxd (%d)" % t)

    switch (__pyx_v_t) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":258 (descr.byteorder == c'<' and not little_endian)): raise ValueError(u"Non-native byte order not supported") if t == NPY_BYTE : f = "b" # <<<<<<<<<<<<<<

      elif t == NPY_UBYTE
      :   f = "B"

      elif t == NPY_SHORT
      :   f = "h"

      case NPY_BYTE:
      __pyx_v_f = __pyx_k_b;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":259 raise ValueError(u"Non-native byte order not supported") if t == NPY_BYTE : f = "b"

      elif t == NPY_UBYTE
      :   f = "B"             # <<<<<<<<<<<<<<

      elif t == NPY_SHORT
      :   f = "h"

      elif t == NPY_USHORT
      :   f = "H"

      case NPY_UBYTE:
      __pyx_v_f = __pyx_k_B;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":260 if t == NPY_BYTE : f = "b"

      elif t == NPY_UBYTE
      :   f = "B"

      elif t == NPY_SHORT
      :   f = "h"             # <<<<<<<<<<<<<<

      elif t == NPY_USHORT
      :   f = "H"

      elif t == NPY_INT
      :   f = "i"

      case NPY_SHORT:
      __pyx_v_f = __pyx_k_h;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":261 elif t == NPY_UBYTE : f = "B"

      elif t == NPY_SHORT
      :   f = "h"

      elif t == NPY_USHORT
      :   f = "H"             # <<<<<<<<<<<<<<

      elif t == NPY_INT
      :   f = "i"

      elif t == NPY_UINT
      :   f = "I"

      case NPY_USHORT:
      __pyx_v_f = __pyx_k_H;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":262 elif t == NPY_SHORT : f = "h"

      elif t == NPY_USHORT
      :   f = "H"

      elif t == NPY_INT
      :   f = "i"             # <<<<<<<<<<<<<<

      elif t == NPY_UINT
      :   f = "I"

      elif t == NPY_LONG
      :   f = "l"

      case NPY_INT:
      __pyx_v_f = __pyx_k_i;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":263 elif t == NPY_USHORT : f = "H"

      elif t == NPY_INT
      :   f = "i"

      elif t == NPY_UINT
      :   f = "I"             # <<<<<<<<<<<<<<

      elif t == NPY_LONG
      :   f = "l"

      elif t == NPY_ULONG
      :   f = "L"

      case NPY_UINT:
      __pyx_v_f = __pyx_k_I;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":264 elif t == NPY_INT : f = "i"

      elif t == NPY_UINT
      :   f = "I"

      elif t == NPY_LONG
      :   f = "l"             # <<<<<<<<<<<<<<

      elif t == NPY_ULONG
      :   f = "L"

      elif t == NPY_LONGLONG
      :   f = "q"

      case NPY_LONG:
      __pyx_v_f = __pyx_k_l;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":265 elif t == NPY_UINT : f = "I"

      elif t == NPY_LONG
      :   f = "l"

      elif t == NPY_ULONG
      :   f = "L"             # <<<<<<<<<<<<<<

      elif t == NPY_LONGLONG
      :   f = "q"

      elif t == NPY_ULONGLONG
      :   f = "Q"

      case NPY_ULONG:
      __pyx_v_f = __pyx_k_L;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":266 elif t == NPY_LONG : f = "l"

      elif t == NPY_ULONG
      :   f = "L"

      elif t == NPY_LONGLONG
      :   f = "q"             # <<<<<<<<<<<<<<

      elif t == NPY_ULONGLONG
      :   f = "Q"

      elif t == NPY_FLOAT
      :   f = "f"

      case NPY_LONGLONG:
      __pyx_v_f = __pyx_k_q;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":267 elif t == NPY_ULONG : f = "L"

      elif t == NPY_LONGLONG
      :   f = "q"

      elif t == NPY_ULONGLONG
      :   f = "Q"             # <<<<<<<<<<<<<<

      elif t == NPY_FLOAT
      :   f = "f"

      elif t == NPY_DOUBLE
      :   f = "d"

      case NPY_ULONGLONG:
      __pyx_v_f = __pyx_k_Q;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":268 elif t == NPY_LONGLONG : f = "q"

      elif t == NPY_ULONGLONG
      :   f = "Q"

      elif t == NPY_FLOAT
      :   f = "f"             # <<<<<<<<<<<<<<

      elif t == NPY_DOUBLE
      :   f = "d"

      elif t == NPY_LONGDOUBLE
      :   f = "g"

      case NPY_FLOAT:
      __pyx_v_f = __pyx_k_f;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":269 elif t == NPY_ULONGLONG : f = "Q"

      elif t == NPY_FLOAT
      :   f = "f"

      elif t == NPY_DOUBLE
      :   f = "d"             # <<<<<<<<<<<<<<

      elif t == NPY_LONGDOUBLE
      :   f = "g"

      elif t == NPY_CFLOAT
      :   f = "Zf"

      case NPY_DOUBLE:
      __pyx_v_f = __pyx_k_d;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":270 elif t == NPY_FLOAT : f = "f"

      elif t == NPY_DOUBLE
      :   f = "d"

      elif t == NPY_LONGDOUBLE
      :   f = "g"             # <<<<<<<<<<<<<<

      elif t == NPY_CFLOAT
      :   f = "Zf"

      elif t == NPY_CDOUBLE
      :   f = "Zd"

      case NPY_LONGDOUBLE:
      __pyx_v_f = __pyx_k_g;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":271 elif t == NPY_DOUBLE : f = "d"

      elif t == NPY_LONGDOUBLE
      :   f = "g"

      elif t == NPY_CFLOAT
      :   f = "Zf"             # <<<<<<<<<<<<<<

      elif t == NPY_CDOUBLE
      :   f = "Zd"

      elif t == NPY_CLONGDOUBLE: f = "Zg"

      case NPY_CFLOAT:
      __pyx_v_f = __pyx_k_Zf;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":272 elif t == NPY_LONGDOUBLE : f = "g"

      elif t == NPY_CFLOAT
      :   f = "Zf"

      elif t == NPY_CDOUBLE
      :   f = "Zd"             # <<<<<<<<<<<<<<

      elif t == NPY_CLONGDOUBLE: f = "Zg"
      elif t == NPY_OBJECT
      :   f = "O"

      case NPY_CDOUBLE:
      __pyx_v_f = __pyx_k_Zd;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":273 elif t == NPY_CFLOAT : f = "Zf"

      elif t == NPY_CDOUBLE
      :   f = "Zd"

      elif t == NPY_CLONGDOUBLE: f = "Zg"             # <<<<<<<<<<<<<<
      elif t == NPY_OBJECT
      :   f = "O"

      else:

      case NPY_CLONGDOUBLE:
      __pyx_v_f = __pyx_k_Zg;
      break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":274 elif t == NPY_CDOUBLE : f = "Zd"

      elif t == NPY_CLONGDOUBLE: f = "Zg"
      elif t == NPY_OBJECT
      :   f = "O"             # <<<<<<<<<<<<<<

      else:
          raise ValueError(u"unknown dtype code in numpy.pxd (%d)" % t)

      case NPY_OBJECT:
      __pyx_v_f = __pyx_k_O;
      break;
      default:

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":276 elif t == NPY_OBJECT : f = "O"

      else:
          raise ValueError(u"unknown dtype code in numpy.pxd (%d)" % t)             # <<<<<<<<<<<<<<
      info.format = f
      return

      __pyx_t_3 = __Pyx_PyInt_From_int(__pyx_v_t); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 276; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_6 = PyUnicode_Format(__pyx_kp_u_unknown_dtype_code_in_numpy_pxd, __pyx_t_3); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 276; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_6);
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_t_3 = PyTuple_New(1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 276; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_t_6);
      __Pyx_GIVEREF(__pyx_t_6);
      __pyx_t_6 = 0;
      __pyx_t_6 = __Pyx_PyObject_Call(__pyx_builtin_ValueError, __pyx_t_3, NULL); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 276; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_6);
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __Pyx_Raise(__pyx_t_6, 0, 0, 0);
      __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0;
      {__pyx_filename = __pyx_f[1]; __pyx_lineno = 276; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      break;
    }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":277 else: raise ValueError(u"unknown dtype code in numpy.pxd (%d)" % t) info.format = f # <<<<<<<<<<<<<< return else:

    __pyx_v_info->format = __pyx_v_f;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":278 raise ValueError(u"unknown dtype code in numpy.pxd (%d)" % t) info.format = f return # <<<<<<<<<<<<<< else: info.format = stdlib.malloc(_buffer_format_string_len)

    __pyx_r = 0;
    goto __pyx_L0;
  }

#

else

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":280 return else: info.format = stdlib.malloc(_buffer_format_string_len) # <<<<<<<<<<<<<< info.format[0] = c'^' # Native data types, manual alignment offset = 0

    __pyx_v_info->format = ((char *)malloc(255));

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":281 else: info.format = stdlib.malloc(_buffer_format_string_len) info.format[0] = c'^' # Native data types, manual alignment # <<<<<<<<<<<<<< offset = 0 f = _util_dtypestring(descr, info.format + 1,

    (__pyx_v_info->format[0]) = '^';

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":282 info.format = stdlib.malloc(_buffer_format_string_len) info.format[0] = c'^' # Native data types, manual alignment offset = 0 # <<<<<<<<<<<<<< f = _util_dtypestring(descr, info.format + 1, info.format + _buffer_format_string_len,

    __pyx_v_offset = 0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":283 info.format[0] = c'^' # Native data types, manual alignment offset = 0 f = _util_dtypestring(descr, info.format + 1, # <<<<<<<<<<<<<< info.format + _buffer_format_string_len, &offset)

    __pyx_t_7 = __pyx_f_5numpy__util_dtypestring(__pyx_v_descr, (__pyx_v_info->format + 1), (__pyx_v_info->format + 255), (&__pyx_v_offset)); if (unlikely(__pyx_t_7 == NULL)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 283; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __pyx_v_f = __pyx_t_7;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":286 info.format + _buffer_format_string_len, &offset) f[0] = c'\0' # Terminate format string # <<<<<<<<<<<<<<

def __releasebuffer__(ndarray self, Py_buffer* info):

    (__pyx_v_f[0]) = '\x00';
  }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":194

experimental exception made for getbuffer and releasebuffer

-- the details of this may change.

  def __getbuffer__(ndarray self, Py_buffer* info, int flags)
  :   # <<<<<<<<<<<<<<

This implementation of getbuffer is geared towards Cython

requirements, and does not yet fullfill the PEP.

#

function exit code

  __pyx_r = 0;
  goto __pyx_L0;
  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_3);
  __Pyx_XDECREF(__pyx_t_6);
  __Pyx_AddTraceback("numpy.ndarray.__getbuffer__", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = -1;
  if (__pyx_v_info != NULL && __pyx_v_info->obj != NULL) {
    __Pyx_GOTREF(__pyx_v_info->obj);
    __Pyx_DECREF(__pyx_v_info->obj); __pyx_v_info->obj = NULL;
  }
  goto __pyx_L2;
  __pyx_L0:;
  if (__pyx_v_info != NULL && __pyx_v_info->obj == Py_None) {
    __Pyx_GOTREF(Py_None);
    __Pyx_DECREF(Py_None); __pyx_v_info->obj = NULL;
  }
  __pyx_L2:;
  __Pyx_XDECREF((PyObject *)__pyx_v_descr);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":288 f[0] = c'\0' # Terminate format string

    def __releasebuffer__(ndarray self, Py_buffer* info)
    :   # <<<<<<<<<<<<<<

        if PyArray_HASFIELDS(self):
            stdlib.free(info.format)

#

Python wrapper

static CYTHON_UNUSED void __pyx_pw_5numpy_7ndarray_3__releasebuffer__(PyObject *__pyx_v_self, Py_buffer *__pyx_v_info); /*proto*/
static CYTHON_UNUSED void __pyx_pw_5numpy_7ndarray_3__releasebuffer__(PyObject *__pyx_v_self, Py_buffer *__pyx_v_info) {
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("__releasebuffer__ (wrapper)", 0);
  __pyx_pf_5numpy_7ndarray_2__releasebuffer__(((PyArrayObject *)__pyx_v_self), ((Py_buffer *)__pyx_v_info));

#

function exit code

  __Pyx_RefNannyFinishContext();
}

static void __pyx_pf_5numpy_7ndarray_2__releasebuffer__(PyArrayObject *__pyx_v_self, Py_buffer *__pyx_v_info) {
  __Pyx_RefNannyDeclarations

#

  int __pyx_t_1;
  __Pyx_RefNannySetupContext("__releasebuffer__", 0);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":289

  def __releasebuffer__(ndarray self, Py_buffer* info):
      if PyArray_HASFIELDS(self)
      :   # <<<<<<<<<<<<<<

          stdlib.free(info.format)
      if sizeof(npy_intp) != sizeof(Py_ssize_t):

  __pyx_t_1 = (PyArray_HASFIELDS(__pyx_v_self) != 0);
  if (__pyx_t_1) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":290 def releasebuffer(ndarray self, Py_buffer* info): if PyArray_HASFIELDS(self): stdlib.free(info.format) # <<<<<<<<<<<<<< if sizeof(npy_intp) != sizeof(Py_ssize_t): stdlib.free(info.strides)

    free(__pyx_v_info->format);
    goto __pyx_L3;
  }
  __pyx_L3:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":291 if PyArray_HASFIELDS(self): stdlib.free(info.format) if sizeof(npy_intp) != sizeof(Py_ssize_t) : # <<<<<<<<<<<<<<

          stdlib.free(info.strides)

info.shape was stored after info.strides in the same block

  __pyx_t_1 = (((sizeof(npy_intp)) != (sizeof(Py_ssize_t))) != 0);
  if (__pyx_t_1) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":292 stdlib.free(info.format) if sizeof(npy_intp) != sizeof(Py_ssize_t): stdlib.free(info.strides) # <<<<<<<<<<<<<<

info.shape was stored after info.strides in the same block

    free(__pyx_v_info->strides);
    goto __pyx_L4;
  }
  __pyx_L4:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":288 f[0] = c'\0' # Terminate format string

  def __releasebuffer__(ndarray self, Py_buffer* info)
  :   # <<<<<<<<<<<<<<

      if PyArray_HASFIELDS(self):
          stdlib.free(info.format)

#

function exit code

  __Pyx_RefNannyFinishContext();
}

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":768 ctypedef npy_cdouble complex_t

cdef inline object PyArray_MultiIterNew1(a): <<<<<<<<<<<<<<

return PyArray_MultiIterNew(1, a)

static CYTHON_INLINE PyObject *__pyx_f_5numpy_PyArray_MultiIterNew1(PyObject *__pyx_v_a) {
  PyObject *__pyx_r = NULL;
  __Pyx_RefNannyDeclarations
  PyObject *__pyx_t_1 = NULL;

  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;

#

  int __pyx_clineno = 0;
  __Pyx_RefNannySetupContext("PyArray_MultiIterNew1", 0);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":769

cdef inline object PyArray_MultiIterNew1(a): return PyArray_MultiIterNew(1, a) # <<<<<<<<<<<<<<

cdef inline object PyArray_MultiIterNew2(a, b):

  __Pyx_XDECREF(__pyx_r);
  __pyx_t_1 = PyArray_MultiIterNew(1, ((void *)__pyx_v_a)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 769; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_r = __pyx_t_1;
  __pyx_t_1 = 0;
  goto __pyx_L0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":768 ctypedef npy_cdouble complex_t

cdef inline object PyArray_MultiIterNew1(a): <<<<<<<<<<<<<<

return PyArray_MultiIterNew(1, a)

#

function exit code

  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_1);
  __Pyx_AddTraceback("numpy.PyArray_MultiIterNew1", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = 0;
  __pyx_L0:;
  __Pyx_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":771 return PyArray_MultiIterNew(1, a)

cdef inline object PyArray_MultiIterNew2(a, b): <<<<<<<<<<<<<<

return PyArray_MultiIterNew(2, a, b)

static CYTHON_INLINE PyObject *__pyx_f_5numpy_PyArray_MultiIterNew2(PyObject *__pyx_v_a, PyObject *__pyx_v_b) {
  PyObject *__pyx_r = NULL;
  __Pyx_RefNannyDeclarations
  PyObject *__pyx_t_1 = NULL;

  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;

#

  int __pyx_clineno = 0;
  __Pyx_RefNannySetupContext("PyArray_MultiIterNew2", 0);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":772

cdef inline object PyArray_MultiIterNew2(a, b): return PyArray_MultiIterNew(2, a, b) # <<<<<<<<<<<<<<

cdef inline object PyArray_MultiIterNew3(a, b, c):

  __Pyx_XDECREF(__pyx_r);
  __pyx_t_1 = PyArray_MultiIterNew(2, ((void *)__pyx_v_a), ((void *)__pyx_v_b)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 772; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_r = __pyx_t_1;
  __pyx_t_1 = 0;
  goto __pyx_L0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":771 return PyArray_MultiIterNew(1, a)

cdef inline object PyArray_MultiIterNew2(a, b): <<<<<<<<<<<<<<

return PyArray_MultiIterNew(2, a, b)

#

function exit code

  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_1);
  __Pyx_AddTraceback("numpy.PyArray_MultiIterNew2", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = 0;
  __pyx_L0:;
  __Pyx_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":774 return PyArray_MultiIterNew(2, a, b)

cdef inline object PyArray_MultiIterNew3(a, b, c): <<<<<<<<<<<<<<

return PyArray_MultiIterNew(3, a, b, c)

static CYTHON_INLINE PyObject *__pyx_f_5numpy_PyArray_MultiIterNew3(PyObject *__pyx_v_a, PyObject *__pyx_v_b, PyObject *__pyx_v_c) {
  PyObject *__pyx_r = NULL;
  __Pyx_RefNannyDeclarations
  PyObject *__pyx_t_1 = NULL;

  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;

#

  int __pyx_clineno = 0;
  __Pyx_RefNannySetupContext("PyArray_MultiIterNew3", 0);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":775

cdef inline object PyArray_MultiIterNew3(a, b, c): return PyArray_MultiIterNew(3, a, b, c) # <<<<<<<<<<<<<<

cdef inline object PyArray_MultiIterNew4(a, b, c, d):

  __Pyx_XDECREF(__pyx_r);
  __pyx_t_1 = PyArray_MultiIterNew(3, ((void *)__pyx_v_a), ((void *)__pyx_v_b), ((void *)__pyx_v_c)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 775; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_r = __pyx_t_1;
  __pyx_t_1 = 0;
  goto __pyx_L0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":774 return PyArray_MultiIterNew(2, a, b)

cdef inline object PyArray_MultiIterNew3(a, b, c): <<<<<<<<<<<<<<

return PyArray_MultiIterNew(3, a, b, c)

#

function exit code

  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_1);
  __Pyx_AddTraceback("numpy.PyArray_MultiIterNew3", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = 0;
  __pyx_L0:;
  __Pyx_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":777 return PyArray_MultiIterNew(3, a, b, c)

cdef inline object PyArray_MultiIterNew4(a, b, c, d): <<<<<<<<<<<<<<

return PyArray_MultiIterNew(4, a, b, c, d)

static CYTHON_INLINE PyObject *__pyx_f_5numpy_PyArray_MultiIterNew4(PyObject *__pyx_v_a, PyObject *__pyx_v_b, PyObject *__pyx_v_c, PyObject *__pyx_v_d) {
  PyObject *__pyx_r = NULL;
  __Pyx_RefNannyDeclarations
  PyObject *__pyx_t_1 = NULL;

  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;

#

  int __pyx_clineno = 0;
  __Pyx_RefNannySetupContext("PyArray_MultiIterNew4", 0);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":778

cdef inline object PyArray_MultiIterNew4(a, b, c, d): return PyArray_MultiIterNew(4, a, b, c, d) # <<<<<<<<<<<<<<

cdef inline object PyArray_MultiIterNew5(a, b, c, d, e):

  __Pyx_XDECREF(__pyx_r);
  __pyx_t_1 = PyArray_MultiIterNew(4, ((void *)__pyx_v_a), ((void *)__pyx_v_b), ((void *)__pyx_v_c), ((void *)__pyx_v_d)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 778; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_r = __pyx_t_1;
  __pyx_t_1 = 0;
  goto __pyx_L0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":777 return PyArray_MultiIterNew(3, a, b, c)

cdef inline object PyArray_MultiIterNew4(a, b, c, d): <<<<<<<<<<<<<<

return PyArray_MultiIterNew(4, a, b, c, d)

#

function exit code

  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_1);
  __Pyx_AddTraceback("numpy.PyArray_MultiIterNew4", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = 0;
  __pyx_L0:;
  __Pyx_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":780 return PyArray_MultiIterNew(4, a, b, c, d)

cdef inline object PyArray_MultiIterNew5(a, b, c, d, e): <<<<<<<<<<<<<<

return PyArray_MultiIterNew(5, a, b, c, d, e)

static CYTHON_INLINE PyObject *__pyx_f_5numpy_PyArray_MultiIterNew5(PyObject *__pyx_v_a, PyObject *__pyx_v_b, PyObject *__pyx_v_c, PyObject *__pyx_v_d, PyObject *__pyx_v_e) {
  PyObject *__pyx_r = NULL;
  __Pyx_RefNannyDeclarations
  PyObject *__pyx_t_1 = NULL;

  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;

#

  int __pyx_clineno = 0;
  __Pyx_RefNannySetupContext("PyArray_MultiIterNew5", 0);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":781

cdef inline object PyArray_MultiIterNew5(a, b, c, d, e): return PyArray_MultiIterNew(5, a, b, c, d, e) # <<<<<<<<<<<<<<

cdef inline char _util_dtypestring(dtype descr, char f, char end, int offset) except NULL:

  __Pyx_XDECREF(__pyx_r);
  __pyx_t_1 = PyArray_MultiIterNew(5, ((void *)__pyx_v_a), ((void *)__pyx_v_b), ((void *)__pyx_v_c), ((void *)__pyx_v_d), ((void *)__pyx_v_e)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 781; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_r = __pyx_t_1;
  __pyx_t_1 = 0;
  goto __pyx_L0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":780 return PyArray_MultiIterNew(4, a, b, c, d)

cdef inline object PyArray_MultiIterNew5(a, b, c, d, e): <<<<<<<<<<<<<<

return PyArray_MultiIterNew(5, a, b, c, d, e)

#

function exit code

  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_1);
  __Pyx_AddTraceback("numpy.PyArray_MultiIterNew5", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = 0;
  __pyx_L0:;
  __Pyx_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":783 return PyArray_MultiIterNew(5, a, b, c, d, e)

cdef inline char _util_dtypestring(dtype descr, char f, char end, int offset) except NULL: <<<<<<<<<<<<<<

Recursive utility function used in getbuffer to get format

string. The new location in the format string is returned.

static CYTHON_INLINE char *__pyx_f_5numpy__util_dtypestring(PyArray_Descr *__pyx_v_descr, char *__pyx_v_f, char *__pyx_v_end, int *__pyx_v_offset) {
  PyArray_Descr *__pyx_v_child = 0;

  int __pyx_v_endian_detector;

#

  int __pyx_v_little_endian;
  PyObject *__pyx_v_fields = 0;
  PyObject *__pyx_v_childname = NULL;
  PyObject *__pyx_v_new_offset = NULL;
  PyObject *__pyx_v_t = NULL;
  char *__pyx_r;
  __Pyx_RefNannyDeclarations
  PyObject *__pyx_t_1 = NULL;
  Py_ssize_t __pyx_t_2;
  PyObject *__pyx_t_3 = NULL;
  PyObject *__pyx_t_4 = NULL;

#

  int __pyx_t_5;

#

  int __pyx_t_6;

#

  int __pyx_t_7;
  long __pyx_t_8;
  char *__pyx_t_9;

#

  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;

#

  int __pyx_clineno = 0;
  __Pyx_RefNannySetupContext("_util_dtypestring", 0);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":790 cdef int delta_offset cdef tuple i cdef int endian_detector = 1 # <<<<<<<<<<<<<< cdef bint little_endian = ((&endian_detector)[0] != 0) cdef tuple fields

  __pyx_v_endian_detector = 1;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":791 cdef tuple i cdef int endian_detector = 1 cdef bint little_endian = ((&endian_detector)[0] != 0) # <<<<<<<<<<<<<< cdef tuple fields

  __pyx_v_little_endian = ((((char *)(&__pyx_v_endian_detector))[0]) != 0);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":794 cdef tuple fields

for childname in descr.names: <<<<<<<<<<<<<<

fields = descr.fields[childname] child, new_offset = fields

  if (unlikely(__pyx_v_descr->names == Py_None)) {
    PyErr_SetString(PyExc_TypeError, "'NoneType' object is not iterable");
    {__pyx_filename = __pyx_f[1]; __pyx_lineno = 794; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  }
  __pyx_t_1 = __pyx_v_descr->names; __Pyx_INCREF(__pyx_t_1); __pyx_t_2 = 0;
  for (;;) {
    if (__pyx_t_2 >= PyTuple_GET_SIZE(__pyx_t_1)) break;
    #if CYTHON_COMPILING_IN_CPYTHON
    __pyx_t_3 = PyTuple_GET_ITEM(__pyx_t_1, __pyx_t_2); __Pyx_INCREF(__pyx_t_3); __pyx_t_2++; if (unlikely(0 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 794; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    #else
    __pyx_t_3 = PySequence_ITEM(__pyx_t_1, __pyx_t_2); __pyx_t_2++; if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 794; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    #endif
    __Pyx_XDECREF_SET(__pyx_v_childname, __pyx_t_3);
    __pyx_t_3 = 0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":795

for childname in descr.names: fields = descr.fields[childname] # <<<<<<<<<<<<<< child, new_offset = fields

    __pyx_t_3 = PyObject_GetItem(__pyx_v_descr->fields, __pyx_v_childname); if (unlikely(__pyx_t_3 == NULL)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 795; __pyx_clineno = __LINE__; goto __pyx_L1_error;};
    __Pyx_GOTREF(__pyx_t_3);
    if (!(likely(PyTuple_CheckExact(__pyx_t_3))||((__pyx_t_3) == Py_None)||(PyErr_Format(PyExc_TypeError, "Expected %.16s, got %.200s", "tuple", Py_TYPE(__pyx_t_3)->tp_name), 0))) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 795; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_XDECREF_SET(__pyx_v_fields, ((PyObject*)__pyx_t_3));
    __pyx_t_3 = 0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":796 for childname in descr.names: fields = descr.fields[childname] child, new_offset = fields # <<<<<<<<<<<<<<

if (end - f) - <int>(new_offset - offset[0]) < 15:

    if (likely(__pyx_v_fields != Py_None)) {
      PyObject* sequence = __pyx_v_fields;
      #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[1]; __pyx_lineno = 796; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      }
      #if CYTHON_COMPILING_IN_CPYTHON
      __pyx_t_3 = PyTuple_GET_ITEM(sequence, 0); 
      __pyx_t_4 = PyTuple_GET_ITEM(sequence, 1); 
      __Pyx_INCREF(__pyx_t_3);
      __Pyx_INCREF(__pyx_t_4);
      #else
      __pyx_t_3 = PySequence_ITEM(sequence, 0); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 796; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_4 = PySequence_ITEM(sequence, 1); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 796; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_4);
      #endif
    } else {
      __Pyx_RaiseNoneNotIterableError(); {__pyx_filename = __pyx_f[1]; __pyx_lineno = 796; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    }
    if (!(likely(((__pyx_t_3) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_3, __pyx_ptype_5numpy_dtype))))) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 796; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_XDECREF_SET(__pyx_v_child, ((PyArray_Descr *)__pyx_t_3));
    __pyx_t_3 = 0;
    __Pyx_XDECREF_SET(__pyx_v_new_offset, __pyx_t_4);
    __pyx_t_4 = 0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":798 child, new_offset = fields

if (end - f) - <int>(new_offset - offset[0]) < 15
:   # <<<<<<<<<<<<<<

    raise RuntimeError(u"Format string allocated too short, see comment in numpy.pxd")

    __pyx_t_4 = __Pyx_PyInt_From_int((__pyx_v_offset[0])); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 798; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_4);
    __pyx_t_3 = PyNumber_Subtract(__pyx_v_new_offset, __pyx_t_4); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 798; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_3);
    __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
    __pyx_t_5 = __Pyx_PyInt_As_int(__pyx_t_3); if (unlikely((__pyx_t_5 == (int)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 798; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
    __pyx_t_6 = ((((__pyx_v_end - __pyx_v_f) - ((int)__pyx_t_5)) < 15) != 0);
    if (__pyx_t_6) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":799

if (end - f) - (new_offset - offset[0]) < 15: raise RuntimeError(u"Format string allocated too short, see comment in numpy.pxd") # <<<<<<<<<<<<<<

if ((child.byteorder == c'>' and little_endian) or

      __pyx_t_3 = __Pyx_PyObject_Call(__pyx_builtin_RuntimeError, __pyx_tuple__4, NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 799; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __Pyx_Raise(__pyx_t_3, 0, 0, 0);
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      {__pyx_filename = __pyx_f[1]; __pyx_lineno = 799; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":801 raise RuntimeError(u"Format string allocated too short, see comment in numpy.pxd")

if ((child.byteorder == c'>' and little_endian) or             # <<<<<<<<<<<<<<
    (child.byteorder == c'<' and not little_endian)):
    raise ValueError(u"Non-native byte order not supported")

    __pyx_t_7 = ((__pyx_v_child->byteorder == '>') != 0);
    if (!__pyx_t_7) {
      goto __pyx_L8_next_or;
    } else {
      goto __pyx_L9_next_and;
    }
    __pyx_L9_next_and:;
    __pyx_t_7 = (__pyx_v_little_endian != 0);
    if (!__pyx_t_7) {
      goto __pyx_L8_next_or;
    } else {
      __pyx_t_6 = __pyx_t_7;
      goto __pyx_L7_bool_binop_done;
    }
    __pyx_L8_next_or:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":802

if ((child.byteorder == c'>' and little_endian) or
    (child.byteorder == c'<' and not little_endian))
    :   # <<<<<<<<<<<<<<

    raise ValueError(u"Non-native byte order not supported")

One could encode it in the format string and have Cython

    __pyx_t_7 = ((__pyx_v_child->byteorder == '<') != 0);
    if (__pyx_t_7) {
      goto __pyx_L10_next_and;
    } else {
      __pyx_t_6 = __pyx_t_7;
      goto __pyx_L7_bool_binop_done;
    }
    __pyx_L10_next_and:;
    __pyx_t_7 = ((!(__pyx_v_little_endian != 0)) != 0);
    __pyx_t_6 = __pyx_t_7;
    __pyx_L7_bool_binop_done:;
    if (__pyx_t_6) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":803 if ((child.byteorder == c'>' and little_endian) or (child.byteorder == c'<' and not little_endian)): raise ValueError(u"Non-native byte order not supported") # <<<<<<<<<<<<<<

One could encode it in the format string and have Cython

      __pyx_t_3 = __Pyx_PyObject_Call(__pyx_builtin_ValueError, __pyx_tuple__5, NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 803; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __Pyx_Raise(__pyx_t_3, 0, 0, 0);
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      {__pyx_filename = __pyx_f[1]; __pyx_lineno = 803; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":813

Output padding bytes

while offset[0] < new_offset
:   # <<<<<<<<<<<<<<

    f[0] = 120 # "x"; pad byte
    f += 1

    while (1) {
      __pyx_t_3 = __Pyx_PyInt_From_int((__pyx_v_offset[0])); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 813; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_4 = PyObject_RichCompare(__pyx_t_3, __pyx_v_new_offset, Py_LT); __Pyx_XGOTREF(__pyx_t_4); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 813; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_4); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 813; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      if (!__pyx_t_6) break;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":814

Output padding bytes

while offset[0] < new_offset: f[0] = 120 # "x"; pad byte # <<<<<<<<<<<<<< f += 1 offset[0] += 1

      (__pyx_v_f[0]) = 120;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":815 while offset[0] < new_offset: f[0] = 120 # "x"; pad byte f += 1 # <<<<<<<<<<<<<< offset[0] += 1

      __pyx_v_f = (__pyx_v_f + 1);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":816 f[0] = 120 # "x"; pad byte f += 1 offset[0] += 1 # <<<<<<<<<<<<<<

offset[0] += child.itemsize

      __pyx_t_8 = 0;
      (__pyx_v_offset[__pyx_t_8]) = ((__pyx_v_offset[__pyx_t_8]) + 1);
    }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":818 offset[0] += 1

offset[0] += child.itemsize             # <<<<<<<<<<<<<<

if not PyDataType_HASFIELDS(child):

    __pyx_t_8 = 0;
    (__pyx_v_offset[__pyx_t_8]) = ((__pyx_v_offset[__pyx_t_8]) + __pyx_v_child->elsize);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":820 offset[0] += child.itemsize

if not PyDataType_HASFIELDS(child)
:   # <<<<<<<<<<<<<<

    t = child.type_num
    if end - f < 5:

    __pyx_t_6 = ((!(PyDataType_HASFIELDS(__pyx_v_child) != 0)) != 0);
    if (__pyx_t_6) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":821

if not PyDataType_HASFIELDS(child): t = child.type_num # <<<<<<<<<<<<<< if end - f < 5: raise RuntimeError(u"Format string allocated too short.")

      __pyx_t_4 = __Pyx_PyInt_From_int(__pyx_v_child->type_num); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 821; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_4);
      __Pyx_XDECREF_SET(__pyx_v_t, __pyx_t_4);
      __pyx_t_4 = 0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":822 if not PyDataType_HASFIELDS(child): t = child.type_num if end - f < 5 : # <<<<<<<<<<<<<<

      raise RuntimeError(u"Format string allocated too short.")

      __pyx_t_6 = (((__pyx_v_end - __pyx_v_f) < 5) != 0);
      if (__pyx_t_6) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":823 t = child.type_num if end - f < 5: raise RuntimeError(u"Format string allocated too short.") # <<<<<<<<<<<<<<

# Until ticket #99 is fixed, use integers to avoid warnings

        __pyx_t_4 = __Pyx_PyObject_Call(__pyx_builtin_RuntimeError, __pyx_tuple__6, NULL); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 823; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_4);
        __Pyx_Raise(__pyx_t_4, 0, 0, 0);
        __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
        {__pyx_filename = __pyx_f[1]; __pyx_lineno = 823; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":826

  # Until ticket #99 is fixed, use integers to avoid warnings
  if   t == NPY_BYTE
  :   f[0] =  98 #"b"             # <<<<<<<<<<<<<<

  elif t == NPY_UBYTE
  :   f[0] =  66 #"B"

  elif t == NPY_SHORT
  :   f[0] = 104 #"h"

      __pyx_t_4 = PyInt_FromLong(NPY_BYTE); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 826; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_4);
      __pyx_t_3 = PyObject_RichCompare(__pyx_v_t, __pyx_t_4, Py_EQ); __Pyx_XGOTREF(__pyx_t_3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 826; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 826; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 98;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":827 # Until ticket #99 is fixed, use integers to avoid warnings if t == NPY_BYTE : f[0] = 98 #"b"

  elif t == NPY_UBYTE
  :   f[0] =  66 #"B"             # <<<<<<<<<<<<<<

  elif t == NPY_SHORT
  :   f[0] = 104 #"h"

  elif t == NPY_USHORT
  :   f[0] =  72 #"H"

      __pyx_t_3 = PyInt_FromLong(NPY_UBYTE); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 827; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_4 = PyObject_RichCompare(__pyx_v_t, __pyx_t_3, Py_EQ); __Pyx_XGOTREF(__pyx_t_4); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 827; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_4); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 827; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 66;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":828 if t == NPY_BYTE : f[0] = 98 #"b"

  elif t == NPY_UBYTE
  :   f[0] =  66 #"B"

  elif t == NPY_SHORT
  :   f[0] = 104 #"h"             # <<<<<<<<<<<<<<

  elif t == NPY_USHORT
  :   f[0] =  72 #"H"

  elif t == NPY_INT
  :   f[0] = 105 #"i"

      __pyx_t_4 = PyInt_FromLong(NPY_SHORT); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 828; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_4);
      __pyx_t_3 = PyObject_RichCompare(__pyx_v_t, __pyx_t_4, Py_EQ); __Pyx_XGOTREF(__pyx_t_3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 828; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 828; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 104;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":829 elif t == NPY_UBYTE : f[0] = 66 #"B"

  elif t == NPY_SHORT
  :   f[0] = 104 #"h"

  elif t == NPY_USHORT
  :   f[0] =  72 #"H"             # <<<<<<<<<<<<<<

  elif t == NPY_INT
  :   f[0] = 105 #"i"

  elif t == NPY_UINT
  :   f[0] =  73 #"I"

      __pyx_t_3 = PyInt_FromLong(NPY_USHORT); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 829; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_4 = PyObject_RichCompare(__pyx_v_t, __pyx_t_3, Py_EQ); __Pyx_XGOTREF(__pyx_t_4); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 829; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_4); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 829; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 72;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":830 elif t == NPY_SHORT : f[0] = 104 #"h"

  elif t == NPY_USHORT
  :   f[0] =  72 #"H"

  elif t == NPY_INT
  :   f[0] = 105 #"i"             # <<<<<<<<<<<<<<

  elif t == NPY_UINT
  :   f[0] =  73 #"I"

  elif t == NPY_LONG
  :   f[0] = 108 #"l"

      __pyx_t_4 = PyInt_FromLong(NPY_INT); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 830; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_4);
      __pyx_t_3 = PyObject_RichCompare(__pyx_v_t, __pyx_t_4, Py_EQ); __Pyx_XGOTREF(__pyx_t_3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 830; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 830; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 105;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":831 elif t == NPY_USHORT : f[0] = 72 #"H"

  elif t == NPY_INT
  :   f[0] = 105 #"i"

  elif t == NPY_UINT
  :   f[0] =  73 #"I"             # <<<<<<<<<<<<<<

  elif t == NPY_LONG
  :   f[0] = 108 #"l"

  elif t == NPY_ULONG
  :   f[0] = 76  #"L"

      __pyx_t_3 = PyInt_FromLong(NPY_UINT); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 831; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_4 = PyObject_RichCompare(__pyx_v_t, __pyx_t_3, Py_EQ); __Pyx_XGOTREF(__pyx_t_4); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 831; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_4); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 831; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 73;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":832 elif t == NPY_INT : f[0] = 105 #"i"

  elif t == NPY_UINT
  :   f[0] =  73 #"I"

  elif t == NPY_LONG
  :   f[0] = 108 #"l"             # <<<<<<<<<<<<<<

  elif t == NPY_ULONG
  :   f[0] = 76  #"L"

  elif t == NPY_LONGLONG
  :   f[0] = 113 #"q"

      __pyx_t_4 = PyInt_FromLong(NPY_LONG); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 832; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_4);
      __pyx_t_3 = PyObject_RichCompare(__pyx_v_t, __pyx_t_4, Py_EQ); __Pyx_XGOTREF(__pyx_t_3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 832; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 832; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 108;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":833 elif t == NPY_UINT : f[0] = 73 #"I"

  elif t == NPY_LONG
  :   f[0] = 108 #"l"

  elif t == NPY_ULONG
  :   f[0] = 76  #"L"             # <<<<<<<<<<<<<<

  elif t == NPY_LONGLONG
  :   f[0] = 113 #"q"

  elif t == NPY_ULONGLONG
  :   f[0] = 81  #"Q"

      __pyx_t_3 = PyInt_FromLong(NPY_ULONG); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 833; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_4 = PyObject_RichCompare(__pyx_v_t, __pyx_t_3, Py_EQ); __Pyx_XGOTREF(__pyx_t_4); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 833; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_4); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 833; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 76;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":834 elif t == NPY_LONG : f[0] = 108 #"l"

  elif t == NPY_ULONG
  :   f[0] = 76  #"L"

  elif t == NPY_LONGLONG
  :   f[0] = 113 #"q"             # <<<<<<<<<<<<<<

  elif t == NPY_ULONGLONG
  :   f[0] = 81  #"Q"

  elif t == NPY_FLOAT
  :   f[0] = 102 #"f"

      __pyx_t_4 = PyInt_FromLong(NPY_LONGLONG); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 834; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_4);
      __pyx_t_3 = PyObject_RichCompare(__pyx_v_t, __pyx_t_4, Py_EQ); __Pyx_XGOTREF(__pyx_t_3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 834; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 834; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 113;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":835 elif t == NPY_ULONG : f[0] = 76 #"L"

  elif t == NPY_LONGLONG
  :   f[0] = 113 #"q"

  elif t == NPY_ULONGLONG
  :   f[0] = 81  #"Q"             # <<<<<<<<<<<<<<

  elif t == NPY_FLOAT
  :   f[0] = 102 #"f"

  elif t == NPY_DOUBLE
  :   f[0] = 100 #"d"

      __pyx_t_3 = PyInt_FromLong(NPY_ULONGLONG); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 835; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_4 = PyObject_RichCompare(__pyx_v_t, __pyx_t_3, Py_EQ); __Pyx_XGOTREF(__pyx_t_4); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 835; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_4); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 835; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 81;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":836 elif t == NPY_LONGLONG : f[0] = 113 #"q"

  elif t == NPY_ULONGLONG
  :   f[0] = 81  #"Q"

  elif t == NPY_FLOAT
  :   f[0] = 102 #"f"             # <<<<<<<<<<<<<<

  elif t == NPY_DOUBLE
  :   f[0] = 100 #"d"

  elif t == NPY_LONGDOUBLE
  :   f[0] = 103 #"g"

      __pyx_t_4 = PyInt_FromLong(NPY_FLOAT); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 836; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_4);
      __pyx_t_3 = PyObject_RichCompare(__pyx_v_t, __pyx_t_4, Py_EQ); __Pyx_XGOTREF(__pyx_t_3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 836; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 836; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 102;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":837 elif t == NPY_ULONGLONG : f[0] = 81 #"Q"

  elif t == NPY_FLOAT
  :   f[0] = 102 #"f"

  elif t == NPY_DOUBLE
  :   f[0] = 100 #"d"             # <<<<<<<<<<<<<<

  elif t == NPY_LONGDOUBLE
  :   f[0] = 103 #"g"

  elif t == NPY_CFLOAT
  :   f[0] = 90; f[1] = 102; f += 1 # Zf

      __pyx_t_3 = PyInt_FromLong(NPY_DOUBLE); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 837; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_4 = PyObject_RichCompare(__pyx_v_t, __pyx_t_3, Py_EQ); __Pyx_XGOTREF(__pyx_t_4); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 837; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_4); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 837; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 100;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":838 elif t == NPY_FLOAT : f[0] = 102 #"f"

  elif t == NPY_DOUBLE
  :   f[0] = 100 #"d"

  elif t == NPY_LONGDOUBLE
  :   f[0] = 103 #"g"             # <<<<<<<<<<<<<<

  elif t == NPY_CFLOAT
  :   f[0] = 90; f[1] = 102; f += 1 # Zf

  elif t == NPY_CDOUBLE
  :   f[0] = 90; f[1] = 100; f += 1 # Zd

      __pyx_t_4 = PyInt_FromLong(NPY_LONGDOUBLE); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 838; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_4);
      __pyx_t_3 = PyObject_RichCompare(__pyx_v_t, __pyx_t_4, Py_EQ); __Pyx_XGOTREF(__pyx_t_3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 838; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 838; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 103;
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":839 elif t == NPY_DOUBLE : f[0] = 100 #"d"

  elif t == NPY_LONGDOUBLE
  :   f[0] = 103 #"g"

  elif t == NPY_CFLOAT
  :   f[0] = 90; f[1] = 102; f += 1 # Zf             # <<<<<<<<<<<<<<

  elif t == NPY_CDOUBLE
  :   f[0] = 90; f[1] = 100; f += 1 # Zd

  elif t == NPY_CLONGDOUBLE: f[0] = 90; f[1] = 103; f += 1 # Zg

      __pyx_t_3 = PyInt_FromLong(NPY_CFLOAT); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 839; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_4 = PyObject_RichCompare(__pyx_v_t, __pyx_t_3, Py_EQ); __Pyx_XGOTREF(__pyx_t_4); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 839; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_4); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 839; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 90;
        (__pyx_v_f[1]) = 102;
        __pyx_v_f = (__pyx_v_f + 1);
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":840 elif t == NPY_LONGDOUBLE : f[0] = 103 #"g"

  elif t == NPY_CFLOAT
  :   f[0] = 90; f[1] = 102; f += 1 # Zf

  elif t == NPY_CDOUBLE
  :   f[0] = 90; f[1] = 100; f += 1 # Zd             # <<<<<<<<<<<<<<

  elif t == NPY_CLONGDOUBLE: f[0] = 90; f[1] = 103; f += 1 # Zg
  elif t == NPY_OBJECT
  :   f[0] = 79 #"O"

      __pyx_t_4 = PyInt_FromLong(NPY_CDOUBLE); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 840; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_4);
      __pyx_t_3 = PyObject_RichCompare(__pyx_v_t, __pyx_t_4, Py_EQ); __Pyx_XGOTREF(__pyx_t_3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 840; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 840; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 90;
        (__pyx_v_f[1]) = 100;
        __pyx_v_f = (__pyx_v_f + 1);
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":841 elif t == NPY_CFLOAT : f[0] = 90; f[1] = 102; f += 1 # Zf

  elif t == NPY_CDOUBLE
  :   f[0] = 90; f[1] = 100; f += 1 # Zd

  elif t == NPY_CLONGDOUBLE: f[0] = 90; f[1] = 103; f += 1 # Zg             # <<<<<<<<<<<<<<
  elif t == NPY_OBJECT
  :   f[0] = 79 #"O"

  else:

      __pyx_t_3 = PyInt_FromLong(NPY_CLONGDOUBLE); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 841; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_3);
      __pyx_t_4 = PyObject_RichCompare(__pyx_v_t, __pyx_t_3, Py_EQ); __Pyx_XGOTREF(__pyx_t_4); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 841; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_4); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 841; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 90;
        (__pyx_v_f[1]) = 103;
        __pyx_v_f = (__pyx_v_f + 1);
        goto __pyx_L15;
      }

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":842 elif t == NPY_CDOUBLE : f[0] = 90; f[1] = 100; f += 1 # Zd

  elif t == NPY_CLONGDOUBLE: f[0] = 90; f[1] = 103; f += 1 # Zg
  elif t == NPY_OBJECT
  :   f[0] = 79 #"O"             # <<<<<<<<<<<<<<

  else:
      raise ValueError(u"unknown dtype code in numpy.pxd (%d)" % t)

      __pyx_t_4 = PyInt_FromLong(NPY_OBJECT); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 842; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_GOTREF(__pyx_t_4);
      __pyx_t_3 = PyObject_RichCompare(__pyx_v_t, __pyx_t_4, Py_EQ); __Pyx_XGOTREF(__pyx_t_3); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 842; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
      __pyx_t_6 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_6 < 0)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 842; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
      if (__pyx_t_6) {
        (__pyx_v_f[0]) = 79;
        goto __pyx_L15;
      }

#

else

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":844 elif t == NPY_OBJECT : f[0] = 79 #"O"

else:
    raise ValueError(u"unknown dtype code in numpy.pxd (%d)" % t)             # <<<<<<<<<<<<<<
f += 1

else:

        __pyx_t_3 = PyUnicode_Format(__pyx_kp_u_unknown_dtype_code_in_numpy_pxd, __pyx_v_t); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 844; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_3);
        __pyx_t_4 = PyTuple_New(1); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 844; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_4);
        PyTuple_SET_ITEM(__pyx_t_4, 0, __pyx_t_3);
        __Pyx_GIVEREF(__pyx_t_3);
        __pyx_t_3 = 0;
        __pyx_t_3 = __Pyx_PyObject_Call(__pyx_builtin_ValueError, __pyx_t_4, NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 844; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_3);
        __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;
        __Pyx_Raise(__pyx_t_3, 0, 0, 0);
        __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;
        {__pyx_filename = __pyx_f[1]; __pyx_lineno = 844; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      }
      __pyx_L15:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":845 else: raise ValueError(u"unknown dtype code in numpy.pxd (%d)" % t) f += 1 # <<<<<<<<<<<<<< else:

Cython ignores struct boundary information ("T{...}"),

      __pyx_v_f = (__pyx_v_f + 1);
      goto __pyx_L13;
    }

#

else

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":849

Cython ignores struct boundary information ("T{...}"),

so don't output it

  f = _util_dtypestring(child, f, end, offset)             # <<<<<<<<<<<<<<
 return f

      __pyx_t_9 = __pyx_f_5numpy__util_dtypestring(__pyx_v_child, __pyx_v_f, __pyx_v_end, __pyx_v_offset); if (unlikely(__pyx_t_9 == NULL)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 849; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
      __pyx_v_f = __pyx_t_9;
    }
    __pyx_L13:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":794 cdef tuple fields

for childname in descr.names: <<<<<<<<<<<<<<

fields = descr.fields[childname] child, new_offset = fields

  }
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":850

so don't output it

      f = _util_dtypestring(child, f, end, offset)

return f # <<<<<<<<<<<<<<

  __pyx_r = __pyx_v_f;
  goto __pyx_L0;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":783 return PyArray_MultiIterNew(5, a, b, c, d, e)

cdef inline char _util_dtypestring(dtype descr, char f, char end, int offset) except NULL: <<<<<<<<<<<<<<

Recursive utility function used in getbuffer to get format

string. The new location in the format string is returned.

#

function exit code

  __pyx_L1_error:;
  __Pyx_XDECREF(__pyx_t_1);
  __Pyx_XDECREF(__pyx_t_3);
  __Pyx_XDECREF(__pyx_t_4);
  __Pyx_AddTraceback("numpy._util_dtypestring", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = NULL;
  __pyx_L0:;
  __Pyx_XDECREF((PyObject *)__pyx_v_child);
  __Pyx_XDECREF(__pyx_v_fields);
  __Pyx_XDECREF(__pyx_v_childname);
  __Pyx_XDECREF(__pyx_v_new_offset);
  __Pyx_XDECREF(__pyx_v_t);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":966

cdef inline void set_array_base(ndarray arr, object base): <<<<<<<<<<<<<<

cdef PyObject* baseptr if base is None:

static CYTHON_INLINE void __pyx_f_5numpy_set_array_base(PyArrayObject *__pyx_v_arr, PyObject *__pyx_v_base) {
  PyObject *__pyx_v_baseptr;
  __Pyx_RefNannyDeclarations

  int __pyx_t_1;

#

  int __pyx_t_2;
  __Pyx_RefNannySetupContext("set_array_base", 0);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":968
cdef inline void set_array_base(ndarray arr, object base):
cdef PyObject* baseptr
if base is None: <<<<<<<<<<<<<<

baseptr = NULL else:

  __pyx_t_1 = (__pyx_v_base == Py_None);
  __pyx_t_2 = (__pyx_t_1 != 0);
  if (__pyx_t_2) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":969 cdef PyObject* baseptr if base is None: baseptr = NULL # <<<<<<<<<<<<<< else: Py_INCREF(base) # important to do this before decref below!

    __pyx_v_baseptr = NULL;
    goto __pyx_L3;
  }

#

else

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":971 baseptr = NULL else: Py_INCREF(base) # important to do this before decref below! # <<<<<<<<<<<<<< baseptr = base Py_XDECREF(arr.base)

    Py_INCREF(__pyx_v_base);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":972 else: Py_INCREF(base) # important to do this before decref below! baseptr = base # <<<<<<<<<<<<<< Py_XDECREF(arr.base) arr.base = baseptr

    __pyx_v_baseptr = ((PyObject *)__pyx_v_base);
  }
  __pyx_L3:;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":973 Py_INCREF(base) # important to do this before decref below! baseptr = base Py_XDECREF(arr.base) # <<<<<<<<<<<<<< arr.base = baseptr

  Py_XDECREF(__pyx_v_arr->base);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":974 baseptr = base Py_XDECREF(arr.base) arr.base = baseptr # <<<<<<<<<<<<<<

cdef inline object get_array_base(ndarray arr):

  __pyx_v_arr->base = __pyx_v_baseptr;

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":966

cdef inline void set_array_base(ndarray arr, object base): <<<<<<<<<<<<<<

cdef PyObject* baseptr if base is None:

#

function exit code

  __Pyx_RefNannyFinishContext();
}

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":976 arr.base = baseptr

cdef inline object get_array_base(ndarray arr): <<<<<<<<<<<<<<

if arr.base is NULL: return None

static CYTHON_INLINE PyObject *__pyx_f_5numpy_get_array_base(PyArrayObject *__pyx_v_arr) {
  PyObject *__pyx_r = NULL;
  __Pyx_RefNannyDeclarations

  int __pyx_t_1;
  __Pyx_RefNannySetupContext("get_array_base", 0);

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":977

cdef inline object get_array_base(ndarray arr):
if arr.base is NULL: <<<<<<<<<<<<<<

return None else:

  __pyx_t_1 = ((__pyx_v_arr->base == NULL) != 0);
  if (__pyx_t_1) {

#

"../../../../../../../usr/local/lib/python2.7/site-packages/Cython/Includes/numpy/init.pxd":978 cdef inline object get_array_base(ndarray arr): if arr.base is NULL: return None # <<<<<<<<<<<<<< else: return