/* Generated by Cython 0.9.8.1.1 on Tue Sep 23 12:23:37 2008 */

#define PY_SSIZE_T_CLEAN
#include "Python.h"
#include "structmember.h"
#ifndef PY_LONG_LONG
  #define PY_LONG_LONG LONG_LONG
#endif
#ifndef DL_EXPORT
  #define DL_EXPORT(t) t
#endif
#if PY_VERSION_HEX < 0x02040000
  #define METH_COEXIST 0
#endif
#if PY_VERSION_HEX < 0x02050000
  typedef int Py_ssize_t;
  #define PY_SSIZE_T_MAX INT_MAX
  #define PY_SSIZE_T_MIN INT_MIN
  #define PyInt_FromSsize_t(z) PyInt_FromLong(z)
  #define PyInt_AsSsize_t(o)   PyInt_AsLong(o)
  #define PyNumber_Index(o)    PyNumber_Int(o)
  #define PyIndex_Check(o)     PyNumber_Check(o)
#endif
#if PY_VERSION_HEX < 0x02060000
  #define Py_REFCNT(ob) (((PyObject*)(ob))->ob_refcnt)
  #define Py_TYPE(ob)   (((PyObject*)(ob))->ob_type)
  #define Py_SIZE(ob)   (((PyVarObject*)(ob))->ob_size)
  #define PyVarObject_HEAD_INIT(type, size) \
          PyObject_HEAD_INIT(type) size,
  #define PyType_Modified(t)

  typedef struct {
       void *buf;
       Py_ssize_t len;
       int readonly;
       const char *format;
       int ndim;
       Py_ssize_t *shape;
       Py_ssize_t *strides;
       Py_ssize_t *suboffsets;
       Py_ssize_t itemsize;
       void *internal;
  } Py_buffer;

  #define PyBUF_SIMPLE 0
  #define PyBUF_WRITABLE 0x0001
  #define PyBUF_LOCK 0x0002
  #define PyBUF_FORMAT 0x0004
  #define PyBUF_ND 0x0008
  #define PyBUF_STRIDES (0x0010 | PyBUF_ND)
  #define PyBUF_C_CONTIGUOUS (0x0020 | PyBUF_STRIDES)
  #define PyBUF_F_CONTIGUOUS (0x0040 | PyBUF_STRIDES)
  #define PyBUF_ANY_CONTIGUOUS (0x0080 | PyBUF_STRIDES)
  #define PyBUF_INDIRECT (0x0100 | PyBUF_STRIDES)

#endif
#if PY_MAJOR_VERSION < 3
  #define __Pyx_BUILTIN_MODULE_NAME "__builtin__"
#else
  #define __Pyx_BUILTIN_MODULE_NAME "builtins"
#endif
#if PY_MAJOR_VERSION >= 3
  #define Py_TPFLAGS_CHECKTYPES 0
  #define Py_TPFLAGS_HAVE_INDEX 0
#endif
#if PY_MAJOR_VERSION >= 3
  #define PyBaseString_Type            PyUnicode_Type
  #define PyString_Type                PyBytes_Type
  #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 __Pyx_PyNumber_Divide(x,y)         PyNumber_TrueDivide(x,y)
#else
  #define __Pyx_PyNumber_Divide(x,y)         PyNumber_Divide(x,y)
  #define PyBytes_Type                 PyString_Type
#endif
#if PY_MAJOR_VERSION >= 3
  #define PyMethod_New(func, self, klass) PyInstanceMethod_New(func)
#endif
#if !defined(WIN32) && !defined(MS_WINDOWS)
  #ifndef __stdcall
    #define __stdcall
  #endif
  #ifndef __cdecl
    #define __cdecl
  #endif
#else
  #define _USE_MATH_DEFINES
#endif
#ifdef __cplusplus
#define __PYX_EXTERN_C extern "C"
#else
#define __PYX_EXTERN_C extern
#endif
#include <math.h>
#define __PYX_HAVE_API__yt__lagos__Interpolators
#include "numpy/arrayobject.h"


#ifdef __GNUC__
#define INLINE __inline__
#elif _WIN32
#define INLINE __inline
#else
#define INLINE 
#endif

typedef struct {PyObject **p; char *s; long n; char is_unicode; char intern; char is_identifier;} __Pyx_StringTabEntry; /*proto*/



static int __pyx_skip_dispatch = 0;


/* Type Conversion Predeclarations */

#if PY_MAJOR_VERSION < 3
#define __Pyx_PyBytes_FromString PyString_FromString
#define __Pyx_PyBytes_AsString   PyString_AsString
#else
#define __Pyx_PyBytes_FromString PyBytes_FromString
#define __Pyx_PyBytes_AsString   PyBytes_AsString
#endif

#define __Pyx_PyBool_FromLong(b) ((b) ? (Py_INCREF(Py_True), Py_True) : (Py_INCREF(Py_False), Py_False))
static INLINE int __Pyx_PyObject_IsTrue(PyObject* x);
static INLINE PY_LONG_LONG __pyx_PyInt_AsLongLong(PyObject* x);
static INLINE unsigned PY_LONG_LONG __pyx_PyInt_AsUnsignedLongLong(PyObject* x);
static INLINE Py_ssize_t __pyx_PyIndex_AsSsize_t(PyObject* b);

#define __pyx_PyInt_AsLong(x) (PyInt_CheckExact(x) ? PyInt_AS_LONG(x) : PyInt_AsLong(x))
#define __pyx_PyFloat_AsDouble(x) (PyFloat_CheckExact(x) ? PyFloat_AS_DOUBLE(x) : PyFloat_AsDouble(x))

static INLINE unsigned char __pyx_PyInt_unsigned_char(PyObject* x);
static INLINE unsigned short __pyx_PyInt_unsigned_short(PyObject* x);
static INLINE char __pyx_PyInt_char(PyObject* x);
static INLINE short __pyx_PyInt_short(PyObject* x);
static INLINE int __pyx_PyInt_int(PyObject* x);
static INLINE long __pyx_PyInt_long(PyObject* x);
static INLINE signed char __pyx_PyInt_signed_char(PyObject* x);
static INLINE signed short __pyx_PyInt_signed_short(PyObject* x);
static INLINE signed int __pyx_PyInt_signed_int(PyObject* x);
static INLINE signed long __pyx_PyInt_signed_long(PyObject* x);
static INLINE long double __pyx_PyInt_long_double(PyObject* x);
#ifdef __GNUC__
/* Test for GCC > 2.95 */
#if __GNUC__ > 2 ||               (__GNUC__ == 2 && (__GNUC_MINOR__ > 95)) 
#define likely(x)   __builtin_expect(!!(x), 1)
#define unlikely(x) __builtin_expect(!!(x), 0)
#else /* __GNUC__ > 2 ... */
#define likely(x)   (x)
#define unlikely(x) (x)
#endif /* __GNUC__ > 2 ... */
#else /* __GNUC__ */
#define likely(x)   (x)
#define unlikely(x) (x)
#endif /* __GNUC__ */
    
static PyObject *__pyx_m;
static PyObject *__pyx_b;
static PyObject *__pyx_empty_tuple;
static int __pyx_lineno;
static int __pyx_clineno = 0;
static const char * __pyx_cfilenm= __FILE__;
static const char *__pyx_filename;
static const char **__pyx_f;

static char __pyx_mdoc[] = "\nSimple interpolators\n\nAuthor: Matthew Turk <matthewturk@gmail.com>\nAffiliation: KIPAC/SLAC/Stanford\nHomepage: http://yt.enzotools.org/\nLicense:\n  Copyright (C) 2008 Matthew Turk.  All Rights Reserved.\n\n  This file is part of yt.\n\n  yt is free software; you can redistribute it and/or modify\n  it under the terms of the GNU General Public License as published by\n  the Free Software Foundation; either version 3 of the License, or\n  (at your option) any later version.\n\n  This program is distributed in the hope that it will be useful,\n  but WITHOUT ANY WARRANTY; without even the implied warranty of\n  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the\n  GNU General Public License for more details.\n\n  You should have received a copy of the GNU General Public License\n  along with this program.  If not, see <http://www.gnu.org/licenses/>.\n";
static INLINE void __Pyx_SafeReleaseBuffer(PyObject* obj, Py_buffer* info);
static INLINE void __Pyx_ZeroBuffer(Py_buffer* buf); /*proto*/
static INLINE const char* __Pyx_ConsumeWhitespace(const char* ts); /*proto*/
static INLINE const char* __Pyx_BufferTypestringCheckEndian(const char* ts); /*proto*/
static void __Pyx_BufferNdimError(Py_buffer* buffer, int expected_ndim); /*proto*/
static const char* __Pyx_BufferTypestringCheck_item_nn___pyx_t_5numpy_float64_t(const char* ts); /*proto*/

static int __Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t(PyObject* obj, Py_buffer* buf, int flags, int nd); /*proto*/
static const char* __Pyx_BufferTypestringCheck_item_nn___pyx_t_5numpy_int32_t(const char* ts); /*proto*/

static int __Pyx_GetBuffer_nn___pyx_t_5numpy_int32_t(PyObject* obj, Py_buffer* buf, int flags, int nd); /*proto*/
#define __Pyx_BufPtrStrided1d(buf, i0, s0) ((char*)buf + i0 * s0)
#define __Pyx_BufPtrStrided2d(buf, i0, s0, i1, s1) ((char*)buf + i0 * s0 + i1 * s1)
static const char* __Pyx_BufferTypestringCheck_item_nn___pyx_t_5numpy_int_t(const char* ts); /*proto*/

static int __Pyx_GetBuffer_nn___pyx_t_5numpy_int_t(PyObject* obj, Py_buffer* buf, int flags, int nd); /*proto*/
#define __Pyx_BufPtrStrided3d(buf, i0, s0, i1, s1, i2, s2) ((char*)buf + i0 * s0 + i1 * s1 + i2 * s2)

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

static INLINE void __Pyx_RaiseArgtupleTooLong(Py_ssize_t num_expected, Py_ssize_t num_found); /*proto*/
#if (PY_MAJOR_VERSION < 3) && !(Py_TPFLAGS_DEFAULT & Py_TPFLAGS_HAVE_NEWBUFFER)
static int __Pyx_GetBuffer(PyObject *obj, Py_buffer *view, int flags);
static void __Pyx_ReleaseBuffer(PyObject *obj, Py_buffer *view);
#else
#define __Pyx_GetBuffer PyObject_GetBuffer
#define __Pyx_ReleaseBuffer PyObject_ReleaseBuffer
#endif

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

static PyObject *__Pyx_Import(PyObject *name, PyObject *from_list); /*proto*/

static PyObject *__Pyx_GetName(PyObject *dict, PyObject *name); /*proto*/

static void __Pyx_Raise(PyObject *type, PyObject *value, PyObject *tb); /*proto*/

static PyTypeObject *__Pyx_ImportType(char *module_name, char *class_name, long size);  /*proto*/

static PyObject *__Pyx_ImportModule(char *name); /*proto*/

static void __Pyx_AddTraceback(const char *funcname); /*proto*/

static int __Pyx_InitStrings(__Pyx_StringTabEntry *t); /*proto*/

/* Type declarations */

typedef npy_int8 __pyx_t_5numpy_int8_t;

typedef npy_int16 __pyx_t_5numpy_int16_t;

typedef npy_int32 __pyx_t_5numpy_int32_t;

typedef npy_int64 __pyx_t_5numpy_int64_t;

typedef npy_uint8 __pyx_t_5numpy_uint8_t;

typedef npy_uint16 __pyx_t_5numpy_uint16_t;

typedef npy_uint32 __pyx_t_5numpy_uint32_t;

typedef npy_uint64 __pyx_t_5numpy_uint64_t;

typedef npy_float32 __pyx_t_5numpy_float32_t;

typedef npy_float64 __pyx_t_5numpy_float64_t;

typedef npy_long __pyx_t_5numpy_int_t;

typedef npy_longlong __pyx_t_5numpy_long_t;

typedef npy_ulong __pyx_t_5numpy_uint_t;

typedef npy_ulonglong __pyx_t_5numpy_ulong_t;

typedef npy_double __pyx_t_5numpy_float_t;

typedef npy_double __pyx_t_5numpy_double_t;

typedef npy_longdouble __pyx_t_5numpy_longdouble_t;
/* Module declarations from numpy */

/* Module declarations from numpy */

static PyTypeObject *__pyx_ptype_5numpy_ndarray = 0;
/* Module declarations from cython */

/* Module declarations from yt.lagos.Interpolators */



/* Implementation of yt.lagos.Interpolators */
static char __pyx_k_numpy[] = "numpy";
static PyObject *__pyx_kp_numpy;
static char __pyx_k_np[] = "np";
static PyObject *__pyx_kp_np;
static char __pyx_k___getbuffer__[] = "__getbuffer__";
static PyObject *__pyx_kp___getbuffer__;
static char __pyx_k_RuntimeError[] = "RuntimeError";
static PyObject *__pyx_kp_RuntimeError;
static char __pyx_k_ValueError[] = "ValueError";
static PyObject *__pyx_kp_ValueError;
static PyObject *__pyx_kp_1;
static PyObject *__pyx_kp_16;
static PyObject *__pyx_builtin_RuntimeError;
static PyObject *__pyx_builtin_ValueError;
static char __pyx_k_1[] = "Py_intptr_t and Py_ssize_t differs in size, numpy.pxd does not support this";
static char __pyx_k_2[] = "b";
static char __pyx_k_3[] = "B";
static char __pyx_k_4[] = "h";
static char __pyx_k_5[] = "H";
static char __pyx_k_6[] = "i";
static char __pyx_k_7[] = "I";
static char __pyx_k_8[] = "l";
static char __pyx_k_9[] = "L";
static char __pyx_k_10[] = "q";
static char __pyx_k_11[] = "Q";
static char __pyx_k_12[] = "f";
static char __pyx_k_13[] = "d";
static char __pyx_k_14[] = "g";
static char __pyx_k_15[] = "O";
static char __pyx_k_16[] = "only objects, int and float dtypes supported for ndarray buffer access so far (dtype is %d)";

/* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":31
 * 
 * @cython.boundscheck(False)
 * def UnilinearlyInterpolate(np.ndarray[np.float64_t, ndim=1] table,             # <<<<<<<<<<<<<<
 *                            np.ndarray[np.float64_t, ndim=1] x_vals,
 *                            np.ndarray[np.float64_t, ndim=1] x_bins,
 */

static PyObject *__pyx_pf_2yt_5lagos_13Interpolators_UnilinearlyInterpolate(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/
static PyObject *__pyx_pf_2yt_5lagos_13Interpolators_UnilinearlyInterpolate(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) {
  PyArrayObject *__pyx_v_table = 0;
  PyArrayObject *__pyx_v_x_vals = 0;
  PyArrayObject *__pyx_v_x_bins = 0;
  PyArrayObject *__pyx_v_x_is = 0;
  PyArrayObject *__pyx_v_output = 0;
  double __pyx_v_x;
  double __pyx_v_xp;
  double __pyx_v_xm;
  int __pyx_v_i;
  int __pyx_v_x_i;
  Py_buffer __pyx_bstruct_x_is;
  Py_ssize_t __pyx_bstride_0_x_is = 0;
  Py_ssize_t __pyx_bshape_0_x_is = 0;
  Py_buffer __pyx_bstruct_output;
  Py_ssize_t __pyx_bstride_0_output = 0;
  Py_ssize_t __pyx_bshape_0_output = 0;
  Py_buffer __pyx_bstruct_x_vals;
  Py_ssize_t __pyx_bstride_0_x_vals = 0;
  Py_ssize_t __pyx_bshape_0_x_vals = 0;
  Py_buffer __pyx_bstruct_table;
  Py_ssize_t __pyx_bstride_0_table = 0;
  Py_ssize_t __pyx_bshape_0_table = 0;
  Py_buffer __pyx_bstruct_x_bins;
  Py_ssize_t __pyx_bstride_0_x_bins = 0;
  Py_ssize_t __pyx_bshape_0_x_bins = 0;
  PyObject *__pyx_r;
  npy_intp __pyx_1;
  __pyx_t_5numpy_int32_t __pyx_2;
  __pyx_t_5numpy_float64_t __pyx_3;
  __pyx_t_5numpy_float64_t __pyx_4;
  __pyx_t_5numpy_float64_t __pyx_5;
  int __pyx_t_1;
  int __pyx_t_2;
  int __pyx_t_3;
  long __pyx_t_4;
  int __pyx_t_5;
  long __pyx_t_6;
  long __pyx_t_7;
  int __pyx_t_8;
  int __pyx_t_9;
  long __pyx_t_10;
  int __pyx_t_11;
  static char *__pyx_argnames[] = {"table","x_vals","x_bins","x_is","output",0};
  __pyx_self = __pyx_self;
  if (likely(!__pyx_kwds) && likely(PyTuple_GET_SIZE(__pyx_args) == 5)) {
    __pyx_v_table = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 0));
    __pyx_v_x_vals = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 1));
    __pyx_v_x_bins = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 2));
    __pyx_v_x_is = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 3));
    __pyx_v_output = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 4));
  }
  else {
    if (unlikely(!PyArg_ParseTupleAndKeywords(__pyx_args, __pyx_kwds, "OOOOO", __pyx_argnames, &__pyx_v_table, &__pyx_v_x_vals, &__pyx_v_x_bins, &__pyx_v_x_is, &__pyx_v_output))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
  }
  goto __pyx_L4;
  __pyx_L3_error:;
  __Pyx_AddTraceback("yt.lagos.Interpolators.UnilinearlyInterpolate");
  return NULL;
  __pyx_L4:;
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_table), __pyx_ptype_5numpy_ndarray, 1, "table", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x_vals), __pyx_ptype_5numpy_ndarray, 1, "x_vals", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x_bins), __pyx_ptype_5numpy_ndarray, 1, "x_bins", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x_is), __pyx_ptype_5numpy_ndarray, 1, "x_is", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_output), __pyx_ptype_5numpy_ndarray, 1, "output", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_table, &__pyx_bstruct_table, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_table = __pyx_bstruct_table.strides[0];
  __pyx_bshape_0_table = __pyx_bstruct_table.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_x_vals, &__pyx_bstruct_x_vals, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_x_vals = __pyx_bstruct_x_vals.strides[0];
  __pyx_bshape_0_x_vals = __pyx_bstruct_x_vals.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_x_bins, &__pyx_bstruct_x_bins, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_x_bins = __pyx_bstruct_x_bins.strides[0];
  __pyx_bshape_0_x_bins = __pyx_bstruct_x_bins.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_int32_t((PyObject*)__pyx_v_x_is, &__pyx_bstruct_x_is, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_x_is = __pyx_bstruct_x_is.strides[0];
  __pyx_bshape_0_x_is = __pyx_bstruct_x_is.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_output, &__pyx_bstruct_output, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_output = __pyx_bstruct_output.strides[0];
  __pyx_bshape_0_output = __pyx_bstruct_output.shape[0];

  /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":38
 *     cdef double x, xp, xm
 *     cdef int i, x_i, y_i
 *     for i in range(x_vals.shape[0]):             # <<<<<<<<<<<<<<
 *         x_i = x_is[i]
 *         x = x_vals[i]
 */
  __pyx_1 = (__pyx_v_x_vals->dimensions[0]);
  for (__pyx_v_i = 0; __pyx_v_i < __pyx_1; __pyx_v_i+=1) {

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":39
 *     cdef int i, x_i, y_i
 *     for i in range(x_vals.shape[0]):
 *         x_i = x_is[i]             # <<<<<<<<<<<<<<
 *         x = x_vals[i]
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 */
    __pyx_t_1 = __pyx_v_i;
    if (__pyx_t_1 < 0) __pyx_t_1 += __pyx_bshape_0_x_is;
    __pyx_2 = *((__pyx_t_5numpy_int32_t *)((__pyx_t_5numpy_int32_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_is.buf, __pyx_t_1, __pyx_bstride_0_x_is)));
    __pyx_v_x_i = __pyx_2;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":40
 *     for i in range(x_vals.shape[0]):
 *         x_i = x_is[i]
 *         x = x_vals[i]             # <<<<<<<<<<<<<<
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 */
    __pyx_t_2 = __pyx_v_i;
    if (__pyx_t_2 < 0) __pyx_t_2 += __pyx_bshape_0_x_vals;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_vals.buf, __pyx_t_2, __pyx_bstride_0_x_vals)));
    __pyx_v_x = __pyx_3;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":41
 *         x_i = x_is[i]
 *         x = x_vals[i]
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])             # <<<<<<<<<<<<<<
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 *         output[i]  = table[x_i  ] * (xm) \
 */
    __pyx_t_3 = __pyx_v_x_i;
    if (__pyx_t_3 < 0) __pyx_t_3 += __pyx_bshape_0_x_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_3, __pyx_bstride_0_x_bins)));
    __pyx_t_4 = (__pyx_v_x_i + 1);
    if (__pyx_t_4 < 0) __pyx_t_4 += __pyx_bshape_0_x_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_4, __pyx_bstride_0_x_bins)));
    __pyx_t_5 = __pyx_v_x_i;
    if (__pyx_t_5 < 0) __pyx_t_5 += __pyx_bshape_0_x_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_5, __pyx_bstride_0_x_bins)));
    __pyx_v_xp = ((__pyx_v_x - __pyx_3) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":42
 *         x = x_vals[i]
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])             # <<<<<<<<<<<<<<
 *         output[i]  = table[x_i  ] * (xm) \
 *                    + table[x_i+1] * (xp)
 */
    __pyx_t_6 = (__pyx_v_x_i + 1);
    if (__pyx_t_6 < 0) __pyx_t_6 += __pyx_bshape_0_x_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_6, __pyx_bstride_0_x_bins)));
    __pyx_t_7 = (__pyx_v_x_i + 1);
    if (__pyx_t_7 < 0) __pyx_t_7 += __pyx_bshape_0_x_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_7, __pyx_bstride_0_x_bins)));
    __pyx_t_8 = __pyx_v_x_i;
    if (__pyx_t_8 < 0) __pyx_t_8 += __pyx_bshape_0_x_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_8, __pyx_bstride_0_x_bins)));
    __pyx_v_xm = ((__pyx_3 - __pyx_v_x) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":43
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 *         output[i]  = table[x_i  ] * (xm) \             # <<<<<<<<<<<<<<
 *                    + table[x_i+1] * (xp)
 * 
 */
    __pyx_t_9 = __pyx_v_x_i;
    if (__pyx_t_9 < 0) __pyx_t_9 += __pyx_bshape_0_table;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_table.buf, __pyx_t_9, __pyx_bstride_0_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":44
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 *         output[i]  = table[x_i  ] * (xm) \
 *                    + table[x_i+1] * (xp)             # <<<<<<<<<<<<<<
 * 
 * @cython.boundscheck(False)
 */
    __pyx_t_10 = (__pyx_v_x_i + 1);
    if (__pyx_t_10 < 0) __pyx_t_10 += __pyx_bshape_0_table;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_table.buf, __pyx_t_10, __pyx_bstride_0_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":43
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 *         output[i]  = table[x_i  ] * (xm) \             # <<<<<<<<<<<<<<
 *                    + table[x_i+1] * (xp)
 * 
 */
    __pyx_t_11 = __pyx_v_i;
    if (__pyx_t_11 < 0) __pyx_t_11 += __pyx_bshape_0_output;
    *((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_output.buf, __pyx_t_11, __pyx_bstride_0_output)) = ((__pyx_3 * __pyx_v_xm) + (__pyx_4 * __pyx_v_xp));
  }

  __pyx_r = Py_None; Py_INCREF(Py_None);
  goto __pyx_L0;
  __pyx_L1_error:;
  { PyObject *__pyx_type, *__pyx_value, *__pyx_tb;
    PyErr_Fetch(&__pyx_type, &__pyx_value, &__pyx_tb);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_is, &__pyx_bstruct_x_is);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_output, &__pyx_bstruct_output);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_vals, &__pyx_bstruct_x_vals);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_table, &__pyx_bstruct_table);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_bins, &__pyx_bstruct_x_bins);
  PyErr_Restore(__pyx_type, __pyx_value, __pyx_tb);}
  __Pyx_AddTraceback("yt.lagos.Interpolators.UnilinearlyInterpolate");
  __pyx_r = NULL;
  goto __pyx_L2;
  __pyx_L0:;
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_is, &__pyx_bstruct_x_is);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_output, &__pyx_bstruct_output);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_vals, &__pyx_bstruct_x_vals);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_table, &__pyx_bstruct_table);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_bins, &__pyx_bstruct_x_bins);
  __pyx_L2:;
  return __pyx_r;
}

/* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":47
 * 
 * @cython.boundscheck(False)
 * def BilinearlyInterpolate(np.ndarray[np.float64_t, ndim=2] table,             # <<<<<<<<<<<<<<
 *                           np.ndarray[np.float64_t, ndim=1] x_vals,
 *                           np.ndarray[np.float64_t, ndim=1] y_vals,
 */

static PyObject *__pyx_pf_2yt_5lagos_13Interpolators_BilinearlyInterpolate(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/
static PyObject *__pyx_pf_2yt_5lagos_13Interpolators_BilinearlyInterpolate(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) {
  PyArrayObject *__pyx_v_table = 0;
  PyArrayObject *__pyx_v_x_vals = 0;
  PyArrayObject *__pyx_v_y_vals = 0;
  PyArrayObject *__pyx_v_x_bins = 0;
  PyArrayObject *__pyx_v_y_bins = 0;
  PyArrayObject *__pyx_v_x_is = 0;
  PyArrayObject *__pyx_v_y_is = 0;
  PyArrayObject *__pyx_v_output = 0;
  double __pyx_v_x;
  double __pyx_v_xp;
  double __pyx_v_xm;
  double __pyx_v_y;
  double __pyx_v_yp;
  double __pyx_v_ym;
  int __pyx_v_i;
  int __pyx_v_x_i;
  int __pyx_v_y_i;
  Py_buffer __pyx_bstruct_y_vals;
  Py_ssize_t __pyx_bstride_0_y_vals = 0;
  Py_ssize_t __pyx_bshape_0_y_vals = 0;
  Py_buffer __pyx_bstruct_output;
  Py_ssize_t __pyx_bstride_0_output = 0;
  Py_ssize_t __pyx_bshape_0_output = 0;
  Py_buffer __pyx_bstruct_y_is;
  Py_ssize_t __pyx_bstride_0_y_is = 0;
  Py_ssize_t __pyx_bshape_0_y_is = 0;
  Py_buffer __pyx_bstruct_y_bins;
  Py_ssize_t __pyx_bstride_0_y_bins = 0;
  Py_ssize_t __pyx_bshape_0_y_bins = 0;
  Py_buffer __pyx_bstruct_x_is;
  Py_ssize_t __pyx_bstride_0_x_is = 0;
  Py_ssize_t __pyx_bshape_0_x_is = 0;
  Py_buffer __pyx_bstruct_x_vals;
  Py_ssize_t __pyx_bstride_0_x_vals = 0;
  Py_ssize_t __pyx_bshape_0_x_vals = 0;
  Py_buffer __pyx_bstruct_table;
  Py_ssize_t __pyx_bstride_0_table = 0;
  Py_ssize_t __pyx_bstride_1_table = 0;
  Py_ssize_t __pyx_bshape_0_table = 0;
  Py_ssize_t __pyx_bshape_1_table = 0;
  Py_buffer __pyx_bstruct_x_bins;
  Py_ssize_t __pyx_bstride_0_x_bins = 0;
  Py_ssize_t __pyx_bshape_0_x_bins = 0;
  PyObject *__pyx_r;
  npy_intp __pyx_1;
  __pyx_t_5numpy_int32_t __pyx_2;
  __pyx_t_5numpy_float64_t __pyx_3;
  __pyx_t_5numpy_float64_t __pyx_4;
  __pyx_t_5numpy_float64_t __pyx_5;
  __pyx_t_5numpy_float64_t __pyx_6;
  int __pyx_t_1;
  int __pyx_t_2;
  int __pyx_t_3;
  int __pyx_t_4;
  int __pyx_t_5;
  long __pyx_t_6;
  int __pyx_t_7;
  int __pyx_t_8;
  long __pyx_t_9;
  int __pyx_t_10;
  long __pyx_t_11;
  long __pyx_t_12;
  int __pyx_t_13;
  long __pyx_t_14;
  long __pyx_t_15;
  int __pyx_t_16;
  int __pyx_t_17;
  int __pyx_t_18;
  long __pyx_t_19;
  int __pyx_t_20;
  int __pyx_t_21;
  long __pyx_t_22;
  long __pyx_t_23;
  long __pyx_t_24;
  int __pyx_t_25;
  static char *__pyx_argnames[] = {"table","x_vals","y_vals","x_bins","y_bins","x_is","y_is","output",0};
  __pyx_self = __pyx_self;
  if (likely(!__pyx_kwds) && likely(PyTuple_GET_SIZE(__pyx_args) == 8)) {
    __pyx_v_table = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 0));
    __pyx_v_x_vals = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 1));
    __pyx_v_y_vals = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 2));
    __pyx_v_x_bins = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 3));
    __pyx_v_y_bins = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 4));
    __pyx_v_x_is = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 5));
    __pyx_v_y_is = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 6));
    __pyx_v_output = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 7));
  }
  else {
    if (unlikely(!PyArg_ParseTupleAndKeywords(__pyx_args, __pyx_kwds, "OOOOOOOO", __pyx_argnames, &__pyx_v_table, &__pyx_v_x_vals, &__pyx_v_y_vals, &__pyx_v_x_bins, &__pyx_v_y_bins, &__pyx_v_x_is, &__pyx_v_y_is, &__pyx_v_output))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
  }
  goto __pyx_L4;
  __pyx_L3_error:;
  __Pyx_AddTraceback("yt.lagos.Interpolators.BilinearlyInterpolate");
  return NULL;
  __pyx_L4:;
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_table), __pyx_ptype_5numpy_ndarray, 1, "table", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x_vals), __pyx_ptype_5numpy_ndarray, 1, "x_vals", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 48; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_y_vals), __pyx_ptype_5numpy_ndarray, 1, "y_vals", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 49; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x_bins), __pyx_ptype_5numpy_ndarray, 1, "x_bins", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 50; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_y_bins), __pyx_ptype_5numpy_ndarray, 1, "y_bins", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 51; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x_is), __pyx_ptype_5numpy_ndarray, 1, "x_is", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 52; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_y_is), __pyx_ptype_5numpy_ndarray, 1, "y_is", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 53; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_output), __pyx_ptype_5numpy_ndarray, 1, "output", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 54; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_table, &__pyx_bstruct_table, PyBUF_FORMAT| PyBUF_STRIDES, 2) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_table = __pyx_bstruct_table.strides[0]; __pyx_bstride_1_table = __pyx_bstruct_table.strides[1];
  __pyx_bshape_0_table = __pyx_bstruct_table.shape[0]; __pyx_bshape_1_table = __pyx_bstruct_table.shape[1];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_x_vals, &__pyx_bstruct_x_vals, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_x_vals = __pyx_bstruct_x_vals.strides[0];
  __pyx_bshape_0_x_vals = __pyx_bstruct_x_vals.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_y_vals, &__pyx_bstruct_y_vals, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_y_vals = __pyx_bstruct_y_vals.strides[0];
  __pyx_bshape_0_y_vals = __pyx_bstruct_y_vals.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_x_bins, &__pyx_bstruct_x_bins, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_x_bins = __pyx_bstruct_x_bins.strides[0];
  __pyx_bshape_0_x_bins = __pyx_bstruct_x_bins.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_y_bins, &__pyx_bstruct_y_bins, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_y_bins = __pyx_bstruct_y_bins.strides[0];
  __pyx_bshape_0_y_bins = __pyx_bstruct_y_bins.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_int32_t((PyObject*)__pyx_v_x_is, &__pyx_bstruct_x_is, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_x_is = __pyx_bstruct_x_is.strides[0];
  __pyx_bshape_0_x_is = __pyx_bstruct_x_is.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_int32_t((PyObject*)__pyx_v_y_is, &__pyx_bstruct_y_is, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_y_is = __pyx_bstruct_y_is.strides[0];
  __pyx_bshape_0_y_is = __pyx_bstruct_y_is.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_output, &__pyx_bstruct_output, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_output = __pyx_bstruct_output.strides[0];
  __pyx_bshape_0_output = __pyx_bstruct_output.shape[0];

  /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":58
 *     cdef double y, yp, ym
 *     cdef int i, x_i, y_i
 *     for i in range(x_vals.shape[0]):             # <<<<<<<<<<<<<<
 *         x_i = x_is[i]
 *         y_i = y_is[i]
 */
  __pyx_1 = (__pyx_v_x_vals->dimensions[0]);
  for (__pyx_v_i = 0; __pyx_v_i < __pyx_1; __pyx_v_i+=1) {

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":59
 *     cdef int i, x_i, y_i
 *     for i in range(x_vals.shape[0]):
 *         x_i = x_is[i]             # <<<<<<<<<<<<<<
 *         y_i = y_is[i]
 *         x = x_vals[i]
 */
    __pyx_t_1 = __pyx_v_i;
    if (__pyx_t_1 < 0) __pyx_t_1 += __pyx_bshape_0_x_is;
    __pyx_2 = *((__pyx_t_5numpy_int32_t *)((__pyx_t_5numpy_int32_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_is.buf, __pyx_t_1, __pyx_bstride_0_x_is)));
    __pyx_v_x_i = __pyx_2;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":60
 *     for i in range(x_vals.shape[0]):
 *         x_i = x_is[i]
 *         y_i = y_is[i]             # <<<<<<<<<<<<<<
 *         x = x_vals[i]
 *         y = y_vals[i]
 */
    __pyx_t_2 = __pyx_v_i;
    if (__pyx_t_2 < 0) __pyx_t_2 += __pyx_bshape_0_y_is;
    __pyx_2 = *((__pyx_t_5numpy_int32_t *)((__pyx_t_5numpy_int32_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_is.buf, __pyx_t_2, __pyx_bstride_0_y_is)));
    __pyx_v_y_i = __pyx_2;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":61
 *         x_i = x_is[i]
 *         y_i = y_is[i]
 *         x = x_vals[i]             # <<<<<<<<<<<<<<
 *         y = y_vals[i]
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 */
    __pyx_t_3 = __pyx_v_i;
    if (__pyx_t_3 < 0) __pyx_t_3 += __pyx_bshape_0_x_vals;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_vals.buf, __pyx_t_3, __pyx_bstride_0_x_vals)));
    __pyx_v_x = __pyx_3;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":62
 *         y_i = y_is[i]
 *         x = x_vals[i]
 *         y = y_vals[i]             # <<<<<<<<<<<<<<
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 *         yp = (y - y_bins[y_i]) / (y_bins[y_i+1] - y_bins[y_i])
 */
    __pyx_t_4 = __pyx_v_i;
    if (__pyx_t_4 < 0) __pyx_t_4 += __pyx_bshape_0_y_vals;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_vals.buf, __pyx_t_4, __pyx_bstride_0_y_vals)));
    __pyx_v_y = __pyx_3;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":63
 *         x = x_vals[i]
 *         y = y_vals[i]
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])             # <<<<<<<<<<<<<<
 *         yp = (y - y_bins[y_i]) / (y_bins[y_i+1] - y_bins[y_i])
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 */
    __pyx_t_5 = __pyx_v_x_i;
    if (__pyx_t_5 < 0) __pyx_t_5 += __pyx_bshape_0_x_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_5, __pyx_bstride_0_x_bins)));
    __pyx_t_6 = (__pyx_v_x_i + 1);
    if (__pyx_t_6 < 0) __pyx_t_6 += __pyx_bshape_0_x_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_6, __pyx_bstride_0_x_bins)));
    __pyx_t_7 = __pyx_v_x_i;
    if (__pyx_t_7 < 0) __pyx_t_7 += __pyx_bshape_0_x_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_7, __pyx_bstride_0_x_bins)));
    __pyx_v_xp = ((__pyx_v_x - __pyx_3) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":64
 *         y = y_vals[i]
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 *         yp = (y - y_bins[y_i]) / (y_bins[y_i+1] - y_bins[y_i])             # <<<<<<<<<<<<<<
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])
 */
    __pyx_t_8 = __pyx_v_y_i;
    if (__pyx_t_8 < 0) __pyx_t_8 += __pyx_bshape_0_y_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_8, __pyx_bstride_0_y_bins)));
    __pyx_t_9 = (__pyx_v_y_i + 1);
    if (__pyx_t_9 < 0) __pyx_t_9 += __pyx_bshape_0_y_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_9, __pyx_bstride_0_y_bins)));
    __pyx_t_10 = __pyx_v_y_i;
    if (__pyx_t_10 < 0) __pyx_t_10 += __pyx_bshape_0_y_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_10, __pyx_bstride_0_y_bins)));
    __pyx_v_yp = ((__pyx_v_y - __pyx_3) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":65
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 *         yp = (y - y_bins[y_i]) / (y_bins[y_i+1] - y_bins[y_i])
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])             # <<<<<<<<<<<<<<
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])
 *         output[i]  = table[x_i  , y_i  ] * (xm*ym) \
 */
    __pyx_t_11 = (__pyx_v_x_i + 1);
    if (__pyx_t_11 < 0) __pyx_t_11 += __pyx_bshape_0_x_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_11, __pyx_bstride_0_x_bins)));
    __pyx_t_12 = (__pyx_v_x_i + 1);
    if (__pyx_t_12 < 0) __pyx_t_12 += __pyx_bshape_0_x_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_12, __pyx_bstride_0_x_bins)));
    __pyx_t_13 = __pyx_v_x_i;
    if (__pyx_t_13 < 0) __pyx_t_13 += __pyx_bshape_0_x_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_13, __pyx_bstride_0_x_bins)));
    __pyx_v_xm = ((__pyx_3 - __pyx_v_x) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":66
 *         yp = (y - y_bins[y_i]) / (y_bins[y_i+1] - y_bins[y_i])
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])             # <<<<<<<<<<<<<<
 *         output[i]  = table[x_i  , y_i  ] * (xm*ym) \
 *                    + table[x_i+1, y_i  ] * (xp*ym) \
 */
    __pyx_t_14 = (__pyx_v_y_i + 1);
    if (__pyx_t_14 < 0) __pyx_t_14 += __pyx_bshape_0_y_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_14, __pyx_bstride_0_y_bins)));
    __pyx_t_15 = (__pyx_v_y_i + 1);
    if (__pyx_t_15 < 0) __pyx_t_15 += __pyx_bshape_0_y_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_15, __pyx_bstride_0_y_bins)));
    __pyx_t_16 = __pyx_v_y_i;
    if (__pyx_t_16 < 0) __pyx_t_16 += __pyx_bshape_0_y_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_16, __pyx_bstride_0_y_bins)));
    __pyx_v_ym = ((__pyx_3 - __pyx_v_y) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":67
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])
 *         output[i]  = table[x_i  , y_i  ] * (xm*ym) \             # <<<<<<<<<<<<<<
 *                    + table[x_i+1, y_i  ] * (xp*ym) \
 *                    + table[x_i  , y_i+1] * (xm*yp) \
 */
    __pyx_t_17 = __pyx_v_x_i;
    __pyx_t_18 = __pyx_v_y_i;
    if (__pyx_t_17 < 0) __pyx_t_17 += __pyx_bshape_0_table;
    if (__pyx_t_18 < 0) __pyx_t_18 += __pyx_bshape_1_table;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided2d(__pyx_bstruct_table.buf, __pyx_t_17, __pyx_bstride_0_table, __pyx_t_18, __pyx_bstride_1_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":68
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])
 *         output[i]  = table[x_i  , y_i  ] * (xm*ym) \
 *                    + table[x_i+1, y_i  ] * (xp*ym) \             # <<<<<<<<<<<<<<
 *                    + table[x_i  , y_i+1] * (xm*yp) \
 *                    + table[x_i+1, y_i+1] * (xp*yp)
 */
    __pyx_t_19 = (__pyx_v_x_i + 1);
    __pyx_t_20 = __pyx_v_y_i;
    if (__pyx_t_19 < 0) __pyx_t_19 += __pyx_bshape_0_table;
    if (__pyx_t_20 < 0) __pyx_t_20 += __pyx_bshape_1_table;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided2d(__pyx_bstruct_table.buf, __pyx_t_19, __pyx_bstride_0_table, __pyx_t_20, __pyx_bstride_1_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":69
 *         output[i]  = table[x_i  , y_i  ] * (xm*ym) \
 *                    + table[x_i+1, y_i  ] * (xp*ym) \
 *                    + table[x_i  , y_i+1] * (xm*yp) \             # <<<<<<<<<<<<<<
 *                    + table[x_i+1, y_i+1] * (xp*yp)
 * 
 */
    __pyx_t_21 = __pyx_v_x_i;
    __pyx_t_22 = (__pyx_v_y_i + 1);
    if (__pyx_t_21 < 0) __pyx_t_21 += __pyx_bshape_0_table;
    if (__pyx_t_22 < 0) __pyx_t_22 += __pyx_bshape_1_table;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided2d(__pyx_bstruct_table.buf, __pyx_t_21, __pyx_bstride_0_table, __pyx_t_22, __pyx_bstride_1_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":70
 *                    + table[x_i+1, y_i  ] * (xp*ym) \
 *                    + table[x_i  , y_i+1] * (xm*yp) \
 *                    + table[x_i+1, y_i+1] * (xp*yp)             # <<<<<<<<<<<<<<
 * 
 * @cython.boundscheck(False)
 */
    __pyx_t_23 = (__pyx_v_x_i + 1);
    __pyx_t_24 = (__pyx_v_y_i + 1);
    if (__pyx_t_23 < 0) __pyx_t_23 += __pyx_bshape_0_table;
    if (__pyx_t_24 < 0) __pyx_t_24 += __pyx_bshape_1_table;
    __pyx_6 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided2d(__pyx_bstruct_table.buf, __pyx_t_23, __pyx_bstride_0_table, __pyx_t_24, __pyx_bstride_1_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":67
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])
 *         output[i]  = table[x_i  , y_i  ] * (xm*ym) \             # <<<<<<<<<<<<<<
 *                    + table[x_i+1, y_i  ] * (xp*ym) \
 *                    + table[x_i  , y_i+1] * (xm*yp) \
 */
    __pyx_t_25 = __pyx_v_i;
    if (__pyx_t_25 < 0) __pyx_t_25 += __pyx_bshape_0_output;
    *((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_output.buf, __pyx_t_25, __pyx_bstride_0_output)) = ((((__pyx_3 * (__pyx_v_xm * __pyx_v_ym)) + (__pyx_4 * (__pyx_v_xp * __pyx_v_ym))) + (__pyx_5 * (__pyx_v_xm * __pyx_v_yp))) + (__pyx_6 * (__pyx_v_xp * __pyx_v_yp)));
  }

  __pyx_r = Py_None; Py_INCREF(Py_None);
  goto __pyx_L0;
  __pyx_L1_error:;
  { PyObject *__pyx_type, *__pyx_value, *__pyx_tb;
    PyErr_Fetch(&__pyx_type, &__pyx_value, &__pyx_tb);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_vals, &__pyx_bstruct_y_vals);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_output, &__pyx_bstruct_output);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_is, &__pyx_bstruct_y_is);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_bins, &__pyx_bstruct_y_bins);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_is, &__pyx_bstruct_x_is);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_vals, &__pyx_bstruct_x_vals);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_table, &__pyx_bstruct_table);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_bins, &__pyx_bstruct_x_bins);
  PyErr_Restore(__pyx_type, __pyx_value, __pyx_tb);}
  __Pyx_AddTraceback("yt.lagos.Interpolators.BilinearlyInterpolate");
  __pyx_r = NULL;
  goto __pyx_L2;
  __pyx_L0:;
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_vals, &__pyx_bstruct_y_vals);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_output, &__pyx_bstruct_output);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_is, &__pyx_bstruct_y_is);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_bins, &__pyx_bstruct_y_bins);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_is, &__pyx_bstruct_x_is);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_vals, &__pyx_bstruct_x_vals);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_table, &__pyx_bstruct_table);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_bins, &__pyx_bstruct_x_bins);
  __pyx_L2:;
  return __pyx_r;
}

/* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":73
 * 
 * @cython.boundscheck(False)
 * def TrilinearlyInterpolate(np.ndarray[np.float64_t, ndim=3] table,             # <<<<<<<<<<<<<<
 *                            np.ndarray[np.float64_t, ndim=1] x_vals,
 *                            np.ndarray[np.float64_t, ndim=1] y_vals,
 */

static PyObject *__pyx_pf_2yt_5lagos_13Interpolators_TrilinearlyInterpolate(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/
static PyObject *__pyx_pf_2yt_5lagos_13Interpolators_TrilinearlyInterpolate(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) {
  PyArrayObject *__pyx_v_table = 0;
  PyArrayObject *__pyx_v_x_vals = 0;
  PyArrayObject *__pyx_v_y_vals = 0;
  PyArrayObject *__pyx_v_z_vals = 0;
  PyArrayObject *__pyx_v_x_bins = 0;
  PyArrayObject *__pyx_v_y_bins = 0;
  PyArrayObject *__pyx_v_z_bins = 0;
  PyArrayObject *__pyx_v_x_is = 0;
  PyArrayObject *__pyx_v_y_is = 0;
  PyArrayObject *__pyx_v_z_is = 0;
  PyArrayObject *__pyx_v_output = 0;
  double __pyx_v_x;
  double __pyx_v_xp;
  double __pyx_v_xm;
  double __pyx_v_y;
  double __pyx_v_yp;
  double __pyx_v_ym;
  double __pyx_v_z;
  double __pyx_v_zp;
  double __pyx_v_zm;
  int __pyx_v_i;
  int __pyx_v_x_i;
  int __pyx_v_y_i;
  int __pyx_v_z_i;
  Py_buffer __pyx_bstruct_z_is;
  Py_ssize_t __pyx_bstride_0_z_is = 0;
  Py_ssize_t __pyx_bshape_0_z_is = 0;
  Py_buffer __pyx_bstruct_x_vals;
  Py_ssize_t __pyx_bstride_0_x_vals = 0;
  Py_ssize_t __pyx_bshape_0_x_vals = 0;
  Py_buffer __pyx_bstruct_y_is;
  Py_ssize_t __pyx_bstride_0_y_is = 0;
  Py_ssize_t __pyx_bshape_0_y_is = 0;
  Py_buffer __pyx_bstruct_y_bins;
  Py_ssize_t __pyx_bstride_0_y_bins = 0;
  Py_ssize_t __pyx_bshape_0_y_bins = 0;
  Py_buffer __pyx_bstruct_z_bins;
  Py_ssize_t __pyx_bstride_0_z_bins = 0;
  Py_ssize_t __pyx_bshape_0_z_bins = 0;
  Py_buffer __pyx_bstruct_x_is;
  Py_ssize_t __pyx_bstride_0_x_is = 0;
  Py_ssize_t __pyx_bshape_0_x_is = 0;
  Py_buffer __pyx_bstruct_output;
  Py_ssize_t __pyx_bstride_0_output = 0;
  Py_ssize_t __pyx_bshape_0_output = 0;
  Py_buffer __pyx_bstruct_table;
  Py_ssize_t __pyx_bstride_0_table = 0;
  Py_ssize_t __pyx_bstride_1_table = 0;
  Py_ssize_t __pyx_bstride_2_table = 0;
  Py_ssize_t __pyx_bshape_0_table = 0;
  Py_ssize_t __pyx_bshape_1_table = 0;
  Py_ssize_t __pyx_bshape_2_table = 0;
  Py_buffer __pyx_bstruct_x_bins;
  Py_ssize_t __pyx_bstride_0_x_bins = 0;
  Py_ssize_t __pyx_bshape_0_x_bins = 0;
  Py_buffer __pyx_bstruct_y_vals;
  Py_ssize_t __pyx_bstride_0_y_vals = 0;
  Py_ssize_t __pyx_bshape_0_y_vals = 0;
  Py_buffer __pyx_bstruct_z_vals;
  Py_ssize_t __pyx_bstride_0_z_vals = 0;
  Py_ssize_t __pyx_bshape_0_z_vals = 0;
  PyObject *__pyx_r;
  npy_intp __pyx_1;
  __pyx_t_5numpy_int_t __pyx_2;
  __pyx_t_5numpy_float64_t __pyx_3;
  __pyx_t_5numpy_float64_t __pyx_4;
  __pyx_t_5numpy_float64_t __pyx_5;
  __pyx_t_5numpy_float64_t __pyx_6;
  __pyx_t_5numpy_float64_t __pyx_7;
  __pyx_t_5numpy_float64_t __pyx_8;
  __pyx_t_5numpy_float64_t __pyx_9;
  __pyx_t_5numpy_float64_t __pyx_10;
  int __pyx_t_1;
  int __pyx_t_2;
  int __pyx_t_3;
  int __pyx_t_4;
  int __pyx_t_5;
  int __pyx_t_6;
  int __pyx_t_7;
  long __pyx_t_8;
  int __pyx_t_9;
  int __pyx_t_10;
  long __pyx_t_11;
  int __pyx_t_12;
  int __pyx_t_13;
  long __pyx_t_14;
  int __pyx_t_15;
  long __pyx_t_16;
  long __pyx_t_17;
  int __pyx_t_18;
  long __pyx_t_19;
  long __pyx_t_20;
  int __pyx_t_21;
  long __pyx_t_22;
  long __pyx_t_23;
  int __pyx_t_24;
  int __pyx_t_25;
  int __pyx_t_26;
  int __pyx_t_27;
  long __pyx_t_28;
  int __pyx_t_29;
  int __pyx_t_30;
  int __pyx_t_31;
  long __pyx_t_32;
  int __pyx_t_33;
  int __pyx_t_34;
  int __pyx_t_35;
  long __pyx_t_36;
  long __pyx_t_37;
  int __pyx_t_38;
  long __pyx_t_39;
  int __pyx_t_40;
  long __pyx_t_41;
  long __pyx_t_42;
  long __pyx_t_43;
  long __pyx_t_44;
  int __pyx_t_45;
  long __pyx_t_46;
  long __pyx_t_47;
  long __pyx_t_48;
  int __pyx_t_49;
  static char *__pyx_argnames[] = {"table","x_vals","y_vals","z_vals","x_bins","y_bins","z_bins","x_is","y_is","z_is","output",0};
  __pyx_self = __pyx_self;
  if (likely(!__pyx_kwds) && likely(PyTuple_GET_SIZE(__pyx_args) == 11)) {
    __pyx_v_table = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 0));
    __pyx_v_x_vals = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 1));
    __pyx_v_y_vals = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 2));
    __pyx_v_z_vals = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 3));
    __pyx_v_x_bins = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 4));
    __pyx_v_y_bins = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 5));
    __pyx_v_z_bins = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 6));
    __pyx_v_x_is = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 7));
    __pyx_v_y_is = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 8));
    __pyx_v_z_is = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 9));
    __pyx_v_output = ((PyArrayObject *)PyTuple_GET_ITEM(__pyx_args, 10));
  }
  else {
    if (unlikely(!PyArg_ParseTupleAndKeywords(__pyx_args, __pyx_kwds, "OOOOOOOOOOO", __pyx_argnames, &__pyx_v_table, &__pyx_v_x_vals, &__pyx_v_y_vals, &__pyx_v_z_vals, &__pyx_v_x_bins, &__pyx_v_y_bins, &__pyx_v_z_bins, &__pyx_v_x_is, &__pyx_v_y_is, &__pyx_v_z_is, &__pyx_v_output))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
  }
  goto __pyx_L4;
  __pyx_L3_error:;
  __Pyx_AddTraceback("yt.lagos.Interpolators.TrilinearlyInterpolate");
  return NULL;
  __pyx_L4:;
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_table), __pyx_ptype_5numpy_ndarray, 1, "table", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x_vals), __pyx_ptype_5numpy_ndarray, 1, "x_vals", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 74; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_y_vals), __pyx_ptype_5numpy_ndarray, 1, "y_vals", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 75; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_z_vals), __pyx_ptype_5numpy_ndarray, 1, "z_vals", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 76; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x_bins), __pyx_ptype_5numpy_ndarray, 1, "x_bins", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 77; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_y_bins), __pyx_ptype_5numpy_ndarray, 1, "y_bins", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_z_bins), __pyx_ptype_5numpy_ndarray, 1, "z_bins", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x_is), __pyx_ptype_5numpy_ndarray, 1, "x_is", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_y_is), __pyx_ptype_5numpy_ndarray, 1, "y_is", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 81; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_z_is), __pyx_ptype_5numpy_ndarray, 1, "z_is", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 82; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_output), __pyx_ptype_5numpy_ndarray, 1, "output", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 83; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_table, &__pyx_bstruct_table, PyBUF_FORMAT| PyBUF_STRIDES, 3) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_table = __pyx_bstruct_table.strides[0]; __pyx_bstride_1_table = __pyx_bstruct_table.strides[1]; __pyx_bstride_2_table = __pyx_bstruct_table.strides[2];
  __pyx_bshape_0_table = __pyx_bstruct_table.shape[0]; __pyx_bshape_1_table = __pyx_bstruct_table.shape[1]; __pyx_bshape_2_table = __pyx_bstruct_table.shape[2];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_x_vals, &__pyx_bstruct_x_vals, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_x_vals = __pyx_bstruct_x_vals.strides[0];
  __pyx_bshape_0_x_vals = __pyx_bstruct_x_vals.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_y_vals, &__pyx_bstruct_y_vals, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_y_vals = __pyx_bstruct_y_vals.strides[0];
  __pyx_bshape_0_y_vals = __pyx_bstruct_y_vals.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_z_vals, &__pyx_bstruct_z_vals, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_z_vals = __pyx_bstruct_z_vals.strides[0];
  __pyx_bshape_0_z_vals = __pyx_bstruct_z_vals.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_x_bins, &__pyx_bstruct_x_bins, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_x_bins = __pyx_bstruct_x_bins.strides[0];
  __pyx_bshape_0_x_bins = __pyx_bstruct_x_bins.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_y_bins, &__pyx_bstruct_y_bins, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_y_bins = __pyx_bstruct_y_bins.strides[0];
  __pyx_bshape_0_y_bins = __pyx_bstruct_y_bins.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_z_bins, &__pyx_bstruct_z_bins, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_z_bins = __pyx_bstruct_z_bins.strides[0];
  __pyx_bshape_0_z_bins = __pyx_bstruct_z_bins.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_int_t((PyObject*)__pyx_v_x_is, &__pyx_bstruct_x_is, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_x_is = __pyx_bstruct_x_is.strides[0];
  __pyx_bshape_0_x_is = __pyx_bstruct_x_is.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_int_t((PyObject*)__pyx_v_y_is, &__pyx_bstruct_y_is, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_y_is = __pyx_bstruct_y_is.strides[0];
  __pyx_bshape_0_y_is = __pyx_bstruct_y_is.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_int_t((PyObject*)__pyx_v_z_is, &__pyx_bstruct_z_is, PyBUF_FORMAT| PyBUF_STRIDES, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_z_is = __pyx_bstruct_z_is.strides[0];
  __pyx_bshape_0_z_is = __pyx_bstruct_z_is.shape[0];
  if (unlikely(__Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t((PyObject*)__pyx_v_output, &__pyx_bstruct_output, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_bstride_0_output = __pyx_bstruct_output.strides[0];
  __pyx_bshape_0_output = __pyx_bstruct_output.shape[0];

  /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":88
 *     cdef double z, zp, zm
 *     cdef int i, x_i, y_i, z_i
 *     for i in range(x_vals.shape[0]):             # <<<<<<<<<<<<<<
 *         x_i = x_is[i]
 *         y_i = y_is[i]
 */
  __pyx_1 = (__pyx_v_x_vals->dimensions[0]);
  for (__pyx_v_i = 0; __pyx_v_i < __pyx_1; __pyx_v_i+=1) {

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":89
 *     cdef int i, x_i, y_i, z_i
 *     for i in range(x_vals.shape[0]):
 *         x_i = x_is[i]             # <<<<<<<<<<<<<<
 *         y_i = y_is[i]
 *         z_i = z_is[i]
 */
    __pyx_t_1 = __pyx_v_i;
    if (__pyx_t_1 < 0) __pyx_t_1 += __pyx_bshape_0_x_is;
    __pyx_2 = *((__pyx_t_5numpy_int_t *)((__pyx_t_5numpy_int_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_is.buf, __pyx_t_1, __pyx_bstride_0_x_is)));
    __pyx_v_x_i = __pyx_2;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":90
 *     for i in range(x_vals.shape[0]):
 *         x_i = x_is[i]
 *         y_i = y_is[i]             # <<<<<<<<<<<<<<
 *         z_i = z_is[i]
 *         x = x_vals[i]
 */
    __pyx_t_2 = __pyx_v_i;
    if (__pyx_t_2 < 0) __pyx_t_2 += __pyx_bshape_0_y_is;
    __pyx_2 = *((__pyx_t_5numpy_int_t *)((__pyx_t_5numpy_int_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_is.buf, __pyx_t_2, __pyx_bstride_0_y_is)));
    __pyx_v_y_i = __pyx_2;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":91
 *         x_i = x_is[i]
 *         y_i = y_is[i]
 *         z_i = z_is[i]             # <<<<<<<<<<<<<<
 *         x = x_vals[i]
 *         y = y_vals[i]
 */
    __pyx_t_3 = __pyx_v_i;
    if (__pyx_t_3 < 0) __pyx_t_3 += __pyx_bshape_0_z_is;
    __pyx_2 = *((__pyx_t_5numpy_int_t *)((__pyx_t_5numpy_int_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_z_is.buf, __pyx_t_3, __pyx_bstride_0_z_is)));
    __pyx_v_z_i = __pyx_2;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":92
 *         y_i = y_is[i]
 *         z_i = z_is[i]
 *         x = x_vals[i]             # <<<<<<<<<<<<<<
 *         y = y_vals[i]
 *         z = z_vals[i]
 */
    __pyx_t_4 = __pyx_v_i;
    if (__pyx_t_4 < 0) __pyx_t_4 += __pyx_bshape_0_x_vals;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_vals.buf, __pyx_t_4, __pyx_bstride_0_x_vals)));
    __pyx_v_x = __pyx_3;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":93
 *         z_i = z_is[i]
 *         x = x_vals[i]
 *         y = y_vals[i]             # <<<<<<<<<<<<<<
 *         z = z_vals[i]
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 */
    __pyx_t_5 = __pyx_v_i;
    if (__pyx_t_5 < 0) __pyx_t_5 += __pyx_bshape_0_y_vals;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_vals.buf, __pyx_t_5, __pyx_bstride_0_y_vals)));
    __pyx_v_y = __pyx_3;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":94
 *         x = x_vals[i]
 *         y = y_vals[i]
 *         z = z_vals[i]             # <<<<<<<<<<<<<<
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 *         yp = (y - y_bins[y_i]) / (y_bins[y_i+1] - y_bins[y_i])
 */
    __pyx_t_6 = __pyx_v_i;
    if (__pyx_t_6 < 0) __pyx_t_6 += __pyx_bshape_0_z_vals;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_z_vals.buf, __pyx_t_6, __pyx_bstride_0_z_vals)));
    __pyx_v_z = __pyx_3;

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":95
 *         y = y_vals[i]
 *         z = z_vals[i]
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])             # <<<<<<<<<<<<<<
 *         yp = (y - y_bins[y_i]) / (y_bins[y_i+1] - y_bins[y_i])
 *         zp = (z - z_bins[z_i]) / (z_bins[z_i+1] - z_bins[z_i])
 */
    __pyx_t_7 = __pyx_v_x_i;
    if (__pyx_t_7 < 0) __pyx_t_7 += __pyx_bshape_0_x_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_7, __pyx_bstride_0_x_bins)));
    __pyx_t_8 = (__pyx_v_x_i + 1);
    if (__pyx_t_8 < 0) __pyx_t_8 += __pyx_bshape_0_x_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_8, __pyx_bstride_0_x_bins)));
    __pyx_t_9 = __pyx_v_x_i;
    if (__pyx_t_9 < 0) __pyx_t_9 += __pyx_bshape_0_x_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_9, __pyx_bstride_0_x_bins)));
    __pyx_v_xp = ((__pyx_v_x - __pyx_3) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":96
 *         z = z_vals[i]
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 *         yp = (y - y_bins[y_i]) / (y_bins[y_i+1] - y_bins[y_i])             # <<<<<<<<<<<<<<
 *         zp = (z - z_bins[z_i]) / (z_bins[z_i+1] - z_bins[z_i])
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 */
    __pyx_t_10 = __pyx_v_y_i;
    if (__pyx_t_10 < 0) __pyx_t_10 += __pyx_bshape_0_y_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_10, __pyx_bstride_0_y_bins)));
    __pyx_t_11 = (__pyx_v_y_i + 1);
    if (__pyx_t_11 < 0) __pyx_t_11 += __pyx_bshape_0_y_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_11, __pyx_bstride_0_y_bins)));
    __pyx_t_12 = __pyx_v_y_i;
    if (__pyx_t_12 < 0) __pyx_t_12 += __pyx_bshape_0_y_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_12, __pyx_bstride_0_y_bins)));
    __pyx_v_yp = ((__pyx_v_y - __pyx_3) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":97
 *         xp = (x - x_bins[x_i]) / (x_bins[x_i+1] - x_bins[x_i])
 *         yp = (y - y_bins[y_i]) / (y_bins[y_i+1] - y_bins[y_i])
 *         zp = (z - z_bins[z_i]) / (z_bins[z_i+1] - z_bins[z_i])             # <<<<<<<<<<<<<<
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])
 */
    __pyx_t_13 = __pyx_v_z_i;
    if (__pyx_t_13 < 0) __pyx_t_13 += __pyx_bshape_0_z_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_z_bins.buf, __pyx_t_13, __pyx_bstride_0_z_bins)));
    __pyx_t_14 = (__pyx_v_z_i + 1);
    if (__pyx_t_14 < 0) __pyx_t_14 += __pyx_bshape_0_z_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_z_bins.buf, __pyx_t_14, __pyx_bstride_0_z_bins)));
    __pyx_t_15 = __pyx_v_z_i;
    if (__pyx_t_15 < 0) __pyx_t_15 += __pyx_bshape_0_z_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_z_bins.buf, __pyx_t_15, __pyx_bstride_0_z_bins)));
    __pyx_v_zp = ((__pyx_v_z - __pyx_3) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":98
 *         yp = (y - y_bins[y_i]) / (y_bins[y_i+1] - y_bins[y_i])
 *         zp = (z - z_bins[z_i]) / (z_bins[z_i+1] - z_bins[z_i])
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])             # <<<<<<<<<<<<<<
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])
 *         zm = (z_bins[z_i+1] - z) / (z_bins[z_i+1] - z_bins[z_i])
 */
    __pyx_t_16 = (__pyx_v_x_i + 1);
    if (__pyx_t_16 < 0) __pyx_t_16 += __pyx_bshape_0_x_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_16, __pyx_bstride_0_x_bins)));
    __pyx_t_17 = (__pyx_v_x_i + 1);
    if (__pyx_t_17 < 0) __pyx_t_17 += __pyx_bshape_0_x_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_17, __pyx_bstride_0_x_bins)));
    __pyx_t_18 = __pyx_v_x_i;
    if (__pyx_t_18 < 0) __pyx_t_18 += __pyx_bshape_0_x_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_x_bins.buf, __pyx_t_18, __pyx_bstride_0_x_bins)));
    __pyx_v_xm = ((__pyx_3 - __pyx_v_x) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":99
 *         zp = (z - z_bins[z_i]) / (z_bins[z_i+1] - z_bins[z_i])
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])             # <<<<<<<<<<<<<<
 *         zm = (z_bins[z_i+1] - z) / (z_bins[z_i+1] - z_bins[z_i])
 *         output[i]  = table[x_i  ,y_i  ,z_i  ] * (xm*ym*zm) \
 */
    __pyx_t_19 = (__pyx_v_y_i + 1);
    if (__pyx_t_19 < 0) __pyx_t_19 += __pyx_bshape_0_y_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_19, __pyx_bstride_0_y_bins)));
    __pyx_t_20 = (__pyx_v_y_i + 1);
    if (__pyx_t_20 < 0) __pyx_t_20 += __pyx_bshape_0_y_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_20, __pyx_bstride_0_y_bins)));
    __pyx_t_21 = __pyx_v_y_i;
    if (__pyx_t_21 < 0) __pyx_t_21 += __pyx_bshape_0_y_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_y_bins.buf, __pyx_t_21, __pyx_bstride_0_y_bins)));
    __pyx_v_ym = ((__pyx_3 - __pyx_v_y) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":100
 *         xm = (x_bins[x_i+1] - x) / (x_bins[x_i+1] - x_bins[x_i])
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])
 *         zm = (z_bins[z_i+1] - z) / (z_bins[z_i+1] - z_bins[z_i])             # <<<<<<<<<<<<<<
 *         output[i]  = table[x_i  ,y_i  ,z_i  ] * (xm*ym*zm) \
 *                    + table[x_i+1,y_i  ,z_i  ] * (xp*ym*zm) \
 */
    __pyx_t_22 = (__pyx_v_z_i + 1);
    if (__pyx_t_22 < 0) __pyx_t_22 += __pyx_bshape_0_z_bins;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_z_bins.buf, __pyx_t_22, __pyx_bstride_0_z_bins)));
    __pyx_t_23 = (__pyx_v_z_i + 1);
    if (__pyx_t_23 < 0) __pyx_t_23 += __pyx_bshape_0_z_bins;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_z_bins.buf, __pyx_t_23, __pyx_bstride_0_z_bins)));
    __pyx_t_24 = __pyx_v_z_i;
    if (__pyx_t_24 < 0) __pyx_t_24 += __pyx_bshape_0_z_bins;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_z_bins.buf, __pyx_t_24, __pyx_bstride_0_z_bins)));
    __pyx_v_zm = ((__pyx_3 - __pyx_v_z) / (__pyx_4 - __pyx_5));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":101
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])
 *         zm = (z_bins[z_i+1] - z) / (z_bins[z_i+1] - z_bins[z_i])
 *         output[i]  = table[x_i  ,y_i  ,z_i  ] * (xm*ym*zm) \             # <<<<<<<<<<<<<<
 *                    + table[x_i+1,y_i  ,z_i  ] * (xp*ym*zm) \
 *                    + table[x_i  ,y_i+1,z_i  ] * (xm*yp*zm) \
 */
    __pyx_t_25 = __pyx_v_x_i;
    __pyx_t_26 = __pyx_v_y_i;
    __pyx_t_27 = __pyx_v_z_i;
    if (__pyx_t_25 < 0) __pyx_t_25 += __pyx_bshape_0_table;
    if (__pyx_t_26 < 0) __pyx_t_26 += __pyx_bshape_1_table;
    if (__pyx_t_27 < 0) __pyx_t_27 += __pyx_bshape_2_table;
    __pyx_3 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided3d(__pyx_bstruct_table.buf, __pyx_t_25, __pyx_bstride_0_table, __pyx_t_26, __pyx_bstride_1_table, __pyx_t_27, __pyx_bstride_2_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":102
 *         zm = (z_bins[z_i+1] - z) / (z_bins[z_i+1] - z_bins[z_i])
 *         output[i]  = table[x_i  ,y_i  ,z_i  ] * (xm*ym*zm) \
 *                    + table[x_i+1,y_i  ,z_i  ] * (xp*ym*zm) \             # <<<<<<<<<<<<<<
 *                    + table[x_i  ,y_i+1,z_i  ] * (xm*yp*zm) \
 *                    + table[x_i  ,y_i  ,z_i+1] * (xm*ym*zp) \
 */
    __pyx_t_28 = (__pyx_v_x_i + 1);
    __pyx_t_29 = __pyx_v_y_i;
    __pyx_t_30 = __pyx_v_z_i;
    if (__pyx_t_28 < 0) __pyx_t_28 += __pyx_bshape_0_table;
    if (__pyx_t_29 < 0) __pyx_t_29 += __pyx_bshape_1_table;
    if (__pyx_t_30 < 0) __pyx_t_30 += __pyx_bshape_2_table;
    __pyx_4 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided3d(__pyx_bstruct_table.buf, __pyx_t_28, __pyx_bstride_0_table, __pyx_t_29, __pyx_bstride_1_table, __pyx_t_30, __pyx_bstride_2_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":103
 *         output[i]  = table[x_i  ,y_i  ,z_i  ] * (xm*ym*zm) \
 *                    + table[x_i+1,y_i  ,z_i  ] * (xp*ym*zm) \
 *                    + table[x_i  ,y_i+1,z_i  ] * (xm*yp*zm) \             # <<<<<<<<<<<<<<
 *                    + table[x_i  ,y_i  ,z_i+1] * (xm*ym*zp) \
 *                    + table[x_i+1,y_i  ,z_i+1] * (xp*ym*zp) \
 */
    __pyx_t_31 = __pyx_v_x_i;
    __pyx_t_32 = (__pyx_v_y_i + 1);
    __pyx_t_33 = __pyx_v_z_i;
    if (__pyx_t_31 < 0) __pyx_t_31 += __pyx_bshape_0_table;
    if (__pyx_t_32 < 0) __pyx_t_32 += __pyx_bshape_1_table;
    if (__pyx_t_33 < 0) __pyx_t_33 += __pyx_bshape_2_table;
    __pyx_5 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided3d(__pyx_bstruct_table.buf, __pyx_t_31, __pyx_bstride_0_table, __pyx_t_32, __pyx_bstride_1_table, __pyx_t_33, __pyx_bstride_2_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":104
 *                    + table[x_i+1,y_i  ,z_i  ] * (xp*ym*zm) \
 *                    + table[x_i  ,y_i+1,z_i  ] * (xm*yp*zm) \
 *                    + table[x_i  ,y_i  ,z_i+1] * (xm*ym*zp) \             # <<<<<<<<<<<<<<
 *                    + table[x_i+1,y_i  ,z_i+1] * (xp*ym*zp) \
 *                    + table[x_i  ,y_i+1,z_i+1] * (xm*yp*zp) \
 */
    __pyx_t_34 = __pyx_v_x_i;
    __pyx_t_35 = __pyx_v_y_i;
    __pyx_t_36 = (__pyx_v_z_i + 1);
    if (__pyx_t_34 < 0) __pyx_t_34 += __pyx_bshape_0_table;
    if (__pyx_t_35 < 0) __pyx_t_35 += __pyx_bshape_1_table;
    if (__pyx_t_36 < 0) __pyx_t_36 += __pyx_bshape_2_table;
    __pyx_6 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided3d(__pyx_bstruct_table.buf, __pyx_t_34, __pyx_bstride_0_table, __pyx_t_35, __pyx_bstride_1_table, __pyx_t_36, __pyx_bstride_2_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":105
 *                    + table[x_i  ,y_i+1,z_i  ] * (xm*yp*zm) \
 *                    + table[x_i  ,y_i  ,z_i+1] * (xm*ym*zp) \
 *                    + table[x_i+1,y_i  ,z_i+1] * (xp*ym*zp) \             # <<<<<<<<<<<<<<
 *                    + table[x_i  ,y_i+1,z_i+1] * (xm*yp*zp) \
 *                    + table[x_i+1,y_i+1,z_i  ] * (xp*yp*zm) \
 */
    __pyx_t_37 = (__pyx_v_x_i + 1);
    __pyx_t_38 = __pyx_v_y_i;
    __pyx_t_39 = (__pyx_v_z_i + 1);
    if (__pyx_t_37 < 0) __pyx_t_37 += __pyx_bshape_0_table;
    if (__pyx_t_38 < 0) __pyx_t_38 += __pyx_bshape_1_table;
    if (__pyx_t_39 < 0) __pyx_t_39 += __pyx_bshape_2_table;
    __pyx_7 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided3d(__pyx_bstruct_table.buf, __pyx_t_37, __pyx_bstride_0_table, __pyx_t_38, __pyx_bstride_1_table, __pyx_t_39, __pyx_bstride_2_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":106
 *                    + table[x_i  ,y_i  ,z_i+1] * (xm*ym*zp) \
 *                    + table[x_i+1,y_i  ,z_i+1] * (xp*ym*zp) \
 *                    + table[x_i  ,y_i+1,z_i+1] * (xm*yp*zp) \             # <<<<<<<<<<<<<<
 *                    + table[x_i+1,y_i+1,z_i  ] * (xp*yp*zm) \
 *                    + table[x_i+1,y_i+1,z_i+1] * (xp*yp*zp)
 */
    __pyx_t_40 = __pyx_v_x_i;
    __pyx_t_41 = (__pyx_v_y_i + 1);
    __pyx_t_42 = (__pyx_v_z_i + 1);
    if (__pyx_t_40 < 0) __pyx_t_40 += __pyx_bshape_0_table;
    if (__pyx_t_41 < 0) __pyx_t_41 += __pyx_bshape_1_table;
    if (__pyx_t_42 < 0) __pyx_t_42 += __pyx_bshape_2_table;
    __pyx_8 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided3d(__pyx_bstruct_table.buf, __pyx_t_40, __pyx_bstride_0_table, __pyx_t_41, __pyx_bstride_1_table, __pyx_t_42, __pyx_bstride_2_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":107
 *                    + table[x_i+1,y_i  ,z_i+1] * (xp*ym*zp) \
 *                    + table[x_i  ,y_i+1,z_i+1] * (xm*yp*zp) \
 *                    + table[x_i+1,y_i+1,z_i  ] * (xp*yp*zm) \             # <<<<<<<<<<<<<<
 *                    + table[x_i+1,y_i+1,z_i+1] * (xp*yp*zp)
 */
    __pyx_t_43 = (__pyx_v_x_i + 1);
    __pyx_t_44 = (__pyx_v_y_i + 1);
    __pyx_t_45 = __pyx_v_z_i;
    if (__pyx_t_43 < 0) __pyx_t_43 += __pyx_bshape_0_table;
    if (__pyx_t_44 < 0) __pyx_t_44 += __pyx_bshape_1_table;
    if (__pyx_t_45 < 0) __pyx_t_45 += __pyx_bshape_2_table;
    __pyx_9 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided3d(__pyx_bstruct_table.buf, __pyx_t_43, __pyx_bstride_0_table, __pyx_t_44, __pyx_bstride_1_table, __pyx_t_45, __pyx_bstride_2_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":108
 *                    + table[x_i  ,y_i+1,z_i+1] * (xm*yp*zp) \
 *                    + table[x_i+1,y_i+1,z_i  ] * (xp*yp*zm) \
 *                    + table[x_i+1,y_i+1,z_i+1] * (xp*yp*zp)             # <<<<<<<<<<<<<<
 */
    __pyx_t_46 = (__pyx_v_x_i + 1);
    __pyx_t_47 = (__pyx_v_y_i + 1);
    __pyx_t_48 = (__pyx_v_z_i + 1);
    if (__pyx_t_46 < 0) __pyx_t_46 += __pyx_bshape_0_table;
    if (__pyx_t_47 < 0) __pyx_t_47 += __pyx_bshape_1_table;
    if (__pyx_t_48 < 0) __pyx_t_48 += __pyx_bshape_2_table;
    __pyx_10 = *((__pyx_t_5numpy_float64_t *)((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided3d(__pyx_bstruct_table.buf, __pyx_t_46, __pyx_bstride_0_table, __pyx_t_47, __pyx_bstride_1_table, __pyx_t_48, __pyx_bstride_2_table)));

    /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":101
 *         ym = (y_bins[y_i+1] - y) / (y_bins[y_i+1] - y_bins[y_i])
 *         zm = (z_bins[z_i+1] - z) / (z_bins[z_i+1] - z_bins[z_i])
 *         output[i]  = table[x_i  ,y_i  ,z_i  ] * (xm*ym*zm) \             # <<<<<<<<<<<<<<
 *                    + table[x_i+1,y_i  ,z_i  ] * (xp*ym*zm) \
 *                    + table[x_i  ,y_i+1,z_i  ] * (xm*yp*zm) \
 */
    __pyx_t_49 = __pyx_v_i;
    if (__pyx_t_49 < 0) __pyx_t_49 += __pyx_bshape_0_output;
    *((__pyx_t_5numpy_float64_t *)__Pyx_BufPtrStrided1d(__pyx_bstruct_output.buf, __pyx_t_49, __pyx_bstride_0_output)) = ((((((((__pyx_3 * ((__pyx_v_xm * __pyx_v_ym) * __pyx_v_zm)) + (__pyx_4 * ((__pyx_v_xp * __pyx_v_ym) * __pyx_v_zm))) + (__pyx_5 * ((__pyx_v_xm * __pyx_v_yp) * __pyx_v_zm))) + (__pyx_6 * ((__pyx_v_xm * __pyx_v_ym) * __pyx_v_zp))) + (__pyx_7 * ((__pyx_v_xp * __pyx_v_ym) * __pyx_v_zp))) + (__pyx_8 * ((__pyx_v_xm * __pyx_v_yp) * __pyx_v_zp))) + (__pyx_9 * ((__pyx_v_xp * __pyx_v_yp) * __pyx_v_zm))) + (__pyx_10 * ((__pyx_v_xp * __pyx_v_yp) * __pyx_v_zp)));
  }

  __pyx_r = Py_None; Py_INCREF(Py_None);
  goto __pyx_L0;
  __pyx_L1_error:;
  { PyObject *__pyx_type, *__pyx_value, *__pyx_tb;
    PyErr_Fetch(&__pyx_type, &__pyx_value, &__pyx_tb);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_z_is, &__pyx_bstruct_z_is);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_vals, &__pyx_bstruct_x_vals);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_is, &__pyx_bstruct_y_is);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_bins, &__pyx_bstruct_y_bins);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_z_bins, &__pyx_bstruct_z_bins);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_is, &__pyx_bstruct_x_is);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_output, &__pyx_bstruct_output);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_table, &__pyx_bstruct_table);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_bins, &__pyx_bstruct_x_bins);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_vals, &__pyx_bstruct_y_vals);
    __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_z_vals, &__pyx_bstruct_z_vals);
  PyErr_Restore(__pyx_type, __pyx_value, __pyx_tb);}
  __Pyx_AddTraceback("yt.lagos.Interpolators.TrilinearlyInterpolate");
  __pyx_r = NULL;
  goto __pyx_L2;
  __pyx_L0:;
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_z_is, &__pyx_bstruct_z_is);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_vals, &__pyx_bstruct_x_vals);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_is, &__pyx_bstruct_y_is);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_bins, &__pyx_bstruct_y_bins);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_z_bins, &__pyx_bstruct_z_bins);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_is, &__pyx_bstruct_x_is);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_output, &__pyx_bstruct_output);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_table, &__pyx_bstruct_table);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_x_bins, &__pyx_bstruct_x_bins);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_y_vals, &__pyx_bstruct_y_vals);
  __Pyx_SafeReleaseBuffer((PyObject*)__pyx_v_z_vals, &__pyx_bstruct_z_vals);
  __pyx_L2:;
  return __pyx_r;
}

/* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":36
 *         # 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 (specifically,
 */

static int __pyx_pf_5numpy_7ndarray___getbuffer__(PyObject *__pyx_v_self, Py_buffer *__pyx_v_info, int __pyx_v_flags); /*proto*/
static int __pyx_pf_5numpy_7ndarray___getbuffer__(PyObject *__pyx_v_self, Py_buffer *__pyx_v_info, int __pyx_v_flags) {
  int __pyx_v_t;
  char *__pyx_v_f;
  int __pyx_r;
  int __pyx_1;
  PyObject *__pyx_2 = 0;
  PyObject *__pyx_3 = 0;

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":42
 *             # so the flags are not even checked).
 * 
 *             if sizeof(npy_intp) != sizeof(Py_ssize_t):             # <<<<<<<<<<<<<<
 *                 raise RuntimeError("Py_intptr_t and Py_ssize_t differs in size, numpy.pxd does not support this")
 * 
 */
  __pyx_1 = ((sizeof(npy_intp)) != (sizeof(Py_ssize_t)));
  if (__pyx_1) {

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":43
 * 
 *             if sizeof(npy_intp) != sizeof(Py_ssize_t):
 *                 raise RuntimeError("Py_intptr_t and Py_ssize_t differs in size, numpy.pxd does not support this")             # <<<<<<<<<<<<<<
 * 
 *             info.buf = PyArray_DATA(self)
 */
    __pyx_2 = PyTuple_New(1); if (unlikely(!__pyx_2)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 43; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    Py_INCREF(__pyx_kp_1);
    PyTuple_SET_ITEM(__pyx_2, 0, __pyx_kp_1);
    __pyx_3 = PyObject_Call(__pyx_builtin_RuntimeError, ((PyObject *)__pyx_2), NULL); if (unlikely(!__pyx_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 43; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    Py_DECREF(((PyObject *)__pyx_2)); __pyx_2 = 0;
    __Pyx_Raise(__pyx_3, 0, 0);
    Py_DECREF(__pyx_3); __pyx_3 = 0;
    {__pyx_filename = __pyx_f[1]; __pyx_lineno = 43; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    goto __pyx_L5;
  }
  __pyx_L5:;

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":45
 *                 raise RuntimeError("Py_intptr_t and Py_ssize_t differs in size, numpy.pxd does not support this")
 * 
 *             info.buf = PyArray_DATA(self)             # <<<<<<<<<<<<<<
 *             info.ndim = PyArray_NDIM(self)
 *             info.strides = <Py_ssize_t*>PyArray_STRIDES(self)
 */
  __pyx_v_info->buf = PyArray_DATA(((PyArrayObject *)__pyx_v_self));

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":46
 * 
 *             info.buf = PyArray_DATA(self)
 *             info.ndim = PyArray_NDIM(self)             # <<<<<<<<<<<<<<
 *             info.strides = <Py_ssize_t*>PyArray_STRIDES(self)
 *             info.shape = <Py_ssize_t*>PyArray_DIMS(self)
 */
  __pyx_v_info->ndim = PyArray_NDIM(((PyArrayObject *)__pyx_v_self));

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":47
 *             info.buf = PyArray_DATA(self)
 *             info.ndim = PyArray_NDIM(self)
 *             info.strides = <Py_ssize_t*>PyArray_STRIDES(self)             # <<<<<<<<<<<<<<
 *             info.shape = <Py_ssize_t*>PyArray_DIMS(self)
 *             info.suboffsets = NULL
 */
  __pyx_v_info->strides = ((Py_ssize_t *)PyArray_STRIDES(((PyArrayObject *)__pyx_v_self)));

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":48
 *             info.ndim = PyArray_NDIM(self)
 *             info.strides = <Py_ssize_t*>PyArray_STRIDES(self)
 *             info.shape = <Py_ssize_t*>PyArray_DIMS(self)             # <<<<<<<<<<<<<<
 *             info.suboffsets = NULL
 *             info.itemsize = PyArray_ITEMSIZE(self)
 */
  __pyx_v_info->shape = ((Py_ssize_t *)PyArray_DIMS(((PyArrayObject *)__pyx_v_self)));

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":49
 *             info.strides = <Py_ssize_t*>PyArray_STRIDES(self)
 *             info.shape = <Py_ssize_t*>PyArray_DIMS(self)
 *             info.suboffsets = NULL             # <<<<<<<<<<<<<<
 *             info.itemsize = PyArray_ITEMSIZE(self)
 *             info.readonly = not PyArray_ISWRITEABLE(self)
 */
  __pyx_v_info->suboffsets = NULL;

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":50
 *             info.shape = <Py_ssize_t*>PyArray_DIMS(self)
 *             info.suboffsets = NULL
 *             info.itemsize = PyArray_ITEMSIZE(self)             # <<<<<<<<<<<<<<
 *             info.readonly = not PyArray_ISWRITEABLE(self)
 * 
 */
  __pyx_v_info->itemsize = PyArray_ITEMSIZE(((PyArrayObject *)__pyx_v_self));

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":51
 *             info.suboffsets = NULL
 *             info.itemsize = PyArray_ITEMSIZE(self)
 *             info.readonly = not PyArray_ISWRITEABLE(self)             # <<<<<<<<<<<<<<
 * 
 *             # Formats that are not tested and working in Cython are not
 */
  __pyx_v_info->readonly = (!PyArray_ISWRITEABLE(((PyArrayObject *)__pyx_v_self)));

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":55
 *             # Formats that are not tested and working in Cython are not
 *             # made available from this pxd file yet.
 *             cdef int t = PyArray_TYPE(self)             # <<<<<<<<<<<<<<
 *             cdef char* f = NULL
 *             if   t == NPY_BYTE:       f = "b"
 */
  __pyx_v_t = PyArray_TYPE(((PyArrayObject *)__pyx_v_self));

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":56
 *             # made available from this pxd file yet.
 *             cdef int t = PyArray_TYPE(self)
 *             cdef char* f = NULL             # <<<<<<<<<<<<<<
 *             if   t == NPY_BYTE:       f = "b"
 *             elif t == NPY_UBYTE:      f = "B"
 */
  __pyx_v_f = NULL;

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":57
 *             cdef int t = PyArray_TYPE(self)
 *             cdef char* f = NULL
 *             if   t == NPY_BYTE:       f = "b"             # <<<<<<<<<<<<<<
 *             elif t == NPY_UBYTE:      f = "B"
 *             elif t == NPY_SHORT:      f = "h"
 */
  switch (__pyx_v_t) {
    case NPY_BYTE:
    __pyx_v_f = __pyx_k_2;
    break;
    case NPY_UBYTE:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":58
 *             cdef char* f = NULL
 *             if   t == NPY_BYTE:       f = "b"
 *             elif t == NPY_UBYTE:      f = "B"             # <<<<<<<<<<<<<<
 *             elif t == NPY_SHORT:      f = "h"
 *             elif t == NPY_USHORT:     f = "H"
 */
    __pyx_v_f = __pyx_k_3;
    break;
    case NPY_SHORT:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":59
 *             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"
 */
    __pyx_v_f = __pyx_k_4;
    break;
    case NPY_USHORT:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":60
 *             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"
 */
    __pyx_v_f = __pyx_k_5;
    break;
    case NPY_INT:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":61
 *             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"
 */
    __pyx_v_f = __pyx_k_6;
    break;
    case NPY_UINT:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":62
 *             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"
 */
    __pyx_v_f = __pyx_k_7;
    break;
    case NPY_LONG:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":63
 *             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"
 */
    __pyx_v_f = __pyx_k_8;
    break;
    case NPY_ULONG:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":64
 *             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"
 */
    __pyx_v_f = __pyx_k_9;
    break;
    case NPY_LONGLONG:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":65
 *             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"
 */
    __pyx_v_f = __pyx_k_10;
    break;
    case NPY_ULONGLONG:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":66
 *             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"
 */
    __pyx_v_f = __pyx_k_11;
    break;
    case NPY_FLOAT:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":67
 *             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"
 */
    __pyx_v_f = __pyx_k_12;
    break;
    case NPY_DOUBLE:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":68
 *             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_OBJECT:     f = "O"
 */
    __pyx_v_f = __pyx_k_13;
    break;
    case NPY_LONGDOUBLE:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":69
 *             elif t == NPY_FLOAT:      f = "f"
 *             elif t == NPY_DOUBLE:     f = "d"
 *             elif t == NPY_LONGDOUBLE: f = "g"             # <<<<<<<<<<<<<<
 *             elif t == NPY_OBJECT:     f = "O"
 * 
 */
    __pyx_v_f = __pyx_k_14;
    break;
    case NPY_OBJECT:

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":70
 *             elif t == NPY_DOUBLE:     f = "d"
 *             elif t == NPY_LONGDOUBLE: f = "g"
 *             elif t == NPY_OBJECT:     f = "O"             # <<<<<<<<<<<<<<
 * 
 *             if f == NULL:
 */
    __pyx_v_f = __pyx_k_15;
    break;
  }

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":72
 *             elif t == NPY_OBJECT:     f = "O"
 * 
 *             if f == NULL:             # <<<<<<<<<<<<<<
 *                 raise ValueError("only objects, int and float dtypes supported for ndarray buffer access so far (dtype is %d)" % t)
 *             info.format = f
 */
  __pyx_1 = (__pyx_v_f == NULL);
  if (__pyx_1) {

    /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":73
 * 
 *             if f == NULL:
 *                 raise ValueError("only objects, int and float dtypes supported for ndarray buffer access so far (dtype is %d)" % t)             # <<<<<<<<<<<<<<
 *             info.format = f
 * 
 */
    __pyx_2 = PyInt_FromLong(__pyx_v_t); if (unlikely(!__pyx_2)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __pyx_3 = PyNumber_Remainder(__pyx_kp_16, __pyx_2); if (unlikely(!__pyx_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    Py_DECREF(__pyx_2); __pyx_2 = 0;
    __pyx_2 = PyTuple_New(1); if (unlikely(!__pyx_2)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    PyTuple_SET_ITEM(__pyx_2, 0, __pyx_3);
    __pyx_3 = 0;
    __pyx_3 = PyObject_Call(__pyx_builtin_ValueError, ((PyObject *)__pyx_2), NULL); if (unlikely(!__pyx_3)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    Py_DECREF(((PyObject *)__pyx_2)); __pyx_2 = 0;
    __Pyx_Raise(__pyx_3, 0, 0);
    Py_DECREF(__pyx_3); __pyx_3 = 0;
    {__pyx_filename = __pyx_f[1]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    goto __pyx_L6;
  }
  __pyx_L6:;

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":74
 *             if f == NULL:
 *                 raise ValueError("only objects, int and float dtypes supported for ndarray buffer access so far (dtype is %d)" % t)
 *             info.format = f             # <<<<<<<<<<<<<<
 * 
 * 
 */
  __pyx_v_info->format = __pyx_v_f;

  __pyx_r = 0;
  goto __pyx_L0;
  __pyx_L1_error:;
  Py_XDECREF(__pyx_2);
  Py_XDECREF(__pyx_3);
  __Pyx_AddTraceback("numpy.ndarray.__getbuffer__");
  __pyx_r = -1;
  __pyx_L0:;
  return __pyx_r;
}

static struct PyMethodDef __pyx_methods[] = {
  {"UnilinearlyInterpolate", (PyCFunction)__pyx_pf_2yt_5lagos_13Interpolators_UnilinearlyInterpolate, METH_VARARGS|METH_KEYWORDS, 0},
  {"BilinearlyInterpolate", (PyCFunction)__pyx_pf_2yt_5lagos_13Interpolators_BilinearlyInterpolate, METH_VARARGS|METH_KEYWORDS, 0},
  {"TrilinearlyInterpolate", (PyCFunction)__pyx_pf_2yt_5lagos_13Interpolators_TrilinearlyInterpolate, METH_VARARGS|METH_KEYWORDS, 0},
  {0, 0, 0, 0}
};

static void __pyx_init_filenames(void); /*proto*/

#if PY_MAJOR_VERSION >= 3
static struct PyModuleDef __pyx_moduledef = {
    PyModuleDef_HEAD_INIT,
    "Interpolators",
    __pyx_mdoc, /* m_doc */
    -1, /* m_size */
    __pyx_methods /* m_methods */,
    NULL, /* m_reload */
    NULL, /* m_traverse */
    NULL, /* m_clear */
    NULL /* m_free */
};
#endif

static __Pyx_StringTabEntry __pyx_string_tab[] = {
  {&__pyx_kp_numpy, __pyx_k_numpy, sizeof(__pyx_k_numpy), 1, 1, 1},
  {&__pyx_kp_np, __pyx_k_np, sizeof(__pyx_k_np), 0, 1, 1},
  {&__pyx_kp___getbuffer__, __pyx_k___getbuffer__, sizeof(__pyx_k___getbuffer__), 0, 1, 1},
  {&__pyx_kp_RuntimeError, __pyx_k_RuntimeError, sizeof(__pyx_k_RuntimeError), 1, 1, 1},
  {&__pyx_kp_ValueError, __pyx_k_ValueError, sizeof(__pyx_k_ValueError), 1, 1, 1},
  {&__pyx_kp_1, __pyx_k_1, sizeof(__pyx_k_1), 0, 0, 0},
  {&__pyx_kp_16, __pyx_k_16, sizeof(__pyx_k_16), 0, 0, 0},
  {0, 0, 0, 0, 0, 0}
};
static int __Pyx_InitCachedBuiltins(void) {
  __pyx_builtin_RuntimeError = __Pyx_GetName(__pyx_b, __pyx_kp_RuntimeError); if (!__pyx_builtin_RuntimeError) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 43; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_builtin_ValueError = __Pyx_GetName(__pyx_b, __pyx_kp_ValueError); if (!__pyx_builtin_ValueError) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  return 0;
  __pyx_L1_error:;
  return -1;
}

static int __Pyx_InitGlobals(void) {
  if (__Pyx_InitStrings(__pyx_string_tab) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};
  return 0;
  __pyx_L1_error:;
  return -1;
}

#if PY_MAJOR_VERSION < 3
PyMODINIT_FUNC initInterpolators(void); /*proto*/
PyMODINIT_FUNC initInterpolators(void)
#else
PyMODINIT_FUNC PyInit_Interpolators(void); /*proto*/
PyMODINIT_FUNC PyInit_Interpolators(void)
#endif
{
  PyObject *__pyx_1 = 0;
  __pyx_empty_tuple = PyTuple_New(0); if (unlikely(!__pyx_empty_tuple)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  /*--- Libary function declarations ---*/
  __pyx_init_filenames();
  /*--- Initialize various global constants etc. ---*/
  if (unlikely(__Pyx_InitGlobals() < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  /*--- Module creation code ---*/
  #if PY_MAJOR_VERSION < 3
  __pyx_m = Py_InitModule4("Interpolators", __pyx_methods, __pyx_mdoc, 0, PYTHON_API_VERSION);
  #else
  __pyx_m = PyModule_Create(&__pyx_moduledef);
  #endif
  if (!__pyx_m) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};
  #if PY_MAJOR_VERSION < 3
  Py_INCREF(__pyx_m);
  #endif
  __pyx_b = PyImport_AddModule(__Pyx_BUILTIN_MODULE_NAME);
  if (!__pyx_b) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};
  if (PyObject_SetAttrString(__pyx_m, "__builtins__", __pyx_b) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};
  /*--- Builtin init code ---*/
  if (unlikely(__Pyx_InitCachedBuiltins() < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_skip_dispatch = 0;
  /*--- Global init code ---*/
  /*--- Function export code ---*/
  /*--- Type init code ---*/
  /*--- Type import code ---*/
  __pyx_ptype_5numpy_ndarray = __Pyx_ImportType("numpy", "ndarray", sizeof(PyArrayObject)); if (unlikely(!__pyx_ptype_5numpy_ndarray)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  /*--- Function import code ---*/
  /*--- Execution code ---*/

  /* "/Users/matthewturk/Development/yt/trunk/yt/lagos/Interpolators.pyx":26
 * """
 * 
 * import numpy as np             # <<<<<<<<<<<<<<
 * cimport numpy as np
 * cimport cython
 */
  __pyx_1 = __Pyx_Import(__pyx_kp_numpy, 0); if (unlikely(!__pyx_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (PyObject_SetAttr(__pyx_m, __pyx_kp_np, __pyx_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  Py_DECREF(__pyx_1); __pyx_1 = 0;

  /* "/Library/Frameworks/Python.framework/Versions/2.5/lib/python2.5/site-packages/Cython-0.9.8.1.1-py2.5-macosx-10.3-fat.egg/Cython/Includes/numpy.pxd":36
 *         # 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 (specifically,
 */
  #if PY_MAJOR_VERSION < 3
  return;
  #else
  return __pyx_m;
  #endif
  __pyx_L1_error:;
  Py_XDECREF(__pyx_1);
  __Pyx_AddTraceback("yt.lagos.Interpolators");
  #if PY_MAJOR_VERSION >= 3
  return NULL;
  #endif
}

static const char *__pyx_filenames[] = {
  "Interpolators.pyx",
  "numpy.pxd",
};

/* Runtime support code */

static void __pyx_init_filenames(void) {
  __pyx_f = __pyx_filenames;
}

static INLINE void __Pyx_SafeReleaseBuffer(PyObject* obj, Py_buffer* info) {
  if (info->buf == NULL) return;
  if (info->suboffsets == __Pyx_minusones) info->suboffsets = NULL;
  __Pyx_ReleaseBuffer(obj, info);
}

static INLINE void __Pyx_ZeroBuffer(Py_buffer* buf) {
  buf->buf = NULL;
  buf->strides = __Pyx_zeros;
  buf->shape = __Pyx_zeros;
  buf->suboffsets = __Pyx_minusones;
}

static INLINE const char* __Pyx_ConsumeWhitespace(const char* ts) {
  while (1) {
    switch (*ts) {
      case 10:
      case 13:
      case ' ':
        ++ts;
      default:
        return ts;
    }
  }
}

static INLINE const char* __Pyx_BufferTypestringCheckEndian(const char* ts) {
  int num = 1;
  int little_endian = ((char*)&num)[0];
  int ok = 1;
  switch (*ts) {
    case '@':
    case '=':
      ++ts; break;
    case '<':
      if (little_endian) ++ts;
      else ok = 0;
      break;
    case '>':
    case '!':
      if (!little_endian) ++ts;
      else ok = 0;
      break;
  }
  if (!ok) {
    PyErr_Format(PyExc_ValueError, "Buffer has wrong endianness (rejecting on '%s')", ts);
    return NULL;
  }
  return ts;
}

static void __Pyx_BufferNdimError(Py_buffer* buffer, int expected_ndim) {
  PyErr_Format(PyExc_ValueError,
               "Buffer has wrong number of dimensions (expected %d, got %d)",
               expected_ndim, buffer->ndim);
}


static const char* __Pyx_BufferTypestringCheck_item_nn___pyx_t_5numpy_float64_t(const char* ts) {
  int ok;
  if (*ts == '1') ++ts;
  switch (*ts) {
    case 'f': ok = (sizeof(__pyx_t_5numpy_float64_t) == sizeof(float) && (__pyx_t_5numpy_float64_t)-1 < 0); break;
    case 'd': ok = (sizeof(__pyx_t_5numpy_float64_t) == sizeof(double) && (__pyx_t_5numpy_float64_t)-1 < 0); break;
    case 'g': ok = (sizeof(__pyx_t_5numpy_float64_t) == sizeof(long double) && (__pyx_t_5numpy_float64_t)-1 < 0); break;  default: ok = 0;
  }
  if (!ok) {
      PyErr_Format(PyExc_ValueError, "Buffer datatype mismatch (rejecting on '%s')", ts);
      return NULL;
  } else return ts + 1;
  
}

static int __Pyx_GetBuffer_nn___pyx_t_5numpy_float64_t(PyObject* obj, Py_buffer* buf, int flags, int nd) {
  const char* ts;
  if (obj == Py_None) {
    __Pyx_ZeroBuffer(buf);
    return 0;
  }
  buf->buf = NULL;
  if (__Pyx_GetBuffer(obj, buf, flags) == -1) goto fail;
  if (buf->ndim != nd) {
    __Pyx_BufferNdimError(buf, nd);
    goto fail;
  }
  ts = buf->format;
  ts = __Pyx_ConsumeWhitespace(ts);
  ts = __Pyx_BufferTypestringCheckEndian(ts);
  if (!ts) goto fail;
  ts = __Pyx_ConsumeWhitespace(ts);
  ts = __Pyx_BufferTypestringCheck_item_nn___pyx_t_5numpy_float64_t(ts);
  if (!ts) goto fail;
  ts = __Pyx_ConsumeWhitespace(ts);
  if (*ts != 0) {
    PyErr_Format(PyExc_ValueError,
      "Expected non-struct buffer data type (expected end, got '%s')", ts);
    goto fail;
  }
  if (buf->suboffsets == NULL) buf->suboffsets = __Pyx_minusones;
  return 0;
fail:;
  __Pyx_ZeroBuffer(buf);
  return -1;
}
static const char* __Pyx_BufferTypestringCheck_item_nn___pyx_t_5numpy_int32_t(const char* ts) {
  int ok;
  if (*ts == '1') ++ts;
  switch (*ts) {
    case 'b': ok = (sizeof(__pyx_t_5numpy_int32_t) == sizeof(char) && (__pyx_t_5numpy_int32_t)-1 < 0); break;
    case 'h': ok = (sizeof(__pyx_t_5numpy_int32_t) == sizeof(short) && (__pyx_t_5numpy_int32_t)-1 < 0); break;
    case 'i': ok = (sizeof(__pyx_t_5numpy_int32_t) == sizeof(int) && (__pyx_t_5numpy_int32_t)-1 < 0); break;
    case 'l': ok = (sizeof(__pyx_t_5numpy_int32_t) == sizeof(long) && (__pyx_t_5numpy_int32_t)-1 < 0); break;
    case 'q': ok = (sizeof(__pyx_t_5numpy_int32_t) == sizeof(long long) && (__pyx_t_5numpy_int32_t)-1 < 0); break;  default: ok = 0;
  }
  if (!ok) {
      PyErr_Format(PyExc_ValueError, "Buffer datatype mismatch (rejecting on '%s')", ts);
      return NULL;
  } else return ts + 1;
  
}

static int __Pyx_GetBuffer_nn___pyx_t_5numpy_int32_t(PyObject* obj, Py_buffer* buf, int flags, int nd) {
  const char* ts;
  if (obj == Py_None) {
    __Pyx_ZeroBuffer(buf);
    return 0;
  }
  buf->buf = NULL;
  if (__Pyx_GetBuffer(obj, buf, flags) == -1) goto fail;
  if (buf->ndim != nd) {
    __Pyx_BufferNdimError(buf, nd);
    goto fail;
  }
  ts = buf->format;
  ts = __Pyx_ConsumeWhitespace(ts);
  ts = __Pyx_BufferTypestringCheckEndian(ts);
  if (!ts) goto fail;
  ts = __Pyx_ConsumeWhitespace(ts);
  ts = __Pyx_BufferTypestringCheck_item_nn___pyx_t_5numpy_int32_t(ts);
  if (!ts) goto fail;
  ts = __Pyx_ConsumeWhitespace(ts);
  if (*ts != 0) {
    PyErr_Format(PyExc_ValueError,
      "Expected non-struct buffer data type (expected end, got '%s')", ts);
    goto fail;
  }
  if (buf->suboffsets == NULL) buf->suboffsets = __Pyx_minusones;
  return 0;
fail:;
  __Pyx_ZeroBuffer(buf);
  return -1;
}
static const char* __Pyx_BufferTypestringCheck_item_nn___pyx_t_5numpy_int_t(const char* ts) {
  int ok;
  if (*ts == '1') ++ts;
  switch (*ts) {
    case 'b': ok = (sizeof(__pyx_t_5numpy_int_t) == sizeof(char) && (__pyx_t_5numpy_int_t)-1 < 0); break;
    case 'h': ok = (sizeof(__pyx_t_5numpy_int_t) == sizeof(short) && (__pyx_t_5numpy_int_t)-1 < 0); break;
    case 'i': ok = (sizeof(__pyx_t_5numpy_int_t) == sizeof(int) && (__pyx_t_5numpy_int_t)-1 < 0); break;
    case 'l': ok = (sizeof(__pyx_t_5numpy_int_t) == sizeof(long) && (__pyx_t_5numpy_int_t)-1 < 0); break;
    case 'q': ok = (sizeof(__pyx_t_5numpy_int_t) == sizeof(long long) && (__pyx_t_5numpy_int_t)-1 < 0); break;  default: ok = 0;
  }
  if (!ok) {
      PyErr_Format(PyExc_ValueError, "Buffer datatype mismatch (rejecting on '%s')", ts);
      return NULL;
  } else return ts + 1;
  
}

static int __Pyx_GetBuffer_nn___pyx_t_5numpy_int_t(PyObject* obj, Py_buffer* buf, int flags, int nd) {
  const char* ts;
  if (obj == Py_None) {
    __Pyx_ZeroBuffer(buf);
    return 0;
  }
  buf->buf = NULL;
  if (__Pyx_GetBuffer(obj, buf, flags) == -1) goto fail;
  if (buf->ndim != nd) {
    __Pyx_BufferNdimError(buf, nd);
    goto fail;
  }
  ts = buf->format;
  ts = __Pyx_ConsumeWhitespace(ts);
  ts = __Pyx_BufferTypestringCheckEndian(ts);
  if (!ts) goto fail;
  ts = __Pyx_ConsumeWhitespace(ts);
  ts = __Pyx_BufferTypestringCheck_item_nn___pyx_t_5numpy_int_t(ts);
  if (!ts) goto fail;
  ts = __Pyx_ConsumeWhitespace(ts);
  if (*ts != 0) {
    PyErr_Format(PyExc_ValueError,
      "Expected non-struct buffer data type (expected end, got '%s')", ts);
    goto fail;
  }
  if (buf->suboffsets == NULL) buf->suboffsets = __Pyx_minusones;
  return 0;
fail:;
  __Pyx_ZeroBuffer(buf);
  return -1;
}
static int __Pyx_ArgTypeTest(PyObject *obj, PyTypeObject *type, int none_allowed, char *name, int exact) {
    if (!type) {
        PyErr_Format(PyExc_SystemError, "Missing type object");
        return 0;
    }
    if (none_allowed && obj == Py_None) return 1;
    else if (exact) {
        if (Py_TYPE(obj) == type) return 1;
    }
    else {
        if (PyObject_TypeCheck(obj, type)) return 1;
    }
    PyErr_Format(PyExc_TypeError,
        "Argument '%s' has incorrect type (expected %s, got %s)",
        name, type->tp_name, Py_TYPE(obj)->tp_name);
    return 0;
}

static INLINE void __Pyx_RaiseArgtupleTooLong(
    Py_ssize_t num_expected,
    Py_ssize_t num_found)
{
    const char* error_message =
    #if PY_VERSION_HEX < 0x02050000
        "function takes at most %d positional arguments (%d given)";
    #else
        "function takes at most %zd positional arguments (%zd given)";
    #endif
    PyErr_Format(PyExc_TypeError, error_message, num_expected, num_found);
}

#if (PY_MAJOR_VERSION < 3) && !(Py_TPFLAGS_DEFAULT & Py_TPFLAGS_HAVE_NEWBUFFER)
static int __Pyx_GetBuffer(PyObject *obj, Py_buffer *view, int flags) {
  if (PyObject_TypeCheck(obj, __pyx_ptype_5numpy_ndarray)) return __pyx_pf_5numpy_7ndarray___getbuffer__(obj, view, flags);
  else {
  PyErr_Format(PyExc_TypeError, "'%100s' does not have the buffer interface", Py_TYPE(obj)->tp_name);
  return -1;
    }
}

static void __Pyx_ReleaseBuffer(PyObject *obj, Py_buffer *view) {

}

#endif

static PyObject *__Pyx_Import(PyObject *name, PyObject *from_list) {
    PyObject *__import__ = 0;
    PyObject *empty_list = 0;
    PyObject *module = 0;
    PyObject *global_dict = 0;
    PyObject *empty_dict = 0;
    PyObject *list;
    __import__ = PyObject_GetAttrString(__pyx_b, "__import__");
    if (!__import__)
        goto bad;
    if (from_list)
        list = from_list;
    else {
        empty_list = PyList_New(0);
        if (!empty_list)
            goto bad;
        list = empty_list;
    }
    global_dict = PyModule_GetDict(__pyx_m);
    if (!global_dict)
        goto bad;
    empty_dict = PyDict_New();
    if (!empty_dict)
        goto bad;
    module = PyObject_CallFunction(__import__, "OOOO",
        name, global_dict, empty_dict, list);
bad:
    Py_XDECREF(empty_list);
    Py_XDECREF(__import__);
    Py_XDECREF(empty_dict);
    return module;
}

static PyObject *__Pyx_GetName(PyObject *dict, PyObject *name) {
    PyObject *result;
    result = PyObject_GetAttr(dict, name);
    if (!result)
        PyErr_SetObject(PyExc_NameError, name);
    return result;
}

static void __Pyx_Raise(PyObject *type, PyObject *value, PyObject *tb) {
    Py_XINCREF(type);
    Py_XINCREF(value);
    Py_XINCREF(tb);
    /* First, check the traceback argument, replacing None with NULL. */
    if (tb == Py_None) {
        Py_DECREF(tb);
        tb = 0;
    }
    else if (tb != NULL && !PyTraceBack_Check(tb)) {
        PyErr_SetString(PyExc_TypeError,
            "raise: arg 3 must be a traceback or None");
        goto raise_error;
    }
    /* Next, replace a missing value with None */
    if (value == NULL) {
        value = Py_None;
        Py_INCREF(value);
    }
    #if PY_VERSION_HEX < 0x02050000
    if (!PyClass_Check(type))
    #else
    if (!PyType_Check(type))
    #endif
    {
        /* Raising an instance.  The value should be a dummy. */
        if (value != Py_None) {
            PyErr_SetString(PyExc_TypeError,
                "instance exception may not have a separate value");
            goto raise_error;
        }
        /* Normalize to raise <class>, <instance> */
        Py_DECREF(value);
        value = type;
        #if PY_VERSION_HEX < 0x02050000
            if (PyInstance_Check(type)) {
                type = (PyObject*) ((PyInstanceObject*)type)->in_class;
                Py_INCREF(type);
            }
            else {
                type = 0;
                PyErr_SetString(PyExc_TypeError,
                    "raise: exception must be an old-style class or instance");
                goto raise_error;
            }
        #else
            type = (PyObject*) Py_TYPE(type);
            Py_INCREF(type);
            if (!PyType_IsSubtype((PyTypeObject *)type, (PyTypeObject *)PyExc_BaseException)) {
                PyErr_SetString(PyExc_TypeError,
                    "raise: exception class must be a subclass of BaseException");
                goto raise_error;
            }
        #endif
    }
    PyErr_Restore(type, value, tb);
    return;
raise_error:
    Py_XDECREF(value);
    Py_XDECREF(type);
    Py_XDECREF(tb);
    return;
}

#ifndef __PYX_HAVE_RT_ImportType
#define __PYX_HAVE_RT_ImportType
static PyTypeObject *__Pyx_ImportType(char *module_name, char *class_name,
    long size)
{
    PyObject *py_module = 0;
    PyObject *result = 0;
    PyObject *py_name = 0;

    #if PY_MAJOR_VERSION < 3
    py_name = PyString_FromString(module_name);
    #else
    py_name = PyUnicode_FromString(module_name);
    #endif
    if (!py_name)
        goto bad;

    py_module = __Pyx_ImportModule(module_name);
    if (!py_module)
        goto bad;
    result = PyObject_GetAttrString(py_module, class_name);
    if (!result)
        goto bad;
    if (!PyType_Check(result)) {
        PyErr_Format(PyExc_TypeError, 
            "%s.%s is not a type object",
            module_name, class_name);
        goto bad;
    }
    if (((PyTypeObject *)result)->tp_basicsize != size) {
        PyErr_Format(PyExc_ValueError, 
            "%s.%s does not appear to be the correct type object",
            module_name, class_name);
        goto bad;
    }
    return (PyTypeObject *)result;
bad:
    Py_XDECREF(py_name);
    Py_XDECREF(result);
    return 0;
}
#endif

#ifndef __PYX_HAVE_RT_ImportModule
#define __PYX_HAVE_RT_ImportModule
static PyObject *__Pyx_ImportModule(char *name) {
    PyObject *py_name = 0;
    PyObject *py_module = 0;

    #if PY_MAJOR_VERSION < 3
    py_name = PyString_FromString(name);
    #else
    py_name = PyUnicode_FromString(name);
    #endif
    if (!py_name)
        goto bad;
    py_module = PyImport_Import(py_name);
    Py_DECREF(py_name);
    return py_module;
bad:
    Py_XDECREF(py_name);
    return 0;
}
#endif

#include "compile.h"
#include "frameobject.h"
#include "traceback.h"

static void __Pyx_AddTraceback(const char *funcname) {
    PyObject *py_srcfile = 0;
    PyObject *py_funcname = 0;
    PyObject *py_globals = 0;
    PyObject *empty_string = 0;
    PyCodeObject *py_code = 0;
    PyFrameObject *py_frame = 0;

    #if PY_MAJOR_VERSION < 3
    py_srcfile = PyString_FromString(__pyx_filename);
    #else
    py_srcfile = PyUnicode_FromString(__pyx_filename);
    #endif
    if (!py_srcfile) goto bad;
    if (__pyx_clineno) {
        #if PY_MAJOR_VERSION < 3
        py_funcname = PyString_FromFormat( "%s (%s:%d)", funcname, __pyx_cfilenm, __pyx_clineno);
        #else
        py_funcname = PyUnicode_FromFormat( "%s (%s:%d)", funcname, __pyx_cfilenm, __pyx_clineno);
        #endif
    }
    else {
        #if PY_MAJOR_VERSION < 3
        py_funcname = PyString_FromString(funcname);
        #else
        py_funcname = PyUnicode_FromString(funcname);
        #endif
    }
    if (!py_funcname) goto bad;
    py_globals = PyModule_GetDict(__pyx_m);
    if (!py_globals) goto bad;
    #if PY_MAJOR_VERSION < 3
    empty_string = PyString_FromStringAndSize("", 0);
    #else
    empty_string = PyBytes_FromStringAndSize("", 0);
    #endif
    if (!empty_string) goto bad;
    py_code = PyCode_New(
        0,            /*int argcount,*/
        #if PY_MAJOR_VERSION >= 3
        0,            /*int kwonlyargcount,*/
        #endif
        0,            /*int nlocals,*/
        0,            /*int stacksize,*/
        0,            /*int flags,*/
        empty_string, /*PyObject *code,*/
        __pyx_empty_tuple,  /*PyObject *consts,*/
        __pyx_empty_tuple,  /*PyObject *names,*/
        __pyx_empty_tuple,  /*PyObject *varnames,*/
        __pyx_empty_tuple,  /*PyObject *freevars,*/
        __pyx_empty_tuple,  /*PyObject *cellvars,*/
        py_srcfile,   /*PyObject *filename,*/
        py_funcname,  /*PyObject *name,*/
        __pyx_lineno,   /*int firstlineno,*/
        empty_string  /*PyObject *lnotab*/
    );
    if (!py_code) goto bad;
    py_frame = PyFrame_New(
        PyThreadState_Get(), /*PyThreadState *tstate,*/
        py_code,             /*PyCodeObject *code,*/
        py_globals,          /*PyObject *globals,*/
        0                    /*PyObject *locals*/
    );
    if (!py_frame) goto bad;
    py_frame->f_lineno = __pyx_lineno;
    PyTraceBack_Here(py_frame);
bad:
    Py_XDECREF(py_srcfile);
    Py_XDECREF(py_funcname);
    Py_XDECREF(empty_string);
    Py_XDECREF(py_code);
    Py_XDECREF(py_frame);
}

static int __Pyx_InitStrings(__Pyx_StringTabEntry *t) {
    while (t->p) {
        #if PY_MAJOR_VERSION < 3
        if (t->is_unicode && (!t->is_identifier)) {
            *t->p = PyUnicode_DecodeUTF8(t->s, t->n - 1, NULL);
        } else if (t->intern) {
            *t->p = PyString_InternFromString(t->s);
        } else {
            *t->p = PyString_FromStringAndSize(t->s, t->n - 1);
        }
        #else  /* Python 3+ has unicode identifiers */
        if (t->is_identifier || (t->is_unicode && t->intern)) {
            *t->p = PyUnicode_InternFromString(t->s);
        } else if (t->is_unicode) {
            *t->p = PyUnicode_FromStringAndSize(t->s, t->n - 1);
        } else {
            *t->p = PyBytes_FromStringAndSize(t->s, t->n - 1);
        }
        #endif
        if (!*t->p)
            return -1;
        ++t;
    }
    return 0;
}

/* Type Conversion Functions */

static INLINE Py_ssize_t __pyx_PyIndex_AsSsize_t(PyObject* b) {
  Py_ssize_t ival;
  PyObject* x = PyNumber_Index(b);
  if (!x) return -1;
  ival = PyInt_AsSsize_t(x);
  Py_DECREF(x);
  return ival;
}

static INLINE int __Pyx_PyObject_IsTrue(PyObject* x) {
   if (x == Py_True) return 1;
   else if (x == Py_False) return 0;
   else return PyObject_IsTrue(x);
}

static INLINE PY_LONG_LONG __pyx_PyInt_AsLongLong(PyObject* x) {
    if (PyInt_CheckExact(x)) {
        return PyInt_AS_LONG(x);
    }
    else if (PyLong_CheckExact(x)) {
        return PyLong_AsLongLong(x);
    }
    else {
        PY_LONG_LONG val;
        PyObject* tmp = PyNumber_Int(x); if (!tmp) return (PY_LONG_LONG)-1;
        val = __pyx_PyInt_AsLongLong(tmp);
        Py_DECREF(tmp);
        return val;
    }
}

static INLINE unsigned PY_LONG_LONG __pyx_PyInt_AsUnsignedLongLong(PyObject* x) {
    if (PyInt_CheckExact(x)) {
        long val = PyInt_AS_LONG(x);
        if (unlikely(val < 0)) {
            PyErr_SetString(PyExc_TypeError, "Negative assignment to unsigned type.");
            return (unsigned PY_LONG_LONG)-1;
        }
        return val;
    }
    else if (PyLong_CheckExact(x)) {
        return PyLong_AsUnsignedLongLong(x);
    }
    else {
        PY_LONG_LONG val;
        PyObject* tmp = PyNumber_Int(x); if (!tmp) return (PY_LONG_LONG)-1;
        val = __pyx_PyInt_AsUnsignedLongLong(tmp);
        Py_DECREF(tmp);
        return val;
    }
}


static INLINE unsigned char __pyx_PyInt_unsigned_char(PyObject* x) {
    if (sizeof(unsigned char) < sizeof(long)) {
        long long_val = __pyx_PyInt_AsLong(x);
        unsigned char val = (unsigned char)long_val;
        if (unlikely((val != long_val)  || (long_val < 0))) {
            PyErr_SetString(PyExc_OverflowError, "value too large to convert to unsigned char");
            return (unsigned char)-1;
        }
        return val;
    }
    else {
        return __pyx_PyInt_AsLong(x);
    }
}

static INLINE unsigned short __pyx_PyInt_unsigned_short(PyObject* x) {
    if (sizeof(unsigned short) < sizeof(long)) {
        long long_val = __pyx_PyInt_AsLong(x);
        unsigned short val = (unsigned short)long_val;
        if (unlikely((val != long_val)  || (long_val < 0))) {
            PyErr_SetString(PyExc_OverflowError, "value too large to convert to unsigned short");
            return (unsigned short)-1;
        }
        return val;
    }
    else {
        return __pyx_PyInt_AsLong(x);
    }
}

static INLINE char __pyx_PyInt_char(PyObject* x) {
    if (sizeof(char) < sizeof(long)) {
        long long_val = __pyx_PyInt_AsLong(x);
        char val = (char)long_val;
        if (unlikely((val != long_val) )) {
            PyErr_SetString(PyExc_OverflowError, "value too large to convert to char");
            return (char)-1;
        }
        return val;
    }
    else {
        return __pyx_PyInt_AsLong(x);
    }
}

static INLINE short __pyx_PyInt_short(PyObject* x) {
    if (sizeof(short) < sizeof(long)) {
        long long_val = __pyx_PyInt_AsLong(x);
        short val = (short)long_val;
        if (unlikely((val != long_val) )) {
            PyErr_SetString(PyExc_OverflowError, "value too large to convert to short");
            return (short)-1;
        }
        return val;
    }
    else {
        return __pyx_PyInt_AsLong(x);
    }
}

static INLINE int __pyx_PyInt_int(PyObject* x) {
    if (sizeof(int) < sizeof(long)) {
        long long_val = __pyx_PyInt_AsLong(x);
        int val = (int)long_val;
        if (unlikely((val != long_val) )) {
            PyErr_SetString(PyExc_OverflowError, "value too large to convert to int");
            return (int)-1;
        }
        return val;
    }
    else {
        return __pyx_PyInt_AsLong(x);
    }
}

static INLINE long __pyx_PyInt_long(PyObject* x) {
    if (sizeof(long) < sizeof(long)) {
        long long_val = __pyx_PyInt_AsLong(x);
        long val = (long)long_val;
        if (unlikely((val != long_val) )) {
            PyErr_SetString(PyExc_OverflowError, "value too large to convert to long");
            return (long)-1;
        }
        return val;
    }
    else {
        return __pyx_PyInt_AsLong(x);
    }
}

static INLINE signed char __pyx_PyInt_signed_char(PyObject* x) {
    if (sizeof(signed char) < sizeof(long)) {
        long long_val = __pyx_PyInt_AsLong(x);
        signed char val = (signed char)long_val;
        if (unlikely((val != long_val) )) {
            PyErr_SetString(PyExc_OverflowError, "value too large to convert to signed char");
            return (signed char)-1;
        }
        return val;
    }
    else {
        return __pyx_PyInt_AsLong(x);
    }
}

static INLINE signed short __pyx_PyInt_signed_short(PyObject* x) {
    if (sizeof(signed short) < sizeof(long)) {
        long long_val = __pyx_PyInt_AsLong(x);
        signed short val = (signed short)long_val;
        if (unlikely((val != long_val) )) {
            PyErr_SetString(PyExc_OverflowError, "value too large to convert to signed short");
            return (signed short)-1;
        }
        return val;
    }
    else {
        return __pyx_PyInt_AsLong(x);
    }
}

static INLINE signed int __pyx_PyInt_signed_int(PyObject* x) {
    if (sizeof(signed int) < sizeof(long)) {
        long long_val = __pyx_PyInt_AsLong(x);
        signed int val = (signed int)long_val;
        if (unlikely((val != long_val) )) {
            PyErr_SetString(PyExc_OverflowError, "value too large to convert to signed int");
            return (signed int)-1;
        }
        return val;
    }
    else {
        return __pyx_PyInt_AsLong(x);
    }
}

static INLINE signed long __pyx_PyInt_signed_long(PyObject* x) {
    if (sizeof(signed long) < sizeof(long)) {
        long long_val = __pyx_PyInt_AsLong(x);
        signed long val = (signed long)long_val;
        if (unlikely((val != long_val) )) {
            PyErr_SetString(PyExc_OverflowError, "value too large to convert to signed long");
            return (signed long)-1;
        }
        return val;
    }
    else {
        return __pyx_PyInt_AsLong(x);
    }
}

static INLINE long double __pyx_PyInt_long_double(PyObject* x) {
    if (sizeof(long double) < sizeof(long)) {
        long long_val = __pyx_PyInt_AsLong(x);
        long double val = (long double)long_val;
        if (unlikely((val != long_val) )) {
            PyErr_SetString(PyExc_OverflowError, "value too large to convert to long double");
            return (long double)-1;
        }
        return val;
    }
    else {
        return __pyx_PyInt_AsLong(x);
    }
}

