Centre Blaise Pascal » Bench4GPU

#!/usr/bin/env python

#

# Ising2D model in serial mode

#

# CC BY-NC-SA 2011 : <emmanuel.quemener@ens-lyon.fr> 

import sys

import numpy

from PIL import Image

from math import exp

from random import random

import time

import getopt

import matplotlib.pyplot as plt

import Metropolis

def Metropolis(sigma,J,B,T,iterations): 

    start=time.time()

    SizeX,SizeY=sigma.shape

    for p in xrange(0,iterations):

        # Random access coordonate

        X,Y=numpy.random.randint(SizeX),numpy.random.randint(SizeY)

        DeltaE=J*sigma[X,Y]*(2*(sigma[X,(Y+1)%SizeY]+

                                sigma[X,(Y-1)%SizeY]+

                                sigma[(X-1)%SizeX,Y]+

                                sigma[(X+1)%SizeX,Y])+B)

        if DeltaE < 0. or random() < exp(-DeltaE/T):

            sigma[X,Y]=-sigma[X,Y]

    duration=time.time()-start

    return(duration)

/* Generated by Cython 0.15.1 on Sun May 26 16:17:32 2013 */

#define PY_SSIZE_T_CLEAN

#include "Python.h"

#ifndef Py_PYTHON_H

    #error Python headers needed to compile C extensions, please install development version of Python.

#else

#include <stddef.h> /* For offsetof */

#ifndef offsetof

#define offsetof(type, member) ( (size_t) & ((type*)0) -> member )

#endif

#if !defined(WIN32) && !defined(MS_WINDOWS)

  #ifndef __stdcall

    #define __stdcall

  #endif

  #ifndef __cdecl

    #define __cdecl

  #endif

  #ifndef __fastcall

    #define __fastcall

  #endif

#endif

#ifndef DL_IMPORT

  #define DL_IMPORT(t) t

#endif

#ifndef DL_EXPORT

  #define DL_EXPORT(t) t

#endif

#ifndef PY_LONG_LONG

  #define PY_LONG_LONG LONG_LONG

#endif

#if PY_VERSION_HEX < 0x02040000

  #define METH_COEXIST 0

  #define PyDict_CheckExact(op) (Py_TYPE(op) == &PyDict_Type)

  #define PyDict_Contains(d,o)   PySequence_Contains(d,o)

#endif

#if PY_VERSION_HEX < 0x02050000

  typedef int Py_ssize_t;

  #define PY_SSIZE_T_MAX INT_MAX

  #define PY_SSIZE_T_MIN INT_MIN

  #define PY_FORMAT_SIZE_T ""

  #define PyInt_FromSsize_t(z) PyInt_FromLong(z)

  #define PyInt_AsSsize_t(o)   __Pyx_PyInt_AsInt(o)

  #define PyNumber_Index(o)    PyNumber_Int(o)

  #define PyIndex_Check(o)     PyNumber_Check(o)

  #define PyErr_WarnEx(category, message, stacklevel) PyErr_Warn(category, message)

#endif

#if PY_VERSION_HEX < 0x02060000

  #define Py_REFCNT(ob) (((PyObject*)(ob))->ob_refcnt)

  #define Py_TYPE(ob)   (((PyObject*)(ob))->ob_type)

  #define Py_SIZE(ob)   (((PyVarObject*)(ob))->ob_size)

  #define PyVarObject_HEAD_INIT(type, size) \

          PyObject_HEAD_INIT(type) size,

  #define PyType_Modified(t)

  typedef struct {

     void *buf;

     PyObject *obj;

     Py_ssize_t len;

     Py_ssize_t itemsize;

     int readonly;

     int ndim;

     char *format;

     Py_ssize_t *shape;

     Py_ssize_t *strides;

     Py_ssize_t *suboffsets;

     void *internal;

  } Py_buffer;

  #define PyBUF_SIMPLE 0

  #define PyBUF_WRITABLE 0x0001

  #define PyBUF_FORMAT 0x0004

  #define PyBUF_ND 0x0008

  #define PyBUF_STRIDES (0x0010 | PyBUF_ND)

  #define PyBUF_C_CONTIGUOUS (0x0020 | PyBUF_STRIDES)

  #define PyBUF_F_CONTIGUOUS (0x0040 | PyBUF_STRIDES)

  #define PyBUF_ANY_CONTIGUOUS (0x0080 | PyBUF_STRIDES)

  #define PyBUF_INDIRECT (0x0100 | PyBUF_STRIDES)

#endif

#if PY_MAJOR_VERSION < 3

  #define __Pyx_BUILTIN_MODULE_NAME "__builtin__"

#else

  #define __Pyx_BUILTIN_MODULE_NAME "builtins"

#endif

#if PY_MAJOR_VERSION >= 3

  #define Py_TPFLAGS_CHECKTYPES 0

  #define Py_TPFLAGS_HAVE_INDEX 0

#endif

#if (PY_VERSION_HEX < 0x02060000) || (PY_MAJOR_VERSION >= 3)

  #define Py_TPFLAGS_HAVE_NEWBUFFER 0

#endif

#if PY_MAJOR_VERSION >= 3

  #define PyBaseString_Type            PyUnicode_Type

  #define PyStringObject               PyUnicodeObject

  #define PyString_Type                PyUnicode_Type

  #define PyString_Check               PyUnicode_Check

  #define PyString_CheckExact          PyUnicode_CheckExact

#endif

#if PY_VERSION_HEX < 0x02060000

  #define PyBytesObject                PyStringObject

  #define PyBytes_Type                 PyString_Type

  #define PyBytes_Check                PyString_Check

  #define PyBytes_CheckExact           PyString_CheckExact

  #define PyBytes_FromString           PyString_FromString

  #define PyBytes_FromStringAndSize    PyString_FromStringAndSize

  #define PyBytes_FromFormat           PyString_FromFormat

  #define PyBytes_DecodeEscape         PyString_DecodeEscape

  #define PyBytes_AsString             PyString_AsString

  #define PyBytes_AsStringAndSize      PyString_AsStringAndSize

  #define PyBytes_Size                 PyString_Size

  #define PyBytes_AS_STRING            PyString_AS_STRING

  #define PyBytes_GET_SIZE             PyString_GET_SIZE

  #define PyBytes_Repr                 PyString_Repr

  #define PyBytes_Concat               PyString_Concat

  #define PyBytes_ConcatAndDel         PyString_ConcatAndDel

#endif

#if PY_VERSION_HEX < 0x02060000

  #define PySet_Check(obj)             PyObject_TypeCheck(obj, &PySet_Type)

  #define PyFrozenSet_Check(obj)       PyObject_TypeCheck(obj, &PyFrozenSet_Type)

#endif

#ifndef PySet_CheckExact

  #define PySet_CheckExact(obj)        (Py_TYPE(obj) == &PySet_Type)

#endif

#define __Pyx_TypeCheck(obj, type) PyObject_TypeCheck(obj, (PyTypeObject *)type)

#if PY_MAJOR_VERSION >= 3

  #define PyIntObject                  PyLongObject

  #define PyInt_Type                   PyLong_Type

  #define PyInt_Check(op)              PyLong_Check(op)

  #define PyInt_CheckExact(op)         PyLong_CheckExact(op)

  #define PyInt_FromString             PyLong_FromString

  #define PyInt_FromUnicode            PyLong_FromUnicode

  #define PyInt_FromLong               PyLong_FromLong

  #define PyInt_FromSize_t             PyLong_FromSize_t

  #define PyInt_FromSsize_t            PyLong_FromSsize_t

  #define PyInt_AsLong                 PyLong_AsLong

  #define PyInt_AS_LONG                PyLong_AS_LONG

  #define PyInt_AsSsize_t              PyLong_AsSsize_t

  #define PyInt_AsUnsignedLongMask     PyLong_AsUnsignedLongMask

  #define PyInt_AsUnsignedLongLongMask PyLong_AsUnsignedLongLongMask

#endif

#if PY_MAJOR_VERSION >= 3

  #define PyBoolObject                 PyLongObject

#endif

#if PY_VERSION_HEX < 0x03020000

  typedef long Py_hash_t;

  #define __Pyx_PyInt_FromHash_t PyInt_FromLong

  #define __Pyx_PyInt_AsHash_t   PyInt_AsLong

#else

  #define __Pyx_PyInt_FromHash_t PyInt_FromSsize_t

  #define __Pyx_PyInt_AsHash_t   PyInt_AsSsize_t

#endif

#if PY_MAJOR_VERSION >= 3

  #define __Pyx_PyNumber_Divide(x,y)         PyNumber_TrueDivide(x,y)

  #define __Pyx_PyNumber_InPlaceDivide(x,y)  PyNumber_InPlaceTrueDivide(x,y)

#else

  #define __Pyx_PyNumber_Divide(x,y)         PyNumber_Divide(x,y)

  #define __Pyx_PyNumber_InPlaceDivide(x,y)  PyNumber_InPlaceDivide(x,y)

#endif

#if (PY_MAJOR_VERSION < 3) || (PY_VERSION_HEX >= 0x03010300)

  #define __Pyx_PySequence_GetSlice(obj, a, b) PySequence_GetSlice(obj, a, b)

  #define __Pyx_PySequence_SetSlice(obj, a, b, value) PySequence_SetSlice(obj, a, b, value)

  #define __Pyx_PySequence_DelSlice(obj, a, b) PySequence_DelSlice(obj, a, b)

#else

  #define __Pyx_PySequence_GetSlice(obj, a, b) (unlikely(!(obj)) ? \

        (PyErr_SetString(PyExc_SystemError, "null argument to internal routine"), (PyObject*)0) : \

        (likely((obj)->ob_type->tp_as_mapping) ? (PySequence_GetSlice(obj, a, b)) : \

            (PyErr_Format(PyExc_TypeError, "'%.200s' object is unsliceable", (obj)->ob_type->tp_name), (PyObject*)0)))

  #define __Pyx_PySequence_SetSlice(obj, a, b, value) (unlikely(!(obj)) ? \

        (PyErr_SetString(PyExc_SystemError, "null argument to internal routine"), -1) : \

        (likely((obj)->ob_type->tp_as_mapping) ? (PySequence_SetSlice(obj, a, b, value)) : \

            (PyErr_Format(PyExc_TypeError, "'%.200s' object doesn't support slice assignment", (obj)->ob_type->tp_name), -1)))

  #define __Pyx_PySequence_DelSlice(obj, a, b) (unlikely(!(obj)) ? \

        (PyErr_SetString(PyExc_SystemError, "null argument to internal routine"), -1) : \

        (likely((obj)->ob_type->tp_as_mapping) ? (PySequence_DelSlice(obj, a, b)) : \

            (PyErr_Format(PyExc_TypeError, "'%.200s' object doesn't support slice deletion", (obj)->ob_type->tp_name), -1)))

#endif

#if PY_MAJOR_VERSION >= 3

  #define PyMethod_New(func, self, klass) ((self) ? PyMethod_New(func, self) : PyInstanceMethod_New(func))

#endif

#if PY_VERSION_HEX < 0x02050000

  #define __Pyx_GetAttrString(o,n)   PyObject_GetAttrString((o),((char *)(n)))

  #define __Pyx_SetAttrString(o,n,a) PyObject_SetAttrString((o),((char *)(n)),(a))

  #define __Pyx_DelAttrString(o,n)   PyObject_DelAttrString((o),((char *)(n)))

#else

  #define __Pyx_GetAttrString(o,n)   PyObject_GetAttrString((o),(n))

  #define __Pyx_SetAttrString(o,n,a) PyObject_SetAttrString((o),(n),(a))

  #define __Pyx_DelAttrString(o,n)   PyObject_DelAttrString((o),(n))

#endif

#if PY_VERSION_HEX < 0x02050000

  #define __Pyx_NAMESTR(n) ((char *)(n))

  #define __Pyx_DOCSTR(n)  ((char *)(n))

#else

  #define __Pyx_NAMESTR(n) (n)

  #define __Pyx_DOCSTR(n)  (n)

#endif

#ifndef __PYX_EXTERN_C

  #ifdef __cplusplus

    #define __PYX_EXTERN_C extern "C"

  #else

    #define __PYX_EXTERN_C extern

  #endif

#endif

#if defined(WIN32) || defined(MS_WINDOWS)

#define _USE_MATH_DEFINES

#endif

#include <math.h>

#define __PYX_HAVE__Metropolis

#define __PYX_HAVE_API__Metropolis

#ifdef _OPENMP

#include <omp.h>

#endif /* _OPENMP */

#ifdef PYREX_WITHOUT_ASSERTIONS

#define CYTHON_WITHOUT_ASSERTIONS

#endif

/* inline attribute */

#ifndef CYTHON_INLINE

  #if defined(__GNUC__)

    #define CYTHON_INLINE __inline__

  #elif defined(_MSC_VER)

    #define CYTHON_INLINE __inline

  #elif defined (__STDC_VERSION__) && __STDC_VERSION__ >= 199901L

    #define CYTHON_INLINE inline

  #else

    #define CYTHON_INLINE

  #endif

#endif

/* unused attribute */

#ifndef CYTHON_UNUSED

# if defined(__GNUC__)

#   if !(defined(__cplusplus)) || (__GNUC__ > 3 || (__GNUC__ == 3 && __GNUC_MINOR__ >= 4))

#     define CYTHON_UNUSED __attribute__ ((__unused__))

#   else

#     define CYTHON_UNUSED

#   endif

# elif defined(__ICC) || defined(__INTEL_COMPILER)

#   define CYTHON_UNUSED __attribute__ ((__unused__))

# else

#   define CYTHON_UNUSED

# endif

#endif

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

/* Type Conversion Predeclarations */

#define __Pyx_PyBytes_FromUString(s) PyBytes_FromString((char*)s)

#define __Pyx_PyBytes_AsUString(s)   ((unsigned char*) PyBytes_AsString(s))

#define __Pyx_Owned_Py_None(b) (Py_INCREF(Py_None), Py_None)

#define __Pyx_PyBool_FromLong(b) ((b) ? (Py_INCREF(Py_True), Py_True) : (Py_INCREF(Py_False), Py_False))

static CYTHON_INLINE int __Pyx_PyObject_IsTrue(PyObject*);

static CYTHON_INLINE PyObject* __Pyx_PyNumber_Int(PyObject* x);

static CYTHON_INLINE Py_ssize_t __Pyx_PyIndex_AsSsize_t(PyObject*);

static CYTHON_INLINE PyObject * __Pyx_PyInt_FromSize_t(size_t);

static CYTHON_INLINE size_t __Pyx_PyInt_AsSize_t(PyObject*);

#define __pyx_PyFloat_AsDouble(x) (PyFloat_CheckExact(x) ? PyFloat_AS_DOUBLE(x) : PyFloat_AsDouble(x))

#ifdef __GNUC__

  /* Test for GCC > 2.95 */

  #if __GNUC__ > 2 || (__GNUC__ == 2 && (__GNUC_MINOR__ > 95))

    #define likely(x)   __builtin_expect(!!(x), 1)

    #define unlikely(x) __builtin_expect(!!(x), 0)

  #else /* __GNUC__ > 2 ... */

    #define likely(x)   (x)

    #define unlikely(x) (x)

  #endif /* __GNUC__ > 2 ... */

#else /* __GNUC__ */

  #define likely(x)   (x)

  #define unlikely(x) (x)

#endif /* __GNUC__ */

static PyObject *__pyx_m;

static PyObject *__pyx_b;

static PyObject *__pyx_empty_tuple;

static PyObject *__pyx_empty_bytes;

static int __pyx_lineno;

static int __pyx_clineno = 0;

static const char * __pyx_cfilenm= __FILE__;

static const char *__pyx_filename;

static const char *__pyx_f[] = {

  "Metropolis.pyx",

};

/*--- Type declarations ---*/

#ifndef CYTHON_REFNANNY

  #define CYTHON_REFNANNY 0

#endif

#if CYTHON_REFNANNY

  typedef struct {

    void (*INCREF)(void*, PyObject*, int);

    void (*DECREF)(void*, PyObject*, int);

    void (*GOTREF)(void*, PyObject*, int);

    void (*GIVEREF)(void*, PyObject*, int);

    void* (*SetupContext)(const char*, int, const char*);

    void (*FinishContext)(void**);

  } __Pyx_RefNannyAPIStruct;

  static __Pyx_RefNannyAPIStruct *__Pyx_RefNanny = NULL;

  static __Pyx_RefNannyAPIStruct *__Pyx_RefNannyImportAPI(const char *modname); /*proto*/

  #define __Pyx_RefNannyDeclarations void *__pyx_refnanny = NULL;

  #define __Pyx_RefNannySetupContext(name)           __pyx_refnanny = __Pyx_RefNanny->SetupContext((name), __LINE__, __FILE__)

  #define __Pyx_RefNannyFinishContext()           __Pyx_RefNanny->FinishContext(&__pyx_refnanny)

  #define __Pyx_INCREF(r)  __Pyx_RefNanny->INCREF(__pyx_refnanny, (PyObject *)(r), __LINE__)

  #define __Pyx_DECREF(r)  __Pyx_RefNanny->DECREF(__pyx_refnanny, (PyObject *)(r), __LINE__)

  #define __Pyx_GOTREF(r)  __Pyx_RefNanny->GOTREF(__pyx_refnanny, (PyObject *)(r), __LINE__)

  #define __Pyx_GIVEREF(r) __Pyx_RefNanny->GIVEREF(__pyx_refnanny, (PyObject *)(r), __LINE__)

  #define __Pyx_XINCREF(r)  do { if((r) != NULL) {__Pyx_INCREF(r); }} while(0)

  #define __Pyx_XDECREF(r)  do { if((r) != NULL) {__Pyx_DECREF(r); }} while(0)

  #define __Pyx_XGOTREF(r)  do { if((r) != NULL) {__Pyx_GOTREF(r); }} while(0)

  #define __Pyx_XGIVEREF(r) do { if((r) != NULL) {__Pyx_GIVEREF(r);}} while(0)

#else

  #define __Pyx_RefNannyDeclarations

  #define __Pyx_RefNannySetupContext(name)

  #define __Pyx_RefNannyFinishContext()

  #define __Pyx_INCREF(r) Py_INCREF(r)

  #define __Pyx_DECREF(r) Py_DECREF(r)

  #define __Pyx_GOTREF(r)

  #define __Pyx_GIVEREF(r)

  #define __Pyx_XINCREF(r) Py_XINCREF(r)

  #define __Pyx_XDECREF(r) Py_XDECREF(r)

  #define __Pyx_XGOTREF(r)

  #define __Pyx_XGIVEREF(r)

#endif /* CYTHON_REFNANNY */

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

static void __Pyx_RaiseArgtupleInvalid(const char* func_name, int exact,

    Py_ssize_t num_min, Py_ssize_t num_max, Py_ssize_t num_found); /*proto*/

static void __Pyx_RaiseDoubleKeywordsError(

    const char* func_name, PyObject* kw_name); /*proto*/

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

static CYTHON_INLINE void __Pyx_RaiseNeedMoreValuesError(Py_ssize_t index);

static CYTHON_INLINE void __Pyx_RaiseTooManyValuesError(Py_ssize_t expected);

static int __Pyx_IternextUnpackEndCheck(PyObject *retval, Py_ssize_t expected); /*proto*/

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

static CYTHON_INLINE unsigned char __Pyx_PyInt_AsUnsignedChar(PyObject *);

static CYTHON_INLINE unsigned short __Pyx_PyInt_AsUnsignedShort(PyObject *);

static CYTHON_INLINE unsigned int __Pyx_PyInt_AsUnsignedInt(PyObject *);

static CYTHON_INLINE char __Pyx_PyInt_AsChar(PyObject *);

static CYTHON_INLINE short __Pyx_PyInt_AsShort(PyObject *);

static CYTHON_INLINE int __Pyx_PyInt_AsInt(PyObject *);

static CYTHON_INLINE signed char __Pyx_PyInt_AsSignedChar(PyObject *);

static CYTHON_INLINE signed short __Pyx_PyInt_AsSignedShort(PyObject *);

static CYTHON_INLINE signed int __Pyx_PyInt_AsSignedInt(PyObject *);

static CYTHON_INLINE int __Pyx_PyInt_AsLongDouble(PyObject *);

static CYTHON_INLINE unsigned long __Pyx_PyInt_AsUnsignedLong(PyObject *);

static CYTHON_INLINE unsigned PY_LONG_LONG __Pyx_PyInt_AsUnsignedLongLong(PyObject *);

static CYTHON_INLINE long __Pyx_PyInt_AsLong(PyObject *);

static CYTHON_INLINE PY_LONG_LONG __Pyx_PyInt_AsLongLong(PyObject *);

static CYTHON_INLINE signed long __Pyx_PyInt_AsSignedLong(PyObject *);

static CYTHON_INLINE signed PY_LONG_LONG __Pyx_PyInt_AsSignedLongLong(PyObject *);

static int __Pyx_check_binary_version(void);

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

                               int __pyx_lineno, const char *__pyx_filename); /*proto*/

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

/* Module declarations from 'Metropolis' */

#define __Pyx_MODULE_NAME "Metropolis"

int __pyx_module_is_main_Metropolis = 0;

/* Implementation of 'Metropolis' */

static PyObject *__pyx_builtin_xrange;

static char __pyx_k_1[] = "matplotlib.pyplot";

static char __pyx_k_2[] = "*";

static char __pyx_k__B[] = "B";

static char __pyx_k__J[] = "J";

static char __pyx_k__T[] = "T";

static char __pyx_k__PIL[] = "PIL";

static char __pyx_k__exp[] = "exp";

static char __pyx_k__plt[] = "plt";

static char __pyx_k__sys[] = "sys";

static char __pyx_k__math[] = "math";

static char __pyx_k__time[] = "time";

static char __pyx_k__Image[] = "Image";

static char __pyx_k__numpy[] = "numpy";

static char __pyx_k__range[] = "range";

static char __pyx_k__shape[] = "shape";

static char __pyx_k__sigma[] = "sigma";

static char __pyx_k__getopt[] = "getopt";

static char __pyx_k__random[] = "random";

static char __pyx_k__xrange[] = "xrange";

static char __pyx_k__randint[] = "randint";

static char __pyx_k____main__[] = "__main__";

static char __pyx_k____test__[] = "__test__";

static char __pyx_k__Metropolis[] = "Metropolis";

static char __pyx_k__iterations[] = "iterations";

static PyObject *__pyx_n_s_1;

static PyObject *__pyx_n_s_2;

static PyObject *__pyx_n_s__B;

static PyObject *__pyx_n_s__Image;

static PyObject *__pyx_n_s__J;

static PyObject *__pyx_n_s__Metropolis;

static PyObject *__pyx_n_s__PIL;

static PyObject *__pyx_n_s__T;

static PyObject *__pyx_n_s____main__;

static PyObject *__pyx_n_s____test__;

static PyObject *__pyx_n_s__exp;

static PyObject *__pyx_n_s__getopt;

static PyObject *__pyx_n_s__iterations;

static PyObject *__pyx_n_s__math;

static PyObject *__pyx_n_s__numpy;

static PyObject *__pyx_n_s__plt;

static PyObject *__pyx_n_s__randint;

static PyObject *__pyx_n_s__random;

static PyObject *__pyx_n_s__range;

static PyObject *__pyx_n_s__shape;

static PyObject *__pyx_n_s__sigma;

static PyObject *__pyx_n_s__sys;

static PyObject *__pyx_n_s__time;

static PyObject *__pyx_n_s__xrange;

static PyObject *__pyx_int_0;

static PyObject *__pyx_int_1;

static PyObject *__pyx_int_2;

/* "Metropolis.pyx":17

 * import Metropolis

 * 

 * def Metropolis(sigma,J,B,T,iterations):             # <<<<<<<<<<<<<<

 *     start=time.time()

 * 

 */

static PyObject *__pyx_pf_10Metropolis_Metropolis(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/

static PyMethodDef __pyx_mdef_10Metropolis_Metropolis = {__Pyx_NAMESTR("Metropolis"), (PyCFunction)__pyx_pf_10Metropolis_Metropolis, METH_VARARGS|METH_KEYWORDS, __Pyx_DOCSTR(0)};

static PyObject *__pyx_pf_10Metropolis_Metropolis(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) {

  PyObject *__pyx_v_sigma = 0;

  PyObject *__pyx_v_J = 0;

  PyObject *__pyx_v_B = 0;

  PyObject *__pyx_v_T = 0;

  PyObject *__pyx_v_iterations = 0;

  PyObject *__pyx_v_start = NULL;

  PyObject *__pyx_v_SizeX = NULL;

  PyObject *__pyx_v_SizeY = NULL;

  PyObject *__pyx_v_p = NULL;

  PyObject *__pyx_v_X = NULL;

  PyObject *__pyx_v_Y = NULL;

  PyObject *__pyx_v_DeltaE = NULL;

  PyObject *__pyx_v_duration = NULL;

  PyObject *__pyx_r = NULL;

  __Pyx_RefNannyDeclarations

  PyObject *__pyx_t_1 = NULL;

  PyObject *__pyx_t_2 = NULL;

  PyObject *__pyx_t_3 = NULL;

  PyObject *__pyx_t_4 = NULL;

  PyObject *(*__pyx_t_5)(PyObject *);

  Py_ssize_t __pyx_t_6;

  PyObject *(*__pyx_t_7)(PyObject *);

  PyObject *__pyx_t_8 = NULL;

  int __pyx_t_9;

  int __pyx_t_10;

  int __pyx_t_11;

  int __pyx_lineno = 0;

  const char *__pyx_filename = NULL;

  int __pyx_clineno = 0;

  static PyObject **__pyx_pyargnames[] = {&__pyx_n_s__sigma,&__pyx_n_s__J,&__pyx_n_s__B,&__pyx_n_s__T,&__pyx_n_s__iterations,0};

  __Pyx_RefNannySetupContext("Metropolis");

  __pyx_self = __pyx_self;

  {

    PyObject* values[5] = {0,0,0,0,0};

    if (unlikely(__pyx_kwds)) {

      Py_ssize_t kw_args;

      switch (PyTuple_GET_SIZE(__pyx_args)) {

        case  5: values[4] = PyTuple_GET_ITEM(__pyx_args, 4);

        case  4: values[3] = PyTuple_GET_ITEM(__pyx_args, 3);

        case  3: values[2] = PyTuple_GET_ITEM(__pyx_args, 2);

        case  2: values[1] = PyTuple_GET_ITEM(__pyx_args, 1);

        case  1: values[0] = PyTuple_GET_ITEM(__pyx_args, 0);

        case  0: break;

        default: goto __pyx_L5_argtuple_error;

      }

      kw_args = PyDict_Size(__pyx_kwds);

      switch (PyTuple_GET_SIZE(__pyx_args)) {

        case  0:

        values[0] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__sigma);

        if (likely(values[0])) kw_args--;

        else goto __pyx_L5_argtuple_error;

        case  1:

        values[1] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__J);

        if (likely(values[1])) kw_args--;

        else {

          __Pyx_RaiseArgtupleInvalid("Metropolis", 1, 5, 5, 1); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L3_error;}

        }

        case  2:

        values[2] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__B);

        if (likely(values[2])) kw_args--;

        else {

          __Pyx_RaiseArgtupleInvalid("Metropolis", 1, 5, 5, 2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L3_error;}

        }

        case  3:

        values[3] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__T);

        if (likely(values[3])) kw_args--;

        else {

          __Pyx_RaiseArgtupleInvalid("Metropolis", 1, 5, 5, 3); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L3_error;}

        }

        case  4:

        values[4] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__iterations);

        if (likely(values[4])) kw_args--;

        else {

          __Pyx_RaiseArgtupleInvalid("Metropolis", 1, 5, 5, 4); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L3_error;}

        }

      }

      if (unlikely(kw_args > 0)) {

        if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_pyargnames, 0, values, PyTuple_GET_SIZE(__pyx_args), "Metropolis") < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L3_error;}

      }

    } else if (PyTuple_GET_SIZE(__pyx_args) != 5) {

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

    }

    __pyx_v_sigma = values[0];

    __pyx_v_J = values[1];

    __pyx_v_B = values[2];

    __pyx_v_T = values[3];

    __pyx_v_iterations = values[4];

  }

  goto __pyx_L4_argument_unpacking_done;

  __pyx_L5_argtuple_error:;

  __Pyx_RaiseArgtupleInvalid("Metropolis", 1, 5, 5, PyTuple_GET_SIZE(__pyx_args)); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L3_error;}

  __pyx_L3_error:;

  __Pyx_AddTraceback("Metropolis.Metropolis", __pyx_clineno, __pyx_lineno, __pyx_filename);

  __Pyx_RefNannyFinishContext();

  return NULL;

  __pyx_L4_argument_unpacking_done:;

  /* "Metropolis.pyx":18

 * 

 * def Metropolis(sigma,J,B,T,iterations):

 *     start=time.time()             # <<<<<<<<<<<<<<

 * 

 *     SizeX,SizeY=sigma.shape

 */

  __pyx_t_1 = __Pyx_GetName(__pyx_m, __pyx_n_s__time); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 18; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  __pyx_t_2 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__time); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 18; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_2);

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  __pyx_t_1 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 18; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

  __pyx_v_start = __pyx_t_1;

  __pyx_t_1 = 0;

  /* "Metropolis.pyx":20

 *     start=time.time()

 * 

 *     SizeX,SizeY=sigma.shape             # <<<<<<<<<<<<<<

 * 

 *     for p in xrange(0,iterations):

 */

  __pyx_t_1 = PyObject_GetAttr(__pyx_v_sigma, __pyx_n_s__shape); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 20; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  if ((likely(PyTuple_CheckExact(__pyx_t_1))) || (PyList_CheckExact(__pyx_t_1))) {

    PyObject* sequence = __pyx_t_1;

    if (likely(PyTuple_CheckExact(sequence))) {

      if (unlikely(PyTuple_GET_SIZE(sequence) != 2)) {

        if (PyTuple_GET_SIZE(sequence) > 2) __Pyx_RaiseTooManyValuesError(2);

        else __Pyx_RaiseNeedMoreValuesError(PyTuple_GET_SIZE(sequence));

        {__pyx_filename = __pyx_f[0]; __pyx_lineno = 20; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      }

      __pyx_t_2 = PyTuple_GET_ITEM(sequence, 0); 

      __pyx_t_3 = PyTuple_GET_ITEM(sequence, 1); 

    } else {

      if (unlikely(PyList_GET_SIZE(sequence) != 2)) {

        if (PyList_GET_SIZE(sequence) > 2) __Pyx_RaiseTooManyValuesError(2);

        else __Pyx_RaiseNeedMoreValuesError(PyList_GET_SIZE(sequence));

        {__pyx_filename = __pyx_f[0]; __pyx_lineno = 20; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      }

      __pyx_t_2 = PyList_GET_ITEM(sequence, 0); 

      __pyx_t_3 = PyList_GET_ITEM(sequence, 1); 

    }

    __Pyx_INCREF(__pyx_t_2);

    __Pyx_INCREF(__pyx_t_3);

    __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  } else {

    Py_ssize_t index = -1;

    __pyx_t_4 = PyObject_GetIter(__pyx_t_1); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 20; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_4);

    __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

    __pyx_t_5 = Py_TYPE(__pyx_t_4)->tp_iternext;

    index = 0; __pyx_t_2 = __pyx_t_5(__pyx_t_4); if (unlikely(!__pyx_t_2)) goto __pyx_L6_unpacking_failed;

    __Pyx_GOTREF(__pyx_t_2);

    index = 1; __pyx_t_3 = __pyx_t_5(__pyx_t_4); if (unlikely(!__pyx_t_3)) goto __pyx_L6_unpacking_failed;

    __Pyx_GOTREF(__pyx_t_3);

    if (__Pyx_IternextUnpackEndCheck(__pyx_t_5(__pyx_t_4), 2) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 20; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;

    goto __pyx_L7_unpacking_done;

    __pyx_L6_unpacking_failed:;

    __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;

    if (PyErr_Occurred() && PyErr_ExceptionMatches(PyExc_StopIteration)) PyErr_Clear();

    if (!PyErr_Occurred()) __Pyx_RaiseNeedMoreValuesError(index);

    {__pyx_filename = __pyx_f[0]; __pyx_lineno = 20; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __pyx_L7_unpacking_done:;

  }

  __pyx_v_SizeX = __pyx_t_2;

  __pyx_t_2 = 0;

  __pyx_v_SizeY = __pyx_t_3;

  __pyx_t_3 = 0;

  /* "Metropolis.pyx":22

 *     SizeX,SizeY=sigma.shape

 * 

 *     for p in xrange(0,iterations):             # <<<<<<<<<<<<<<

 *         # Random access coordonate

 *         X,Y=numpy.random.randint(SizeX),numpy.random.randint(SizeY)

 */

  __pyx_t_1 = PyTuple_New(2); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 22; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(((PyObject *)__pyx_t_1));

  __Pyx_INCREF(__pyx_int_0);

  PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_int_0);

  __Pyx_GIVEREF(__pyx_int_0);

  __Pyx_INCREF(__pyx_v_iterations);

  PyTuple_SET_ITEM(__pyx_t_1, 1, __pyx_v_iterations);

  __Pyx_GIVEREF(__pyx_v_iterations);

  __pyx_t_3 = PyObject_Call(__pyx_builtin_xrange, ((PyObject *)__pyx_t_1), NULL); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 22; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_3);

  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;

  if (PyList_CheckExact(__pyx_t_3) || PyTuple_CheckExact(__pyx_t_3)) {

    __pyx_t_1 = __pyx_t_3; __Pyx_INCREF(__pyx_t_1); __pyx_t_6 = 0;

    __pyx_t_7 = NULL;

  } else {

    __pyx_t_6 = -1; __pyx_t_1 = PyObject_GetIter(__pyx_t_3); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 22; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_1);

    __pyx_t_7 = Py_TYPE(__pyx_t_1)->tp_iternext;

  }

  __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

  for (;;) {

    if (PyList_CheckExact(__pyx_t_1)) {

      if (__pyx_t_6 >= PyList_GET_SIZE(__pyx_t_1)) break;

      __pyx_t_3 = PyList_GET_ITEM(__pyx_t_1, __pyx_t_6); __Pyx_INCREF(__pyx_t_3); __pyx_t_6++;

    } else if (PyTuple_CheckExact(__pyx_t_1)) {

      if (__pyx_t_6 >= PyTuple_GET_SIZE(__pyx_t_1)) break;

      __pyx_t_3 = PyTuple_GET_ITEM(__pyx_t_1, __pyx_t_6); __Pyx_INCREF(__pyx_t_3); __pyx_t_6++;

    } else {

      __pyx_t_3 = __pyx_t_7(__pyx_t_1);

      if (unlikely(!__pyx_t_3)) {

        if (PyErr_Occurred()) {

          if (likely(PyErr_ExceptionMatches(PyExc_StopIteration))) PyErr_Clear();

          else {__pyx_filename = __pyx_f[0]; __pyx_lineno = 22; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

        }

        break;

      }

      __Pyx_GOTREF(__pyx_t_3);

    }

    __Pyx_XDECREF(__pyx_v_p);

    __pyx_v_p = __pyx_t_3;

    __pyx_t_3 = 0;

    /* "Metropolis.pyx":24

 *     for p in xrange(0,iterations):

 *         # Random access coordonate

 *         X,Y=numpy.random.randint(SizeX),numpy.random.randint(SizeY)             # <<<<<<<<<<<<<<

 * 

 *         DeltaE=J*sigma[X,Y]*(2*(sigma[X,(Y+1)%SizeY]+

 */

    __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_3);

    __pyx_t_2 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__random); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_2);

    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

    __pyx_t_3 = PyObject_GetAttr(__pyx_t_2, __pyx_n_s__randint); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_3);

    __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

    __pyx_t_2 = PyTuple_New(1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(((PyObject *)__pyx_t_2));

    __Pyx_INCREF(__pyx_v_SizeX);

    PyTuple_SET_ITEM(__pyx_t_2, 0, __pyx_v_SizeX);

    __Pyx_GIVEREF(__pyx_v_SizeX);

    __pyx_t_4 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_t_2), NULL); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_4);

    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

    __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0;

    __pyx_t_2 = __Pyx_GetName(__pyx_m, __pyx_n_s__numpy); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_2);

    __pyx_t_3 = PyObject_GetAttr(__pyx_t_2, __pyx_n_s__random); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_3);

    __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

    __pyx_t_2 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__randint); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_2);

    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

    __pyx_t_3 = PyTuple_New(1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(((PyObject *)__pyx_t_3));

    __Pyx_INCREF(__pyx_v_SizeY);

    PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_v_SizeY);

    __Pyx_GIVEREF(__pyx_v_SizeY);

    __pyx_t_8 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_t_3), NULL); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_8);

    __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

    __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0;

    __Pyx_XDECREF(__pyx_v_X);

    __pyx_v_X = __pyx_t_4;

    __pyx_t_4 = 0;

    __Pyx_XDECREF(__pyx_v_Y);

    __pyx_v_Y = __pyx_t_8;

    __pyx_t_8 = 0;

    /* "Metropolis.pyx":26

 *         X,Y=numpy.random.randint(SizeX),numpy.random.randint(SizeY)

 * 

 *         DeltaE=J*sigma[X,Y]*(2*(sigma[X,(Y+1)%SizeY]+             # <<<<<<<<<<<<<<

 *                                 sigma[X,(Y-1)%SizeY]+

 *                                 sigma[(X-1)%SizeX,Y]+

 */

    __pyx_t_8 = PyTuple_New(2); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(((PyObject *)__pyx_t_8));

    __Pyx_INCREF(__pyx_v_X);

    PyTuple_SET_ITEM(__pyx_t_8, 0, __pyx_v_X);

    __Pyx_GIVEREF(__pyx_v_X);

    __Pyx_INCREF(__pyx_v_Y);

    PyTuple_SET_ITEM(__pyx_t_8, 1, __pyx_v_Y);

    __Pyx_GIVEREF(__pyx_v_Y);

    __pyx_t_4 = PyObject_GetItem(__pyx_v_sigma, ((PyObject *)__pyx_t_8)); if (!__pyx_t_4) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_4);

    __Pyx_DECREF(((PyObject *)__pyx_t_8)); __pyx_t_8 = 0;

    __pyx_t_8 = PyNumber_Multiply(__pyx_v_J, __pyx_t_4); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_8);

    __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;

    /* "Metropolis.pyx":27

 * 

 *         DeltaE=J*sigma[X,Y]*(2*(sigma[X,(Y+1)%SizeY]+

 *                                 sigma[X,(Y-1)%SizeY]+             # <<<<<<<<<<<<<<

 *                                 sigma[(X-1)%SizeX,Y]+

 *                                 sigma[(X+1)%SizeX,Y])+B)

 */

    __pyx_t_4 = PyNumber_Add(__pyx_v_Y, __pyx_int_1); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_4);

    /* "Metropolis.pyx":26

 *         X,Y=numpy.random.randint(SizeX),numpy.random.randint(SizeY)

 * 

 *         DeltaE=J*sigma[X,Y]*(2*(sigma[X,(Y+1)%SizeY]+             # <<<<<<<<<<<<<<

 *                                 sigma[X,(Y-1)%SizeY]+

 *                                 sigma[(X-1)%SizeX,Y]+

 */

    __pyx_t_3 = PyNumber_Remainder(__pyx_t_4, __pyx_v_SizeY); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_3);

    __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;

    __pyx_t_4 = PyTuple_New(2); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(((PyObject *)__pyx_t_4));

    __Pyx_INCREF(__pyx_v_X);

    PyTuple_SET_ITEM(__pyx_t_4, 0, __pyx_v_X);

    __Pyx_GIVEREF(__pyx_v_X);

    PyTuple_SET_ITEM(__pyx_t_4, 1, __pyx_t_3);

    __Pyx_GIVEREF(__pyx_t_3);

    __pyx_t_3 = 0;

    __pyx_t_3 = PyObject_GetItem(__pyx_v_sigma, ((PyObject *)__pyx_t_4)); if (!__pyx_t_3) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_3);

    __Pyx_DECREF(((PyObject *)__pyx_t_4)); __pyx_t_4 = 0;

    /* "Metropolis.pyx":27

 * 

 *         DeltaE=J*sigma[X,Y]*(2*(sigma[X,(Y+1)%SizeY]+

 *                                 sigma[X,(Y-1)%SizeY]+             # <<<<<<<<<<<<<<

 *                                 sigma[(X-1)%SizeX,Y]+

 *                                 sigma[(X+1)%SizeX,Y])+B)

 */

    __pyx_t_4 = PyNumber_Subtract(__pyx_v_Y, __pyx_int_1); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_4);

    __pyx_t_2 = PyNumber_Remainder(__pyx_t_4, __pyx_v_SizeY); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_2);

    __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;

    __pyx_t_4 = PyTuple_New(2); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(((PyObject *)__pyx_t_4));

    __Pyx_INCREF(__pyx_v_X);

    PyTuple_SET_ITEM(__pyx_t_4, 0, __pyx_v_X);

    __Pyx_GIVEREF(__pyx_v_X);

    PyTuple_SET_ITEM(__pyx_t_4, 1, __pyx_t_2);

    __Pyx_GIVEREF(__pyx_t_2);

    __pyx_t_2 = 0;

    __pyx_t_2 = PyObject_GetItem(__pyx_v_sigma, ((PyObject *)__pyx_t_4)); if (!__pyx_t_2) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_2);

    __Pyx_DECREF(((PyObject *)__pyx_t_4)); __pyx_t_4 = 0;

    __pyx_t_4 = PyNumber_Add(__pyx_t_3, __pyx_t_2); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_4);

    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

    __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

    /* "Metropolis.pyx":28

 *         DeltaE=J*sigma[X,Y]*(2*(sigma[X,(Y+1)%SizeY]+

 *                                 sigma[X,(Y-1)%SizeY]+

 *                                 sigma[(X-1)%SizeX,Y]+             # <<<<<<<<<<<<<<

 *                                 sigma[(X+1)%SizeX,Y])+B)

 * 

 */

    __pyx_t_2 = PyNumber_Subtract(__pyx_v_X, __pyx_int_1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_2);

    __pyx_t_3 = PyNumber_Remainder(__pyx_t_2, __pyx_v_SizeX); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_3);

    __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

    __pyx_t_2 = PyTuple_New(2); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(((PyObject *)__pyx_t_2));

    PyTuple_SET_ITEM(__pyx_t_2, 0, __pyx_t_3);

    __Pyx_GIVEREF(__pyx_t_3);

    __Pyx_INCREF(__pyx_v_Y);

    PyTuple_SET_ITEM(__pyx_t_2, 1, __pyx_v_Y);

    __Pyx_GIVEREF(__pyx_v_Y);

    __pyx_t_3 = 0;

    __pyx_t_3 = PyObject_GetItem(__pyx_v_sigma, ((PyObject *)__pyx_t_2)); if (!__pyx_t_3) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_3);

    __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0;

    __pyx_t_2 = PyNumber_Add(__pyx_t_4, __pyx_t_3); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_2);

    __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;

    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

    /* "Metropolis.pyx":29

 *                                 sigma[X,(Y-1)%SizeY]+

 *                                 sigma[(X-1)%SizeX,Y]+

 *                                 sigma[(X+1)%SizeX,Y])+B)             # <<<<<<<<<<<<<<

 * 

 *         if DeltaE < 0. or random() < exp(-DeltaE/T):

 */

    __pyx_t_3 = PyNumber_Add(__pyx_v_X, __pyx_int_1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_3);

    __pyx_t_4 = PyNumber_Remainder(__pyx_t_3, __pyx_v_SizeX); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_4);

    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

    __pyx_t_3 = PyTuple_New(2); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(((PyObject *)__pyx_t_3));

    PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_t_4);

    __Pyx_GIVEREF(__pyx_t_4);

    __Pyx_INCREF(__pyx_v_Y);

    PyTuple_SET_ITEM(__pyx_t_3, 1, __pyx_v_Y);

    __Pyx_GIVEREF(__pyx_v_Y);

    __pyx_t_4 = 0;

    __pyx_t_4 = PyObject_GetItem(__pyx_v_sigma, ((PyObject *)__pyx_t_3)); if (!__pyx_t_4) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_4);

    __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0;

    __pyx_t_3 = PyNumber_Add(__pyx_t_2, __pyx_t_4); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_3);

    __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

    __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;

    __pyx_t_4 = PyNumber_Multiply(__pyx_int_2, __pyx_t_3); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_4);

    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

    __pyx_t_3 = PyNumber_Add(__pyx_t_4, __pyx_v_B); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_3);

    __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;

    __pyx_t_4 = PyNumber_Multiply(__pyx_t_8, __pyx_t_3); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_4);

    __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0;

    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

    __Pyx_XDECREF(__pyx_v_DeltaE);

    __pyx_v_DeltaE = __pyx_t_4;

    __pyx_t_4 = 0;

    /* "Metropolis.pyx":31

 *                                 sigma[(X+1)%SizeX,Y])+B)

 * 

 *         if DeltaE < 0. or random() < exp(-DeltaE/T):             # <<<<<<<<<<<<<<

 *             sigma[X,Y]=-sigma[X,Y]

 *     duration=time.time()-start

 */

    __pyx_t_4 = PyFloat_FromDouble(0.); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_4);

    __pyx_t_3 = PyObject_RichCompare(__pyx_v_DeltaE, __pyx_t_4, Py_LT); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_GOTREF(__pyx_t_3);

    __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;

    __pyx_t_9 = __Pyx_PyObject_IsTrue(__pyx_t_3); if (unlikely(__pyx_t_9 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

    __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

    if (!__pyx_t_9) {

      __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__random); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(__pyx_t_3);

      __pyx_t_4 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_4)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(__pyx_t_4);

      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

      __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__exp); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(__pyx_t_3);

      __pyx_t_8 = PyNumber_Negative(__pyx_v_DeltaE); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(__pyx_t_8);

      __pyx_t_2 = __Pyx_PyNumber_Divide(__pyx_t_8, __pyx_v_T); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(__pyx_t_2);

      __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0;

      __pyx_t_8 = PyTuple_New(1); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(((PyObject *)__pyx_t_8));

      PyTuple_SET_ITEM(__pyx_t_8, 0, __pyx_t_2);

      __Pyx_GIVEREF(__pyx_t_2);

      __pyx_t_2 = 0;

      __pyx_t_2 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_t_8), NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(__pyx_t_2);

      __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0;

      __Pyx_DECREF(((PyObject *)__pyx_t_8)); __pyx_t_8 = 0;

      __pyx_t_8 = PyObject_RichCompare(__pyx_t_4, __pyx_t_2, Py_LT); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(__pyx_t_8);

      __Pyx_DECREF(__pyx_t_4); __pyx_t_4 = 0;

      __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

      __pyx_t_10 = __Pyx_PyObject_IsTrue(__pyx_t_8); if (unlikely(__pyx_t_10 < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0;

      __pyx_t_11 = __pyx_t_10;

    } else {

      __pyx_t_11 = __pyx_t_9;

    }

    if (__pyx_t_11) {

      /* "Metropolis.pyx":32

 * 

 *         if DeltaE < 0. or random() < exp(-DeltaE/T):

 *             sigma[X,Y]=-sigma[X,Y]             # <<<<<<<<<<<<<<

 *     duration=time.time()-start

 *     return(duration)

 */

      __pyx_t_8 = PyTuple_New(2); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(((PyObject *)__pyx_t_8));

      __Pyx_INCREF(__pyx_v_X);

      PyTuple_SET_ITEM(__pyx_t_8, 0, __pyx_v_X);

      __Pyx_GIVEREF(__pyx_v_X);

      __Pyx_INCREF(__pyx_v_Y);

      PyTuple_SET_ITEM(__pyx_t_8, 1, __pyx_v_Y);

      __Pyx_GIVEREF(__pyx_v_Y);

      __pyx_t_2 = PyObject_GetItem(__pyx_v_sigma, ((PyObject *)__pyx_t_8)); if (!__pyx_t_2) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(__pyx_t_2);

      __Pyx_DECREF(((PyObject *)__pyx_t_8)); __pyx_t_8 = 0;

      __pyx_t_8 = PyNumber_Negative(__pyx_t_2); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(__pyx_t_8);

      __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

      __pyx_t_2 = PyTuple_New(2); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_GOTREF(((PyObject *)__pyx_t_2));

      __Pyx_INCREF(__pyx_v_X);

      PyTuple_SET_ITEM(__pyx_t_2, 0, __pyx_v_X);

      __Pyx_GIVEREF(__pyx_v_X);

      __Pyx_INCREF(__pyx_v_Y);

      PyTuple_SET_ITEM(__pyx_t_2, 1, __pyx_v_Y);

      __Pyx_GIVEREF(__pyx_v_Y);

      if (PyObject_SetItem(__pyx_v_sigma, ((PyObject *)__pyx_t_2), __pyx_t_8) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

      __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0;

      __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0;

      goto __pyx_L10;

    }

    __pyx_L10:;

  }

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  /* "Metropolis.pyx":33

 *         if DeltaE < 0. or random() < exp(-DeltaE/T):

 *             sigma[X,Y]=-sigma[X,Y]

 *     duration=time.time()-start             # <<<<<<<<<<<<<<

 *     return(duration)

 */

  __pyx_t_1 = __Pyx_GetName(__pyx_m, __pyx_n_s__time); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  __pyx_t_8 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__time); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_8);

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  __pyx_t_1 = PyObject_Call(__pyx_t_8, ((PyObject *)__pyx_empty_tuple), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0;

  __pyx_t_8 = PyNumber_Subtract(__pyx_t_1, __pyx_v_start); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_8);

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  __pyx_v_duration = __pyx_t_8;

  __pyx_t_8 = 0;

  /* "Metropolis.pyx":34

 *             sigma[X,Y]=-sigma[X,Y]

 *     duration=time.time()-start

 *     return(duration)             # <<<<<<<<<<<<<<

 */

  __Pyx_XDECREF(__pyx_r);

  __Pyx_INCREF(__pyx_v_duration);

  __pyx_r = __pyx_v_duration;

  goto __pyx_L0;

  __pyx_r = Py_None; __Pyx_INCREF(Py_None);

  goto __pyx_L0;

  __pyx_L1_error:;

  __Pyx_XDECREF(__pyx_t_1);

  __Pyx_XDECREF(__pyx_t_2);

  __Pyx_XDECREF(__pyx_t_3);

  __Pyx_XDECREF(__pyx_t_4);

  __Pyx_XDECREF(__pyx_t_8);

  __Pyx_AddTraceback("Metropolis.Metropolis", __pyx_clineno, __pyx_lineno, __pyx_filename);

  __pyx_r = NULL;

  __pyx_L0:;

  __Pyx_XDECREF(__pyx_v_start);

  __Pyx_XDECREF(__pyx_v_SizeX);

  __Pyx_XDECREF(__pyx_v_SizeY);

  __Pyx_XDECREF(__pyx_v_p);

  __Pyx_XDECREF(__pyx_v_X);

  __Pyx_XDECREF(__pyx_v_Y);

  __Pyx_XDECREF(__pyx_v_DeltaE);

  __Pyx_XDECREF(__pyx_v_duration);

  __Pyx_XGIVEREF(__pyx_r);

  __Pyx_RefNannyFinishContext();

  return __pyx_r;

}

static PyMethodDef __pyx_methods[] = {

  {0, 0, 0, 0}

};

#if PY_MAJOR_VERSION >= 3

static struct PyModuleDef __pyx_moduledef = {

    PyModuleDef_HEAD_INIT,

    __Pyx_NAMESTR("Metropolis"),

    0, /* m_doc */

    -1, /* m_size */

    __pyx_methods /* m_methods */,

    NULL, /* m_reload */

    NULL, /* m_traverse */

    NULL, /* m_clear */

    NULL /* m_free */

};

#endif

static __Pyx_StringTabEntry __pyx_string_tab[] = {

  {&__pyx_n_s_1, __pyx_k_1, sizeof(__pyx_k_1), 0, 0, 1, 1},

  {&__pyx_n_s_2, __pyx_k_2, sizeof(__pyx_k_2), 0, 0, 1, 1},

  {&__pyx_n_s__B, __pyx_k__B, sizeof(__pyx_k__B), 0, 0, 1, 1},

  {&__pyx_n_s__Image, __pyx_k__Image, sizeof(__pyx_k__Image), 0, 0, 1, 1},

  {&__pyx_n_s__J, __pyx_k__J, sizeof(__pyx_k__J), 0, 0, 1, 1},

  {&__pyx_n_s__Metropolis, __pyx_k__Metropolis, sizeof(__pyx_k__Metropolis), 0, 0, 1, 1},

  {&__pyx_n_s__PIL, __pyx_k__PIL, sizeof(__pyx_k__PIL), 0, 0, 1, 1},

  {&__pyx_n_s__T, __pyx_k__T, sizeof(__pyx_k__T), 0, 0, 1, 1},

  {&__pyx_n_s____main__, __pyx_k____main__, sizeof(__pyx_k____main__), 0, 0, 1, 1},

  {&__pyx_n_s____test__, __pyx_k____test__, sizeof(__pyx_k____test__), 0, 0, 1, 1},

  {&__pyx_n_s__exp, __pyx_k__exp, sizeof(__pyx_k__exp), 0, 0, 1, 1},

  {&__pyx_n_s__getopt, __pyx_k__getopt, sizeof(__pyx_k__getopt), 0, 0, 1, 1},

  {&__pyx_n_s__iterations, __pyx_k__iterations, sizeof(__pyx_k__iterations), 0, 0, 1, 1},

  {&__pyx_n_s__math, __pyx_k__math, sizeof(__pyx_k__math), 0, 0, 1, 1},

  {&__pyx_n_s__numpy, __pyx_k__numpy, sizeof(__pyx_k__numpy), 0, 0, 1, 1},

  {&__pyx_n_s__plt, __pyx_k__plt, sizeof(__pyx_k__plt), 0, 0, 1, 1},

  {&__pyx_n_s__randint, __pyx_k__randint, sizeof(__pyx_k__randint), 0, 0, 1, 1},

  {&__pyx_n_s__random, __pyx_k__random, sizeof(__pyx_k__random), 0, 0, 1, 1},

  {&__pyx_n_s__range, __pyx_k__range, sizeof(__pyx_k__range), 0, 0, 1, 1},

  {&__pyx_n_s__shape, __pyx_k__shape, sizeof(__pyx_k__shape), 0, 0, 1, 1},

  {&__pyx_n_s__sigma, __pyx_k__sigma, sizeof(__pyx_k__sigma), 0, 0, 1, 1},

  {&__pyx_n_s__sys, __pyx_k__sys, sizeof(__pyx_k__sys), 0, 0, 1, 1},

  {&__pyx_n_s__time, __pyx_k__time, sizeof(__pyx_k__time), 0, 0, 1, 1},

  {&__pyx_n_s__xrange, __pyx_k__xrange, sizeof(__pyx_k__xrange), 0, 0, 1, 1},

  {0, 0, 0, 0, 0, 0, 0}

};

static int __Pyx_InitCachedBuiltins(void) {

  #if PY_MAJOR_VERSION >= 3

  __pyx_builtin_xrange = __Pyx_GetName(__pyx_b, __pyx_n_s__range); if (!__pyx_builtin_xrange) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 22; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  #else

  __pyx_builtin_xrange = __Pyx_GetName(__pyx_b, __pyx_n_s__xrange); if (!__pyx_builtin_xrange) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 22; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  #endif

  return 0;

  __pyx_L1_error:;

  return -1;

}

static int __Pyx_InitCachedConstants(void) {

  __Pyx_RefNannyDeclarations

  __Pyx_RefNannySetupContext("__Pyx_InitCachedConstants");

  __Pyx_RefNannyFinishContext();

  return 0;

}

static int __Pyx_InitGlobals(void) {

  if (__Pyx_InitStrings(__pyx_string_tab) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};

  __pyx_int_0 = PyInt_FromLong(0); if (unlikely(!__pyx_int_0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};

  __pyx_int_1 = PyInt_FromLong(1); if (unlikely(!__pyx_int_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};

  __pyx_int_2 = PyInt_FromLong(2); if (unlikely(!__pyx_int_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};

  return 0;

  __pyx_L1_error:;

  return -1;

}

#if PY_MAJOR_VERSION < 3

PyMODINIT_FUNC initMetropolis(void); /*proto*/

PyMODINIT_FUNC initMetropolis(void)

#else

PyMODINIT_FUNC PyInit_Metropolis(void); /*proto*/

PyMODINIT_FUNC PyInit_Metropolis(void)

#endif

{

  PyObject *__pyx_t_1 = NULL;

  PyObject *__pyx_t_2 = NULL;

  __Pyx_RefNannyDeclarations

  #if CYTHON_REFNANNY

  __Pyx_RefNanny = __Pyx_RefNannyImportAPI("refnanny");

  if (!__Pyx_RefNanny) {

      PyErr_Clear();

      __Pyx_RefNanny = __Pyx_RefNannyImportAPI("Cython.Runtime.refnanny");

      if (!__Pyx_RefNanny)

          Py_FatalError("failed to import 'refnanny' module");

  }

  #endif

  __Pyx_RefNannySetupContext("PyMODINIT_FUNC PyInit_Metropolis(void)");

  if ( __Pyx_check_binary_version() < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __pyx_empty_tuple = PyTuple_New(0); if (unlikely(!__pyx_empty_tuple)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __pyx_empty_bytes = PyBytes_FromStringAndSize("", 0); if (unlikely(!__pyx_empty_bytes)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  #ifdef __pyx_binding_PyCFunctionType_USED

  if (__pyx_binding_PyCFunctionType_init() < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  #endif

  /*--- Library function declarations ---*/

  /*--- Threads initialization code ---*/

  #if defined(__PYX_FORCE_INIT_THREADS) && __PYX_FORCE_INIT_THREADS

  #ifdef WITH_THREAD /* Python build with threading support? */

  PyEval_InitThreads();

  #endif

  #endif

  /*--- Module creation code ---*/

  #if PY_MAJOR_VERSION < 3

  __pyx_m = Py_InitModule4(__Pyx_NAMESTR("Metropolis"), __pyx_methods, 0, 0, PYTHON_API_VERSION);

  #else

  __pyx_m = PyModule_Create(&__pyx_moduledef);

  #endif

  if (!__pyx_m) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};

  #if PY_MAJOR_VERSION < 3

  Py_INCREF(__pyx_m);

  #endif

  __pyx_b = PyImport_AddModule(__Pyx_NAMESTR(__Pyx_BUILTIN_MODULE_NAME));

  if (!__pyx_b) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};

  if (__Pyx_SetAttrString(__pyx_m, "__builtins__", __pyx_b) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};

  /*--- Initialize various global constants etc. ---*/

  if (unlikely(__Pyx_InitGlobals() < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  if (__pyx_module_is_main_Metropolis) {

    if (__Pyx_SetAttrString(__pyx_m, "__name__", __pyx_n_s____main__) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;};

  }

  /*--- Builtin init code ---*/

  if (unlikely(__Pyx_InitCachedBuiltins() < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  /*--- Constants init code ---*/

  if (unlikely(__Pyx_InitCachedConstants() < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  /*--- Global init code ---*/

  /*--- Variable export code ---*/

  /*--- Function export code ---*/

  /*--- Type init code ---*/

  /*--- Type import code ---*/

  /*--- Variable import code ---*/

  /*--- Function import code ---*/

  /*--- Execution code ---*/

  /* "Metropolis.pyx":7

 * # CC BY-NC-SA 2011 : <emmanuel.quemener@ens-lyon.fr>

 * 

 * import sys             # <<<<<<<<<<<<<<

 * import numpy

 * from PIL import Image

 */

  __pyx_t_1 = __Pyx_Import(((PyObject *)__pyx_n_s__sys), 0, -1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 7; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  if (PyObject_SetAttr(__pyx_m, __pyx_n_s__sys, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 7; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  /* "Metropolis.pyx":8

 * 

 * import sys

 * import numpy             # <<<<<<<<<<<<<<

 * from PIL import Image

 * from math import exp

 */

  __pyx_t_1 = __Pyx_Import(((PyObject *)__pyx_n_s__numpy), 0, -1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 8; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  if (PyObject_SetAttr(__pyx_m, __pyx_n_s__numpy, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 8; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  /* "Metropolis.pyx":9

 * import sys

 * import numpy

 * from PIL import Image             # <<<<<<<<<<<<<<

 * from math import exp

 * from random import random

 */

  __pyx_t_1 = PyList_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 9; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(((PyObject *)__pyx_t_1));

  __Pyx_INCREF(((PyObject *)__pyx_n_s__Image));

  PyList_SET_ITEM(__pyx_t_1, 0, ((PyObject *)__pyx_n_s__Image));

  __Pyx_GIVEREF(((PyObject *)__pyx_n_s__Image));

  __pyx_t_2 = __Pyx_Import(((PyObject *)__pyx_n_s__PIL), ((PyObject *)__pyx_t_1), -1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 9; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_2);

  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;

  __pyx_t_1 = PyObject_GetAttr(__pyx_t_2, __pyx_n_s__Image); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 9; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  if (PyObject_SetAttr(__pyx_m, __pyx_n_s__Image, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 9; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

  /* "Metropolis.pyx":10

 * import numpy

 * from PIL import Image

 * from math import exp             # <<<<<<<<<<<<<<

 * from random import random

 * import time

 */

  __pyx_t_2 = PyList_New(1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 10; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(((PyObject *)__pyx_t_2));

  __Pyx_INCREF(((PyObject *)__pyx_n_s__exp));

  PyList_SET_ITEM(__pyx_t_2, 0, ((PyObject *)__pyx_n_s__exp));

  __Pyx_GIVEREF(((PyObject *)__pyx_n_s__exp));

  __pyx_t_1 = __Pyx_Import(((PyObject *)__pyx_n_s__math), ((PyObject *)__pyx_t_2), -1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 10; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0;

  __pyx_t_2 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__exp); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 10; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_2);

  if (PyObject_SetAttr(__pyx_m, __pyx_n_s__exp, __pyx_t_2) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 10; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  /* "Metropolis.pyx":11

 * from PIL import Image

 * from math import exp

 * from random import random             # <<<<<<<<<<<<<<

 * import time

 * import getopt

 */

  __pyx_t_1 = PyList_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 11; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(((PyObject *)__pyx_t_1));

  __Pyx_INCREF(((PyObject *)__pyx_n_s__random));

  PyList_SET_ITEM(__pyx_t_1, 0, ((PyObject *)__pyx_n_s__random));

  __Pyx_GIVEREF(((PyObject *)__pyx_n_s__random));

  __pyx_t_2 = __Pyx_Import(((PyObject *)__pyx_n_s__random), ((PyObject *)__pyx_t_1), -1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 11; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_2);

  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;

  __pyx_t_1 = PyObject_GetAttr(__pyx_t_2, __pyx_n_s__random); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 11; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  if (PyObject_SetAttr(__pyx_m, __pyx_n_s__random, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 11; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

  /* "Metropolis.pyx":12

 * from math import exp

 * from random import random

 * import time             # <<<<<<<<<<<<<<

 * import getopt

 * import matplotlib.pyplot as plt

 */

  __pyx_t_2 = __Pyx_Import(((PyObject *)__pyx_n_s__time), 0, -1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_2);

  if (PyObject_SetAttr(__pyx_m, __pyx_n_s__time, __pyx_t_2) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

  /* "Metropolis.pyx":13

 * from random import random

 * import time

 * import getopt             # <<<<<<<<<<<<<<

 * import matplotlib.pyplot as plt

 * import Metropolis

 */

  __pyx_t_2 = __Pyx_Import(((PyObject *)__pyx_n_s__getopt), 0, -1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 13; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_2);

  if (PyObject_SetAttr(__pyx_m, __pyx_n_s__getopt, __pyx_t_2) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 13; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;

  /* "Metropolis.pyx":14

 * import time

 * import getopt

 * import matplotlib.pyplot as plt             # <<<<<<<<<<<<<<

 * import Metropolis

 * 

 */

  __pyx_t_2 = PyList_New(1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 14; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(((PyObject *)__pyx_t_2));

  __Pyx_INCREF(((PyObject *)__pyx_n_s_2));

  PyList_SET_ITEM(__pyx_t_2, 0, ((PyObject *)__pyx_n_s_2));

  __Pyx_GIVEREF(((PyObject *)__pyx_n_s_2));

  __pyx_t_1 = __Pyx_Import(((PyObject *)__pyx_n_s_1), ((PyObject *)__pyx_t_2), -1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 14; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0;

  if (PyObject_SetAttr(__pyx_m, __pyx_n_s__plt, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 14; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  /* "Metropolis.pyx":15

 * import getopt

 * import matplotlib.pyplot as plt

 * import Metropolis             # <<<<<<<<<<<<<<

 * 

 * def Metropolis(sigma,J,B,T,iterations):

 */

  __pyx_t_1 = __Pyx_Import(((PyObject *)__pyx_n_s__Metropolis), 0, -1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 15; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  if (PyObject_SetAttr(__pyx_m, __pyx_n_s__Metropolis, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 15; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  /* "Metropolis.pyx":17

 * import Metropolis

 * 

 * def Metropolis(sigma,J,B,T,iterations):             # <<<<<<<<<<<<<<

 *     start=time.time()

 * 

 */

  __pyx_t_1 = PyCFunction_NewEx(&__pyx_mdef_10Metropolis_Metropolis, NULL, __pyx_n_s__Metropolis); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(__pyx_t_1);

  if (PyObject_SetAttr(__pyx_m, __pyx_n_s__Metropolis, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

  /* "Metropolis.pyx":1

 * #!/usr/bin/env python             # <<<<<<<<<<<<<<

 * #

 * # Ising2D model in serial mode

 */

  __pyx_t_1 = PyDict_New(); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_GOTREF(((PyObject *)__pyx_t_1));

  if (PyObject_SetAttr(__pyx_m, __pyx_n_s____test__, ((PyObject *)__pyx_t_1)) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}

  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;

  goto __pyx_L0;

  __pyx_L1_error:;

  __Pyx_XDECREF(__pyx_t_1);

  __Pyx_XDECREF(__pyx_t_2);

  if (__pyx_m) {

    __Pyx_AddTraceback("init Metropolis", __pyx_clineno, __pyx_lineno, __pyx_filename);

    Py_DECREF(__pyx_m); __pyx_m = 0;

  } else if (!PyErr_Occurred()) {

    PyErr_SetString(PyExc_ImportError, "init Metropolis");

  }

  __pyx_L0:;

  __Pyx_RefNannyFinishContext();

  #if PY_MAJOR_VERSION < 3

  return;

  #else

  return __pyx_m;

  #endif

}

/* Runtime support code */

#if CYTHON_REFNANNY

static __Pyx_RefNannyAPIStruct *__Pyx_RefNannyImportAPI(const char *modname) {

    PyObject *m = NULL, *p = NULL;

    void *r = NULL;

    m = PyImport_ImportModule((char *)modname);

    if (!m) goto end;

    p = PyObject_GetAttrString(m, (char *)"RefNannyAPI");

    if (!p) goto end;

    r = PyLong_AsVoidPtr(p);

end:

    Py_XDECREF(p);

    Py_XDECREF(m);

    return (__Pyx_RefNannyAPIStruct *)r;

}

#endif /* CYTHON_REFNANNY */

static PyObject *__Pyx_GetName(PyObject *dict, PyObject *name) {

    PyObject *result;

    result = PyObject_GetAttr(dict, name);

    if (!result) {

        if (dict != __pyx_b) {

            PyErr_Clear();

            result = PyObject_GetAttr(__pyx_b, name);

        }

        if (!result) {

            PyErr_SetObject(PyExc_NameError, name);

        }

    }

    return result;

}

static void __Pyx_RaiseArgtupleInvalid(

    const char* func_name,

    int exact,

    Py_ssize_t num_min,

    Py_ssize_t num_max,

    Py_ssize_t num_found)

{

    Py_ssize_t num_expected;

    const char *more_or_less;

    if (num_found < num_min) {

        num_expected = num_min;

        more_or_less = "at least";

    } else {

        num_expected = num_max;

        more_or_less = "at most";

    }

    if (exact) {

        more_or_less = "exactly";

    }

    PyErr_Format(PyExc_TypeError,

                 "%s() takes %s %"PY_FORMAT_SIZE_T"d positional argument%s (%"PY_FORMAT_SIZE_T"d given)",

                 func_name, more_or_less, num_expected,

                 (num_expected == 1) ? "" : "s", num_found);

}

static void __Pyx_RaiseDoubleKeywordsError(

    const char* func_name,

    PyObject* kw_name)

{

    PyErr_Format(PyExc_TypeError,

        #if PY_MAJOR_VERSION >= 3

        "%s() got multiple values for keyword argument '%U'", func_name, kw_name);

        #else

        "%s() got multiple values for keyword argument '%s'", func_name,

        PyString_AS_STRING(kw_name));

        #endif

}

static int __Pyx_ParseOptionalKeywords(

    PyObject *kwds,

    PyObject **argnames[],

    PyObject *kwds2,

    PyObject *values[],

    Py_ssize_t num_pos_args,

    const char* function_name)

{