�
��]c@`
s%dZd
dl
mZmZmZ
d
dlZd
dlZddlmZ
ej	Z
ejjZ
i
gd6gd6gd6gd	6gd
6gd6gd6gd
6gd6gd6aiai
dd6Zi
dd6Zi
dd
6Zi
dd6Zi
dd	6Zi
dd6Zi
dd6Zi
dd
6Zidd6dd6dd6Zidd6dd6Zd ed!<d"ed#<d$ed%<d&ed'<d'g
td(<d)ed(<d*ed(<d+ed,<d-ed.<d/ed0<d1ed2<d3ed4<d5ed6<d7ed8<d9ed:<d;ed<<d=ed><d?ed@<dAedB<dCedD<dEedF<dGedH<dIedJ<dKedL<dMedN<dOedP<dQedR<dSg
tdT<dUedT<d,g
tdS<dVedS<dWedX<dYedZ<d4g
td[<d\ed[<d]ed^<d6g
td_<d`ed_<d:g
tda<dbeda<dcedd<deedf<d<g
tdg<dhedg<d@g
tdi<djedi<d>g
tdk<dledk<dBg
tdm<dnedm<dBg
tdo<dpedo<dPg
tdq<dredq<dPg
tds<dteds<dudPgtdv<dwedv<dxedy<dzed{<d#d{gtd|<d}ed|<d~ed<d�ed�<d#d{gtdu<d�edu<d#d{gtd�<d�ed�<d�ed�<d�ed�<d�ed�<d�ed�<d�ed�<d�ed�<d�ed�<d�ed�<d�ed�<dudPdBgtd�<d�ed�<dBddugtd�<d�ed�<d�g
td�<d�ed�<d�d4gtd�<d�ed�<d�g
td�<d�ed�<d�ed�<d�ed�<d6g
td�<d�ed�<d�d:gtd�<d�ed�<d�ed�<d�g
td�<d�ed�<d<d:d�gtd�<d�ed�<d@g
td�<d�ed�<d>d�gtd�<d�ed�<dXdqgtd�<d�ed�<dqd.gtd�<d�ed�<dqd4gtd�<d�ed�<dZdqgtd�<d�ed�<d[dqgtd�<d�ed�<d^dqgtd�<d�ed�<d_dqd6gtd�<d�ed�<dfdqgtd�<d�ed�<dddqgtd�<d�ed�<dkdsd>gtd�<d�ed�<didsd@gtd�<d�ed�<dDdPdRgtd�<d�ed�<d��Zdd��
Zd��ZdS(�s�


C declarations, CPP macros, and C functions for f2py2e.
Only required declarations/macros/functions will be used.

Copyright 1999,2000 Pearu Peterson all rights reserved,
Pearu Peterson <pearu@ioc.ee>
Permission to use, modify, and distribute this software is given under the
terms of the NumPy License.

NO WARRANTY IS EXPRESSED OR IMPLIED.  USE AT YOUR OWN RISK.
$Date: 2005/05/06 11:42:34 $
Pearu Peterson

i(tdivisiontabsolute_importtprint_functionNi
(
t__version__t	includes0tincludesttypedefsttypedefs_generatedtuserincludest	cppmacrostcfuncst	callbackstf90modhookstcommonhookss/*need_includes0*/s/*need_includes*/s/*need_userincludes*/s/*need_typedefs*/s/*need_typedefs_generated*/s/*need_cppmacros*/s/*need_cfuncs*/s/*need_callbacks*/s/*need_f90modhooks*/s/*initf90modhooksstatic*/tinitf90modhooksstatics/*initf90modhooksdynamic*/tinitf90modhooksdynamics/*need_commonhooks*/s/*need_initcommonhooks*/tinitcommonhookss#include <math.h>smath.hs#include <string.h>sstring.hs#include <setjmp.h>ssetjmp.hs#include "Python.h"sPython.hs
arrayobject.hsC#define PY_ARRAY_UNIQUE_SYMBOL PyArray_API
#include "arrayobject.h"s#include "fortranobject.h"s#include <stdarg.h>sstdarg.hs$typedef unsigned char unsigned_char;t
unsigned_chars&typedef unsigned short unsigned_short;tunsigned_shorts$typedef unsigned long unsigned_long;t
unsigned_longs typedef signed char signed_char;tsigned_chars�#ifdef _WIN32
typedef __int64 long_long;
#else
typedef long long long_long;
typedef unsigned long long unsigned_long_long;
#endif
t	long_longsf#ifdef _WIN32
typedef __uint64 long_long;
#else
typedef unsigned long long unsigned_long_long;
#endif
tunsigned_long_longs=#ifndef _LONG_DOUBLE
typedef long double long_double;
#endif
tlong_doubles6typedef struct {long double r,i;} complex_long_double;tcomplex_long_doubles*typedef struct {float r,i;} complex_float;t
complex_floats,typedef struct {double r,i;} complex_double;tcomplex_doublestypedef char * string;tstrings
#ifdef DEBUGCFUNCS
#define CFUNCSMESS(mess) fprintf(stderr,"debug-capi:"mess);
#define CFUNCSMESSPY(mess,obj) CFUNCSMESS(mess) \
	PyObject_Print((PyObject *)obj,stderr,Py_PRINT_RAW);\
	fprintf(stderr,"\n");
#else
#define CFUNCSMESS(mess)
#define CFUNCSMESSPY(mess,obj)
#endif
t
CFUNCSMESSso#if defined(PREPEND_FORTRAN)
#if defined(NO_APPEND_FORTRAN)
#if defined(UPPERCASE_FORTRAN)
#define F_FUNC(f,F) _##F
#else
#define F_FUNC(f,F) _##f
#endif
#else
#if defined(UPPERCASE_FORTRAN)
#define F_FUNC(f,F) _##F##_
#else
#define F_FUNC(f,F) _##f##_
#endif
#endif
#else
#if defined(NO_APPEND_FORTRAN)
#if defined(UPPERCASE_FORTRAN)
#define F_FUNC(f,F) F
#else
#define F_FUNC(f,F) f
#endif
#else
#if defined(UPPERCASE_FORTRAN)
#define F_FUNC(f,F) F##_
#else
#define F_FUNC(f,F) f##_
#endif
#endif
#endif
#if defined(UNDERSCORE_G77)
#define F_FUNC_US(f,F) F_FUNC(f##_,F##_)
#else
#define F_FUNC_US(f,F) F_FUNC(f,F)
#endif
tF_FUNCs#if defined(PREPEND_FORTRAN)
#if defined(NO_APPEND_FORTRAN)
#if defined(UPPERCASE_FORTRAN)
#define F_WRAPPEDFUNC(f,F) _F2PYWRAP##F
#else
#define F_WRAPPEDFUNC(f,F) _f2pywrap##f
#endif
#else
#if defined(UPPERCASE_FORTRAN)
#define F_WRAPPEDFUNC(f,F) _F2PYWRAP##F##_
#else
#define F_WRAPPEDFUNC(f,F) _f2pywrap##f##_
#endif
#endif
#else
#if defined(NO_APPEND_FORTRAN)
#if defined(UPPERCASE_FORTRAN)
#define F_WRAPPEDFUNC(f,F) F2PYWRAP##F
#else
#define F_WRAPPEDFUNC(f,F) f2pywrap##f
#endif
#else
#if defined(UPPERCASE_FORTRAN)
#define F_WRAPPEDFUNC(f,F) F2PYWRAP##F##_
#else
#define F_WRAPPEDFUNC(f,F) f2pywrap##f##_
#endif
#endif
#endif
#if defined(UNDERSCORE_G77)
#define F_WRAPPEDFUNC_US(f,F) F_WRAPPEDFUNC(f##_,F##_)
#else
#define F_WRAPPEDFUNC_US(f,F) F_WRAPPEDFUNC(f,F)
#endif
t
F_WRAPPEDFUNCs0#if defined(F90MOD2CCONV1) /*E.g. Compaq Fortran */
#if defined(NO_APPEND_FORTRAN)
#define F_MODFUNCNAME(m,f) $ ## m ## $ ## f
#else
#define F_MODFUNCNAME(m,f) $ ## m ## $ ## f ## _
#endif
#endif

#if defined(F90MOD2CCONV2) /*E.g. IBM XL Fortran, not tested though */
#if defined(NO_APPEND_FORTRAN)
#define F_MODFUNCNAME(m,f)  __ ## m ## _MOD_ ## f
#else
#define F_MODFUNCNAME(m,f)  __ ## m ## _MOD_ ## f ## _
#endif
#endif

#if defined(F90MOD2CCONV3) /*E.g. MIPSPro Compilers */
#if defined(NO_APPEND_FORTRAN)
#define F_MODFUNCNAME(m,f)  f ## .in. ## m
#else
#define F_MODFUNCNAME(m,f)  f ## .in. ## m ## _
#endif
#endif
/*
#if defined(UPPERCASE_FORTRAN)
#define F_MODFUNC(m,M,f,F) F_MODFUNCNAME(M,F)
#else
#define F_MODFUNC(m,M,f,F) F_MODFUNCNAME(m,f)
#endif
*/

#define F_MODFUNC(m,f) (*(f2pymodstruct##m##.##f))
t	F_MODFUNCs�#define SWAP(a,b) (size_t)(a) = ((size_t)(a) ^ (size_t)(b));\
 (size_t)(b) = ((size_t)(a) ^ (size_t)(b));\
 (size_t)(a) = ((size_t)(a) ^ (size_t)(b))
t
SWAPUNSAFEs:#define SWAP(a,b,t) {\
	t *c;\
	c = a;\
	a = b;\
	b = c;}
tSWAPs�#define PRINTPYOBJERR(obj)\
	fprintf(stderr,"#modulename#.error is related to ");\
	PyObject_Print((PyObject *)obj,stderr,Py_PRINT_RAW);\
	fprintf(stderr,"\n");
t
PRINTPYOBJERRs�#ifndef max
#define max(a,b) ((a > b) ? (a) : (b))
#endif
#ifndef min
#define min(a,b) ((a < b) ? (a) : (b))
#endif
#ifndef MAX
#define MAX(a,b) ((a > b) ? (a) : (b))
#endif
#ifndef MIN
#define MIN(a,b) ((a < b) ? (a) : (b))
#endif
tMINMAXt	f2py_sizeslen..sQ#define rank(var) var ## _Rank
#define shape(var,dim) var ## _Dims[dim]
#define old_rank(var) (PyArray_NDIM((PyArrayObject *)(capi_ ## var ## _tmp)))
#define old_shape(var,dim) PyArray_DIM(((PyArrayObject *)(capi_ ## var ## _tmp)),dim)
#define fshape(var,dim) shape(var,rank(var)-dim-1)
#define len(var) shape(var,0)
#define flen(var) fshape(var,0)
#define old_size(var) PyArray_SIZE((PyArrayObject *)(capi_ ## var ## _tmp))
/* #define index(i) capi_i ## i */
#define slen(var) capi_ ## var ## _len
#define size(var, ...) f2py_size((PyArrayObject *)(capi_ ## var ## _tmp), ## __VA_ARGS__, -1)
s�
static int f2py_size(PyArrayObject* var, ...)
{
  npy_int sz = 0;
  npy_int dim;
  npy_int rank;
  va_list argp;
  va_start(argp, var);
  dim = va_arg(argp, npy_int);
  if (dim==-1)
    {
      sz = PyArray_SIZE(var);
    }
  else
    {
      rank = PyArray_NDIM(var);
      if (dim>=1 && dim<=rank)
        sz = PyArray_DIM(var, dim-1);
      else
        fprintf(stderr, "f2py_size: 2nd argument value=%d fails to satisfy 1<=value<=%d. Result will be 0.\n", dim, rank);
    }
  va_end(argp);
  return sz;
}
s/#define pyobj_from_char1(v) (PyInt_FromLong(v))tpyobj_from_char1s0#define pyobj_from_short1(v) (PyInt_FromLong(v))tpyobj_from_short1tpyobj_from_int1s.#define pyobj_from_int1(v) (PyInt_FromLong(v))s0#define pyobj_from_long1(v) (PyLong_FromLong(v))tpyobj_from_long1tpyobj_from_long_long1s�#ifdef HAVE_LONG_LONG
#define pyobj_from_long_long1(v) (PyLong_FromLongLong(v))
#else
#warning HAVE_LONG_LONG is not available. Redefining pyobj_from_long_long.
#define pyobj_from_long_long1(v) (PyLong_FromLong(v))
#endif
tpyobj_from_long_double1s:#define pyobj_from_long_double1(v) (PyFloat_FromDouble(v))s5#define pyobj_from_double1(v) (PyFloat_FromDouble(v))tpyobj_from_double1s4#define pyobj_from_float1(v) (PyFloat_FromDouble(v))tpyobj_from_float1tpyobj_from_complex_long_double1sK#define pyobj_from_complex_long_double1(v) (PyComplex_FromDoubles(v.r,v.i))tpyobj_from_complex_double1sF#define pyobj_from_complex_double1(v) (PyComplex_FromDoubles(v.r,v.i))tpyobj_from_complex_float1sE#define pyobj_from_complex_float1(v) (PyComplex_FromDoubles(v.r,v.i))tpyobj_from_string1s>#define pyobj_from_string1(v) (PyString_FromString((char *)v))tpyobj_from_string1sizesS#define pyobj_from_string1size(v,len) (PyUString_FromStringAndSize((char *)v, len))tTRYPYARRAYTEMPLATEs�/* New SciPy */
#define TRYPYARRAYTEMPLATECHAR case NPY_STRING: *(char *)(PyArray_DATA(arr))=*v; break;
#define TRYPYARRAYTEMPLATELONG case NPY_LONG: *(long *)(PyArray_DATA(arr))=*v; break;
#define TRYPYARRAYTEMPLATEOBJECT case NPY_OBJECT: (PyArray_DESCR(arr)->f->setitem)(pyobj_from_ ## ctype ## 1(*v),PyArray_DATA(arr)); break;

#define TRYPYARRAYTEMPLATE(ctype,typecode) \
        PyArrayObject *arr = NULL;\
        if (!obj) return -2;\
        if (!PyArray_Check(obj)) return -1;\
        if (!(arr=(PyArrayObject *)obj)) {fprintf(stderr,"TRYPYARRAYTEMPLATE:");PRINTPYOBJERR(obj);return 0;}\
        if (PyArray_DESCR(arr)->type==typecode)  {*(ctype *)(PyArray_DATA(arr))=*v; return 1;}\
        switch (PyArray_TYPE(arr)) {\
                case NPY_DOUBLE: *(double *)(PyArray_DATA(arr))=*v; break;\
                case NPY_INT: *(int *)(PyArray_DATA(arr))=*v; break;\
                case NPY_LONG: *(long *)(PyArray_DATA(arr))=*v; break;\
                case NPY_FLOAT: *(float *)(PyArray_DATA(arr))=*v; break;\
                case NPY_CDOUBLE: *(double *)(PyArray_DATA(arr))=*v; break;\
                case NPY_CFLOAT: *(float *)(PyArray_DATA(arr))=*v; break;\
                case NPY_BOOL: *(npy_bool *)(PyArray_DATA(arr))=(*v!=0); break;\
                case NPY_UBYTE: *(unsigned char *)(PyArray_DATA(arr))=*v; break;\
                case NPY_BYTE: *(signed char *)(PyArray_DATA(arr))=*v; break;\
                case NPY_SHORT: *(short *)(PyArray_DATA(arr))=*v; break;\
                case NPY_USHORT: *(npy_ushort *)(PyArray_DATA(arr))=*v; break;\
                case NPY_UINT: *(npy_uint *)(PyArray_DATA(arr))=*v; break;\
                case NPY_ULONG: *(npy_ulong *)(PyArray_DATA(arr))=*v; break;\
                case NPY_LONGLONG: *(npy_longlong *)(PyArray_DATA(arr))=*v; break;\
                case NPY_ULONGLONG: *(npy_ulonglong *)(PyArray_DATA(arr))=*v; break;\
                case NPY_LONGDOUBLE: *(npy_longdouble *)(PyArray_DATA(arr))=*v; break;\
                case NPY_CLONGDOUBLE: *(npy_longdouble *)(PyArray_DATA(arr))=*v; break;\
                case NPY_OBJECT: (PyArray_DESCR(arr)->f->setitem)(pyobj_from_ ## ctype ## 1(*v),PyArray_DATA(arr), arr); break;\
        default: return -2;\
        };\
        return 1
tTRYCOMPLEXPYARRAYTEMPLATEs�	#define TRYCOMPLEXPYARRAYTEMPLATEOBJECT case NPY_OBJECT: (PyArray_DESCR(arr)->f->setitem)(pyobj_from_complex_ ## ctype ## 1((*v)),PyArray_DATA(arr), arr); break;
#define TRYCOMPLEXPYARRAYTEMPLATE(ctype,typecode)\
        PyArrayObject *arr = NULL;\
        if (!obj) return -2;\
        if (!PyArray_Check(obj)) return -1;\
        if (!(arr=(PyArrayObject *)obj)) {fprintf(stderr,"TRYCOMPLEXPYARRAYTEMPLATE:");PRINTPYOBJERR(obj);return 0;}\
        if (PyArray_DESCR(arr)->type==typecode) {\
            *(ctype *)(PyArray_DATA(arr))=(*v).r;\
            *(ctype *)(PyArray_DATA(arr)+sizeof(ctype))=(*v).i;\
            return 1;\
        }\
        switch (PyArray_TYPE(arr)) {\
                case NPY_CDOUBLE: *(double *)(PyArray_DATA(arr))=(*v).r;*(double *)(PyArray_DATA(arr)+sizeof(double))=(*v).i;break;\
                case NPY_CFLOAT: *(float *)(PyArray_DATA(arr))=(*v).r;*(float *)(PyArray_DATA(arr)+sizeof(float))=(*v).i;break;\
                case NPY_DOUBLE: *(double *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_LONG: *(long *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_FLOAT: *(float *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_INT: *(int *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_SHORT: *(short *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_UBYTE: *(unsigned char *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_BYTE: *(signed char *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_BOOL: *(npy_bool *)(PyArray_DATA(arr))=((*v).r!=0 && (*v).i!=0); break;\
                case NPY_USHORT: *(npy_ushort *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_UINT: *(npy_uint *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_ULONG: *(npy_ulong *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_LONGLONG: *(npy_longlong *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_ULONGLONG: *(npy_ulonglong *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_LONGDOUBLE: *(npy_longdouble *)(PyArray_DATA(arr))=(*v).r; break;\
                case NPY_CLONGDOUBLE: *(npy_longdouble *)(PyArray_DATA(arr))=(*v).r;*(npy_longdouble *)(PyArray_DATA(arr)+sizeof(npy_longdouble))=(*v).i;break;\
                case NPY_OBJECT: (PyArray_DESCR(arr)->f->setitem)(pyobj_from_complex_ ## ctype ## 1((*v)),PyArray_DATA(arr), arr); break;\
                default: return -2;\
        };\
        return -1;
tSTRINGCOPYNtGETSTRFROMPYTUPLEs�
#define GETSTRFROMPYTUPLE(tuple,index,str,len) {\
		PyObject *rv_cb_str = PyTuple_GetItem((tuple),(index));\
		if (rv_cb_str == NULL)\
			goto capi_fail;\
		if (PyString_Check(rv_cb_str)) {\
			str[len-1]='\0';\
			STRINGCOPYN((str),PyString_AS_STRING((PyStringObject*)rv_cb_str),(len));\
		} else {\
			PRINTPYOBJERR(rv_cb_str);\
			PyErr_SetString(#modulename#_error,"string object expected");\
			goto capi_fail;\
		}\
	}
s�#define GETSCALARFROMPYTUPLE(tuple,index,var,ctype,mess) {\
		if ((capi_tmp = PyTuple_GetItem((tuple),(index)))==NULL) goto capi_fail;\
		if (!(ctype ## _from_pyobj((var),capi_tmp,mess)))\
			goto capi_fail;\
	}
tGETSCALARFROMPYTUPLEsl
\
#define FAILNULL(p) do {                                            \
    if ((p) == NULL) {                                              \
        PyErr_SetString(PyExc_MemoryError, "NULL pointer found");   \
        goto capi_fail;                                             \
    }                                                               \
} while (0)
tFAILNULLtMEMCOPYsh#define MEMCOPY(to,from,n)\
    do { FAILNULL(to); FAILNULL(from); (void)memcpy(to,from,n); } while (0)
s�#define STRINGMALLOC(str,len)\
	if ((str = (string)malloc(sizeof(char)*(len+1))) == NULL) {\
		PyErr_SetString(PyExc_MemoryError, "out of memory");\
		goto capi_fail;\
	} else {\
		(str)[len] = '\0';\
	}
tSTRINGMALLOCsF#define STRINGFREE(str) do {if (!(str == NULL)) free(str);} while (0)
t
STRINGFREEs&#define STRINGCOPYN(to,from,buf_size)                           \
    do {                                                        \
        int _m = (buf_size);                                    \
        char *_to = (to);                                       \
        char *_from = (from);                                   \
        FAILNULL(_to); FAILNULL(_from);                         \
        (void)strncpy(_to, _from, sizeof(char)*_m);             \
        _to[_m-1] = '\0';                                      \
        /* Padding with spaces instead of nulls */              \
        for (_m -= 2; _m >= 0 && _to[_m] == '\0'; _m--) {      \
            _to[_m] = ' ';                                      \
        }                                                       \
    } while (0)
t
STRINGCOPYsg#define STRINGCOPY(to,from)\
    do { FAILNULL(to); FAILNULL(from); (void)strcpy(to,from); } while (0)
s�#define CHECKGENERIC(check,tcheck,name) \
	if (!(check)) {\
		PyErr_SetString(#modulename#_error,"("tcheck") failed for "name);\
		/*goto capi_fail;*/\
	} else tCHECKGENERICs�#define CHECKARRAY(check,tcheck,name) \
	if (!(check)) {\
		PyErr_SetString(#modulename#_error,"("tcheck") failed for "name);\
		/*goto capi_fail;*/\
	} else t
CHECKARRAYs�#define CHECKSTRING(check,tcheck,name,show,var)\
	if (!(check)) {\
		char errstring[256];\
		sprintf(errstring, "%s: "show, "("tcheck") failed for "name, slen(var), var);\
		PyErr_SetString(#modulename#_error, errstring);\
		/*goto capi_fail;*/\
	} else tCHECKSTRINGs�#define CHECKSCALAR(check,tcheck,name,show,var)\
	if (!(check)) {\
		char errstring[256];\
		sprintf(errstring, "%s: "show, "("tcheck") failed for "name, var);\
		PyErr_SetString(#modulename#_error,errstring);\
		/*goto capi_fail;*/\
	} else tCHECKSCALARs>#define ARRSIZE(dims,rank) (_PyArray_multiply_list(dims,rank))tARRSIZEs�#ifdef OLDPYNUM
#error You need to intall Numeric Python version 13 or higher. Get it from http:/sourceforge.net/project/?group_id=1369
#endif
tOLDPYNUMs�static int calcarrindex(int *i,PyArrayObject *arr) {
	int k,ii = i[0];
	for (k=1; k < PyArray_NDIM(arr); k++)
		ii += (ii*(PyArray_DIM(arr,k) - 1)+i[k]); /* assuming contiguous arr */
	return ii;
}tcalcarrindexs

static int calcarrindextr(int *i,PyArrayObject *arr) {
	int k,ii = i[PyArray_NDIM(arr)-1];
	for (k=1; k < PyArray_NDIM(arr); k++)
		ii += (ii*(PyArray_DIM(arr,PyArray_NDIM(arr)-k-1) - 1)+i[PyArray_NDIM(arr)-k-1]); /* assuming contiguous arr */
	return ii;
}tcalcarrindextrsTstatic struct { int nd;npy_intp *d;int *i,*i_tr,tr; } forcombcache;
static int initforcomb(npy_intp *dims,int nd,int tr) {
  int k;
  if (dims==NULL) return 0;
  if (nd<0) return 0;
  forcombcache.nd = nd;
  forcombcache.d = dims;
  forcombcache.tr = tr;
  if ((forcombcache.i = (int *)malloc(sizeof(int)*nd))==NULL) return 0;
  if ((forcombcache.i_tr = (int *)malloc(sizeof(int)*nd))==NULL) return 0;
  for (k=1;k<nd;k++) {
    forcombcache.i[k] = forcombcache.i_tr[nd-k-1] = 0;
  }
  forcombcache.i[0] = forcombcache.i_tr[nd-1] = -1;
  return 1;
}
static int *nextforcomb(void) {
  int j,*i,*i_tr,k;
  int nd=forcombcache.nd;
  if ((i=forcombcache.i) == NULL) return NULL;
  if ((i_tr=forcombcache.i_tr) == NULL) return NULL;
  if (forcombcache.d == NULL) return NULL;
  i[0]++;
  if (i[0]==forcombcache.d[0]) {
    j=1;
    while ((j<nd) && (i[j]==forcombcache.d[j]-1)) j++;
    if (j==nd) {
      free(i);
      free(i_tr);
      return NULL;
    }
    for (k=0;k<j;k++) i[k] = i_tr[nd-k-1] = 0;
    i[j]++;
    i_tr[nd-j-1]++;
  } else
    i_tr[nd-1]++;
  if (forcombcache.tr) return i_tr;
  return i;
}tforcombttry_pyarr_from_strings`
static int try_pyarr_from_string(PyObject *obj,const string str) {
	PyArrayObject *arr = NULL;
	if (PyArray_Check(obj) && (!((arr = (PyArrayObject *)obj) == NULL)))
		{ STRINGCOPYN(PyArray_DATA(arr),str,PyArray_NBYTES(arr)); }
	return 1;
capi_fail:
	PRINTPYOBJERR(obj);
	PyErr_SetString(#modulename#_error,"try_pyarr_from_string failed");
	return 0;
}
tstring_from_pyobjs&static int string_from_pyobj(string *str,int *len,const string inistr,PyObject *obj,const char *errmess) {
	PyArrayObject *arr = NULL;
	PyObject *tmp = NULL;
#ifdef DEBUGCFUNCS
fprintf(stderr,"string_from_pyobj(str='%s',len=%d,inistr='%s',obj=%p)\n",(char*)str,*len,(char *)inistr,obj);
#endif
	if (obj == Py_None) {
		if (*len == -1)
			*len = strlen(inistr); /* Will this cause problems? */
		STRINGMALLOC(*str,*len);
		STRINGCOPYN(*str,inistr,*len+1);
		return 1;
	}
	if (PyArray_Check(obj)) {
		if ((arr = (PyArrayObject *)obj) == NULL)
			goto capi_fail;
		if (!ISCONTIGUOUS(arr)) {
			PyErr_SetString(PyExc_ValueError,"array object is non-contiguous.");
			goto capi_fail;
		}
		if (*len == -1)
			*len = (PyArray_ITEMSIZE(arr))*PyArray_SIZE(arr);
		STRINGMALLOC(*str,*len);
		STRINGCOPYN(*str,PyArray_DATA(arr),*len+1);
		return 1;
	}
	if (PyString_Check(obj)) {
		tmp = obj;
		Py_INCREF(tmp);
	}
#if PY_VERSION_HEX >= 0x03000000
	else if (PyUnicode_Check(obj)) {
		tmp = PyUnicode_AsASCIIString(obj);
	}
	else {
		PyObject *tmp2;
		tmp2 = PyObject_Str(obj);
		if (tmp2) {
			tmp = PyUnicode_AsASCIIString(tmp2);
			Py_DECREF(tmp2);
		}
		else {
			tmp = NULL;
		}
	}
#else
	else {
		tmp = PyObject_Str(obj);
	}
#endif
	if (tmp == NULL) goto capi_fail;
	if (*len == -1)
		*len = PyString_GET_SIZE(tmp);
	STRINGMALLOC(*str,*len);
	STRINGCOPYN(*str,PyString_AS_STRING(tmp),*len+1);
	Py_DECREF(tmp);
	return 1;
capi_fail:
	Py_XDECREF(tmp);
	{
		PyObject* err = PyErr_Occurred();
		if (err==NULL) err = #modulename#_error;
		PyErr_SetString(err,errmess);
	}
	return 0;
}
tint_from_pyobjtchar_from_pyobjs�static int char_from_pyobj(char* v,PyObject *obj,const char *errmess) {
	int i=0;
	if (int_from_pyobj(&i,obj,errmess)) {
		*v = (char)i;
		return 1;
	}
	return 0;
}
tsigned_char_from_pyobjs�static int signed_char_from_pyobj(signed_char* v,PyObject *obj,const char *errmess) {
	int i=0;
	if (int_from_pyobj(&i,obj,errmess)) {
		*v = (signed_char)i;
		return 1;
	}
	return 0;
}
tshort_from_pyobjs�static int short_from_pyobj(short* v,PyObject *obj,const char *errmess) {
	int i=0;
	if (int_from_pyobj(&i,obj,errmess)) {
		*v = (short)i;
		return 1;
	}
	return 0;
}
s�static int int_from_pyobj(int* v,PyObject *obj,const char *errmess) {
	PyObject* tmp = NULL;
	if (PyInt_Check(obj)) {
		*v = (int)PyInt_AS_LONG(obj);
		return 1;
	}
	tmp = PyNumber_Int(obj);
	if (tmp) {
		*v = PyInt_AS_LONG(tmp);
		Py_DECREF(tmp);
		return 1;
	}
	if (PyComplex_Check(obj))
		tmp = PyObject_GetAttrString(obj,"real");
	else if (PyString_Check(obj) || PyUnicode_Check(obj))
		/*pass*/;
	else if (PySequence_Check(obj))
		tmp = PySequence_GetItem(obj,0);
	if (tmp) {
		PyErr_Clear();
		if (int_from_pyobj(v,tmp,errmess)) {Py_DECREF(tmp); return 1;}
		Py_DECREF(tmp);
	}
	{
		PyObject* err = PyErr_Occurred();
		if (err==NULL) err = #modulename#_error;
		PyErr_SetString(err,errmess);
	}
	return 0;
}
s�static int long_from_pyobj(long* v,PyObject *obj,const char *errmess) {
	PyObject* tmp = NULL;
	if (PyInt_Check(obj)) {
		*v = PyInt_AS_LONG(obj);
		return 1;
	}
	tmp = PyNumber_Int(obj);
	if (tmp) {
		*v = PyInt_AS_LONG(tmp);
		Py_DECREF(tmp);
		return 1;
	}
	if (PyComplex_Check(obj))
		tmp = PyObject_GetAttrString(obj,"real");
	else if (PyString_Check(obj) || PyUnicode_Check(obj))
		/*pass*/;
	else if (PySequence_Check(obj))
		tmp = PySequence_GetItem(obj,0);
	if (tmp) {
		PyErr_Clear();
		if (long_from_pyobj(v,tmp,errmess)) {Py_DECREF(tmp); return 1;}
		Py_DECREF(tmp);
	}
	{
		PyObject* err = PyErr_Occurred();
		if (err==NULL) err = #modulename#_error;
		PyErr_SetString(err,errmess);
	}
	return 0;
}
tlong_from_pyobjtlong_long_from_pyobjsSstatic int long_long_from_pyobj(long_long* v,PyObject *obj,const char *errmess) {
	PyObject* tmp = NULL;
	if (PyLong_Check(obj)) {
		*v = PyLong_AsLongLong(obj);
		return (!PyErr_Occurred());
	}
	if (PyInt_Check(obj)) {
		*v = (long_long)PyInt_AS_LONG(obj);
		return 1;
	}
	tmp = PyNumber_Long(obj);
	if (tmp) {
		*v = PyLong_AsLongLong(tmp);
		Py_DECREF(tmp);
		return (!PyErr_Occurred());
	}
	if (PyComplex_Check(obj))
		tmp = PyObject_GetAttrString(obj,"real");
	else if (PyString_Check(obj) || PyUnicode_Check(obj))
		/*pass*/;
	else if (PySequence_Check(obj))
		tmp = PySequence_GetItem(obj,0);
	if (tmp) {
		PyErr_Clear();
		if (long_long_from_pyobj(v,tmp,errmess)) {Py_DECREF(tmp); return 1;}
		Py_DECREF(tmp);
	}
	{
		PyObject* err = PyErr_Occurred();
		if (err==NULL) err = #modulename#_error;
		PyErr_SetString(err,errmess);
	}
	return 0;
}
tdouble_from_pyobjtlong_double_from_pyobjs�
static int long_double_from_pyobj(long_double* v,PyObject *obj,const char *errmess) {
	double d=0;
	if (PyArray_CheckScalar(obj)){
		if PyArray_IsScalar(obj, LongDouble) {
			PyArray_ScalarAsCtype(obj, v);
			return 1;
		}
		else if (PyArray_Check(obj) && PyArray_TYPE(obj)==NPY_LONGDOUBLE) {
			(*v) = *((npy_longdouble *)PyArray_DATA(obj));
			return 1;
		}
	}
	if (double_from_pyobj(&d,obj,errmess)) {
		*v = (long_double)d;
		return 1;
	}
	return 0;
}
sJstatic int double_from_pyobj(double* v,PyObject *obj,const char *errmess) {
	PyObject* tmp = NULL;
	if (PyFloat_Check(obj)) {
#ifdef __sgi
		*v = PyFloat_AsDouble(obj);
#else
		*v = PyFloat_AS_DOUBLE(obj);
#endif
		return 1;
	}
	tmp = PyNumber_Float(obj);
	if (tmp) {
#ifdef __sgi
		*v = PyFloat_AsDouble(tmp);
#else
		*v = PyFloat_AS_DOUBLE(tmp);
#endif
		Py_DECREF(tmp);
		return 1;
	}
	if (PyComplex_Check(obj))
		tmp = PyObject_GetAttrString(obj,"real");
	else if (PyString_Check(obj) || PyUnicode_Check(obj))
		/*pass*/;
	else if (PySequence_Check(obj))
		tmp = PySequence_GetItem(obj,0);
	if (tmp) {
		PyErr_Clear();
		if (double_from_pyobj(v,tmp,errmess)) {Py_DECREF(tmp); return 1;}
		Py_DECREF(tmp);
	}
	{
		PyObject* err = PyErr_Occurred();
		if (err==NULL) err = #modulename#_error;
		PyErr_SetString(err,errmess);
	}
	return 0;
}
tfloat_from_pyobjs�static int float_from_pyobj(float* v,PyObject *obj,const char *errmess) {
	double d=0.0;
	if (double_from_pyobj(&d,obj,errmess)) {
		*v = (float)d;
		return 1;
	}
	return 0;
}
tcomplex_double_from_pyobjtcomplex_long_double_from_pyobjs[static int complex_long_double_from_pyobj(complex_long_double* v,PyObject *obj,const char *errmess) {
	complex_double cd={0.0,0.0};
	if (PyArray_CheckScalar(obj)){
		if PyArray_IsScalar(obj, CLongDouble) {
			PyArray_ScalarAsCtype(obj, v);
			return 1;
		}
		else if (PyArray_Check(obj) && PyArray_TYPE(obj)==NPY_CLONGDOUBLE) {
			(*v).r = ((npy_clongdouble *)PyArray_DATA(obj))->real;
			(*v).i = ((npy_clongdouble *)PyArray_DATA(obj))->imag;
			return 1;
		}
	}
	if (complex_double_from_pyobj(&cd,obj,errmess)) {
		(*v).r = (long_double)cd.r;
		(*v).i = (long_double)cd.i;
		return 1;
	}
	return 0;
}
sSstatic int complex_double_from_pyobj(complex_double* v,PyObject *obj,const char *errmess) {
	Py_complex c;
	if (PyComplex_Check(obj)) {
		c=PyComplex_AsCComplex(obj);
		(*v).r=c.real, (*v).i=c.imag;
		return 1;
	}
	if (PyArray_IsScalar(obj, ComplexFloating)) {
		if (PyArray_IsScalar(obj, CFloat)) {
			npy_cfloat new;
			PyArray_ScalarAsCtype(obj, &new);
			(*v).r = (double)new.real;
			(*v).i = (double)new.imag;
		}
		else if (PyArray_IsScalar(obj, CLongDouble)) {
			npy_clongdouble new;
			PyArray_ScalarAsCtype(obj, &new);
			(*v).r = (double)new.real;
			(*v).i = (double)new.imag;
		}
		else { /* if (PyArray_IsScalar(obj, CDouble)) */
			PyArray_ScalarAsCtype(obj, v);
		}
		return 1;
	}
	if (PyArray_CheckScalar(obj)) { /* 0-dim array or still array scalar */
		PyObject *arr;
		if (PyArray_Check(obj)) {
			arr = PyArray_Cast((PyArrayObject *)obj, NPY_CDOUBLE);
		}
		else {
			arr = PyArray_FromScalar(obj, PyArray_DescrFromType(NPY_CDOUBLE));
		}
		if (arr==NULL) return 0;
		(*v).r = ((npy_cdouble *)PyArray_DATA(arr))->real;
		(*v).i = ((npy_cdouble *)PyArray_DATA(arr))->imag;
		return 1;
	}
	/* Python does not provide PyNumber_Complex function :-( */
	(*v).i=0.0;
	if (PyFloat_Check(obj)) {
#ifdef __sgi
		(*v).r = PyFloat_AsDouble(obj);
#else
		(*v).r = PyFloat_AS_DOUBLE(obj);
#endif
		return 1;
	}
	if (PyInt_Check(obj)) {
		(*v).r = (double)PyInt_AS_LONG(obj);
		return 1;
	}
	if (PyLong_Check(obj)) {
		(*v).r = PyLong_AsDouble(obj);
		return (!PyErr_Occurred());
	}
	if (PySequence_Check(obj) && !(PyString_Check(obj) || PyUnicode_Check(obj))) {
		PyObject *tmp = PySequence_GetItem(obj,0);
		if (tmp) {
			if (complex_double_from_pyobj(v,tmp,errmess)) {
				Py_DECREF(tmp);
				return 1;
			}
			Py_DECREF(tmp);
		}
	}
	{
		PyObject* err = PyErr_Occurred();
		if (err==NULL)
			err = PyExc_TypeError;
		PyErr_SetString(err,errmess);
	}
	return 0;
}
tcomplex_float_from_pyobjs�static int complex_float_from_pyobj(complex_float* v,PyObject *obj,const char *errmess) {
	complex_double cd={0.0,0.0};
	if (complex_double_from_pyobj(&cd,obj,errmess)) {
		(*v).r = (float)cd.r;
		(*v).i = (float)cd.i;
		return 1;
	}
	return 0;
}
ttry_pyarr_from_charsYstatic int try_pyarr_from_char(PyObject* obj,char* v) {
	TRYPYARRAYTEMPLATE(char,'c');
}
ttry_pyarr_from_signed_charststatic int try_pyarr_from_unsigned_char(PyObject* obj,unsigned_char* v) {
	TRYPYARRAYTEMPLATE(unsigned_char,'b');
}
ttry_pyarr_from_unsigned_charsnstatic int try_pyarr_from_signed_char(PyObject* obj,signed_char* v) {
	TRYPYARRAYTEMPLATE(signed_char,'1');
}
ttry_pyarr_from_shorts\static int try_pyarr_from_short(PyObject* obj,short* v) {
	TRYPYARRAYTEMPLATE(short,'s');
}
ttry_pyarr_from_intsVstatic int try_pyarr_from_int(PyObject* obj,int* v) {
	TRYPYARRAYTEMPLATE(int,'i');
}
ttry_pyarr_from_longsYstatic int try_pyarr_from_long(PyObject* obj,long* v) {
	TRYPYARRAYTEMPLATE(long,'l');
}
ttry_pyarr_from_long_longshstatic int try_pyarr_from_long_long(PyObject* obj,long_long* v) {
	TRYPYARRAYTEMPLATE(long_long,'L');
}
ttry_pyarr_from_floats\static int try_pyarr_from_float(PyObject* obj,float* v) {
	TRYPYARRAYTEMPLATE(float,'f');
}
ttry_pyarr_from_doubles_static int try_pyarr_from_double(PyObject* obj,double* v) {
	TRYPYARRAYTEMPLATE(double,'d');
}
ttry_pyarr_from_complex_floatssstatic int try_pyarr_from_complex_float(PyObject* obj,complex_float* v) {
	TRYCOMPLEXPYARRAYTEMPLATE(float,'F');
}
ttry_pyarr_from_complex_doublesvstatic int try_pyarr_from_complex_double(PyObject* obj,complex_double* v) {
	TRYCOMPLEXPYARRAYTEMPLATE(double,'D');
}
tcreate_cb_arglistst
static int create_cb_arglist(PyObject* fun,PyTupleObject* xa,const int maxnofargs,const int nofoptargs,int *nofargs,PyTupleObject **args,const char *errmess) {
	PyObject *tmp = NULL;
	PyObject *tmp_fun = NULL;
	int tot,opt,ext,siz,i,di=0;
	CFUNCSMESS("create_cb_arglist\n");
	tot=opt=ext=siz=0;
	/* Get the total number of arguments */
	if (PyFunction_Check(fun))
		tmp_fun = fun;
	else {
		di = 1;
		if (PyObject_HasAttrString(fun,"im_func")) {
			tmp_fun = PyObject_GetAttrString(fun,"im_func");
		}
		else if (PyObject_HasAttrString(fun,"__call__")) {
			tmp = PyObject_GetAttrString(fun,"__call__");
			if (PyObject_HasAttrString(tmp,"im_func"))
				tmp_fun = PyObject_GetAttrString(tmp,"im_func");
			else {
				tmp_fun = fun; /* built-in function */
				tot = maxnofargs;
				if (xa != NULL)
					tot += PyTuple_Size((PyObject *)xa);
			}
			Py_XDECREF(tmp);
		}
		else if (PyFortran_Check(fun) || PyFortran_Check1(fun)) {
			tot = maxnofargs;
			if (xa != NULL)
				tot += PyTuple_Size((PyObject *)xa);
			tmp_fun = fun;
		}
		else if (F2PyCapsule_Check(fun)) {
			tot = maxnofargs;
			if (xa != NULL)
				ext = PyTuple_Size((PyObject *)xa);
			if(ext>0) {
				fprintf(stderr,"extra arguments tuple cannot be used with CObject call-back\n");
				goto capi_fail;
			}
			tmp_fun = fun;
		}
	}
if (tmp_fun==NULL) {
fprintf(stderr,"Call-back argument must be function|instance|instance.__call__|f2py-function but got %s.\n",(fun==NULL?"NULL":Py_TYPE(fun)->tp_name));
goto capi_fail;
}
#if PY_VERSION_HEX >= 0x03000000
	if (PyObject_HasAttrString(tmp_fun,"__code__")) {
		if (PyObject_HasAttrString(tmp = PyObject_GetAttrString(tmp_fun,"__code__"),"co_argcount"))
#else
	if (PyObject_HasAttrString(tmp_fun,"func_code")) {
		if (PyObject_HasAttrString(tmp = PyObject_GetAttrString(tmp_fun,"func_code"),"co_argcount"))
#endif
			tot = PyInt_AsLong(PyObject_GetAttrString(tmp,"co_argcount")) - di;
		Py_XDECREF(tmp);
	}
	/* Get the number of optional arguments */
#if PY_VERSION_HEX >= 0x03000000
	if (PyObject_HasAttrString(tmp_fun,"__defaults__")) {
		if (PyTuple_Check(tmp = PyObject_GetAttrString(tmp_fun,"__defaults__")))
#else
	if (PyObject_HasAttrString(tmp_fun,"func_defaults")) {
		if (PyTuple_Check(tmp = PyObject_GetAttrString(tmp_fun,"func_defaults")))
#endif
			opt = PyTuple_Size(tmp);
		Py_XDECREF(tmp);
	}
	/* Get the number of extra arguments */
	if (xa != NULL)
		ext = PyTuple_Size((PyObject *)xa);
	/* Calculate the size of call-backs argument list */
	siz = MIN(maxnofargs+ext,tot);
	*nofargs = MAX(0,siz-ext);
#ifdef DEBUGCFUNCS
	fprintf(stderr,"debug-capi:create_cb_arglist:maxnofargs(-nofoptargs),tot,opt,ext,siz,nofargs=%d(-%d),%d,%d,%d,%d,%d\n",maxnofargs,nofoptargs,tot,opt,ext,siz,*nofargs);
#endif
	if (siz<tot-opt) {
		fprintf(stderr,"create_cb_arglist: Failed to build argument list (siz) with enough arguments (tot-opt) required by user-supplied function (siz,tot,opt=%d,%d,%d).\n",siz,tot,opt);
		goto capi_fail;
	}
	/* Initialize argument list */
	*args = (PyTupleObject *)PyTuple_New(siz);
	for (i=0;i<*nofargs;i++) {
		Py_INCREF(Py_None);
		PyTuple_SET_ITEM((PyObject *)(*args),i,Py_None);
	}
	if (xa != NULL)
		for (i=(*nofargs);i<siz;i++) {
			tmp = PyTuple_GetItem((PyObject *)xa,i-(*nofargs));
			Py_INCREF(tmp);
			PyTuple_SET_ITEM(*args,i,tmp);
		}
	CFUNCSMESS("create_cb_arglist-end\n");
	return 1;
capi_fail:
	if ((PyErr_Occurred())==NULL)
		PyErr_SetString(#modulename#_error,errmess);
	return 0;
}
cC`
skd
dlm
}
x6|j�D](}
d|
}d|||
ft|<qWd}
d|
}d|t|<dS(Ni
(
t
c2capi_mapspyarr_from_p_%s1s>#define %s(v) (PyArray_SimpleNewFromData(0,NULL,%s,(char *)v))Rsh#define %s(v,dims) (PyArray_New(&PyArray_Type, 1, dims, NPY_STRING, NULL, v, 1, NPY_ARRAY_CARRAY, NULL))(t	capi_mapsR_tkeysR	(R_t
kt
m((sD/opt/alt/python27/lib64/python2.7/site-packages/numpy/f2py/cfuncs.pytbuildcfuncs{s






cC`
st|t
�r0x�|D]}t||
�
qWn�t|t�r�|sIdS|tkr^d
}n�|tkrsd}n�|tkr�d}n�|tkr�d}n�|tkr�d}nl|t	kr�d}nW|t
kr�d}nB|tkr�d}n-|tkr
d	}nt
d
t|�
�

dS|t|kr2
dS|
r5i}|tkr�
x�t|D]q}t|d�}t|t�rU
xJ|j�D]9}||kr�
||||||<q�
||||<q�
WqU
qU
WnxQ|j�D]C}x:||D].}|t|kr�
|g
t|t|<q�
q�
Wq�
Wt|j|�

qi}|tkr�x�t|D]q}t||
�}t|t�rRxJ|j�D]9}||kr�||||||<q�||||<q�WqRqRWn||kr�g||<n||j|�

|Snt
dt|�
�

dS(
NRRRRR	R
RRR
sappend_needs: unknown need %s
is2append_needs: expected list or string but got :%s
(t
isinstancetlisttappend_needststrRRRRR	R
RRR
terrmesstreprtoutneedstneedstdictRatappend(tneedtflagt
nttmptnnt
ttnnn((sD/opt/alt/python27/lib64/python2.7/site-packages/numpy/f2py/cfuncs.pyRg�st







	

	

	

	

	

	

	

	

	















 













cC`
s�
i}x�
tj
�D]�
}
g}tjt|
�
}xI
tt|
�
d
kr}
t|
d
tkr�|jt|
d
�

t|
d
=n�d
}x8t|
dD](}|tt|
d
kr�d}Pq�q�W|r�t|
dt|
d
g
t|
<n |jt|
d
�

t|
d
=|rg
d
td�|t|
�krg
t|
gkrg
t|
|�
td�

||}Pntjt|
�
}q5W|gkr�
|
g
}n|||
<qW|S(Nii
cS`
s
||
kS(
N((t
xt
y((sD/opt/alt/python27/lib64/python2.7/site-packages/numpy/f2py/cfuncs.pyt<lambda>�tsQget_needs: no progress in sorting needs, probably circular dependence, skipping.
(	RkRatcopytlenRlRntmaptprintRi(tresRqtouttsaveoutRpRb((sD/opt/alt/python27/lib64/python2.7/site-packages/numpy/f2py/cfuncs.pyt	get_needs�s:












$

%












(t__doc__t
__future__RR
RtsysRzRyRtversiontf2py_versiontstderrtwriteRiRkRlRRRRRR	R
RRR
RdRgR�(((sD/opt/alt/python27/lib64/python2.7/site-packages/numpy/f2py/cfuncs.pyt<module>s4


	



























































$
$
"










































#

#
"

	


	












	

(



D









 
 


$


(









L




















f
	C