Skip to content
Snippets Groups Projects
Normal view Permalink
ancillary.c 990 KiB
Newer Older
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027 2028 2029 2030 2031 2032 2033 2034 2035 2036 2037 2038 2039 2040 2041 2042 2043 2044 2045 2046 2047 2048 2049 2050 2051 2052 2053 2054 2055 2056 2057 2058 2059 2060 2061 2062 2063 2064 2065 2066 2067 2068 2069 2070 2071 2072 2073 2074 2075 2076 2077 2078 2079 2080 2081 2082 2083 2084 2085 2086 2087 2088 

/* PyObject_GenericGetAttr.proto */
#if CYTHON_USE_TYPE_SLOTS && CYTHON_USE_PYTYPE_LOOKUP && PY_VERSION_HEX < 0x03070000
static PyObject* __Pyx_PyObject_GenericGetAttr(PyObject* obj, PyObject* attr_name);
#else
#define __Pyx_PyObject_GenericGetAttr PyObject_GenericGetAttr
#endif

/* SetVTable.proto */
static int __Pyx_SetVtable(PyObject *dict, void *vtable);

/* PyObjectGetAttrStrNoError.proto */
static CYTHON_INLINE PyObject* __Pyx_PyObject_GetAttrStrNoError(PyObject* obj, PyObject* attr_name);

/* SetupReduce.proto */
static int __Pyx_setup_reduce(PyObject* type_obj);

/* TypeImport.proto */
#ifndef __PYX_HAVE_RT_ImportType_proto
#define __PYX_HAVE_RT_ImportType_proto
enum __Pyx_ImportType_CheckSize {
   __Pyx_ImportType_CheckSize_Error = 0,
   __Pyx_ImportType_CheckSize_Warn = 1,
   __Pyx_ImportType_CheckSize_Ignore = 2
};
static PyTypeObject *__Pyx_ImportType(PyObject* module, const char *module_name, const char *class_name, size_t size, enum __Pyx_ImportType_CheckSize check_size);
#endif

/* CLineInTraceback.proto */
#ifdef CYTHON_CLINE_IN_TRACEBACK
#define __Pyx_CLineForTraceback(tstate, c_line)  (((CYTHON_CLINE_IN_TRACEBACK)) ? c_line : 0)
#else
static int __Pyx_CLineForTraceback(PyThreadState *tstate, int c_line);
#endif

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

/* AddTraceback.proto */
static void __Pyx_AddTraceback(const char *funcname, int c_line,
                               int py_line, const char *filename);

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


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

/* MemviewSliceIsContig.proto */
static int __pyx_memviewslice_is_contig(const __Pyx_memviewslice mvs, char order, int ndim);

/* OverlappingSlices.proto */
static int __pyx_slices_overlap(__Pyx_memviewslice *slice1,
                                __Pyx_memviewslice *slice2,
                                int ndim, size_t itemsize);

/* Capsule.proto */
static CYTHON_INLINE PyObject *__pyx_capsule_create(void *p, const char *sig);
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2089 2090 2091 2092 2093 
/* GCCDiagnostics.proto */
#if defined(__GNUC__) && (__GNUC__ > 4 || (__GNUC__ == 4 && __GNUC_MINOR__ >= 6))
#define __Pyx_HAS_GCC_DIAGNOSTIC
#endif
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2094 2095 2096 2097 2098 2099 2100 2101 2102 2103 2104 2105 2106 2107 2108 2109 2110 2111 2112 2113 
/* TypeInfoCompare.proto */
static int __pyx_typeinfo_cmp(__Pyx_TypeInfo *a, __Pyx_TypeInfo *b);

/* MemviewSliceValidateAndInit.proto */
static int __Pyx_ValidateAndInit_memviewslice(
                int *axes_specs,
                int c_or_f_flag,
                int buf_flags,
                int ndim,
                __Pyx_TypeInfo *dtype,
                __Pyx_BufFmt_StackElem stack[],
                __Pyx_memviewslice *memviewslice,
                PyObject *original_obj);

/* ObjectToMemviewSlice.proto */
static CYTHON_INLINE __Pyx_memviewslice __Pyx_PyObject_to_MemoryviewSlice_dc_float(PyObject *, int writable_flag);

/* ObjectToMemviewSlice.proto */
static CYTHON_INLINE __Pyx_memviewslice __Pyx_PyObject_to_MemoryviewSlice_dc_unsigned_char(PyObject *, int writable_flag);
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2114 2115 2116 
/* ObjectToMemviewSlice.proto */
static CYTHON_INLINE __Pyx_memviewslice __Pyx_PyObject_to_MemoryviewSlice_d_dc_int(PyObject *, int writable_flag);
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2117 2118 2119 2120 2121 2122 2123 2124 2125 2126 2127 2128 2129 2130 2131 2132 2133 2134 2135 2136 2137 2138 2139 2140 2141 2142 2143 2144 2145 2146 2147 2148 2149 2150 2151 2152 2153 2154 2155 2156 2157 2158 2159 2160 2161 2162 2163 2164 2165 2166 2167 2168 2169 2170 2171 2172 2173 2174 2175 2176 2177 2178 2179 2180 2181 2182 2183 2184 2185 2186 2187 2188 2189 2190 2191 2192 2193 2194 2195 2196 2197 2198 2199 2200 2201 2202 2203 2204 2205 2206 2207 2208 2209 2210 2211 2212 2213 2214 2215 2216 2217 2218 2219 2220 2221 2222 2223 2224 2225 
/* RealImag.proto */
#if CYTHON_CCOMPLEX
  #ifdef __cplusplus
    #define __Pyx_CREAL(z) ((z).real())
    #define __Pyx_CIMAG(z) ((z).imag())
  #else
    #define __Pyx_CREAL(z) (__real__(z))
    #define __Pyx_CIMAG(z) (__imag__(z))
  #endif
#else
    #define __Pyx_CREAL(z) ((z).real)
    #define __Pyx_CIMAG(z) ((z).imag)
#endif
#if defined(__cplusplus) && CYTHON_CCOMPLEX\
        && (defined(_WIN32) || defined(__clang__) || (defined(__GNUC__) && (__GNUC__ >= 5 || __GNUC__ == 4 && __GNUC_MINOR__ >= 4 )) || __cplusplus >= 201103)
    #define __Pyx_SET_CREAL(z,x) ((z).real(x))
    #define __Pyx_SET_CIMAG(z,y) ((z).imag(y))
#else
    #define __Pyx_SET_CREAL(z,x) __Pyx_CREAL(z) = (x)
    #define __Pyx_SET_CIMAG(z,y) __Pyx_CIMAG(z) = (y)
#endif

/* Arithmetic.proto */
#if CYTHON_CCOMPLEX
    #define __Pyx_c_eq_float(a, b)   ((a)==(b))
    #define __Pyx_c_sum_float(a, b)  ((a)+(b))
    #define __Pyx_c_diff_float(a, b) ((a)-(b))
    #define __Pyx_c_prod_float(a, b) ((a)*(b))
    #define __Pyx_c_quot_float(a, b) ((a)/(b))
    #define __Pyx_c_neg_float(a)     (-(a))
  #ifdef __cplusplus
    #define __Pyx_c_is_zero_float(z) ((z)==(float)0)
    #define __Pyx_c_conj_float(z)    (::std::conj(z))
    #if 1
        #define __Pyx_c_abs_float(z)     (::std::abs(z))
        #define __Pyx_c_pow_float(a, b)  (::std::pow(a, b))
    #endif
  #else
    #define __Pyx_c_is_zero_float(z) ((z)==0)
    #define __Pyx_c_conj_float(z)    (conjf(z))
    #if 1
        #define __Pyx_c_abs_float(z)     (cabsf(z))
        #define __Pyx_c_pow_float(a, b)  (cpowf(a, b))
    #endif
 #endif
#else
    static CYTHON_INLINE int __Pyx_c_eq_float(__pyx_t_float_complex, __pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_sum_float(__pyx_t_float_complex, __pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_diff_float(__pyx_t_float_complex, __pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_prod_float(__pyx_t_float_complex, __pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_quot_float(__pyx_t_float_complex, __pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_neg_float(__pyx_t_float_complex);
    static CYTHON_INLINE int __Pyx_c_is_zero_float(__pyx_t_float_complex);
    static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_conj_float(__pyx_t_float_complex);
    #if 1
        static CYTHON_INLINE float __Pyx_c_abs_float(__pyx_t_float_complex);
        static CYTHON_INLINE __pyx_t_float_complex __Pyx_c_pow_float(__pyx_t_float_complex, __pyx_t_float_complex);
    #endif
#endif

/* Arithmetic.proto */
#if CYTHON_CCOMPLEX
    #define __Pyx_c_eq_double(a, b)   ((a)==(b))
    #define __Pyx_c_sum_double(a, b)  ((a)+(b))
    #define __Pyx_c_diff_double(a, b) ((a)-(b))
    #define __Pyx_c_prod_double(a, b) ((a)*(b))
    #define __Pyx_c_quot_double(a, b) ((a)/(b))
    #define __Pyx_c_neg_double(a)     (-(a))
  #ifdef __cplusplus
    #define __Pyx_c_is_zero_double(z) ((z)==(double)0)
    #define __Pyx_c_conj_double(z)    (::std::conj(z))
    #if 1
        #define __Pyx_c_abs_double(z)     (::std::abs(z))
        #define __Pyx_c_pow_double(a, b)  (::std::pow(a, b))
    #endif
  #else
    #define __Pyx_c_is_zero_double(z) ((z)==0)
    #define __Pyx_c_conj_double(z)    (conj(z))
    #if 1
        #define __Pyx_c_abs_double(z)     (cabs(z))
        #define __Pyx_c_pow_double(a, b)  (cpow(a, b))
    #endif
 #endif
#else
    static CYTHON_INLINE int __Pyx_c_eq_double(__pyx_t_double_complex, __pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_sum_double(__pyx_t_double_complex, __pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_diff_double(__pyx_t_double_complex, __pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_prod_double(__pyx_t_double_complex, __pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_quot_double(__pyx_t_double_complex, __pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_neg_double(__pyx_t_double_complex);
    static CYTHON_INLINE int __Pyx_c_is_zero_double(__pyx_t_double_complex);
    static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_conj_double(__pyx_t_double_complex);
    #if 1
        static CYTHON_INLINE double __Pyx_c_abs_double(__pyx_t_double_complex);
        static CYTHON_INLINE __pyx_t_double_complex __Pyx_c_pow_double(__pyx_t_double_complex, __pyx_t_double_complex);
    #endif
#endif

/* MemviewSliceCopyTemplate.proto */
static __Pyx_memviewslice
__pyx_memoryview_copy_new_contig(const __Pyx_memviewslice *from_mvs,
                                 const char *mode, int ndim,
                                 size_t sizeof_dtype, int contig_flag,
                                 int dtype_is_object);

/* CIntFromPy.proto */
static CYTHON_INLINE int __Pyx_PyInt_As_int(PyObject *);

/* CIntFromPy.proto */
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2226 2227 2228 2229 2230 2231 2232 
static CYTHON_INLINE unsigned char __Pyx_PyInt_As_unsigned_char(PyObject *);

/* CIntToPy.proto */
static CYTHON_INLINE PyObject* __Pyx_PyInt_From_Py_intptr_t(Py_intptr_t value);

/* CIntFromPy.proto */
static CYTHON_INLINE Py_intptr_t __Pyx_PyInt_As_Py_intptr_t(PyObject *);
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2233 2234 2235 2236 

/* CIntToPy.proto */
static CYTHON_INLINE PyObject* __Pyx_PyInt_From_int(int value);
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2237 2238 2239 
/* CIntFromPy.proto */
static CYTHON_INLINE long __Pyx_PyInt_As_long(PyObject *);
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2240 2241 2242 2243 2244 2245 2246 2247 2248 2249 2250 2251 2252 2253 2254 2255 2256 2257 2258 2259 2260 2261 2262 2263 2264 2265 2266 2267 2268 2269 2270 2271 2272 2273 2274 2275 2276 2277 2278 2279 2280 2281 2282 2283 2284 2285 2286 2287 2288 2289 2290 2291 2292 2293 2294 2295 2296 2297 2298 2299 2300 2301 2302 2303 2304 2305 2306 2307 2308 2309 2310 2311 2312 2313 2314 2315 2316 2317 
/* CIntToPy.proto */
static CYTHON_INLINE PyObject* __Pyx_PyInt_From_long(long value);

/* CIntFromPy.proto */
static CYTHON_INLINE char __Pyx_PyInt_As_char(PyObject *);

/* CheckBinaryVersion.proto */
static int __Pyx_check_binary_version(void);

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

static PyObject *__pyx_array_get_memview(struct __pyx_array_obj *__pyx_v_self); /* proto*/
static char *__pyx_memoryview_get_item_pointer(struct __pyx_memoryview_obj *__pyx_v_self, PyObject *__pyx_v_index); /* proto*/
static PyObject *__pyx_memoryview_is_slice(struct __pyx_memoryview_obj *__pyx_v_self, PyObject *__pyx_v_obj); /* proto*/
static PyObject *__pyx_memoryview_setitem_slice_assignment(struct __pyx_memoryview_obj *__pyx_v_self, PyObject *__pyx_v_dst, PyObject *__pyx_v_src); /* proto*/
static PyObject *__pyx_memoryview_setitem_slice_assign_scalar(struct __pyx_memoryview_obj *__pyx_v_self, struct __pyx_memoryview_obj *__pyx_v_dst, PyObject *__pyx_v_value); /* proto*/
static PyObject *__pyx_memoryview_setitem_indexed(struct __pyx_memoryview_obj *__pyx_v_self, PyObject *__pyx_v_index, PyObject *__pyx_v_value); /* proto*/
static PyObject *__pyx_memoryview_convert_item_to_object(struct __pyx_memoryview_obj *__pyx_v_self, char *__pyx_v_itemp); /* proto*/
static PyObject *__pyx_memoryview_assign_item_from_object(struct __pyx_memoryview_obj *__pyx_v_self, char *__pyx_v_itemp, PyObject *__pyx_v_value); /* proto*/
static PyObject *__pyx_memoryviewslice_convert_item_to_object(struct __pyx_memoryviewslice_obj *__pyx_v_self, char *__pyx_v_itemp); /* proto*/
static PyObject *__pyx_memoryviewslice_assign_item_from_object(struct __pyx_memoryviewslice_obj *__pyx_v_self, char *__pyx_v_itemp, PyObject *__pyx_v_value); /* proto*/

/* Module declarations from 'cython.view' */

/* Module declarations from 'cython' */

/* Module declarations from 'cpython.buffer' */

/* Module declarations from 'libc.string' */

/* Module declarations from 'libc.stdio' */

/* Module declarations from '__builtin__' */

/* Module declarations from 'cpython.type' */
static PyTypeObject *__pyx_ptype_7cpython_4type_type = 0;

/* Module declarations from 'cpython' */

/* Module declarations from 'cpython.object' */

/* Module declarations from 'cpython.ref' */

/* Module declarations from 'cpython.mem' */

/* Module declarations from 'numpy' */

/* Module declarations from 'numpy' */
static PyTypeObject *__pyx_ptype_5numpy_dtype = 0;
static PyTypeObject *__pyx_ptype_5numpy_flatiter = 0;
static PyTypeObject *__pyx_ptype_5numpy_broadcast = 0;
static PyTypeObject *__pyx_ptype_5numpy_ndarray = 0;
static PyTypeObject *__pyx_ptype_5numpy_generic = 0;
static PyTypeObject *__pyx_ptype_5numpy_number = 0;
static PyTypeObject *__pyx_ptype_5numpy_integer = 0;
static PyTypeObject *__pyx_ptype_5numpy_signedinteger = 0;
static PyTypeObject *__pyx_ptype_5numpy_unsignedinteger = 0;
static PyTypeObject *__pyx_ptype_5numpy_inexact = 0;
static PyTypeObject *__pyx_ptype_5numpy_floating = 0;
static PyTypeObject *__pyx_ptype_5numpy_complexfloating = 0;
static PyTypeObject *__pyx_ptype_5numpy_flexible = 0;
static PyTypeObject *__pyx_ptype_5numpy_character = 0;
static PyTypeObject *__pyx_ptype_5numpy_ufunc = 0;
static CYTHON_INLINE int __pyx_f_5numpy_import_array(void); /*proto*/

/* Module declarations from 'ancillary_data' */
static PyTypeObject *__pyx_array_type = 0;
static PyTypeObject *__pyx_MemviewEnum_type = 0;
static PyTypeObject *__pyx_memoryview_type = 0;
static PyTypeObject *__pyx_memoryviewslice_type = 0;
static PyObject *generic = 0;
static PyObject *strided = 0;
static PyObject *indirect = 0;
static PyObject *contiguous = 0;
static PyObject *indirect_contiguous = 0;
static int __pyx_memoryview_thread_locks_used;
static PyThread_type_lock __pyx_memoryview_thread_locks[8];
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2318 2319 2320 2321 2322 2323 
__PYX_EXTERN_C DL_IMPORT(void) get_Reynolds_SST(float *, float *, int, char *, char *, float *); /*proto*/
__PYX_EXTERN_C DL_IMPORT(void) get_NDVI_background(float *, float *, int, char *, char *, float *); /*proto*/
__PYX_EXTERN_C DL_IMPORT(void) get_Olson_eco(float *, float *, int, char *, unsigned char *); /*proto*/
__PYX_EXTERN_C DL_IMPORT(void) get_GEOS(float *, float *, int, char *, char *, char *, char *, char *, char *, char *, float *, float *, float *, float *, float *, float *); /*proto*/
__PYX_EXTERN_C DL_IMPORT(float) cithr(int, float, float); /*proto*/
__PYX_EXTERN_C DL_IMPORT(void) check_reg_uniformity(int, int, int, int, int, unsigned char, unsigned char *, float *, float *, unsigned char *, int *, int *, int *, int *); /*proto*/
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2324 2325 2326 2327 2328 2329 2330 2331 2332 2333 2334 2335 2336 2337 2338 2339 2340 2341 2342 2343 2344 2345 2346 2347 2348 2349 2350 2351 2352 2353 2354 2355 2356 2357 2358 
static struct __pyx_array_obj *__pyx_array_new(PyObject *, Py_ssize_t, char *, char *, char *); /*proto*/
static void *__pyx_align_pointer(void *, size_t); /*proto*/
static PyObject *__pyx_memoryview_new(PyObject *, int, int, __Pyx_TypeInfo *); /*proto*/
static CYTHON_INLINE int __pyx_memoryview_check(PyObject *); /*proto*/
static PyObject *_unellipsify(PyObject *, int); /*proto*/
static PyObject *assert_direct_dimensions(Py_ssize_t *, int); /*proto*/
static struct __pyx_memoryview_obj *__pyx_memview_slice(struct __pyx_memoryview_obj *, PyObject *); /*proto*/
static int __pyx_memoryview_slice_memviewslice(__Pyx_memviewslice *, Py_ssize_t, Py_ssize_t, Py_ssize_t, int, int, int *, Py_ssize_t, Py_ssize_t, Py_ssize_t, int, int, int, int); /*proto*/
static char *__pyx_pybuffer_index(Py_buffer *, char *, Py_ssize_t, Py_ssize_t); /*proto*/
static int __pyx_memslice_transpose(__Pyx_memviewslice *); /*proto*/
static PyObject *__pyx_memoryview_fromslice(__Pyx_memviewslice, int, PyObject *(*)(char *), int (*)(char *, PyObject *), int); /*proto*/
static __Pyx_memviewslice *__pyx_memoryview_get_slice_from_memoryview(struct __pyx_memoryview_obj *, __Pyx_memviewslice *); /*proto*/
static void __pyx_memoryview_slice_copy(struct __pyx_memoryview_obj *, __Pyx_memviewslice *); /*proto*/
static PyObject *__pyx_memoryview_copy_object(struct __pyx_memoryview_obj *); /*proto*/
static PyObject *__pyx_memoryview_copy_object_from_slice(struct __pyx_memoryview_obj *, __Pyx_memviewslice *); /*proto*/
static Py_ssize_t abs_py_ssize_t(Py_ssize_t); /*proto*/
static char __pyx_get_best_slice_order(__Pyx_memviewslice *, int); /*proto*/
static void _copy_strided_to_strided(char *, Py_ssize_t *, char *, Py_ssize_t *, Py_ssize_t *, Py_ssize_t *, int, size_t); /*proto*/
static void copy_strided_to_strided(__Pyx_memviewslice *, __Pyx_memviewslice *, int, size_t); /*proto*/
static Py_ssize_t __pyx_memoryview_slice_get_size(__Pyx_memviewslice *, int); /*proto*/
static Py_ssize_t __pyx_fill_contig_strides_array(Py_ssize_t *, Py_ssize_t *, Py_ssize_t, int, char); /*proto*/
static void *__pyx_memoryview_copy_data_to_temp(__Pyx_memviewslice *, __Pyx_memviewslice *, char, int); /*proto*/
static int __pyx_memoryview_err_extents(int, Py_ssize_t, Py_ssize_t); /*proto*/
static int __pyx_memoryview_err_dim(PyObject *, char *, int); /*proto*/
static int __pyx_memoryview_err(PyObject *, char *); /*proto*/
static int __pyx_memoryview_copy_contents(__Pyx_memviewslice, __Pyx_memviewslice, int, int, int); /*proto*/
static void __pyx_memoryview_broadcast_leading(__Pyx_memviewslice *, int, int); /*proto*/
static void __pyx_memoryview_refcount_copying(__Pyx_memviewslice *, int, int, int); /*proto*/
static void __pyx_memoryview_refcount_objects_in_slice_with_gil(char *, Py_ssize_t *, Py_ssize_t *, int, int); /*proto*/
static void __pyx_memoryview_refcount_objects_in_slice(char *, Py_ssize_t *, Py_ssize_t *, int, int); /*proto*/
static void __pyx_memoryview_slice_assign_scalar(__Pyx_memviewslice *, int, size_t, void *, int); /*proto*/
static void __pyx_memoryview__slice_assign_scalar(char *, Py_ssize_t *, Py_ssize_t *, int, size_t, void *); /*proto*/
static PyObject *__pyx_unpickle_Enum__set_state(struct __pyx_MemviewEnum_obj *, PyObject *); /*proto*/
static __Pyx_TypeInfo __Pyx_TypeInfo_float = { "float", NULL, sizeof(float), { 0 }, 0, 'R', 0, 0 };
static __Pyx_TypeInfo __Pyx_TypeInfo_unsigned_char = { "unsigned char", NULL, sizeof(unsigned char), { 0 }, 0, IS_UNSIGNED(unsigned char) ? 'U' : 'I', IS_UNSIGNED(unsigned char), 0 };
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2359 
static __Pyx_TypeInfo __Pyx_TypeInfo_int = { "int", NULL, sizeof(int), { 0 }, 0, IS_UNSIGNED(int) ? 'U' : 'I', IS_UNSIGNED(int), 0 };
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2360 2361 2362 2363 2364 
#define __Pyx_MODULE_NAME "ancillary_data"
extern int __pyx_module_is_main_ancillary_data;
int __pyx_module_is_main_ancillary_data = 0;

/* Implementation of 'ancillary_data' */
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2365 
static PyObject *__pyx_builtin_range;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2366 2367 2368 2369 2370 2371 2372 2373 
static PyObject *__pyx_builtin_ImportError;
static PyObject *__pyx_builtin_ValueError;
static PyObject *__pyx_builtin_MemoryError;
static PyObject *__pyx_builtin_enumerate;
static PyObject *__pyx_builtin_TypeError;
static PyObject *__pyx_builtin_Ellipsis;
static PyObject *__pyx_builtin_id;
static PyObject *__pyx_builtin_IndexError;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2374 
static const char __pyx_k_C[] = "C";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2375 2376 
static const char __pyx_k_O[] = "O";
static const char __pyx_k_c[] = "c";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2377 2378 
static const char __pyx_k_i[] = "i";
static const char __pyx_k_j[] = "j";
Paolo Veglio's avatar
function to determine scene flags created  
Paolo Veglio committed
2379 
static const char __pyx_k_v[] = "v";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2380 2381 2382 
static const char __pyx_k_id[] = "id";
static const char __pyx_k_np[] = "np";
static const char __pyx_k_eco[] = "eco";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2383 
static const char __pyx_k_key[] = "key";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2384 2385 
static const char __pyx_k_lat[] = "lat";
static const char __pyx_k_lon[] = "lon";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2386 
static const char __pyx_k_lsf[] = "lsf";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2387 2388 
static const char __pyx_k_new[] = "__new__";
static const char __pyx_k_obj[] = "obj";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2389 
static const char __pyx_k_res[] = "res";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2390 2391 2392 
static const char __pyx_k_sst[] = "sst";
static const char __pyx_k_tpw[] = "tpw";
static const char __pyx_k_base[] = "base";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2393 
static const char __pyx_k_bt11[] = "bt11";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2394 
static const char __pyx_k_dict[] = "__dict__";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2395 2396 
static const char __pyx_k_eles[] = "eles";
static const char __pyx_k_land[] = "land";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2397 2398 2399 2400 2401 2402 2403 2404 2405 2406 2407 2408 2409 2410 
static const char __pyx_k_main[] = "__main__";
static const char __pyx_k_mode[] = "mode";
static const char __pyx_k_name[] = "name";
static const char __pyx_k_ndim[] = "ndim";
static const char __pyx_k_ndvi[] = "ndvi";
static const char __pyx_k_pack[] = "pack";
static const char __pyx_k_sfct[] = "sfct";
static const char __pyx_k_size[] = "size";
static const char __pyx_k_step[] = "step";
static const char __pyx_k_stop[] = "stop";
static const char __pyx_k_test[] = "__test__";
static const char __pyx_k_ASCII[] = "ASCII";
static const char __pyx_k_DTYPE[] = "DTYPE";
static const char __pyx_k_class[] = "__class__";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2411 2412 
static const char __pyx_k_coast[] = "coast";
static const char __pyx_k_dtype[] = "dtype";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2413 2414 2415 2416 2417 2418 
static const char __pyx_k_error[] = "error";
static const char __pyx_k_flags[] = "flags";
static const char __pyx_k_float[] = "float";
static const char __pyx_k_geos1[] = "geos1";
static const char __pyx_k_geos2[] = "geos2";
static const char __pyx_k_icefr[] = "icefr";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2419 2420 
static const char __pyx_k_int32[] = "int32";
static const char __pyx_k_lines[] = "lines";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2421 2422 
static const char __pyx_k_numpy[] = "numpy";
static const char __pyx_k_ocnfr[] = "ocnfr";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2423 
static const char __pyx_k_order[] = "order";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2424 2425 2426 
static const char __pyx_k_range[] = "range";
static const char __pyx_k_shape[] = "shape";
static const char __pyx_k_start[] = "start";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2427 2428 
static const char __pyx_k_water[] = "water";
static const char __pyx_k_zeros[] = "zeros";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2429 2430 2431 2432 
static const char __pyx_k_eco_mv[] = "eco_mv";
static const char __pyx_k_encode[] = "encode";
static const char __pyx_k_format[] = "format";
static const char __pyx_k_import[] = "__import__";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2433 
static const char __pyx_k_land_r[] = "land_r";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2434 2435 2436 2437 2438 2439 2440 2441 2442 2443 
static const char __pyx_k_name_2[] = "__name__";
static const char __pyx_k_pickle[] = "pickle";
static const char __pyx_k_reduce[] = "__reduce__";
static const char __pyx_k_snowfr[] = "snowfr";
static const char __pyx_k_sst_mv[] = "sst_mv";
static const char __pyx_k_struct[] = "struct";
static const char __pyx_k_tpw_mv[] = "tpw_mv";
static const char __pyx_k_unpack[] = "unpack";
static const char __pyx_k_update[] = "update";
static const char __pyx_k_anc_dir[] = "anc_dir";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2444 
static const char __pyx_k_coast_r[] = "coast_r";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2445 
static const char __pyx_k_fortran[] = "fortran";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2446 
static const char __pyx_k_land_mv[] = "land_mv";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2447 2448 
static const char __pyx_k_memview[] = "memview";
static const char __pyx_k_ndvi_mv[] = "ndvi_mv";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2449 2450 
static const char __pyx_k_satname[] = "satname";
static const char __pyx_k_sec_vza[] = "sec_vza";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2451 
static const char __pyx_k_sfct_mv[] = "sfct_mv";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2452 2453 
static const char __pyx_k_tci_thr[] = "tci_thr";
static const char __pyx_k_water_r[] = "water_r";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2454 
static const char __pyx_k_Ellipsis[] = "Ellipsis";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2455 2456 
static const char __pyx_k_coast_mv[] = "coast_mv";
static const char __pyx_k_eco_type[] = "eco_type";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2457 2458 2459 2460 2461 2462 
static const char __pyx_k_geos_lnd[] = "geos_lnd";
static const char __pyx_k_geos_ocn[] = "geos_ocn";
static const char __pyx_k_getstate[] = "__getstate__";
static const char __pyx_k_icefr_mv[] = "icefr_mv";
static const char __pyx_k_itemsize[] = "itemsize";
static const char __pyx_k_ocnfr_mv[] = "ocnfr_mv";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2463 
static const char __pyx_k_py_cithr[] = "py_cithr";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2464 2465 2466 
static const char __pyx_k_pyx_type[] = "__pyx_type";
static const char __pyx_k_setstate[] = "__setstate__";
static const char __pyx_k_sst_file[] = "sst_file";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2467 
static const char __pyx_k_water_mv[] = "water_mv";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2468 
static const char __pyx_k_TypeError[] = "TypeError";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2469 
static const char __pyx_k_elem_edge[] = "elem_edge";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2470 2471 
static const char __pyx_k_enumerate[] = "enumerate";
static const char __pyx_k_geos_cnst[] = "geos_cnst";
Paolo Veglio's avatar
function to determine scene flags created  
Paolo Veglio committed
2472 2473 
static const char __pyx_k_geos_data[] = "geos_data";
static const char __pyx_k_geos_dict[] = "geos_dict";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2474 
static const char __pyx_k_landicefr[] = "landicefr";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2475 
static const char __pyx_k_line_edge[] = "line_edge";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2476 2477 2478 2479 2480 2481 2482 2483 2484 2485 2486 2487 
static const char __pyx_k_ndvi_file[] = "ndvi_file";
static const char __pyx_k_pyx_state[] = "__pyx_state";
static const char __pyx_k_reduce_ex[] = "__reduce_ex__";
static const char __pyx_k_snowfr_mv[] = "snowfr_mv";
static const char __pyx_k_startTime[] = "startTime";
static const char __pyx_k_IndexError[] = "IndexError";
static const char __pyx_k_ValueError[] = "ValueError";
static const char __pyx_k_pyx_result[] = "__pyx_result";
static const char __pyx_k_pyx_vtable[] = "__pyx_vtable__";
static const char __pyx_k_ImportError[] = "ImportError";
static const char __pyx_k_MemoryError[] = "MemoryError";
static const char __pyx_k_PickleError[] = "PickleError";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2488 
static const char __pyx_k_loc_uniform[] = "loc_uniform";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2489 2490 2491 
static const char __pyx_k_py_get_GEOS[] = "py_get_GEOS";
static const char __pyx_k_C_CONTIGUOUS[] = "C_CONTIGUOUS";
static const char __pyx_k_landicefr_mv[] = "landicefr_mv";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2492 
static const char __pyx_k_py_snow_mask[] = "py_snow_mask";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2493 2494 2495 
static const char __pyx_k_pyx_checksum[] = "__pyx_checksum";
static const char __pyx_k_stringsource[] = "stringsource";
static const char __pyx_k_ancillary_pyx[] = "ancillary.pyx";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2496 
static const char __pyx_k_loc_uniform_r[] = "loc_uniform_r";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2497 2498 2499 
static const char __pyx_k_pyx_getbuffer[] = "__pyx_getbuffer";
static const char __pyx_k_reduce_cython[] = "__reduce_cython__";
static const char __pyx_k_ancillary_data[] = "ancillary_data";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2500 
static const char __pyx_k_loc_uniform_mv[] = "loc_uniform_mv";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2501 2502 2503 2504 2505 2506 
static const char __pyx_k_View_MemoryView[] = "View.MemoryView";
static const char __pyx_k_allocate_buffer[] = "allocate_buffer";
static const char __pyx_k_dtype_is_object[] = "dtype_is_object";
static const char __pyx_k_pyx_PickleError[] = "__pyx_PickleError";
static const char __pyx_k_setstate_cython[] = "__setstate_cython__";
static const char __pyx_k_py_get_Olson_eco[] = "py_get_Olson_eco";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2507 
static const char __pyx_k_scene_uniformity[] = "scene_uniformity";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2508 2509 2510 2511 2512 2513 2514 2515 2516 2517 2518 
static const char __pyx_k_ascontiguousarray[] = "ascontiguousarray";
static const char __pyx_k_pyx_unpickle_Enum[] = "__pyx_unpickle_Enum";
static const char __pyx_k_cline_in_traceback[] = "cline_in_traceback";
static const char __pyx_k_strided_and_direct[] = "<strided and direct>";
static const char __pyx_k_py_get_Reynolds_SST[] = "py_get_Reynolds_SST";
static const char __pyx_k_strided_and_indirect[] = "<strided and indirect>";
static const char __pyx_k_contiguous_and_direct[] = "<contiguous and direct>";
static const char __pyx_k_MemoryView_of_r_object[] = "<MemoryView of %r object>";
static const char __pyx_k_py_get_NDVI_background[] = "py_get_NDVI_background";
static const char __pyx_k_MemoryView_of_r_at_0x_x[] = "<MemoryView of %r at 0x%x>";
static const char __pyx_k_contiguous_and_indirect[] = "<contiguous and indirect>";
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2519 
static const char __pyx_k_py_check_reg_uniformity[] = "py_check_reg_uniformity";
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2520 2521 2522 2523 2524 2525 2526 2527 2528 2529 2530 2531 2532 2533 2534 2535 2536 2537 2538 2539 2540 2541 
static const char __pyx_k_Cannot_index_with_type_s[] = "Cannot index with type '%s'";
static const char __pyx_k_Invalid_shape_in_axis_d_d[] = "Invalid shape in axis %d: %d.";
static const char __pyx_k_itemsize_0_for_cython_array[] = "itemsize <= 0 for cython.array";
static const char __pyx_k_unable_to_allocate_array_data[] = "unable to allocate array data.";
static const char __pyx_k_strided_and_direct_or_indirect[] = "<strided and direct or indirect>";
static const char __pyx_k_numpy_core_multiarray_failed_to[] = "numpy.core.multiarray failed to import";
static const char __pyx_k_Buffer_view_does_not_expose_stri[] = "Buffer view does not expose strides";
static const char __pyx_k_Can_only_create_a_buffer_that_is[] = "Can only create a buffer that is contiguous in memory.";
static const char __pyx_k_Cannot_assign_to_read_only_memor[] = "Cannot assign to read-only memoryview";
static const char __pyx_k_Cannot_create_writable_memory_vi[] = "Cannot create writable memory view from read-only memoryview";
static const char __pyx_k_Empty_shape_tuple_for_cython_arr[] = "Empty shape tuple for cython.array";
static const char __pyx_k_Incompatible_checksums_s_vs_0xb0[] = "Incompatible checksums (%s vs 0xb068931 = (name))";
static const char __pyx_k_Indirect_dimensions_not_supporte[] = "Indirect dimensions not supported";
static const char __pyx_k_Invalid_mode_expected_c_or_fortr[] = "Invalid mode, expected 'c' or 'fortran', got %s";
static const char __pyx_k_Out_of_bounds_on_buffer_access_a[] = "Out of bounds on buffer access (axis %d)";
static const char __pyx_k_Unable_to_convert_item_to_object[] = "Unable to convert item to object";
static const char __pyx_k_got_differing_extents_in_dimensi[] = "got differing extents in dimension %d (got %d and %d)";
static const char __pyx_k_no_default___reduce___due_to_non[] = "no default __reduce__ due to non-trivial __cinit__";
static const char __pyx_k_numpy_core_umath_failed_to_impor[] = "numpy.core.umath failed to import";
static const char __pyx_k_unable_to_allocate_shape_and_str[] = "unable to allocate shape and strides.";
static PyObject *__pyx_n_s_ASCII;
static PyObject *__pyx_kp_s_Buffer_view_does_not_expose_stri;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2542 
static PyObject *__pyx_n_u_C;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2543 2544 2545 2546 2547 2548 2549 2550 2551 2552 2553 2554 2555 2556 2557 2558 2559 2560 2561 2562 2563 2564 2565 2566 2567 2568 2569 2570 2571 2572 
static PyObject *__pyx_n_u_C_CONTIGUOUS;
static PyObject *__pyx_kp_s_Can_only_create_a_buffer_that_is;
static PyObject *__pyx_kp_s_Cannot_assign_to_read_only_memor;
static PyObject *__pyx_kp_s_Cannot_create_writable_memory_vi;
static PyObject *__pyx_kp_s_Cannot_index_with_type_s;
static PyObject *__pyx_n_s_DTYPE;
static PyObject *__pyx_n_s_Ellipsis;
static PyObject *__pyx_kp_s_Empty_shape_tuple_for_cython_arr;
static PyObject *__pyx_n_s_ImportError;
static PyObject *__pyx_kp_s_Incompatible_checksums_s_vs_0xb0;
static PyObject *__pyx_n_s_IndexError;
static PyObject *__pyx_kp_s_Indirect_dimensions_not_supporte;
static PyObject *__pyx_kp_s_Invalid_mode_expected_c_or_fortr;
static PyObject *__pyx_kp_s_Invalid_shape_in_axis_d_d;
static PyObject *__pyx_n_s_MemoryError;
static PyObject *__pyx_kp_s_MemoryView_of_r_at_0x_x;
static PyObject *__pyx_kp_s_MemoryView_of_r_object;
static PyObject *__pyx_n_b_O;
static PyObject *__pyx_kp_s_Out_of_bounds_on_buffer_access_a;
static PyObject *__pyx_n_s_PickleError;
static PyObject *__pyx_n_s_TypeError;
static PyObject *__pyx_kp_s_Unable_to_convert_item_to_object;
static PyObject *__pyx_n_s_ValueError;
static PyObject *__pyx_n_s_View_MemoryView;
static PyObject *__pyx_n_s_allocate_buffer;
static PyObject *__pyx_n_s_anc_dir;
static PyObject *__pyx_n_s_ancillary_data;
static PyObject *__pyx_kp_s_ancillary_pyx;
static PyObject *__pyx_n_s_ascontiguousarray;
static PyObject *__pyx_n_s_base;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2573 
static PyObject *__pyx_n_s_bt11;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2574 2575 2576 2577 
static PyObject *__pyx_n_s_c;
static PyObject *__pyx_n_u_c;
static PyObject *__pyx_n_s_class;
static PyObject *__pyx_n_s_cline_in_traceback;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2578 2579 2580 2581 
static PyObject *__pyx_n_s_coast;
static PyObject *__pyx_n_u_coast;
static PyObject *__pyx_n_s_coast_mv;
static PyObject *__pyx_n_s_coast_r;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2582 2583 2584 
static PyObject *__pyx_kp_s_contiguous_and_direct;
static PyObject *__pyx_kp_s_contiguous_and_indirect;
static PyObject *__pyx_n_s_dict;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2585 
static PyObject *__pyx_n_s_dtype;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2586 2587 2588 
static PyObject *__pyx_n_s_dtype_is_object;
static PyObject *__pyx_n_s_eco;
static PyObject *__pyx_n_s_eco_mv;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2589 2590 2591 
static PyObject *__pyx_n_s_eco_type;
static PyObject *__pyx_n_s_elem_edge;
static PyObject *__pyx_n_s_eles;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2592 2593 2594 2595 2596 2597 2598 2599 2600 2601 2602 
static PyObject *__pyx_n_s_encode;
static PyObject *__pyx_n_s_enumerate;
static PyObject *__pyx_n_s_error;
static PyObject *__pyx_n_s_flags;
static PyObject *__pyx_n_s_float;
static PyObject *__pyx_n_s_format;
static PyObject *__pyx_n_s_fortran;
static PyObject *__pyx_n_u_fortran;
static PyObject *__pyx_n_s_geos1;
static PyObject *__pyx_n_s_geos2;
static PyObject *__pyx_n_s_geos_cnst;
Paolo Veglio's avatar
function to determine scene flags created  
Paolo Veglio committed
2603 2604 
static PyObject *__pyx_n_s_geos_data;
static PyObject *__pyx_n_s_geos_dict;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2605 2606 2607 2608 
static PyObject *__pyx_n_s_geos_lnd;
static PyObject *__pyx_n_s_geos_ocn;
static PyObject *__pyx_n_s_getstate;
static PyObject *__pyx_kp_s_got_differing_extents_in_dimensi;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2609 2610 
static PyObject *__pyx_n_s_i;
static PyObject *__pyx_n_s_icefr;
Paolo Veglio's avatar
function to determine scene flags created  
Paolo Veglio committed
2611 
static PyObject *__pyx_n_u_icefr;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2612 2613 2614 
static PyObject *__pyx_n_s_icefr_mv;
static PyObject *__pyx_n_s_id;
static PyObject *__pyx_n_s_import;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2615 
static PyObject *__pyx_n_s_int32;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2616 2617 
static PyObject *__pyx_n_s_itemsize;
static PyObject *__pyx_kp_s_itemsize_0_for_cython_array;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2618 2619 2620 2621 2622 2623 
static PyObject *__pyx_n_s_j;
static PyObject *__pyx_n_s_key;
static PyObject *__pyx_n_s_land;
static PyObject *__pyx_n_u_land;
static PyObject *__pyx_n_s_land_mv;
static PyObject *__pyx_n_s_land_r;
Paolo Veglio's avatar
function to determine scene flags created  
Paolo Veglio committed
2624 
static PyObject *__pyx_n_u_landicefr;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2625 2626 
static PyObject *__pyx_n_s_landicefr_mv;
static PyObject *__pyx_n_s_lat;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2627 2628 2629 2630 2631 2632 
static PyObject *__pyx_n_s_line_edge;
static PyObject *__pyx_n_s_lines;
static PyObject *__pyx_n_s_loc_uniform;
static PyObject *__pyx_n_u_loc_uniform;
static PyObject *__pyx_n_s_loc_uniform_mv;
static PyObject *__pyx_n_s_loc_uniform_r;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2633 
static PyObject *__pyx_n_s_lon;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2634 
static PyObject *__pyx_n_s_lsf;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2635 2636 2637 2638 2639 2640 2641 2642 2643 2644 2645 2646 2647 2648 2649 2650 
static PyObject *__pyx_n_s_main;
static PyObject *__pyx_n_s_memview;
static PyObject *__pyx_n_s_mode;
static PyObject *__pyx_n_s_name;
static PyObject *__pyx_n_s_name_2;
static PyObject *__pyx_n_s_ndim;
static PyObject *__pyx_n_s_ndvi;
static PyObject *__pyx_n_s_ndvi_file;
static PyObject *__pyx_n_s_ndvi_mv;
static PyObject *__pyx_n_s_new;
static PyObject *__pyx_kp_s_no_default___reduce___due_to_non;
static PyObject *__pyx_n_s_np;
static PyObject *__pyx_n_s_numpy;
static PyObject *__pyx_kp_u_numpy_core_multiarray_failed_to;
static PyObject *__pyx_kp_u_numpy_core_umath_failed_to_impor;
static PyObject *__pyx_n_s_obj;
Paolo Veglio's avatar
function to determine scene flags created  
Paolo Veglio committed
2651 
static PyObject *__pyx_n_u_ocnfr;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2652 
static PyObject *__pyx_n_s_ocnfr_mv;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2653 
static PyObject *__pyx_n_s_order;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2654 2655 
static PyObject *__pyx_n_s_pack;
static PyObject *__pyx_n_s_pickle;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2656 2657 
static PyObject *__pyx_n_s_py_check_reg_uniformity;
static PyObject *__pyx_n_s_py_cithr;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2658 2659 2660 2661 
static PyObject *__pyx_n_s_py_get_GEOS;
static PyObject *__pyx_n_s_py_get_NDVI_background;
static PyObject *__pyx_n_s_py_get_Olson_eco;
static PyObject *__pyx_n_s_py_get_Reynolds_SST;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2662 
static PyObject *__pyx_n_s_py_snow_mask;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2663 2664 2665 2666 2667 2668 2669 2670 2671 2672 2673 2674 
static PyObject *__pyx_n_s_pyx_PickleError;
static PyObject *__pyx_n_s_pyx_checksum;
static PyObject *__pyx_n_s_pyx_getbuffer;
static PyObject *__pyx_n_s_pyx_result;
static PyObject *__pyx_n_s_pyx_state;
static PyObject *__pyx_n_s_pyx_type;
static PyObject *__pyx_n_s_pyx_unpickle_Enum;
static PyObject *__pyx_n_s_pyx_vtable;
static PyObject *__pyx_n_s_range;
static PyObject *__pyx_n_s_reduce;
static PyObject *__pyx_n_s_reduce_cython;
static PyObject *__pyx_n_s_reduce_ex;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2675 2676 2677 2678 
static PyObject *__pyx_n_s_res;
static PyObject *__pyx_n_s_satname;
static PyObject *__pyx_n_s_scene_uniformity;
static PyObject *__pyx_n_s_sec_vza;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2679 2680 
static PyObject *__pyx_n_s_setstate;
static PyObject *__pyx_n_s_setstate_cython;
Paolo Veglio's avatar
function to determine scene flags created  
Paolo Veglio committed
2681 
static PyObject *__pyx_n_u_sfct;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2682 2683 2684 
static PyObject *__pyx_n_s_sfct_mv;
static PyObject *__pyx_n_s_shape;
static PyObject *__pyx_n_s_size;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2685 
static PyObject *__pyx_n_s_snowfr;
Paolo Veglio's avatar
function to determine scene flags created  
Paolo Veglio committed
2686 
static PyObject *__pyx_n_u_snowfr;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2687 2688 2689 2690 2691 2692 2693 2694 2695 2696 2697 2698 2699 
static PyObject *__pyx_n_s_snowfr_mv;
static PyObject *__pyx_n_s_sst;
static PyObject *__pyx_n_s_sst_file;
static PyObject *__pyx_n_s_sst_mv;
static PyObject *__pyx_n_s_start;
static PyObject *__pyx_n_s_startTime;
static PyObject *__pyx_n_s_step;
static PyObject *__pyx_n_s_stop;
static PyObject *__pyx_kp_s_strided_and_direct;
static PyObject *__pyx_kp_s_strided_and_direct_or_indirect;
static PyObject *__pyx_kp_s_strided_and_indirect;
static PyObject *__pyx_kp_s_stringsource;
static PyObject *__pyx_n_s_struct;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2700 
static PyObject *__pyx_n_s_tci_thr;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2701 
static PyObject *__pyx_n_s_test;
Paolo Veglio's avatar
function to determine scene flags created  
Paolo Veglio committed
2702 
static PyObject *__pyx_n_u_tpw;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2703 2704 2705 2706 2707 
static PyObject *__pyx_n_s_tpw_mv;
static PyObject *__pyx_kp_s_unable_to_allocate_array_data;
static PyObject *__pyx_kp_s_unable_to_allocate_shape_and_str;
static PyObject *__pyx_n_s_unpack;
static PyObject *__pyx_n_s_update;
Paolo Veglio's avatar
function to determine scene flags created  
Paolo Veglio committed
2708 
static PyObject *__pyx_n_s_v;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2709 2710 2711 2712 2713 2714 2715 2716 2717 2718 2719 2720 
static PyObject *__pyx_n_s_water;
static PyObject *__pyx_n_u_water;
static PyObject *__pyx_n_s_water_mv;
static PyObject *__pyx_n_s_water_r;
static PyObject *__pyx_n_s_zeros;
static PyObject *__pyx_pf_14ancillary_data_py_get_Reynolds_SST(CYTHON_UNUSED PyObject *__pyx_self, PyArrayObject *__pyx_v_lat, PyArrayObject *__pyx_v_lon, PyObject *__pyx_v_res, char *__pyx_v_anc_dir, char *__pyx_v_sst_file, PyObject *__pyx_v_sst); /* proto */
static PyObject *__pyx_pf_14ancillary_data_2py_get_NDVI_background(CYTHON_UNUSED PyObject *__pyx_self, PyArrayObject *__pyx_v_lat, PyArrayObject *__pyx_v_lon, PyObject *__pyx_v_res, char *__pyx_v_anc_dir, char *__pyx_v_ndvi_file, PyObject *__pyx_v_ndvi); /* proto */
static PyObject *__pyx_pf_14ancillary_data_4py_get_Olson_eco(CYTHON_UNUSED PyObject *__pyx_self, PyArrayObject *__pyx_v_lat, PyArrayObject *__pyx_v_lon, PyObject *__pyx_v_res, char *__pyx_v_anc_dir, PyObject *__pyx_v_eco); /* proto */
static PyObject *__pyx_pf_14ancillary_data_6py_get_GEOS(CYTHON_UNUSED PyObject *__pyx_self, PyArrayObject *__pyx_v_lat, PyArrayObject *__pyx_v_lon, int __pyx_v_res, char *__pyx_v_startTime, char *__pyx_v_anc_dir, char *__pyx_v_geos1, char *__pyx_v_geos2, char *__pyx_v_geos_lnd, char *__pyx_v_geos_ocn, char *__pyx_v_geos_cnst, PyObject *__pyx_v_geos_data); /* proto */
static PyObject *__pyx_pf_14ancillary_data_8py_snow_mask(CYTHON_UNUSED PyObject *__pyx_self, CYTHON_UNUSED char *__pyx_v_satname, CYTHON_UNUSED unsigned char __pyx_v_lsf); /* proto */
static PyObject *__pyx_pf_14ancillary_data_10py_cithr(CYTHON_UNUSED PyObject *__pyx_self, int __pyx_v_key, PyArrayObject *__pyx_v_sec_vza, PyArrayObject *__pyx_v_bt11); /* proto */
static PyObject *__pyx_pf_14ancillary_data_12py_check_reg_uniformity(CYTHON_UNUSED PyObject *__pyx_self, PyArrayObject *__pyx_v_eco_type, PyArrayObject *__pyx_v_eco, PyArrayObject *__pyx_v_snowfr, PyArrayObject *__pyx_v_icefr, PyArrayObject *__pyx_v_lsf); /* proto */
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2721 2722 2723 2724 2725 2726 2727 2728 2729 2730 2731 2732 2733 2734 2735 2736 2737 2738 2739 2740 2741 2742 2743 2744 2745 2746 2747 2748 2749 2750 2751 2752 2753 2754 2755 2756 2757 2758 2759 2760 2761 2762 2763 2764 2765 2766 2767 2768 2769 2770 2771 2772 2773 2774 2775 2776 2777 2778 2779 2780 2781 2782 2783 2784 2785 2786 2787 2788 2789 2790 2791 2792 2793 2794 2795 2796 2797 
static int __pyx_array___pyx_pf_15View_dot_MemoryView_5array___cinit__(struct __pyx_array_obj *__pyx_v_self, PyObject *__pyx_v_shape, Py_ssize_t __pyx_v_itemsize, PyObject *__pyx_v_format, PyObject *__pyx_v_mode, int __pyx_v_allocate_buffer); /* proto */
static int __pyx_array___pyx_pf_15View_dot_MemoryView_5array_2__getbuffer__(struct __pyx_array_obj *__pyx_v_self, Py_buffer *__pyx_v_info, int __pyx_v_flags); /* proto */
static void __pyx_array___pyx_pf_15View_dot_MemoryView_5array_4__dealloc__(struct __pyx_array_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView_5array_7memview___get__(struct __pyx_array_obj *__pyx_v_self); /* proto */
static Py_ssize_t __pyx_array___pyx_pf_15View_dot_MemoryView_5array_6__len__(struct __pyx_array_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_array___pyx_pf_15View_dot_MemoryView_5array_8__getattr__(struct __pyx_array_obj *__pyx_v_self, PyObject *__pyx_v_attr); /* proto */
static PyObject *__pyx_array___pyx_pf_15View_dot_MemoryView_5array_10__getitem__(struct __pyx_array_obj *__pyx_v_self, PyObject *__pyx_v_item); /* proto */
static int __pyx_array___pyx_pf_15View_dot_MemoryView_5array_12__setitem__(struct __pyx_array_obj *__pyx_v_self, PyObject *__pyx_v_item, PyObject *__pyx_v_value); /* proto */
static PyObject *__pyx_pf___pyx_array___reduce_cython__(CYTHON_UNUSED struct __pyx_array_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf___pyx_array_2__setstate_cython__(CYTHON_UNUSED struct __pyx_array_obj *__pyx_v_self, CYTHON_UNUSED PyObject *__pyx_v___pyx_state); /* proto */
static int __pyx_MemviewEnum___pyx_pf_15View_dot_MemoryView_4Enum___init__(struct __pyx_MemviewEnum_obj *__pyx_v_self, PyObject *__pyx_v_name); /* proto */
static PyObject *__pyx_MemviewEnum___pyx_pf_15View_dot_MemoryView_4Enum_2__repr__(struct __pyx_MemviewEnum_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf___pyx_MemviewEnum___reduce_cython__(struct __pyx_MemviewEnum_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf___pyx_MemviewEnum_2__setstate_cython__(struct __pyx_MemviewEnum_obj *__pyx_v_self, PyObject *__pyx_v___pyx_state); /* proto */
static int __pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview___cinit__(struct __pyx_memoryview_obj *__pyx_v_self, PyObject *__pyx_v_obj, int __pyx_v_flags, int __pyx_v_dtype_is_object); /* proto */
static void __pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview_2__dealloc__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview_4__getitem__(struct __pyx_memoryview_obj *__pyx_v_self, PyObject *__pyx_v_index); /* proto */
static int __pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview_6__setitem__(struct __pyx_memoryview_obj *__pyx_v_self, PyObject *__pyx_v_index, PyObject *__pyx_v_value); /* proto */
static int __pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview_8__getbuffer__(struct __pyx_memoryview_obj *__pyx_v_self, Py_buffer *__pyx_v_info, int __pyx_v_flags); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView_10memoryview_1T___get__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView_10memoryview_4base___get__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView_10memoryview_5shape___get__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView_10memoryview_7strides___get__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView_10memoryview_10suboffsets___get__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView_10memoryview_4ndim___get__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView_10memoryview_8itemsize___get__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView_10memoryview_6nbytes___get__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView_10memoryview_4size___get__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static Py_ssize_t __pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview_10__len__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview_12__repr__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview_14__str__(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview_16is_c_contig(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview_18is_f_contig(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview_20copy(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_memoryview___pyx_pf_15View_dot_MemoryView_10memoryview_22copy_fortran(struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf___pyx_memoryview___reduce_cython__(CYTHON_UNUSED struct __pyx_memoryview_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf___pyx_memoryview_2__setstate_cython__(CYTHON_UNUSED struct __pyx_memoryview_obj *__pyx_v_self, CYTHON_UNUSED PyObject *__pyx_v___pyx_state); /* proto */
static void __pyx_memoryviewslice___pyx_pf_15View_dot_MemoryView_16_memoryviewslice___dealloc__(struct __pyx_memoryviewslice_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView_16_memoryviewslice_4base___get__(struct __pyx_memoryviewslice_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf___pyx_memoryviewslice___reduce_cython__(CYTHON_UNUSED struct __pyx_memoryviewslice_obj *__pyx_v_self); /* proto */
static PyObject *__pyx_pf___pyx_memoryviewslice_2__setstate_cython__(CYTHON_UNUSED struct __pyx_memoryviewslice_obj *__pyx_v_self, CYTHON_UNUSED PyObject *__pyx_v___pyx_state); /* proto */
static PyObject *__pyx_pf_15View_dot_MemoryView___pyx_unpickle_Enum(CYTHON_UNUSED PyObject *__pyx_self, PyObject *__pyx_v___pyx_type, long __pyx_v___pyx_checksum, PyObject *__pyx_v___pyx_state); /* proto */
static PyObject *__pyx_tp_new_array(PyTypeObject *t, PyObject *a, PyObject *k); /*proto*/
static PyObject *__pyx_tp_new_Enum(PyTypeObject *t, PyObject *a, PyObject *k); /*proto*/
static PyObject *__pyx_tp_new_memoryview(PyTypeObject *t, PyObject *a, PyObject *k); /*proto*/
static PyObject *__pyx_tp_new__memoryviewslice(PyTypeObject *t, PyObject *a, PyObject *k); /*proto*/
static PyObject *__pyx_int_0;
static PyObject *__pyx_int_1;
static PyObject *__pyx_int_184977713;
static PyObject *__pyx_int_neg_1;
static PyObject *__pyx_tuple_;
static PyObject *__pyx_tuple__2;
static PyObject *__pyx_tuple__3;
static PyObject *__pyx_tuple__4;
static PyObject *__pyx_tuple__5;
static PyObject *__pyx_tuple__6;
static PyObject *__pyx_tuple__7;
static PyObject *__pyx_tuple__8;
static PyObject *__pyx_tuple__9;
static PyObject *__pyx_slice__17;
static PyObject *__pyx_tuple__10;
static PyObject *__pyx_tuple__11;
static PyObject *__pyx_tuple__12;
static PyObject *__pyx_tuple__13;
static PyObject *__pyx_tuple__14;
static PyObject *__pyx_tuple__15;
static PyObject *__pyx_tuple__16;
static PyObject *__pyx_tuple__18;
static PyObject *__pyx_tuple__19;
static PyObject *__pyx_tuple__20;
static PyObject *__pyx_tuple__21;
static PyObject *__pyx_tuple__23;
static PyObject *__pyx_tuple__25;
static PyObject *__pyx_tuple__27;
static PyObject *__pyx_tuple__29;
static PyObject *__pyx_tuple__31;
static PyObject *__pyx_tuple__33;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2798 2799 2800 2801 2802 2803 
static PyObject *__pyx_tuple__35;
static PyObject *__pyx_tuple__36;
static PyObject *__pyx_tuple__37;
static PyObject *__pyx_tuple__38;
static PyObject *__pyx_tuple__39;
static PyObject *__pyx_tuple__40;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2804 2805 2806 2807 
static PyObject *__pyx_codeobj__22;
static PyObject *__pyx_codeobj__24;
static PyObject *__pyx_codeobj__26;
static PyObject *__pyx_codeobj__28;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2808 2809 2810 2811 
static PyObject *__pyx_codeobj__30;
static PyObject *__pyx_codeobj__32;
static PyObject *__pyx_codeobj__34;
static PyObject *__pyx_codeobj__41;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2812 2813 
/* Late includes */
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2814 
/* "ancillary.pyx":28
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2815 2816 2817 
 * @cython.wraparound(False)
 * @cython.initializedcheck(False)
 * def py_get_Reynolds_SST(np.ndarray[float, ndim=1] lat,             # <<<<<<<<<<<<<<
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2818 
 *                         np.ndarray[float, ndim=1] lon, res,
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2819 2820 2821 2822 2823 2824 2825 2826 2827 
 *                         char *anc_dir, char *sst_file, sst):
 */

/* Python wrapper */
static PyObject *__pyx_pw_14ancillary_data_1py_get_Reynolds_SST(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/
static PyMethodDef __pyx_mdef_14ancillary_data_1py_get_Reynolds_SST = {"py_get_Reynolds_SST", (PyCFunction)(void*)(PyCFunctionWithKeywords)__pyx_pw_14ancillary_data_1py_get_Reynolds_SST, METH_VARARGS|METH_KEYWORDS, 0};
static PyObject *__pyx_pw_14ancillary_data_1py_get_Reynolds_SST(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) {
  PyArrayObject *__pyx_v_lat = 0;
  PyArrayObject *__pyx_v_lon = 0;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2828 
  PyObject *__pyx_v_res = 0;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2829 2830 2831 2832 2833 2834 2835 2836 2837 2838 
  char *__pyx_v_anc_dir;
  char *__pyx_v_sst_file;
  PyObject *__pyx_v_sst = 0;
  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;
  int __pyx_clineno = 0;
  PyObject *__pyx_r = 0;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("py_get_Reynolds_SST (wrapper)", 0);
  {
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2839 2840 
    static PyObject **__pyx_pyargnames[] = {&__pyx_n_s_lat,&__pyx_n_s_lon,&__pyx_n_s_res,&__pyx_n_s_anc_dir,&__pyx_n_s_sst_file,&__pyx_n_s_sst,0};
    PyObject* values[6] = {0,0,0,0,0,0};
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2841 2842 2843 2844 
    if (unlikely(__pyx_kwds)) {
      Py_ssize_t kw_args;
      const Py_ssize_t pos_args = PyTuple_GET_SIZE(__pyx_args);
      switch (pos_args) {
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2845 2846 
        case  6: values[5] = PyTuple_GET_ITEM(__pyx_args, 5);
        CYTHON_FALLTHROUGH;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2847 2848 2849 2850 2851 2852 2853 2854 2855 2856 2857 2858 2859 2860 2861 2862 2863 2864 2865 2866 2867 2868 
        case  5: values[4] = PyTuple_GET_ITEM(__pyx_args, 4);
        CYTHON_FALLTHROUGH;
        case  4: values[3] = PyTuple_GET_ITEM(__pyx_args, 3);
        CYTHON_FALLTHROUGH;
        case  3: values[2] = PyTuple_GET_ITEM(__pyx_args, 2);
        CYTHON_FALLTHROUGH;
        case  2: values[1] = PyTuple_GET_ITEM(__pyx_args, 1);
        CYTHON_FALLTHROUGH;
        case  1: values[0] = PyTuple_GET_ITEM(__pyx_args, 0);
        CYTHON_FALLTHROUGH;
        case  0: break;
        default: goto __pyx_L5_argtuple_error;
      }
      kw_args = PyDict_Size(__pyx_kwds);
      switch (pos_args) {
        case  0:
        if (likely((values[0] = __Pyx_PyDict_GetItemStr(__pyx_kwds, __pyx_n_s_lat)) != 0)) kw_args--;
        else goto __pyx_L5_argtuple_error;
        CYTHON_FALLTHROUGH;
        case  1:
        if (likely((values[1] = __Pyx_PyDict_GetItemStr(__pyx_kwds, __pyx_n_s_lon)) != 0)) kw_args--;
        else {
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2869 
          __Pyx_RaiseArgtupleInvalid("py_get_Reynolds_SST", 1, 6, 6, 1); __PYX_ERR(0, 28, __pyx_L3_error)
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2870 2871 2872 
        }
        CYTHON_FALLTHROUGH;
        case  2:
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2873 
        if (likely((values[2] = __Pyx_PyDict_GetItemStr(__pyx_kwds, __pyx_n_s_res)) != 0)) kw_args--;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2874 
        else {
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2875 
          __Pyx_RaiseArgtupleInvalid("py_get_Reynolds_SST", 1, 6, 6, 2); __PYX_ERR(0, 28, __pyx_L3_error)
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2876 2877 2878 
        }
        CYTHON_FALLTHROUGH;
        case  3:
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2879 
        if (likely((values[3] = __Pyx_PyDict_GetItemStr(__pyx_kwds, __pyx_n_s_anc_dir)) != 0)) kw_args--;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2880 
        else {
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2881 
          __Pyx_RaiseArgtupleInvalid("py_get_Reynolds_SST", 1, 6, 6, 3); __PYX_ERR(0, 28, __pyx_L3_error)
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2882 2883 2884 
        }
        CYTHON_FALLTHROUGH;
        case  4:
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2885 2886 2887 2888 2889 2890 2891 
        if (likely((values[4] = __Pyx_PyDict_GetItemStr(__pyx_kwds, __pyx_n_s_sst_file)) != 0)) kw_args--;
        else {
          __Pyx_RaiseArgtupleInvalid("py_get_Reynolds_SST", 1, 6, 6, 4); __PYX_ERR(0, 28, __pyx_L3_error)
        }
        CYTHON_FALLTHROUGH;
        case  5:
        if (likely((values[5] = __Pyx_PyDict_GetItemStr(__pyx_kwds, __pyx_n_s_sst)) != 0)) kw_args--;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2892 
        else {
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2893 
          __Pyx_RaiseArgtupleInvalid("py_get_Reynolds_SST", 1, 6, 6, 5); __PYX_ERR(0, 28, __pyx_L3_error)
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2894 2895 2896 
        }
      }
      if (unlikely(kw_args > 0)) {
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2897 
        if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_pyargnames, 0, values, pos_args, "py_get_Reynolds_SST") < 0)) __PYX_ERR(0, 28, __pyx_L3_error)
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2898 
      }
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2899 
    } else if (PyTuple_GET_SIZE(__pyx_args) != 6) {
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2900 2901 2902 2903 2904 2905 2906 
      goto __pyx_L5_argtuple_error;
    } else {
      values[0] = PyTuple_GET_ITEM(__pyx_args, 0);
      values[1] = PyTuple_GET_ITEM(__pyx_args, 1);
      values[2] = PyTuple_GET_ITEM(__pyx_args, 2);
      values[3] = PyTuple_GET_ITEM(__pyx_args, 3);
      values[4] = PyTuple_GET_ITEM(__pyx_args, 4);
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2907 
      values[5] = PyTuple_GET_ITEM(__pyx_args, 5);
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2908 2909 2910 
    }
    __pyx_v_lat = ((PyArrayObject *)values[0]);
    __pyx_v_lon = ((PyArrayObject *)values[1]);
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2911 2912 2913 2914 
    __pyx_v_res = values[2];
    __pyx_v_anc_dir = __Pyx_PyObject_AsWritableString(values[3]); if (unlikely((!__pyx_v_anc_dir) && PyErr_Occurred())) __PYX_ERR(0, 30, __pyx_L3_error)
    __pyx_v_sst_file = __Pyx_PyObject_AsWritableString(values[4]); if (unlikely((!__pyx_v_sst_file) && PyErr_Occurred())) __PYX_ERR(0, 30, __pyx_L3_error)
    __pyx_v_sst = values[5];
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2915 2916 2917 
  }
  goto __pyx_L4_argument_unpacking_done;
  __pyx_L5_argtuple_error:;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2918 
  __Pyx_RaiseArgtupleInvalid("py_get_Reynolds_SST", 1, 6, 6, PyTuple_GET_SIZE(__pyx_args)); __PYX_ERR(0, 28, __pyx_L3_error)
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2919 2920 2921 2922 2923 
  __pyx_L3_error:;
  __Pyx_AddTraceback("ancillary_data.py_get_Reynolds_SST", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __Pyx_RefNannyFinishContext();
  return NULL;
  __pyx_L4_argument_unpacking_done:;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2924 2925 2926 
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_lat), __pyx_ptype_5numpy_ndarray, 1, "lat", 0))) __PYX_ERR(0, 28, __pyx_L1_error)
  if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_lon), __pyx_ptype_5numpy_ndarray, 1, "lon", 0))) __PYX_ERR(0, 29, __pyx_L1_error)
  __pyx_r = __pyx_pf_14ancillary_data_py_get_Reynolds_SST(__pyx_self, __pyx_v_lat, __pyx_v_lon, __pyx_v_res, __pyx_v_anc_dir, __pyx_v_sst_file, __pyx_v_sst);
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2927 2928 2929 2930 2931 2932 2933 2934 2935 2936 

  /* function exit code */
  goto __pyx_L0;
  __pyx_L1_error:;
  __pyx_r = NULL;
  __pyx_L0:;
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2937 
static PyObject *__pyx_pf_14ancillary_data_py_get_Reynolds_SST(CYTHON_UNUSED PyObject *__pyx_self, PyArrayObject *__pyx_v_lat, PyArrayObject *__pyx_v_lon, PyObject *__pyx_v_res, char *__pyx_v_anc_dir, char *__pyx_v_sst_file, PyObject *__pyx_v_sst) {
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2938 2939 2940 2941 2942 2943 2944 2945 2946 2947 2948 2949 2950 2951 2952 
  __Pyx_memviewslice __pyx_v_sst_mv = { 0, 0, { 0 }, { 0 }, { 0 } };
  __Pyx_LocalBuf_ND __pyx_pybuffernd_lat;
  __Pyx_Buffer __pyx_pybuffer_lat;
  __Pyx_LocalBuf_ND __pyx_pybuffernd_lon;
  __Pyx_Buffer __pyx_pybuffer_lon;
  PyObject *__pyx_r = NULL;
  __Pyx_RefNannyDeclarations
  PyObject *__pyx_t_1 = NULL;
  PyObject *__pyx_t_2 = NULL;
  int __pyx_t_3;
  int __pyx_t_4;
  PyObject *__pyx_t_5 = NULL;
  __Pyx_memviewslice __pyx_t_6 = { 0, 0, { 0 }, { 0 }, { 0 } };
  Py_ssize_t __pyx_t_7;
  Py_ssize_t __pyx_t_8;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2953 2954 
  int __pyx_t_9;
  Py_ssize_t __pyx_t_10;
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2955 2956 2957 2958 2959 2960 2961 2962 2963 2964 2965 2966 2967 2968 2969 
  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;
  int __pyx_clineno = 0;
  __Pyx_RefNannySetupContext("py_get_Reynolds_SST", 0);
  __Pyx_INCREF(__pyx_v_sst);
  __pyx_pybuffer_lat.pybuffer.buf = NULL;
  __pyx_pybuffer_lat.refcount = 0;
  __pyx_pybuffernd_lat.data = NULL;
  __pyx_pybuffernd_lat.rcbuffer = &__pyx_pybuffer_lat;
  __pyx_pybuffer_lon.pybuffer.buf = NULL;
  __pyx_pybuffer_lon.refcount = 0;
  __pyx_pybuffernd_lon.data = NULL;
  __pyx_pybuffernd_lon.rcbuffer = &__pyx_pybuffer_lon;
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2970 
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_lat.rcbuffer->pybuffer, (PyObject*)__pyx_v_lat, &__Pyx_TypeInfo_float, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) __PYX_ERR(0, 28, __pyx_L1_error)
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2971 2972 2973 2974 
  }
  __pyx_pybuffernd_lat.diminfo[0].strides = __pyx_pybuffernd_lat.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_lat.diminfo[0].shape = __pyx_pybuffernd_lat.rcbuffer->pybuffer.shape[0];
  {
    __Pyx_BufFmt_StackElem __pyx_stack[1];
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2975 
    if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_lon.rcbuffer->pybuffer, (PyObject*)__pyx_v_lon, &__Pyx_TypeInfo_float, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) __PYX_ERR(0, 28, __pyx_L1_error)
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2976 2977 2978 
  }
  __pyx_pybuffernd_lon.diminfo[0].strides = __pyx_pybuffernd_lon.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_lon.diminfo[0].shape = __pyx_pybuffernd_lon.rcbuffer->pybuffer.shape[0];
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2979 
  /* "ancillary.pyx":33
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2980 2981 2982 2983 2984 2985 
 * 
 *     # cdef np.ndarray sst = np.zeros((3232*3200, ), order='C', dtype=np.float32)
 *     if not sst.flags['C_CONTIGUOUS']:             # <<<<<<<<<<<<<<
 *         sst = np.ascontiguousarray(sst)
 * 
 */
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2986 
  __pyx_t_1 = __Pyx_PyObject_GetAttrStr(__pyx_v_sst, __pyx_n_s_flags); if (unlikely(!__pyx_t_1)) __PYX_ERR(0, 33, __pyx_L1_error)
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2987 
  __Pyx_GOTREF(__pyx_t_1);
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2988 
  __pyx_t_2 = __Pyx_PyObject_Dict_GetItem(__pyx_t_1, __pyx_n_u_C_CONTIGUOUS); if (unlikely(!__pyx_t_2)) __PYX_ERR(0, 33, __pyx_L1_error)
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2989 2990 
  __Pyx_GOTREF(__pyx_t_2);
  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2991 
  __pyx_t_3 = __Pyx_PyObject_IsTrue(__pyx_t_2); if (unlikely(__pyx_t_3 < 0)) __PYX_ERR(0, 33, __pyx_L1_error)
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2992 2993 2994 2995 
  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
  __pyx_t_4 = ((!__pyx_t_3) != 0);
  if (__pyx_t_4) {
Paolo Veglio's avatar
first setup to write all tests in a clean way  
Paolo Veglio committed
2996 
    /* "ancillary.pyx":34
Paolo Veglio's avatar
temp commit
Paolo Veglio committed
2997 2998 2999 3000 
 *     # cdef np.ndarray sst = np.zeros((3232*3200, ), order='C', dtype=np.float32)
 *     if not sst.flags['C_CONTIGUOUS']:
 *         sst = np.ascontiguousarray(sst)             # <<<<<<<<<<<<<<
 *
Show full blame