/* Generated by Cython 0.16 on Wed Apr 25 13:44:33 2012 */

#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.
#elif PY_VERSION_HEX < 0x02040000
    #error Cython requires Python 2.4+.
#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

#ifndef Py_HUGE_VAL
  #define Py_HUGE_VAL HUGE_VAL
#endif

#ifdef PYPY_VERSION
#define CYTHON_COMPILING_IN_PYPY 1
#define CYTHON_COMPILING_IN_CPYTHON 0
#else
#define CYTHON_COMPILING_IN_PYPY 0
#define CYTHON_COMPILING_IN_CPYTHON 1
#endif

#if CYTHON_COMPILING_IN_PYPY
  #define __Pyx_PyCFunction_Call PyObject_Call
#else
  #define __Pyx_PyCFunction_Call PyCFunction_Call
#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)
  #define __PYX_BUILD_PY_SSIZE_T "i"
#else
  #define __PYX_BUILD_PY_SSIZE_T "n"
#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)
  #define PyBUF_RECORDS (PyBUF_STRIDES | PyBUF_FORMAT | PyBUF_WRITABLE)
  #define PyBUF_FULL (PyBUF_INDIRECT | PyBUF_FORMAT | PyBUF_WRITABLE)

  typedef int (*getbufferproc)(PyObject *, Py_buffer *, int);
  typedef void (*releasebufferproc)(PyObject *, Py_buffer *);
#endif

#if PY_MAJOR_VERSION < 3
  #define __Pyx_BUILTIN_MODULE_NAME "__builtin__"
  #define __Pyx_PyCode_New(a, k, l, s, f, code, c, n, v, fv, cell, fn, name, fline, lnos) \
          PyCode_New(a, l, s, f, code, c, n, v, fv, cell, fn, name, fline, lnos)
#else
  #define __Pyx_BUILTIN_MODULE_NAME "builtins"
  #define __Pyx_PyCode_New(a, k, l, s, f, code, c, n, v, fv, cell, fn, name, fline, lnos) \
          PyCode_New(a, k, l, s, f, code, c, n, v, fv, cell, fn, name, fline, lnos)
#endif

#if PY_MAJOR_VERSION < 3 && PY_MINOR_VERSION < 6
  #define PyUnicode_FromString(s) PyUnicode_Decode(s, strlen(s), "UTF-8", "strict")
#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_VERSION_HEX > 0x03030000 && defined(PyUnicode_GET_LENGTH)
  #define CYTHON_PEP393_ENABLED 1
  #define __Pyx_PyUnicode_GET_LENGTH(u) PyUnicode_GET_LENGTH(u)
  #define __Pyx_PyUnicode_READ_CHAR(u, i) PyUnicode_READ_CHAR(u, i)
#else
  #define CYTHON_PEP393_ENABLED 0
  #define __Pyx_PyUnicode_GET_LENGTH(u) PyUnicode_GET_SIZE(u)
  #define __Pyx_PyUnicode_READ_CHAR(u, i) ((Py_UCS4)(PyUnicode_AS_UNICODE(u)[i]))
#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) || (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

#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

#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__obitools__align___gprofilenws
#define __PYX_HAVE_API__obitools__align___gprofilenws
#include "stdlib.h"
#include "string.h"
#include "math.h"
#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) && !defined(_MSC_VER))
#   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))
#define __pyx_PyFloat_AsFloat(x) ((float) __pyx_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[] = {
  "_gprofilenws.pyx",
  "_dynamic.pxd",
  "_nws.pxd",
  "_profile.pxd",
};

/*--- Type declarations ---*/
struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming;
struct __pyx_obj_8obitools_5align_4_nws_NWS;
struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile;
struct __pyx_obj_8obitools_5align_11_profilenws_ProfileNWS;
struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray;
struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS;
struct __pyx_t_8obitools_5align_8_dynamic_AlignCell;
struct __pyx_t_8obitools_5align_8_dynamic_AlignMatrix;
struct __pyx_opt_args_8obitools_5align_8_dynamic_allocateMatrix;
struct __pyx_t_8obitools_5align_8_dynamic_alignSequence;
struct __pyx_opt_args_8obitools_5align_8_dynamic_allocateSequence;
struct __pyx_t_8obitools_5align_8_dynamic_alignPath;
struct __pyx_opt_args_8obitools_5align_8_dynamic_allocatePath;

/* "_dynamic.pxd":11
 *     void memcpy(void* s1, void* s2, int n)
 * 
 * cdef struct AlignCell :             # <<<<<<<<<<<<<<
 *     double score
 *     int   path
 */
struct __pyx_t_8obitools_5align_8_dynamic_AlignCell {
  double score;
  int path;
};

/* "_dynamic.pxd":15
 *     int   path
 * 
 * cdef struct AlignMatrix :             # <<<<<<<<<<<<<<
 *     AlignCell*  matrix
 *     int*        bestVJump
 */
struct __pyx_t_8obitools_5align_8_dynamic_AlignMatrix {
  struct __pyx_t_8obitools_5align_8_dynamic_AlignCell *matrix;
  int *bestVJump;
  int *bestHJump;
  int msize;
  int vsize;
  int hsize;
};

/* "_dynamic.pxd":25
 * 
 * 
 * cdef AlignMatrix* allocateMatrix(int hsize, int vsize,AlignMatrix *matrix=?)             # <<<<<<<<<<<<<<
 * 
 * cdef void freeMatrix(AlignMatrix* matrix)
 */
struct __pyx_opt_args_8obitools_5align_8_dynamic_allocateMatrix {
  int __pyx_n;
  struct __pyx_t_8obitools_5align_8_dynamic_AlignMatrix *matrix;
};

/* "_dynamic.pxd":32
 * 
 * 
 * cdef struct alignSequence:             # <<<<<<<<<<<<<<
 *     long    length
 *     long    buffsize
 */
struct __pyx_t_8obitools_5align_8_dynamic_alignSequence {
  long length;
  long buffsize;
  int hasQuality;
  char *sequence;
  double *quality;
};

/* "_dynamic.pxd":39
 *     double* quality
 * 
 * cdef alignSequence* allocateSequence(object bioseq, alignSequence* seq=?) except *             # <<<<<<<<<<<<<<
 * 
 * cdef void freeSequence(alignSequence* seq)
 */
struct __pyx_opt_args_8obitools_5align_8_dynamic_allocateSequence {
  int __pyx_n;
  struct __pyx_t_8obitools_5align_8_dynamic_alignSequence *seq;
};

/* "_dynamic.pxd":43
 * cdef void freeSequence(alignSequence* seq)
 * 
 * cdef struct alignPath:             # <<<<<<<<<<<<<<
 *     long length
 *     long buffsize
 */
struct __pyx_t_8obitools_5align_8_dynamic_alignPath {
  long length;
  long buffsize;
  long vStart;
  long hStart;
  long *path;
};

/* "_dynamic.pxd":50
 *     long *path
 * 
 * cdef alignPath* allocatePath(long l1,long l2,alignPath* path=?)             # <<<<<<<<<<<<<<
 * 
 * cdef void reversePath(alignPath* path)
 */
struct __pyx_opt_args_8obitools_5align_8_dynamic_allocatePath {
  int __pyx_n;
  struct __pyx_t_8obitools_5align_8_dynamic_alignPath *path;
};
struct __pyx_t_8obitools_7profile_8_profile_dnaprofile_t;
struct __pyx_opt_args_8obitools_7profile_8_profile_12_MemIntArray_frequency;

/* "obitools/profile/_profile.pxd":20
 * 
 * 
 * cdef struct dnaprofile_t:             # <<<<<<<<<<<<<<
 *     int     length
 *     int     weight
 */
struct __pyx_t_8obitools_7profile_8_profile_dnaprofile_t {
  int length;
  int weight;
  int value;
  double pseudo;
  int* A;
  int* C;
  int* G;
  int* T;
  int* Og;
  int* Eg;
};

/* "obitools/profile/_profile.pxd":44
 *     cdef initialize(self, int_p start,int size)
 *     cdef int normalize(self, int pos)
 *     cpdef double frequency(self,int pos, int weight, double pseudo=?)             # <<<<<<<<<<<<<<
 * 
 * 
 */
struct __pyx_opt_args_8obitools_7profile_8_profile_12_MemIntArray_frequency {
  int __pyx_n;
  double pseudo;
};
struct __pyx_t_8obitools_5align_11_profilenws_alignProfile;
struct __pyx_opt_args_8obitools_5align_11_profilenws_allocateProfile;

/* "_profilenws.pxd":4
 * from obitools.profile._profile cimport *
 * 
 * cdef struct alignProfile:             # <<<<<<<<<<<<<<
 *     long    length
 *     long    buffsize
 */
struct __pyx_t_8obitools_5align_11_profilenws_alignProfile {
  long length;
  long buffsize;
  double *frequency;
};

/* "_profilenws.pxd":9
 *     double* frequency
 * 
 * cdef alignProfile* allocateProfile(object profile, alignProfile* prof=?)             # <<<<<<<<<<<<<<
 * 
 * cdef void freeProfile(alignProfile* prof)
 */
struct __pyx_opt_args_8obitools_5align_11_profilenws_allocateProfile {
  int __pyx_n;
  struct __pyx_t_8obitools_5align_11_profilenws_alignProfile *prof;
};

/* "_dynamic.pxd":63
 * cpdef unsigned char encodeBase(unsigned char lettre)
 * 
 * cdef class DynamicProgramming:             # <<<<<<<<<<<<<<
 *     cdef AlignMatrix* matrix
 * 
 */
struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming {
  PyObject_HEAD
  struct __pyx_vtabstruct_8obitools_5align_8_dynamic_DynamicProgramming *__pyx_vtab;
  struct __pyx_t_8obitools_5align_8_dynamic_AlignMatrix *matrix;
  PyObject *horizontalSeq;
  PyObject *verticalSeq;
  struct __pyx_t_8obitools_5align_8_dynamic_alignSequence *hSeq;
  struct __pyx_t_8obitools_5align_8_dynamic_alignSequence *vSeq;
  struct __pyx_t_8obitools_5align_8_dynamic_alignPath *path;
  double _opengap;
  double _extgap;
  PyObject *alignment;
  int sequenceChanged;
  int scoreChanged;
};


/* "_nws.pxd":3
 * from _dynamic cimport *
 * 
 * cdef class NWS(DynamicProgramming):             # <<<<<<<<<<<<<<
 *     cdef double _match
 *     cdef double _mismatch
 */
struct __pyx_obj_8obitools_5align_4_nws_NWS {
  struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming __pyx_base;
  double _match;
  double _mismatch;
};


/* "obitools/profile/_profile.pxd":47
 * 
 * 
 * cdef class DNAProfile:             # <<<<<<<<<<<<<<
 * 
 *     cdef dnaprofile_t* profile
 */
struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile {
  PyObject_HEAD
  struct __pyx_vtabstruct_8obitools_7profile_8_profile_DNAProfile *__pyx_vtab;
  struct __pyx_t_8obitools_7profile_8_profile_dnaprofile_t *profile;
  struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray *_A;
  struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray *_C;
  struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray *_G;
  struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray *_T;
  struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray *_Og;
  struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray *_Eg;
};


/* "_profilenws.pxd":13
 * cdef void freeProfile(alignProfile* prof)
 * 
 * cdef class ProfileNWS(NWS):             # <<<<<<<<<<<<<<
 * 
 *     cdef alignProfile* hProf
 */
struct __pyx_obj_8obitools_5align_11_profilenws_ProfileNWS {
  struct __pyx_obj_8obitools_5align_4_nws_NWS __pyx_base;
  struct __pyx_t_8obitools_5align_11_profilenws_alignProfile *hProf;
  struct __pyx_t_8obitools_5align_11_profilenws_alignProfile *vProf;
};


/* "obitools/profile/_profile.pxd":37
 * cdef void copyDNAProfile(dnaprofile_t* dest, dnaprofile_t* source)
 * 
 * cdef class _MemIntArray:             # <<<<<<<<<<<<<<
 * 
 *     cdef    int_p   start
 */
struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray {
  PyObject_HEAD
  struct __pyx_vtabstruct_8obitools_7profile_8_profile__MemIntArray *__pyx_vtab;
  int* start;
  int size;
};


/* "obitools/align/_gprofilenws.pxd":3
 * from _profilenws cimport *
 * 
 * cdef class GProfileNWS(ProfileNWS):             # <<<<<<<<<<<<<<
 * 
 *     cdef double matchScore(self,int h, int v)
 */
struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS {
  struct __pyx_obj_8obitools_5align_11_profilenws_ProfileNWS __pyx_base;
  PyObject *alignment1;
  PyObject *alignment2;
};



/* "_dynamic.pxd":63
 * cpdef unsigned char encodeBase(unsigned char lettre)
 * 
 * cdef class DynamicProgramming:             # <<<<<<<<<<<<<<
 *     cdef AlignMatrix* matrix
 * 
 */

struct __pyx_vtabstruct_8obitools_5align_8_dynamic_DynamicProgramming {
  int (*_vlen)(struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *);
  int (*_hlen)(struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *);
  int (*allocate)(struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *);
  double (*doAlignment)(struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *);
  void (*reset)(struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *);
  int (*index)(struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *, int, int);
  int (*_needToCompute)(struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *);
  void (*backtrack)(struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *);
  void (*clean)(struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *);
};
static struct __pyx_vtabstruct_8obitools_5align_8_dynamic_DynamicProgramming *__pyx_vtabptr_8obitools_5align_8_dynamic_DynamicProgramming;


/* "_nws.pxd":3
 * from _dynamic cimport *
 * 
 * cdef class NWS(DynamicProgramming):             # <<<<<<<<<<<<<<
 *     cdef double _match
 *     cdef double _mismatch
 */

struct __pyx_vtabstruct_8obitools_5align_4_nws_NWS {
  struct __pyx_vtabstruct_8obitools_5align_8_dynamic_DynamicProgramming __pyx_base;
  double (*matchScore)(struct __pyx_obj_8obitools_5align_4_nws_NWS *, int, int);
};
static struct __pyx_vtabstruct_8obitools_5align_4_nws_NWS *__pyx_vtabptr_8obitools_5align_4_nws_NWS;


/* "_profilenws.pxd":13
 * cdef void freeProfile(alignProfile* prof)
 * 
 * cdef class ProfileNWS(NWS):             # <<<<<<<<<<<<<<
 * 
 *     cdef alignProfile* hProf
 */

struct __pyx_vtabstruct_8obitools_5align_11_profilenws_ProfileNWS {
  struct __pyx_vtabstruct_8obitools_5align_4_nws_NWS __pyx_base;
};
static struct __pyx_vtabstruct_8obitools_5align_11_profilenws_ProfileNWS *__pyx_vtabptr_8obitools_5align_11_profilenws_ProfileNWS;


/* "obitools/align/_gprofilenws.pyx":11
 * 
 * 
 * cdef class GProfileNWS(ProfileNWS):             # <<<<<<<<<<<<<<
 * 
 *     cdef double matchScore(self,int h, int v):
 */

struct __pyx_vtabstruct_8obitools_5align_12_gprofilenws_GProfileNWS {
  struct __pyx_vtabstruct_8obitools_5align_11_profilenws_ProfileNWS __pyx_base;
};
static struct __pyx_vtabstruct_8obitools_5align_12_gprofilenws_GProfileNWS *__pyx_vtabptr_8obitools_5align_12_gprofilenws_GProfileNWS;


/* "obitools/profile/_profile.pxd":47
 * 
 * 
 * cdef class DNAProfile:             # <<<<<<<<<<<<<<
 * 
 *     cdef dnaprofile_t* profile
 */

struct __pyx_vtabstruct_8obitools_7profile_8_profile_DNAProfile {
  int (*equal)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, int __pyx_skip_dispatch);
  struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *(*add)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, int __pyx_skip_dispatch);
  double (*lproba)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, int __pyx_skip_dispatch);
  double (*proba)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, int __pyx_skip_dispatch);
  void (*_initLetter)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *);
  void (*_initFromString)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, char *);
  double (*fA)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, int, int __pyx_skip_dispatch);
  double (*fC)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, int, int __pyx_skip_dispatch);
  double (*fG)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, int, int __pyx_skip_dispatch);
  double (*fT)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, int, int __pyx_skip_dispatch);
  double (*fOg)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, int, int __pyx_skip_dispatch);
  double (*fEg)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, int, int __pyx_skip_dispatch);
  int (*WP)(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *, int __pyx_skip_dispatch);
};
static struct __pyx_vtabstruct_8obitools_7profile_8_profile_DNAProfile *__pyx_vtabptr_8obitools_7profile_8_profile_DNAProfile;


/* "obitools/profile/_profile.pxd":37
 * cdef void copyDNAProfile(dnaprofile_t* dest, dnaprofile_t* source)
 * 
 * cdef class _MemIntArray:             # <<<<<<<<<<<<<<
 * 
 *     cdef    int_p   start
 */

struct __pyx_vtabstruct_8obitools_7profile_8_profile__MemIntArray {
  PyObject *(*initialize)(struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray *, int*, int);
  int (*normalize)(struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray *, int);
  double (*frequency)(struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray *, int, int, int __pyx_skip_dispatch, struct __pyx_opt_args_8obitools_7profile_8_profile_12_MemIntArray_frequency *__pyx_optional_args);
};
static struct __pyx_vtabstruct_8obitools_7profile_8_profile__MemIntArray *__pyx_vtabptr_8obitools_7profile_8_profile__MemIntArray;
#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;
#ifdef WITH_THREAD
  #define __Pyx_RefNannySetupContext(name, acquire_gil) \
          if (acquire_gil) { \
              PyGILState_STATE __pyx_gilstate_save = PyGILState_Ensure(); \
              __pyx_refnanny = __Pyx_RefNanny->SetupContext((name), __LINE__, __FILE__); \
              PyGILState_Release(__pyx_gilstate_save); \
          } else { \
              __pyx_refnanny = __Pyx_RefNanny->SetupContext((name), __LINE__, __FILE__); \
          }
#else
  #define __Pyx_RefNannySetupContext(name, acquire_gil) \
          __pyx_refnanny = __Pyx_RefNanny->SetupContext((name), __LINE__, __FILE__)
#endif
  #define __Pyx_RefNannyFinishContext() \
          __Pyx_RefNanny->FinishContext(&__pyx_refnanny)
  #define __Pyx_INCREF(r)  __Pyx_RefNanny->INCREF(__pyx_refnanny, (PyObject *)(r), __LINE__)
  #define __Pyx_DECREF(r)  __Pyx_RefNanny->DECREF(__pyx_refnanny, (PyObject *)(r), __LINE__)
  #define __Pyx_GOTREF(r)  __Pyx_RefNanny->GOTREF(__pyx_refnanny, (PyObject *)(r), __LINE__)
  #define __Pyx_GIVEREF(r) __Pyx_RefNanny->GIVEREF(__pyx_refnanny, (PyObject *)(r), __LINE__)
  #define __Pyx_XINCREF(r)  do { if((r) != NULL) {__Pyx_INCREF(r); }} while(0)
  #define __Pyx_XDECREF(r)  do { if((r) != NULL) {__Pyx_DECREF(r); }} while(0)
  #define __Pyx_XGOTREF(r)  do { if((r) != NULL) {__Pyx_GOTREF(r); }} while(0)
  #define __Pyx_XGIVEREF(r) do { if((r) != NULL) {__Pyx_GIVEREF(r);}} while(0)
#else
  #define __Pyx_RefNannyDeclarations
  #define __Pyx_RefNannySetupContext(name, acquire_gil)
  #define __Pyx_RefNannyFinishContext()
  #define __Pyx_INCREF(r) Py_INCREF(r)
  #define __Pyx_DECREF(r) Py_DECREF(r)
  #define __Pyx_GOTREF(r)
  #define __Pyx_GIVEREF(r)
  #define __Pyx_XINCREF(r) Py_XINCREF(r)
  #define __Pyx_XDECREF(r) Py_XDECREF(r)
  #define __Pyx_XGOTREF(r)
  #define __Pyx_XGIVEREF(r)
#endif /* CYTHON_REFNANNY */
#define __Pyx_CLEAR(r)    do { PyObject* tmp = ((PyObject*)(r)); r = NULL; __Pyx_DECREF(tmp);} while(0)
#define __Pyx_XCLEAR(r)   do { if((r) != NULL) {PyObject* tmp = ((PyObject*)(r)); r = NULL; __Pyx_DECREF(tmp);}} while(0)

static PyObject *__Pyx_GetName(PyObject *dict, PyObject *name); /*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 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 CYTHON_INLINE int __Pyx_TypeTest(PyObject *obj, PyTypeObject *type); /*proto*/

#if HAVE_LONG_LONG && defined (__STDC_VERSION__) && __STDC_VERSION__ >= 199901L
#define __Pyx_abs_int(x)     ((sizeof(x) <= sizeof(int)) ? ((unsigned int)abs(x)) :      ((sizeof(x) <= sizeof(long)) ? ((unsigned long)labs(x)) :       ((unsigned PY_LONG_LONG)llabs(x))))
#else
#define __Pyx_abs_int(x)     ((sizeof(x) <= sizeof(int)) ? ((unsigned int)abs(x)) : ((unsigned long)labs(x)))
#endif
#define __Pyx_abs_long(x) __Pyx_abs_int(x)

static CYTHON_INLINE PyObject *__Pyx_GetItemInt_Generic(PyObject *o, PyObject* j) {
    PyObject *r;
    if (!j) return NULL;
    r = PyObject_GetItem(o, j);
    Py_DECREF(j);
    return r;
}
#define __Pyx_GetItemInt_List(o, i, size, to_py_func) (((size) <= sizeof(Py_ssize_t)) ? \
                                                    __Pyx_GetItemInt_List_Fast(o, i) : \
                                                    __Pyx_GetItemInt_Generic(o, to_py_func(i)))
static CYTHON_INLINE PyObject *__Pyx_GetItemInt_List_Fast(PyObject *o, Py_ssize_t i) {
    if (likely(o != Py_None)) {
        if (likely((0 <= i) & (i < PyList_GET_SIZE(o)))) {
            PyObject *r = PyList_GET_ITEM(o, i);
            Py_INCREF(r);
            return r;
        }
        else if ((-PyList_GET_SIZE(o) <= i) & (i < 0)) {
            PyObject *r = PyList_GET_ITEM(o, PyList_GET_SIZE(o) + i);
            Py_INCREF(r);
            return r;
        }
    }
    return __Pyx_GetItemInt_Generic(o, PyInt_FromSsize_t(i));
}
#define __Pyx_GetItemInt_Tuple(o, i, size, to_py_func) (((size) <= sizeof(Py_ssize_t)) ? \
                                                    __Pyx_GetItemInt_Tuple_Fast(o, i) : \
                                                    __Pyx_GetItemInt_Generic(o, to_py_func(i)))
static CYTHON_INLINE PyObject *__Pyx_GetItemInt_Tuple_Fast(PyObject *o, Py_ssize_t i) {
    if (likely(o != Py_None)) {
        if (likely((0 <= i) & (i < PyTuple_GET_SIZE(o)))) {
            PyObject *r = PyTuple_GET_ITEM(o, i);
            Py_INCREF(r);
            return r;
        }
        else if ((-PyTuple_GET_SIZE(o) <= i) & (i < 0)) {
            PyObject *r = PyTuple_GET_ITEM(o, PyTuple_GET_SIZE(o) + i);
            Py_INCREF(r);
            return r;
        }
    }
    return __Pyx_GetItemInt_Generic(o, PyInt_FromSsize_t(i));
}
#define __Pyx_GetItemInt(o, i, size, to_py_func) (((size) <= sizeof(Py_ssize_t)) ? \
                                                    __Pyx_GetItemInt_Fast(o, i) : \
                                                    __Pyx_GetItemInt_Generic(o, to_py_func(i)))
static CYTHON_INLINE PyObject *__Pyx_GetItemInt_Fast(PyObject *o, Py_ssize_t i) {
    if (PyList_CheckExact(o)) {
        Py_ssize_t n = (likely(i >= 0)) ? i : i + PyList_GET_SIZE(o);
        if (likely((n >= 0) & (n < PyList_GET_SIZE(o)))) {
            PyObject *r = PyList_GET_ITEM(o, n);
            Py_INCREF(r);
            return r;
        }
    }
    else if (PyTuple_CheckExact(o)) {
        Py_ssize_t n = (likely(i >= 0)) ? i : i + PyTuple_GET_SIZE(o);
        if (likely((n >= 0) & (n < PyTuple_GET_SIZE(o)))) {
            PyObject *r = PyTuple_GET_ITEM(o, n);
            Py_INCREF(r);
            return r;
        }
    }
    else if (likely(i >= 0)) {
        PySequenceMethods *m = Py_TYPE(o)->tp_as_sequence;
        if (likely(m && m->sq_item)) {
            return m->sq_item(o, i);
        }
    }
    return __Pyx_GetItemInt_Generic(o, PyInt_FromSsize_t(i));
}

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

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

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

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

static void* __Pyx_GetVtable(PyObject *dict); /*proto*/

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

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

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

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


/* Module declarations from 'obitools.align._dynamic' */
static PyTypeObject *__pyx_ptype_8obitools_5align_8_dynamic_DynamicProgramming = 0;

/* Module declarations from 'obitools.align._nws' */
static PyTypeObject *__pyx_ptype_8obitools_5align_4_nws_NWS = 0;

/* Module declarations from 'obitools.profile._profile' */
static PyTypeObject *__pyx_ptype_8obitools_7profile_8_profile__MemIntArray = 0;
static PyTypeObject *__pyx_ptype_8obitools_7profile_8_profile_DNAProfile = 0;

/* Module declarations from 'obitools.align._profilenws' */
static PyTypeObject *__pyx_ptype_8obitools_5align_11_profilenws_ProfileNWS = 0;

/* Module declarations from 'obitools.align._gprofilenws' */
static PyTypeObject *__pyx_ptype_8obitools_5align_12_gprofilenws_GProfileNWS = 0;
#define __Pyx_MODULE_NAME "obitools.align._gprofilenws"
int __pyx_module_is_main_obitools__align___gprofilenws = 0;

/* Implementation of 'obitools.align._gprofilenws' */
static PyObject *__pyx_builtin_xrange;
static PyObject *__pyx_builtin_range;
static PyObject *__pyx_pf_8obitools_5align_12_gprofilenws_11GProfileNWS___call__(struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *__pyx_v_self, PyObject *__pyx_v_pseudocounts); /* proto */
static char __pyx_k_1[] = "\nCreated on 16 Feb. 2011\n\n@author: celine\n";
static char __pyx_k__A[] = "A";
static char __pyx_k__C[] = "C";
static char __pyx_k__G[] = "G";
static char __pyx_k__T[] = "T";
static char __pyx_k__Eg[] = "Eg";
static char __pyx_k__Og[] = "Og";
static char __pyx_k__size[] = "size";
static char __pyx_k__range[] = "range";
static char __pyx_k__pseudo[] = "pseudo";
static char __pyx_k__xrange[] = "xrange";
static char __pyx_k____main__[] = "__main__";
static char __pyx_k____test__[] = "__test__";
static char __pyx_k__pseudocounts[] = "pseudocounts";
static PyObject *__pyx_n_s__A;
static PyObject *__pyx_n_s__C;
static PyObject *__pyx_n_s__Eg;
static PyObject *__pyx_n_s__G;
static PyObject *__pyx_n_s__Og;
static PyObject *__pyx_n_s__T;
static PyObject *__pyx_n_s____main__;
static PyObject *__pyx_n_s____test__;
static PyObject *__pyx_n_s__pseudo;
static PyObject *__pyx_n_s__pseudocounts;
static PyObject *__pyx_n_s__range;
static PyObject *__pyx_n_s__size;
static PyObject *__pyx_n_s__xrange;
static PyObject *__pyx_int_0;
static PyObject *__pyx_int_1;

/* "obitools/align/_gprofilenws.pyx":13
 * cdef class GProfileNWS(ProfileNWS):
 * 
 *     cdef double matchScore(self,int h, int v):             # <<<<<<<<<<<<<<
 *         cdef double pmatch
 *         cdef double* hp = self.hProf.frequency
 */

static double __pyx_f_8obitools_5align_12_gprofilenws_11GProfileNWS_matchScore(struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *__pyx_v_self, int __pyx_v_h, int __pyx_v_v) {
  double __pyx_v_pmatch;
  double *__pyx_v_hp;
  double *__pyx_v_vp;
  int __pyx_v_hl;
  int __pyx_v_vl;
  double __pyx_r;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("matchScore", 0);

  /* "obitools/align/_gprofilenws.pyx":15
 *     cdef double matchScore(self,int h, int v):
 *         cdef double pmatch
 *         cdef double* hp = self.hProf.frequency             # <<<<<<<<<<<<<<
 *         cdef double* vp = self.vProf.frequency
 *         cdef int     hl = self.hProf.length
 */
  __pyx_v_hp = __pyx_v_self->__pyx_base.hProf->frequency;

  /* "obitools/align/_gprofilenws.pyx":16
 *         cdef double pmatch
 *         cdef double* hp = self.hProf.frequency
 *         cdef double* vp = self.vProf.frequency             # <<<<<<<<<<<<<<
 *         cdef int     hl = self.hProf.length
 *         cdef int     vl = self.vProf.length
 */
  __pyx_v_vp = __pyx_v_self->__pyx_base.vProf->frequency;

  /* "obitools/align/_gprofilenws.pyx":17
 *         cdef double* hp = self.hProf.frequency
 *         cdef double* vp = self.vProf.frequency
 *         cdef int     hl = self.hProf.length             # <<<<<<<<<<<<<<
 *         cdef int     vl = self.vProf.length
 * 
 */
  __pyx_v_hl = __pyx_v_self->__pyx_base.hProf->length;

  /* "obitools/align/_gprofilenws.pyx":18
 *         cdef double* vp = self.vProf.frequency
 *         cdef int     hl = self.hProf.length
 *         cdef int     vl = self.vProf.length             # <<<<<<<<<<<<<<
 * 
 *         h-=1
 */
  __pyx_v_vl = __pyx_v_self->__pyx_base.vProf->length;

  /* "obitools/align/_gprofilenws.pyx":20
 *         cdef int     vl = self.vProf.length
 * 
 *         h-=1             # <<<<<<<<<<<<<<
 *         v-=1
 *         pmatch =  hp[h]*vp[v] + \
 */
  __pyx_v_h = (__pyx_v_h - 1);

  /* "obitools/align/_gprofilenws.pyx":21
 * 
 *         h-=1
 *         v-=1             # <<<<<<<<<<<<<<
 *         pmatch =  hp[h]*vp[v] + \
 *                   hp[h+hl]*vp[v+vl] + \
 */
  __pyx_v_v = (__pyx_v_v - 1);

  /* "obitools/align/_gprofilenws.pyx":27
 *                   hp[h+3*hl]*vp[v+3*vl] + \
 *                   hp[h+4*hl]*vp[v+4*vl] + \
 *                   hp[h+5*hl]*vp[v+5*vl]             # <<<<<<<<<<<<<<
 *         return self._match * pmatch + (1-pmatch) * self._mismatch
 * 
 */
  __pyx_v_pmatch = (((((((__pyx_v_hp[__pyx_v_h]) * (__pyx_v_vp[__pyx_v_v])) + ((__pyx_v_hp[(__pyx_v_h + __pyx_v_hl)]) * (__pyx_v_vp[(__pyx_v_v + __pyx_v_vl)]))) + ((__pyx_v_hp[(__pyx_v_h + (2 * __pyx_v_hl))]) * (__pyx_v_vp[(__pyx_v_v + (2 * __pyx_v_vl))]))) + ((__pyx_v_hp[(__pyx_v_h + (3 * __pyx_v_hl))]) * (__pyx_v_vp[(__pyx_v_v + (3 * __pyx_v_vl))]))) + ((__pyx_v_hp[(__pyx_v_h + (4 * __pyx_v_hl))]) * (__pyx_v_vp[(__pyx_v_v + (4 * __pyx_v_vl))]))) + ((__pyx_v_hp[(__pyx_v_h + (5 * __pyx_v_hl))]) * (__pyx_v_vp[(__pyx_v_v + (5 * __pyx_v_vl))])));

  /* "obitools/align/_gprofilenws.pyx":28
 *                   hp[h+4*hl]*vp[v+4*vl] + \
 *                   hp[h+5*hl]*vp[v+5*vl]
 *         return self._match * pmatch + (1-pmatch) * self._mismatch             # <<<<<<<<<<<<<<
 * 
 * 
 */
  __pyx_r = ((__pyx_v_self->__pyx_base.__pyx_base._match * __pyx_v_pmatch) + ((1.0 - __pyx_v_pmatch) * __pyx_v_self->__pyx_base.__pyx_base._mismatch));
  goto __pyx_L0;

  __pyx_r = 0;
  __pyx_L0:;
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

/* Python wrapper */
static PyObject *__pyx_pw_8obitools_5align_12_gprofilenws_11GProfileNWS_1__call__(PyObject *__pyx_v_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/
static PyObject *__pyx_pw_8obitools_5align_12_gprofilenws_11GProfileNWS_1__call__(PyObject *__pyx_v_self, PyObject *__pyx_args, PyObject *__pyx_kwds) {
  PyObject *__pyx_v_pseudocounts = 0;
  static PyObject **__pyx_pyargnames[] = {&__pyx_n_s__pseudocounts,0};
  PyObject *__pyx_r = 0;
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("__call__ (wrapper)", 0);
  {
    PyObject* values[1] = {0};
    values[0] = ((PyObject *)__pyx_int_0);
    if (unlikely(__pyx_kwds)) {
      Py_ssize_t kw_args;
      const Py_ssize_t pos_args = PyTuple_GET_SIZE(__pyx_args);
      switch (pos_args) {
        case  1: values[0] = PyTuple_GET_ITEM(__pyx_args, 0);
        case  0: break;
        default: goto __pyx_L5_argtuple_error;
      }
      kw_args = PyDict_Size(__pyx_kwds);
      switch (pos_args) {
        case  0:
        if (kw_args > 0) {
          PyObject* value = PyDict_GetItem(__pyx_kwds, __pyx_n_s__pseudocounts);
          if (value) { values[0] = value; kw_args--; }
        }
      }
      if (unlikely(kw_args > 0)) {
        if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_pyargnames, 0, values, pos_args, "__call__") < 0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
      }
    } else {
      switch (PyTuple_GET_SIZE(__pyx_args)) {
        case  1: values[0] = PyTuple_GET_ITEM(__pyx_args, 0);
        case  0: break;
        default: goto __pyx_L5_argtuple_error;
      }
    }
    __pyx_v_pseudocounts = values[0];
  }
  goto __pyx_L4_argument_unpacking_done;
  __pyx_L5_argtuple_error:;
  __Pyx_RaiseArgtupleInvalid("__call__", 0, 0, 1, PyTuple_GET_SIZE(__pyx_args)); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L3_error;}
  __pyx_L3_error:;
  __Pyx_AddTraceback("obitools.align._gprofilenws.GProfileNWS.__call__", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __Pyx_RefNannyFinishContext();
  return NULL;
  __pyx_L4_argument_unpacking_done:;
  __pyx_r = __pyx_pf_8obitools_5align_12_gprofilenws_11GProfileNWS___call__(((struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *)__pyx_v_self), __pyx_v_pseudocounts);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}

/* "obitools/align/_gprofilenws.pyx":31
 * 
 * 
 *     def __call__(self,pseudocounts=0):             # <<<<<<<<<<<<<<
 * 
 *         cdef list hgaps=[]
 */

static PyObject *__pyx_pf_8obitools_5align_12_gprofilenws_11GProfileNWS___call__(struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *__pyx_v_self, PyObject *__pyx_v_pseudocounts) {
  CYTHON_UNUSED PyObject *__pyx_v_hgaps = 0;
  CYTHON_UNUSED PyObject *__pyx_v_vgaps = 0;
  int __pyx_v_hp;
  int __pyx_v_vp;
  CYTHON_UNUSED int __pyx_v_lenh;
  CYTHON_UNUSED int __pyx_v_lenv;
  int __pyx_v_v;
  int __pyx_v_p;
  int __pyx_v_i;
  CYTHON_UNUSED double __pyx_v_score;
  struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *__pyx_v_newProfile1 = 0;
  struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *__pyx_v_newProfile2 = 0;
  struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *__pyx_v_horizontalSeq = 0;
  struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *__pyx_v_verticalSeq = 0;
  PyObject *__pyx_v_sum = NULL;
  PyObject *__pyx_v_rp1 = NULL;
  PyObject *__pyx_v_rp2 = NULL;
  CYTHON_UNUSED PyObject *__pyx_v_x = NULL;
  struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *__pyx_v_ali1 = NULL;
  struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *__pyx_v_ali2 = NULL;
  PyObject *__pyx_r = NULL;
  __Pyx_RefNannyDeclarations
  PyObject *__pyx_t_1 = NULL;
  int __pyx_t_2;
  double __pyx_t_3;
  long __pyx_t_4;
  int __pyx_t_5;
  unsigned int __pyx_t_6;
  PyObject *__pyx_t_7 = NULL;
  Py_ssize_t __pyx_t_8;
  PyObject *(*__pyx_t_9)(PyObject *);
  PyObject *__pyx_t_10 = NULL;
  int __pyx_lineno = 0;
  const char *__pyx_filename = NULL;
  int __pyx_clineno = 0;
  __Pyx_RefNannySetupContext("__call__", 0);

  /* "obitools/align/_gprofilenws.pyx":33
 *     def __call__(self,pseudocounts=0):
 * 
 *         cdef list hgaps=[]             # <<<<<<<<<<<<<<
 *         cdef list vgaps=[]
 *         cdef list b
 */
  __pyx_t_1 = PyList_New(0); 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_v_hgaps = __pyx_t_1;
  __pyx_t_1 = 0;

  /* "obitools/align/_gprofilenws.pyx":34
 * 
 *         cdef list hgaps=[]
 *         cdef list vgaps=[]             # <<<<<<<<<<<<<<
 *         cdef list b
 *         cdef int  hp
 */
  __pyx_t_1 = PyList_New(0); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_1);
  __pyx_v_vgaps = __pyx_t_1;
  __pyx_t_1 = 0;

  /* "obitools/align/_gprofilenws.pyx":39
 *         cdef int  vp
 *         cdef int  rp
 *         cdef int  lenh=0             # <<<<<<<<<<<<<<
 *         cdef int  lenv=0
 *         cdef int  h,v,p
 */
  __pyx_v_lenh = 0;

  /* "obitools/align/_gprofilenws.pyx":40
 *         cdef int  rp
 *         cdef int  lenh=0
 *         cdef int  lenv=0             # <<<<<<<<<<<<<<
 *         cdef int  h,v,p
 *         cdef int  i
 */
  __pyx_v_lenv = 0;

  /* "obitools/align/_gprofilenws.pyx":47
 *         cdef DNAProfile newProfile1
 *         cdef DNAProfile newProfile2
 *         cdef DNAProfile horizontalSeq=self.horizontalSeq             # <<<<<<<<<<<<<<
 *         cdef DNAProfile verticalSeq=self.verticalSeq
 * 
 */
  if (!(likely(((__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.horizontalSeq) == Py_None) || likely(__Pyx_TypeTest(__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.horizontalSeq, __pyx_ptype_8obitools_7profile_8_profile_DNAProfile))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_INCREF(__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.horizontalSeq);
  __pyx_v_horizontalSeq = ((struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *)__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.horizontalSeq);

  /* "obitools/align/_gprofilenws.pyx":48
 *         cdef DNAProfile newProfile2
 *         cdef DNAProfile horizontalSeq=self.horizontalSeq
 *         cdef DNAProfile verticalSeq=self.verticalSeq             # <<<<<<<<<<<<<<
 * 
 *         if self._needToCompute():
 */
  if (!(likely(((__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.verticalSeq) == Py_None) || likely(__Pyx_TypeTest(__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.verticalSeq, __pyx_ptype_8obitools_7profile_8_profile_DNAProfile))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 48; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_INCREF(__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.verticalSeq);
  __pyx_v_verticalSeq = ((struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *)__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.verticalSeq);

  /* "obitools/align/_gprofilenws.pyx":50
 *         cdef DNAProfile verticalSeq=self.verticalSeq
 * 
 *         if self._needToCompute():             # <<<<<<<<<<<<<<
 * 
 *             score = self.doAlignment()
 */
  __pyx_t_2 = ((struct __pyx_vtabstruct_8obitools_5align_12_gprofilenws_GProfileNWS *)__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.__pyx_vtab)->__pyx_base.__pyx_base.__pyx_base._needToCompute(((struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *)__pyx_v_self));
  if (__pyx_t_2) {

    /* "obitools/align/_gprofilenws.pyx":52
 *         if self._needToCompute():
 * 
 *             score = self.doAlignment()             # <<<<<<<<<<<<<<
 *             self.backtrack()
 *             sum = 0
 */
    __pyx_t_3 = ((struct __pyx_vtabstruct_8obitools_5align_12_gprofilenws_GProfileNWS *)__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.__pyx_vtab)->__pyx_base.__pyx_base.__pyx_base.doAlignment(((struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *)__pyx_v_self)); if (unlikely(__pyx_t_3 == 0.0 && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 52; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __pyx_v_score = __pyx_t_3;

    /* "obitools/align/_gprofilenws.pyx":53
 * 
 *             score = self.doAlignment()
 *             self.backtrack()             # <<<<<<<<<<<<<<
 *             sum = 0
 *             for p in xrange(self.path.length) :
 */
    ((struct __pyx_vtabstruct_8obitools_5align_12_gprofilenws_GProfileNWS *)__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.__pyx_vtab)->__pyx_base.__pyx_base.__pyx_base.backtrack(((struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming *)__pyx_v_self));

    /* "obitools/align/_gprofilenws.pyx":54
 *             score = self.doAlignment()
 *             self.backtrack()
 *             sum = 0             # <<<<<<<<<<<<<<
 *             for p in xrange(self.path.length) :
 *                 v = self.path.path[p]
 */
    __Pyx_INCREF(__pyx_int_0);
    __pyx_v_sum = __pyx_int_0;

    /* "obitools/align/_gprofilenws.pyx":55
 *             self.backtrack()
 *             sum = 0
 *             for p in xrange(self.path.length) :             # <<<<<<<<<<<<<<
 *                 v = self.path.path[p]
 *                 if v == 0 :
 */
    __pyx_t_4 = __pyx_v_self->__pyx_base.__pyx_base.__pyx_base.path->length;
    for (__pyx_t_5 = 0; __pyx_t_5 < __pyx_t_4; __pyx_t_5+=1) {
      __pyx_v_p = __pyx_t_5;

      /* "obitools/align/_gprofilenws.pyx":56
 *             sum = 0
 *             for p in xrange(self.path.length) :
 *                 v = self.path.path[p]             # <<<<<<<<<<<<<<
 *                 if v == 0 :
 *                     sum += 1
 */
      __pyx_v_v = (__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.path->path[__pyx_v_p]);

      /* "obitools/align/_gprofilenws.pyx":57
 *             for p in xrange(self.path.length) :
 *                 v = self.path.path[p]
 *                 if v == 0 :             # <<<<<<<<<<<<<<
 *                     sum += 1
 *                 else :
 */
      __pyx_t_2 = (__pyx_v_v == 0);
      if (__pyx_t_2) {

        /* "obitools/align/_gprofilenws.pyx":58
 *                 v = self.path.path[p]
 *                 if v == 0 :
 *                     sum += 1             # <<<<<<<<<<<<<<
 *                 else :
 *                     sum += abs(v)
 */
        __pyx_t_1 = PyNumber_InPlaceAdd(__pyx_v_sum, __pyx_int_1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 58; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_v_sum);
        __pyx_v_sum = __pyx_t_1;
        __pyx_t_1 = 0;
        goto __pyx_L6;
      }
      /*else*/ {

        /* "obitools/align/_gprofilenws.pyx":60
 *                     sum += 1
 *                 else :
 *                     sum += abs(v)             # <<<<<<<<<<<<<<
 * 
 *             newProfile1 = DNAProfile(size=sum,pseudo=pseudocounts)
 */
        __pyx_t_6 = __Pyx_abs_int(__pyx_v_v); 
        __pyx_t_1 = PyLong_FromUnsignedLong(__pyx_t_6); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = PyNumber_InPlaceAdd(__pyx_v_sum, __pyx_t_1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_v_sum);
        __pyx_v_sum = __pyx_t_7;
        __pyx_t_7 = 0;
      }
      __pyx_L6:;
    }

    /* "obitools/align/_gprofilenws.pyx":62
 *                     sum += abs(v)
 * 
 *             newProfile1 = DNAProfile(size=sum,pseudo=pseudocounts)             # <<<<<<<<<<<<<<
 *             newProfile1.profile.weight = horizontalSeq.profile.weight
 * 
 */
    __pyx_t_7 = PyDict_New(); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 62; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(((PyObject *)__pyx_t_7));
    if (PyDict_SetItem(__pyx_t_7, ((PyObject *)__pyx_n_s__size), __pyx_v_sum) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 62; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    if (PyDict_SetItem(__pyx_t_7, ((PyObject *)__pyx_n_s__pseudo), __pyx_v_pseudocounts) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 62; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __pyx_t_1 = PyObject_Call(((PyObject *)((PyObject*)__pyx_ptype_8obitools_7profile_8_profile_DNAProfile)), ((PyObject *)__pyx_empty_tuple), ((PyObject *)__pyx_t_7)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 62; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_1);
    __Pyx_DECREF(((PyObject *)__pyx_t_7)); __pyx_t_7 = 0;
    __pyx_v_newProfile1 = ((struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *)__pyx_t_1);
    __pyx_t_1 = 0;

    /* "obitools/align/_gprofilenws.pyx":63
 * 
 *             newProfile1 = DNAProfile(size=sum,pseudo=pseudocounts)
 *             newProfile1.profile.weight = horizontalSeq.profile.weight             # <<<<<<<<<<<<<<
 * 
 *             newProfile2 = DNAProfile(size=sum,pseudo=pseudocounts)
 */
    __pyx_v_newProfile1->profile->weight = __pyx_v_horizontalSeq->profile->weight;

    /* "obitools/align/_gprofilenws.pyx":65
 *             newProfile1.profile.weight = horizontalSeq.profile.weight
 * 
 *             newProfile2 = DNAProfile(size=sum,pseudo=pseudocounts)             # <<<<<<<<<<<<<<
 *             newProfile2.profile.weight = verticalSeq.profile.weight
 * 
 */
    __pyx_t_1 = PyDict_New(); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 65; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(((PyObject *)__pyx_t_1));
    if (PyDict_SetItem(__pyx_t_1, ((PyObject *)__pyx_n_s__size), __pyx_v_sum) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 65; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    if (PyDict_SetItem(__pyx_t_1, ((PyObject *)__pyx_n_s__pseudo), __pyx_v_pseudocounts) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 65; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __pyx_t_7 = PyObject_Call(((PyObject *)((PyObject*)__pyx_ptype_8obitools_7profile_8_profile_DNAProfile)), ((PyObject *)__pyx_empty_tuple), ((PyObject *)__pyx_t_1)); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 65; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_7);
    __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
    __pyx_v_newProfile2 = ((struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *)__pyx_t_7);
    __pyx_t_7 = 0;

    /* "obitools/align/_gprofilenws.pyx":66
 * 
 *             newProfile2 = DNAProfile(size=sum,pseudo=pseudocounts)
 *             newProfile2.profile.weight = verticalSeq.profile.weight             # <<<<<<<<<<<<<<
 * 
 *             hp=horizontalSeq.profile.length-1
 */
    __pyx_v_newProfile2->profile->weight = __pyx_v_verticalSeq->profile->weight;

    /* "obitools/align/_gprofilenws.pyx":68
 *             newProfile2.profile.weight = verticalSeq.profile.weight
 * 
 *             hp=horizontalSeq.profile.length-1             # <<<<<<<<<<<<<<
 *             vp=verticalSeq.profile.length-1
 *             rp1=newProfile1.profile.length-1
 */
    __pyx_v_hp = (__pyx_v_horizontalSeq->profile->length - 1);

    /* "obitools/align/_gprofilenws.pyx":69
 * 
 *             hp=horizontalSeq.profile.length-1
 *             vp=verticalSeq.profile.length-1             # <<<<<<<<<<<<<<
 *             rp1=newProfile1.profile.length-1
 *             rp2=newProfile2.profile.length-1
 */
    __pyx_v_vp = (__pyx_v_verticalSeq->profile->length - 1);

    /* "obitools/align/_gprofilenws.pyx":70
 *             hp=horizontalSeq.profile.length-1
 *             vp=verticalSeq.profile.length-1
 *             rp1=newProfile1.profile.length-1             # <<<<<<<<<<<<<<
 *             rp2=newProfile2.profile.length-1
 * 
 */
    __pyx_t_7 = PyInt_FromLong((__pyx_v_newProfile1->profile->length - 1)); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_7);
    __pyx_v_rp1 = __pyx_t_7;
    __pyx_t_7 = 0;

    /* "obitools/align/_gprofilenws.pyx":71
 *             vp=verticalSeq.profile.length-1
 *             rp1=newProfile1.profile.length-1
 *             rp2=newProfile2.profile.length-1             # <<<<<<<<<<<<<<
 * 
 *             for i in range(self.path.length):
 */
    __pyx_t_7 = PyInt_FromLong((__pyx_v_newProfile2->profile->length - 1)); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 71; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
    __Pyx_GOTREF(__pyx_t_7);
    __pyx_v_rp2 = __pyx_t_7;
    __pyx_t_7 = 0;

    /* "obitools/align/_gprofilenws.pyx":73
 *             rp2=newProfile2.profile.length-1
 * 
 *             for i in range(self.path.length):             # <<<<<<<<<<<<<<
 *                 p=self.path.path[i]
 * 
 */
    __pyx_t_4 = __pyx_v_self->__pyx_base.__pyx_base.__pyx_base.path->length;
    for (__pyx_t_5 = 0; __pyx_t_5 < __pyx_t_4; __pyx_t_5+=1) {
      __pyx_v_i = __pyx_t_5;

      /* "obitools/align/_gprofilenws.pyx":74
 * 
 *             for i in range(self.path.length):
 *                 p=self.path.path[i]             # <<<<<<<<<<<<<<
 * 
 *             for i in range(self.path.length):
 */
      __pyx_v_p = (__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.path->path[__pyx_v_i]);
    }

    /* "obitools/align/_gprofilenws.pyx":76
 *                 p=self.path.path[i]
 * 
 *             for i in range(self.path.length):             # <<<<<<<<<<<<<<
 *                 p=self.path.path[i]
 * 
 */
    __pyx_t_4 = __pyx_v_self->__pyx_base.__pyx_base.__pyx_base.path->length;
    for (__pyx_t_5 = 0; __pyx_t_5 < __pyx_t_4; __pyx_t_5+=1) {
      __pyx_v_i = __pyx_t_5;

      /* "obitools/align/_gprofilenws.pyx":77
 * 
 *             for i in range(self.path.length):
 *                 p=self.path.path[i]             # <<<<<<<<<<<<<<
 * 
 *                 if p==0:
 */
      __pyx_v_p = (__pyx_v_self->__pyx_base.__pyx_base.__pyx_base.path->path[__pyx_v_i]);

      /* "obitools/align/_gprofilenws.pyx":79
 *                 p=self.path.path[i]
 * 
 *                 if p==0:             # <<<<<<<<<<<<<<
 * 
 *                     newProfile1.A[rp1] = horizontalSeq.A[hp]
 */
      __pyx_t_2 = (__pyx_v_p == 0);
      if (__pyx_t_2) {

        /* "obitools/align/_gprofilenws.pyx":81
 *                 if p==0:
 * 
 *                     newProfile1.A[rp1] = horizontalSeq.A[hp]             # <<<<<<<<<<<<<<
 *                     newProfile1.C[rp1] = horizontalSeq.C[hp]
 *                     newProfile1.G[rp1] = horizontalSeq.G[hp]
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__A); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 81; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 81; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__A); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 81; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp1, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 81; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":82
 * 
 *                     newProfile1.A[rp1] = horizontalSeq.A[hp]
 *                     newProfile1.C[rp1] = horizontalSeq.C[hp]             # <<<<<<<<<<<<<<
 *                     newProfile1.G[rp1] = horizontalSeq.G[hp]
 *                     newProfile1.T[rp1] = horizontalSeq.T[hp]
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__C); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 82; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 82; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__C); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 82; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp1, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 82; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":83
 *                     newProfile1.A[rp1] = horizontalSeq.A[hp]
 *                     newProfile1.C[rp1] = horizontalSeq.C[hp]
 *                     newProfile1.G[rp1] = horizontalSeq.G[hp]             # <<<<<<<<<<<<<<
 *                     newProfile1.T[rp1] = horizontalSeq.T[hp]
 *                     newProfile1.Og[rp1] = horizontalSeq.Og[hp]
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__G); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 83; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 83; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__G); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 83; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp1, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 83; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":84
 *                     newProfile1.C[rp1] = horizontalSeq.C[hp]
 *                     newProfile1.G[rp1] = horizontalSeq.G[hp]
 *                     newProfile1.T[rp1] = horizontalSeq.T[hp]             # <<<<<<<<<<<<<<
 *                     newProfile1.Og[rp1] = horizontalSeq.Og[hp]
 *                     newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__T); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 84; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 84; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__T); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 84; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp1, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 84; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":85
 *                     newProfile1.G[rp1] = horizontalSeq.G[hp]
 *                     newProfile1.T[rp1] = horizontalSeq.T[hp]
 *                     newProfile1.Og[rp1] = horizontalSeq.Og[hp]             # <<<<<<<<<<<<<<
 *                     newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]
 * 
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__Og); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 85; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 85; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__Og); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 85; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp1, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 85; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":86
 *                     newProfile1.T[rp1] = horizontalSeq.T[hp]
 *                     newProfile1.Og[rp1] = horizontalSeq.Og[hp]
 *                     newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]             # <<<<<<<<<<<<<<
 * 
 *                     newProfile2.A[rp2] = verticalSeq.A[vp]
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__Eg); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 86; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 86; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__Eg); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 86; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp1, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 86; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":88
 *                     newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]
 * 
 *                     newProfile2.A[rp2] = verticalSeq.A[vp]             # <<<<<<<<<<<<<<
 *                     newProfile2.C[rp2] = verticalSeq.C[vp]
 *                     newProfile2.G[rp2] = verticalSeq.G[vp]
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__A); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__A); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp2, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 88; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":89
 * 
 *                     newProfile2.A[rp2] = verticalSeq.A[vp]
 *                     newProfile2.C[rp2] = verticalSeq.C[vp]             # <<<<<<<<<<<<<<
 *                     newProfile2.G[rp2] = verticalSeq.G[vp]
 *                     newProfile2.T[rp2] = verticalSeq.T[vp]
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__C); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 89; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 89; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__C); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 89; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp2, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 89; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":90
 *                     newProfile2.A[rp2] = verticalSeq.A[vp]
 *                     newProfile2.C[rp2] = verticalSeq.C[vp]
 *                     newProfile2.G[rp2] = verticalSeq.G[vp]             # <<<<<<<<<<<<<<
 *                     newProfile2.T[rp2] = verticalSeq.T[vp]
 *                     newProfile2.Og[rp2] = verticalSeq.Og[vp]
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__G); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 90; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 90; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__G); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 90; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp2, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 90; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":91
 *                     newProfile2.C[rp2] = verticalSeq.C[vp]
 *                     newProfile2.G[rp2] = verticalSeq.G[vp]
 *                     newProfile2.T[rp2] = verticalSeq.T[vp]             # <<<<<<<<<<<<<<
 *                     newProfile2.Og[rp2] = verticalSeq.Og[vp]
 *                     newProfile2.Eg[rp2] = verticalSeq.Eg[vp]
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__T); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 91; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 91; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__T); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 91; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp2, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 91; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":92
 *                     newProfile2.G[rp2] = verticalSeq.G[vp]
 *                     newProfile2.T[rp2] = verticalSeq.T[vp]
 *                     newProfile2.Og[rp2] = verticalSeq.Og[vp]             # <<<<<<<<<<<<<<
 *                     newProfile2.Eg[rp2] = verticalSeq.Eg[vp]
 * 
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__Og); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 92; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 92; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__Og); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 92; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp2, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 92; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":93
 *                     newProfile2.T[rp2] = verticalSeq.T[vp]
 *                     newProfile2.Og[rp2] = verticalSeq.Og[vp]
 *                     newProfile2.Eg[rp2] = verticalSeq.Eg[vp]             # <<<<<<<<<<<<<<
 * 
 *                     hp-=1
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__Eg); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 93; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 93; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__Eg); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 93; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp2, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 93; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":95
 *                     newProfile2.Eg[rp2] = verticalSeq.Eg[vp]
 * 
 *                     hp-=1             # <<<<<<<<<<<<<<
 *                     vp-=1
 *                     rp1-=1
 */
        __pyx_v_hp = (__pyx_v_hp - 1);

        /* "obitools/align/_gprofilenws.pyx":96
 * 
 *                     hp-=1
 *                     vp-=1             # <<<<<<<<<<<<<<
 *                     rp1-=1
 *                     rp2-=1
 */
        __pyx_v_vp = (__pyx_v_vp - 1);

        /* "obitools/align/_gprofilenws.pyx":97
 *                     hp-=1
 *                     vp-=1
 *                     rp1-=1             # <<<<<<<<<<<<<<
 *                     rp2-=1
 * 
 */
        __pyx_t_7 = PyNumber_InPlaceSubtract(__pyx_v_rp1, __pyx_int_1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 97; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_v_rp1);
        __pyx_v_rp1 = __pyx_t_7;
        __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":98
 *                     vp-=1
 *                     rp1-=1
 *                     rp2-=1             # <<<<<<<<<<<<<<
 * 
 *                 elif p>0:
 */
        __pyx_t_7 = PyNumber_InPlaceSubtract(__pyx_v_rp2, __pyx_int_1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 98; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_v_rp2);
        __pyx_v_rp2 = __pyx_t_7;
        __pyx_t_7 = 0;
        goto __pyx_L11;
      }

      /* "obitools/align/_gprofilenws.pyx":100
 *                     rp2-=1
 * 
 *                 elif p>0:             # <<<<<<<<<<<<<<
 * 
 *                     for x in xrange(p-1) :
 */
      __pyx_t_2 = (__pyx_v_p > 0);
      if (__pyx_t_2) {

        /* "obitools/align/_gprofilenws.pyx":102
 *                 elif p>0:
 * 
 *                     for x in xrange(p-1) :             # <<<<<<<<<<<<<<
 * 
 *                         newProfile1.A[rp1] = horizontalSeq.A[hp]
 */
        __pyx_t_7 = PyInt_FromLong((__pyx_v_p - 1)); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 102; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = PyTuple_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 102; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_t_7);
        __Pyx_GIVEREF(__pyx_t_7);
        __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_Call(__pyx_builtin_xrange, ((PyObject *)__pyx_t_1), NULL); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 102; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
        if (PyList_CheckExact(__pyx_t_7) || PyTuple_CheckExact(__pyx_t_7)) {
          __pyx_t_1 = __pyx_t_7; __Pyx_INCREF(__pyx_t_1); __pyx_t_8 = 0;
          __pyx_t_9 = NULL;
        } else {
          __pyx_t_8 = -1; __pyx_t_1 = PyObject_GetIter(__pyx_t_7); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 102; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          __pyx_t_9 = Py_TYPE(__pyx_t_1)->tp_iternext;
        }
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        for (;;) {
          if (!__pyx_t_9 && PyList_CheckExact(__pyx_t_1)) {
            if (__pyx_t_8 >= PyList_GET_SIZE(__pyx_t_1)) break;
            __pyx_t_7 = PyList_GET_ITEM(__pyx_t_1, __pyx_t_8); __Pyx_INCREF(__pyx_t_7); __pyx_t_8++;
          } else if (!__pyx_t_9 && PyTuple_CheckExact(__pyx_t_1)) {
            if (__pyx_t_8 >= PyTuple_GET_SIZE(__pyx_t_1)) break;
            __pyx_t_7 = PyTuple_GET_ITEM(__pyx_t_1, __pyx_t_8); __Pyx_INCREF(__pyx_t_7); __pyx_t_8++;
          } else {
            __pyx_t_7 = __pyx_t_9(__pyx_t_1);
            if (unlikely(!__pyx_t_7)) {
              if (PyErr_Occurred()) {
                if (likely(PyErr_ExceptionMatches(PyExc_StopIteration))) PyErr_Clear();
                else {__pyx_filename = __pyx_f[0]; __pyx_lineno = 102; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
              }
              break;
            }
            __Pyx_GOTREF(__pyx_t_7);
          }
          __Pyx_XDECREF(__pyx_v_x);
          __pyx_v_x = __pyx_t_7;
          __pyx_t_7 = 0;

          /* "obitools/align/_gprofilenws.pyx":104
 *                     for x in xrange(p-1) :
 * 
 *                         newProfile1.A[rp1] = horizontalSeq.A[hp]             # <<<<<<<<<<<<<<
 *                         newProfile1.C[rp1] = horizontalSeq.C[hp]
 *                         newProfile1.G[rp1] = horizontalSeq.G[hp]
 */
          __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__A); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 104; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          __pyx_t_10 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_10) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 104; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
          __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__A); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 104; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp1, __pyx_t_10) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 104; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;

          /* "obitools/align/_gprofilenws.pyx":105
 * 
 *                         newProfile1.A[rp1] = horizontalSeq.A[hp]
 *                         newProfile1.C[rp1] = horizontalSeq.C[hp]             # <<<<<<<<<<<<<<
 *                         newProfile1.G[rp1] = horizontalSeq.G[hp]
 *                         newProfile1.T[rp1] = horizontalSeq.T[hp]
 */
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__C); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 105; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_10, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 105; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__C); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 105; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          if (PyObject_SetItem(__pyx_t_10, __pyx_v_rp1, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 105; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

          /* "obitools/align/_gprofilenws.pyx":106
 *                         newProfile1.A[rp1] = horizontalSeq.A[hp]
 *                         newProfile1.C[rp1] = horizontalSeq.C[hp]
 *                         newProfile1.G[rp1] = horizontalSeq.G[hp]             # <<<<<<<<<<<<<<
 *                         newProfile1.T[rp1] = horizontalSeq.T[hp]
 *                         newProfile1.Og[rp1] = horizontalSeq.Og[hp]
 */
          __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__G); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 106; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          __pyx_t_10 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_10) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 106; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
          __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__G); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 106; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp1, __pyx_t_10) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 106; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;

          /* "obitools/align/_gprofilenws.pyx":107
 *                         newProfile1.C[rp1] = horizontalSeq.C[hp]
 *                         newProfile1.G[rp1] = horizontalSeq.G[hp]
 *                         newProfile1.T[rp1] = horizontalSeq.T[hp]             # <<<<<<<<<<<<<<
 *                         newProfile1.Og[rp1] = horizontalSeq.Og[hp]
 *                         newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]
 */
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__T); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 107; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_10, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 107; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__T); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 107; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          if (PyObject_SetItem(__pyx_t_10, __pyx_v_rp1, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 107; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

          /* "obitools/align/_gprofilenws.pyx":108
 *                         newProfile1.G[rp1] = horizontalSeq.G[hp]
 *                         newProfile1.T[rp1] = horizontalSeq.T[hp]
 *                         newProfile1.Og[rp1] = horizontalSeq.Og[hp]             # <<<<<<<<<<<<<<
 *                         newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]
 * 
 */
          __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__Og); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 108; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          __pyx_t_10 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_10) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 108; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
          __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__Og); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 108; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp1, __pyx_t_10) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 108; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;

          /* "obitools/align/_gprofilenws.pyx":109
 *                         newProfile1.T[rp1] = horizontalSeq.T[hp]
 *                         newProfile1.Og[rp1] = horizontalSeq.Og[hp]
 *                         newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]             # <<<<<<<<<<<<<<
 * 
 *                         newProfile2.Eg[rp2] = verticalSeq.profile.weight
 */
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__Eg); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 109; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_10, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 109; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__Eg); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 109; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          if (PyObject_SetItem(__pyx_t_10, __pyx_v_rp1, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 109; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

          /* "obitools/align/_gprofilenws.pyx":111
 *                         newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]
 * 
 *                         newProfile2.Eg[rp2] = verticalSeq.profile.weight             # <<<<<<<<<<<<<<
 * 
 *                         hp-=1
 */
          __pyx_t_7 = PyInt_FromLong(__pyx_v_verticalSeq->profile->weight); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 111; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__Eg); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 111; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          if (PyObject_SetItem(__pyx_t_10, __pyx_v_rp2, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 111; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

          /* "obitools/align/_gprofilenws.pyx":113
 *                         newProfile2.Eg[rp2] = verticalSeq.profile.weight
 * 
 *                         hp-=1             # <<<<<<<<<<<<<<
 *                         rp1-=1
 *                         rp2-=1
 */
          __pyx_v_hp = (__pyx_v_hp - 1);

          /* "obitools/align/_gprofilenws.pyx":114
 * 
 *                         hp-=1
 *                         rp1-=1             # <<<<<<<<<<<<<<
 *                         rp2-=1
 * 
 */
          __pyx_t_7 = PyNumber_InPlaceSubtract(__pyx_v_rp1, __pyx_int_1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 114; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          __Pyx_DECREF(__pyx_v_rp1);
          __pyx_v_rp1 = __pyx_t_7;
          __pyx_t_7 = 0;

          /* "obitools/align/_gprofilenws.pyx":115
 *                         hp-=1
 *                         rp1-=1
 *                         rp2-=1             # <<<<<<<<<<<<<<
 * 
 *                     newProfile1.A[rp1] = horizontalSeq.A[hp]
 */
          __pyx_t_7 = PyNumber_InPlaceSubtract(__pyx_v_rp2, __pyx_int_1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 115; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          __Pyx_DECREF(__pyx_v_rp2);
          __pyx_v_rp2 = __pyx_t_7;
          __pyx_t_7 = 0;
        }
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":117
 *                         rp2-=1
 * 
 *                     newProfile1.A[rp1] = horizontalSeq.A[hp]             # <<<<<<<<<<<<<<
 *                     newProfile1.C[rp1] = horizontalSeq.C[hp]
 *                     newProfile1.G[rp1] = horizontalSeq.G[hp]
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__A); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 117; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 117; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__A); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 117; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp1, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 117; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":118
 * 
 *                     newProfile1.A[rp1] = horizontalSeq.A[hp]
 *                     newProfile1.C[rp1] = horizontalSeq.C[hp]             # <<<<<<<<<<<<<<
 *                     newProfile1.G[rp1] = horizontalSeq.G[hp]
 *                     newProfile1.T[rp1] = horizontalSeq.T[hp]
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__C); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 118; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 118; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__C); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 118; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp1, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 118; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":119
 *                     newProfile1.A[rp1] = horizontalSeq.A[hp]
 *                     newProfile1.C[rp1] = horizontalSeq.C[hp]
 *                     newProfile1.G[rp1] = horizontalSeq.G[hp]             # <<<<<<<<<<<<<<
 *                     newProfile1.T[rp1] = horizontalSeq.T[hp]
 *                     newProfile1.Og[rp1] = horizontalSeq.Og[hp]
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__G); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 119; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 119; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__G); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 119; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp1, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 119; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":120
 *                     newProfile1.C[rp1] = horizontalSeq.C[hp]
 *                     newProfile1.G[rp1] = horizontalSeq.G[hp]
 *                     newProfile1.T[rp1] = horizontalSeq.T[hp]             # <<<<<<<<<<<<<<
 *                     newProfile1.Og[rp1] = horizontalSeq.Og[hp]
 *                     newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__T); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 120; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 120; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__T); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 120; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp1, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 120; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":121
 *                     newProfile1.G[rp1] = horizontalSeq.G[hp]
 *                     newProfile1.T[rp1] = horizontalSeq.T[hp]
 *                     newProfile1.Og[rp1] = horizontalSeq.Og[hp]             # <<<<<<<<<<<<<<
 *                     newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]
 * 
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__Og); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 121; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 121; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__Og); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 121; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp1, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 121; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":122
 *                     newProfile1.T[rp1] = horizontalSeq.T[hp]
 *                     newProfile1.Og[rp1] = horizontalSeq.Og[hp]
 *                     newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]             # <<<<<<<<<<<<<<
 * 
 *                     newProfile2.Og[rp2] = verticalSeq.profile.weight
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_horizontalSeq), __pyx_n_s__Eg); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 122; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_hp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 122; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__Eg); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 122; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp1, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 122; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":124
 *                     newProfile1.Eg[rp1] = horizontalSeq.Eg[hp]
 * 
 *                     newProfile2.Og[rp2] = verticalSeq.profile.weight             # <<<<<<<<<<<<<<
 * 
 *                     hp-=1
 */
        __pyx_t_1 = PyInt_FromLong(__pyx_v_verticalSeq->profile->weight); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 124; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__Og); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 124; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp2, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 124; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":126
 *                     newProfile2.Og[rp2] = verticalSeq.profile.weight
 * 
 *                     hp-=1             # <<<<<<<<<<<<<<
 *                     rp1-=1
 *                     rp2-=1
 */
        __pyx_v_hp = (__pyx_v_hp - 1);

        /* "obitools/align/_gprofilenws.pyx":127
 * 
 *                     hp-=1
 *                     rp1-=1             # <<<<<<<<<<<<<<
 *                     rp2-=1
 * 
 */
        __pyx_t_1 = PyNumber_InPlaceSubtract(__pyx_v_rp1, __pyx_int_1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 127; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_v_rp1);
        __pyx_v_rp1 = __pyx_t_1;
        __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":128
 *                     hp-=1
 *                     rp1-=1
 *                     rp2-=1             # <<<<<<<<<<<<<<
 * 
 *                 else:
 */
        __pyx_t_1 = PyNumber_InPlaceSubtract(__pyx_v_rp2, __pyx_int_1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 128; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_v_rp2);
        __pyx_v_rp2 = __pyx_t_1;
        __pyx_t_1 = 0;
        goto __pyx_L11;
      }
      /*else*/ {

        /* "obitools/align/_gprofilenws.pyx":132
 *                 else:
 * 
 *                     for x in xrange(abs(p)-1) :             # <<<<<<<<<<<<<<
 * 
 *                         newProfile2.A[rp2] = verticalSeq.A[vp]
 */
        __pyx_t_6 = __Pyx_abs_int(__pyx_v_p); 
        __pyx_t_1 = PyInt_FromLong((__pyx_t_6 - 1)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 132; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = PyTuple_New(1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 132; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        PyTuple_SET_ITEM(__pyx_t_7, 0, __pyx_t_1);
        __Pyx_GIVEREF(__pyx_t_1);
        __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_Call(__pyx_builtin_xrange, ((PyObject *)__pyx_t_7), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 132; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(((PyObject *)__pyx_t_7)); __pyx_t_7 = 0;
        if (PyList_CheckExact(__pyx_t_1) || PyTuple_CheckExact(__pyx_t_1)) {
          __pyx_t_7 = __pyx_t_1; __Pyx_INCREF(__pyx_t_7); __pyx_t_8 = 0;
          __pyx_t_9 = NULL;
        } else {
          __pyx_t_8 = -1; __pyx_t_7 = PyObject_GetIter(__pyx_t_1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 132; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_7);
          __pyx_t_9 = Py_TYPE(__pyx_t_7)->tp_iternext;
        }
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        for (;;) {
          if (!__pyx_t_9 && PyList_CheckExact(__pyx_t_7)) {
            if (__pyx_t_8 >= PyList_GET_SIZE(__pyx_t_7)) break;
            __pyx_t_1 = PyList_GET_ITEM(__pyx_t_7, __pyx_t_8); __Pyx_INCREF(__pyx_t_1); __pyx_t_8++;
          } else if (!__pyx_t_9 && PyTuple_CheckExact(__pyx_t_7)) {
            if (__pyx_t_8 >= PyTuple_GET_SIZE(__pyx_t_7)) break;
            __pyx_t_1 = PyTuple_GET_ITEM(__pyx_t_7, __pyx_t_8); __Pyx_INCREF(__pyx_t_1); __pyx_t_8++;
          } else {
            __pyx_t_1 = __pyx_t_9(__pyx_t_7);
            if (unlikely(!__pyx_t_1)) {
              if (PyErr_Occurred()) {
                if (likely(PyErr_ExceptionMatches(PyExc_StopIteration))) PyErr_Clear();
                else {__pyx_filename = __pyx_f[0]; __pyx_lineno = 132; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
              }
              break;
            }
            __Pyx_GOTREF(__pyx_t_1);
          }
          __Pyx_XDECREF(__pyx_v_x);
          __pyx_v_x = __pyx_t_1;
          __pyx_t_1 = 0;

          /* "obitools/align/_gprofilenws.pyx":134
 *                     for x in xrange(abs(p)-1) :
 * 
 *                         newProfile2.A[rp2] = verticalSeq.A[vp]             # <<<<<<<<<<<<<<
 *                         newProfile2.C[rp2] = verticalSeq.C[vp]
 *                         newProfile2.G[rp2] = verticalSeq.G[vp]
 */
          __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__A); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 134; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          __pyx_t_10 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_10) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 134; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
          __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__A); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 134; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp2, __pyx_t_10) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 134; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;

          /* "obitools/align/_gprofilenws.pyx":135
 * 
 *                         newProfile2.A[rp2] = verticalSeq.A[vp]
 *                         newProfile2.C[rp2] = verticalSeq.C[vp]             # <<<<<<<<<<<<<<
 *                         newProfile2.G[rp2] = verticalSeq.G[vp]
 *                         newProfile2.T[rp2] = verticalSeq.T[vp]
 */
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__C); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 135; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_10, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 135; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__C); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 135; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          if (PyObject_SetItem(__pyx_t_10, __pyx_v_rp2, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 135; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

          /* "obitools/align/_gprofilenws.pyx":136
 *                         newProfile2.A[rp2] = verticalSeq.A[vp]
 *                         newProfile2.C[rp2] = verticalSeq.C[vp]
 *                         newProfile2.G[rp2] = verticalSeq.G[vp]             # <<<<<<<<<<<<<<
 *                         newProfile2.T[rp2] = verticalSeq.T[vp]
 *                         newProfile2.Og[rp2] = verticalSeq.Og[vp]
 */
          __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__G); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 136; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          __pyx_t_10 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_10) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 136; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
          __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__G); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 136; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp2, __pyx_t_10) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 136; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;

          /* "obitools/align/_gprofilenws.pyx":137
 *                         newProfile2.C[rp2] = verticalSeq.C[vp]
 *                         newProfile2.G[rp2] = verticalSeq.G[vp]
 *                         newProfile2.T[rp2] = verticalSeq.T[vp]             # <<<<<<<<<<<<<<
 *                         newProfile2.Og[rp2] = verticalSeq.Og[vp]
 *                         newProfile2.Eg[rp2] = verticalSeq.Eg[vp]
 */
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__T); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 137; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_10, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 137; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__T); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 137; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          if (PyObject_SetItem(__pyx_t_10, __pyx_v_rp2, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 137; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

          /* "obitools/align/_gprofilenws.pyx":138
 *                         newProfile2.G[rp2] = verticalSeq.G[vp]
 *                         newProfile2.T[rp2] = verticalSeq.T[vp]
 *                         newProfile2.Og[rp2] = verticalSeq.Og[vp]             # <<<<<<<<<<<<<<
 *                         newProfile2.Eg[rp2] = verticalSeq.Eg[vp]
 * 
 */
          __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__Og); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 138; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          __pyx_t_10 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_10) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 138; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
          __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__Og); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 138; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp2, __pyx_t_10) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 138; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;

          /* "obitools/align/_gprofilenws.pyx":139
 *                         newProfile2.T[rp2] = verticalSeq.T[vp]
 *                         newProfile2.Og[rp2] = verticalSeq.Og[vp]
 *                         newProfile2.Eg[rp2] = verticalSeq.Eg[vp]             # <<<<<<<<<<<<<<
 * 
 *                         newProfile1.Eg[rp1] = horizontalSeq.profile.weight
 */
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__Eg); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 139; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_10, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 139; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__Eg); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 139; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          if (PyObject_SetItem(__pyx_t_10, __pyx_v_rp2, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 139; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

          /* "obitools/align/_gprofilenws.pyx":141
 *                         newProfile2.Eg[rp2] = verticalSeq.Eg[vp]
 * 
 *                         newProfile1.Eg[rp1] = horizontalSeq.profile.weight             # <<<<<<<<<<<<<<
 * 
 *                         vp-=1
 */
          __pyx_t_1 = PyInt_FromLong(__pyx_v_horizontalSeq->profile->weight); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 141; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          __pyx_t_10 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__Eg); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 141; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_10);
          if (PyObject_SetItem(__pyx_t_10, __pyx_v_rp1, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 141; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_DECREF(__pyx_t_10); __pyx_t_10 = 0;
          __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

          /* "obitools/align/_gprofilenws.pyx":143
 *                         newProfile1.Eg[rp1] = horizontalSeq.profile.weight
 * 
 *                         vp-=1             # <<<<<<<<<<<<<<
 *                         rp1-=1
 *                         rp2-=1
 */
          __pyx_v_vp = (__pyx_v_vp - 1);

          /* "obitools/align/_gprofilenws.pyx":144
 * 
 *                         vp-=1
 *                         rp1-=1             # <<<<<<<<<<<<<<
 *                         rp2-=1
 * 
 */
          __pyx_t_1 = PyNumber_InPlaceSubtract(__pyx_v_rp1, __pyx_int_1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 144; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          __Pyx_DECREF(__pyx_v_rp1);
          __pyx_v_rp1 = __pyx_t_1;
          __pyx_t_1 = 0;

          /* "obitools/align/_gprofilenws.pyx":145
 *                         vp-=1
 *                         rp1-=1
 *                         rp2-=1             # <<<<<<<<<<<<<<
 * 
 *                     newProfile2.A[rp2] = verticalSeq.A[vp]
 */
          __pyx_t_1 = PyNumber_InPlaceSubtract(__pyx_v_rp2, __pyx_int_1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 145; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
          __Pyx_GOTREF(__pyx_t_1);
          __Pyx_DECREF(__pyx_v_rp2);
          __pyx_v_rp2 = __pyx_t_1;
          __pyx_t_1 = 0;
        }
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":147
 *                         rp2-=1
 * 
 *                     newProfile2.A[rp2] = verticalSeq.A[vp]             # <<<<<<<<<<<<<<
 *                     newProfile2.C[rp2] = verticalSeq.C[vp]
 *                     newProfile2.G[rp2] = verticalSeq.G[vp]
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__A); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 147; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 147; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__A); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 147; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp2, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 147; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":148
 * 
 *                     newProfile2.A[rp2] = verticalSeq.A[vp]
 *                     newProfile2.C[rp2] = verticalSeq.C[vp]             # <<<<<<<<<<<<<<
 *                     newProfile2.G[rp2] = verticalSeq.G[vp]
 *                     newProfile2.T[rp2] = verticalSeq.T[vp]
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__C); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 148; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 148; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__C); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 148; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp2, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 148; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":149
 *                     newProfile2.A[rp2] = verticalSeq.A[vp]
 *                     newProfile2.C[rp2] = verticalSeq.C[vp]
 *                     newProfile2.G[rp2] = verticalSeq.G[vp]             # <<<<<<<<<<<<<<
 *                     newProfile2.T[rp2] = verticalSeq.T[vp]
 *                     newProfile2.Og[rp2] = verticalSeq.Og[vp]
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__G); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 149; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 149; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__G); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 149; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp2, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 149; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":150
 *                     newProfile2.C[rp2] = verticalSeq.C[vp]
 *                     newProfile2.G[rp2] = verticalSeq.G[vp]
 *                     newProfile2.T[rp2] = verticalSeq.T[vp]             # <<<<<<<<<<<<<<
 *                     newProfile2.Og[rp2] = verticalSeq.Og[vp]
 *                     newProfile2.Eg[rp2] = verticalSeq.Eg[vp]
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__T); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 150; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 150; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__T); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 150; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp2, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 150; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":151
 *                     newProfile2.G[rp2] = verticalSeq.G[vp]
 *                     newProfile2.T[rp2] = verticalSeq.T[vp]
 *                     newProfile2.Og[rp2] = verticalSeq.Og[vp]             # <<<<<<<<<<<<<<
 *                     newProfile2.Eg[rp2] = verticalSeq.Eg[vp]
 * 
 */
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__Og); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 151; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = __Pyx_GetItemInt(__pyx_t_7, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 151; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __pyx_t_7 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__Og); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 151; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        if (PyObject_SetItem(__pyx_t_7, __pyx_v_rp2, __pyx_t_1) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 151; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;

        /* "obitools/align/_gprofilenws.pyx":152
 *                     newProfile2.T[rp2] = verticalSeq.T[vp]
 *                     newProfile2.Og[rp2] = verticalSeq.Og[vp]
 *                     newProfile2.Eg[rp2] = verticalSeq.Eg[vp]             # <<<<<<<<<<<<<<
 * 
 *                     newProfile1.Og[rp1] = horizontalSeq.profile.weight
 */
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_verticalSeq), __pyx_n_s__Eg); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 152; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        __pyx_t_7 = __Pyx_GetItemInt(__pyx_t_1, __pyx_v_vp, sizeof(int), PyInt_FromLong); if (!__pyx_t_7) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 152; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile2), __pyx_n_s__Eg); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 152; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp2, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 152; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":154
 *                     newProfile2.Eg[rp2] = verticalSeq.Eg[vp]
 * 
 *                     newProfile1.Og[rp1] = horizontalSeq.profile.weight             # <<<<<<<<<<<<<<
 * 
 *                     vp-=1
 */
        __pyx_t_7 = PyInt_FromLong(__pyx_v_horizontalSeq->profile->weight); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 154; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __pyx_t_1 = PyObject_GetAttr(((PyObject *)__pyx_v_newProfile1), __pyx_n_s__Og); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 154; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_1);
        if (PyObject_SetItem(__pyx_t_1, __pyx_v_rp1, __pyx_t_7) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 154; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
        __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":156
 *                     newProfile1.Og[rp1] = horizontalSeq.profile.weight
 * 
 *                     vp-=1             # <<<<<<<<<<<<<<
 *                     rp1-=1
 *                     rp2-=1
 */
        __pyx_v_vp = (__pyx_v_vp - 1);

        /* "obitools/align/_gprofilenws.pyx":157
 * 
 *                     vp-=1
 *                     rp1-=1             # <<<<<<<<<<<<<<
 *                     rp2-=1
 * 
 */
        __pyx_t_7 = PyNumber_InPlaceSubtract(__pyx_v_rp1, __pyx_int_1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 157; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_v_rp1);
        __pyx_v_rp1 = __pyx_t_7;
        __pyx_t_7 = 0;

        /* "obitools/align/_gprofilenws.pyx":158
 *                     vp-=1
 *                     rp1-=1
 *                     rp2-=1             # <<<<<<<<<<<<<<
 * 
 *             self.alignment1 = newProfile1
 */
        __pyx_t_7 = PyNumber_InPlaceSubtract(__pyx_v_rp2, __pyx_int_1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 158; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
        __Pyx_GOTREF(__pyx_t_7);
        __Pyx_DECREF(__pyx_v_rp2);
        __pyx_v_rp2 = __pyx_t_7;
        __pyx_t_7 = 0;
      }
      __pyx_L11:;
    }

    /* "obitools/align/_gprofilenws.pyx":160
 *                     rp2-=1
 * 
 *             self.alignment1 = newProfile1             # <<<<<<<<<<<<<<
 *             self.alignment2 = newProfile2
 * 
 */
    __Pyx_INCREF(((PyObject *)__pyx_v_newProfile1));
    __Pyx_GIVEREF(((PyObject *)__pyx_v_newProfile1));
    __Pyx_GOTREF(__pyx_v_self->alignment1);
    __Pyx_DECREF(__pyx_v_self->alignment1);
    __pyx_v_self->alignment1 = ((PyObject *)__pyx_v_newProfile1);

    /* "obitools/align/_gprofilenws.pyx":161
 * 
 *             self.alignment1 = newProfile1
 *             self.alignment2 = newProfile2             # <<<<<<<<<<<<<<
 * 
 *         ali1=DNAProfile(self.alignment1,pseudo=pseudocounts)
 */
    __Pyx_INCREF(((PyObject *)__pyx_v_newProfile2));
    __Pyx_GIVEREF(((PyObject *)__pyx_v_newProfile2));
    __Pyx_GOTREF(__pyx_v_self->alignment2);
    __Pyx_DECREF(__pyx_v_self->alignment2);
    __pyx_v_self->alignment2 = ((PyObject *)__pyx_v_newProfile2);
    goto __pyx_L3;
  }
  __pyx_L3:;

  /* "obitools/align/_gprofilenws.pyx":163
 *             self.alignment2 = newProfile2
 * 
 *         ali1=DNAProfile(self.alignment1,pseudo=pseudocounts)             # <<<<<<<<<<<<<<
 *         ali2=DNAProfile(self.alignment2,pseudo=pseudocounts)
 * 
 */
  __pyx_t_7 = PyTuple_New(1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 163; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  __Pyx_INCREF(__pyx_v_self->alignment1);
  PyTuple_SET_ITEM(__pyx_t_7, 0, __pyx_v_self->alignment1);
  __Pyx_GIVEREF(__pyx_v_self->alignment1);
  __pyx_t_1 = PyDict_New(); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 163; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_1));
  if (PyDict_SetItem(__pyx_t_1, ((PyObject *)__pyx_n_s__pseudo), __pyx_v_pseudocounts) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 163; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_10 = PyObject_Call(((PyObject *)((PyObject*)__pyx_ptype_8obitools_7profile_8_profile_DNAProfile)), ((PyObject *)__pyx_t_7), ((PyObject *)__pyx_t_1)); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 163; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_10);
  __Pyx_DECREF(((PyObject *)__pyx_t_7)); __pyx_t_7 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
  __pyx_v_ali1 = ((struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *)__pyx_t_10);
  __pyx_t_10 = 0;

  /* "obitools/align/_gprofilenws.pyx":164
 * 
 *         ali1=DNAProfile(self.alignment1,pseudo=pseudocounts)
 *         ali2=DNAProfile(self.alignment2,pseudo=pseudocounts)             # <<<<<<<<<<<<<<
 * 
 *         return ali1, ali2
 */
  __pyx_t_10 = PyTuple_New(1); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 164; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_10);
  __Pyx_INCREF(__pyx_v_self->alignment2);
  PyTuple_SET_ITEM(__pyx_t_10, 0, __pyx_v_self->alignment2);
  __Pyx_GIVEREF(__pyx_v_self->alignment2);
  __pyx_t_1 = PyDict_New(); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 164; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(((PyObject *)__pyx_t_1));
  if (PyDict_SetItem(__pyx_t_1, ((PyObject *)__pyx_n_s__pseudo), __pyx_v_pseudocounts) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 164; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_t_7 = PyObject_Call(((PyObject *)((PyObject*)__pyx_ptype_8obitools_7profile_8_profile_DNAProfile)), ((PyObject *)__pyx_t_10), ((PyObject *)__pyx_t_1)); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 164; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  __Pyx_DECREF(((PyObject *)__pyx_t_10)); __pyx_t_10 = 0;
  __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0;
  __pyx_v_ali2 = ((struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile *)__pyx_t_7);
  __pyx_t_7 = 0;

  /* "obitools/align/_gprofilenws.pyx":166
 *         ali2=DNAProfile(self.alignment2,pseudo=pseudocounts)
 * 
 *         return ali1, ali2             # <<<<<<<<<<<<<<
 * 
 */
  __Pyx_XDECREF(__pyx_r);
  __pyx_t_7 = PyTuple_New(2); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 166; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __Pyx_GOTREF(__pyx_t_7);
  __Pyx_INCREF(((PyObject *)__pyx_v_ali1));
  PyTuple_SET_ITEM(__pyx_t_7, 0, ((PyObject *)__pyx_v_ali1));
  __Pyx_GIVEREF(((PyObject *)__pyx_v_ali1));
  __Pyx_INCREF(((PyObject *)__pyx_v_ali2));
  PyTuple_SET_ITEM(__pyx_t_7, 1, ((PyObject *)__pyx_v_ali2));
  __Pyx_GIVEREF(((PyObject *)__pyx_v_ali2));
  __pyx_r = ((PyObject *)__pyx_t_7);
  __pyx_t_7 = 0;
  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_7);
  __Pyx_XDECREF(__pyx_t_10);
  __Pyx_AddTraceback("obitools.align._gprofilenws.GProfileNWS.__call__", __pyx_clineno, __pyx_lineno, __pyx_filename);
  __pyx_r = NULL;
  __pyx_L0:;
  __Pyx_XDECREF(__pyx_v_hgaps);
  __Pyx_XDECREF(__pyx_v_vgaps);
  __Pyx_XDECREF((PyObject *)__pyx_v_newProfile1);
  __Pyx_XDECREF((PyObject *)__pyx_v_newProfile2);
  __Pyx_XDECREF((PyObject *)__pyx_v_horizontalSeq);
  __Pyx_XDECREF((PyObject *)__pyx_v_verticalSeq);
  __Pyx_XDECREF(__pyx_v_sum);
  __Pyx_XDECREF(__pyx_v_rp1);
  __Pyx_XDECREF(__pyx_v_rp2);
  __Pyx_XDECREF(__pyx_v_x);
  __Pyx_XDECREF((PyObject *)__pyx_v_ali1);
  __Pyx_XDECREF((PyObject *)__pyx_v_ali2);
  __Pyx_XGIVEREF(__pyx_r);
  __Pyx_RefNannyFinishContext();
  return __pyx_r;
}
static struct __pyx_vtabstruct_8obitools_5align_12_gprofilenws_GProfileNWS __pyx_vtable_8obitools_5align_12_gprofilenws_GProfileNWS;

static PyObject *__pyx_tp_new_8obitools_5align_12_gprofilenws_GProfileNWS(PyTypeObject *t, PyObject *a, PyObject *k) {
  struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *p;
  PyObject *o = __pyx_ptype_8obitools_5align_11_profilenws_ProfileNWS->tp_new(t, a, k);
  if (!o) return 0;
  p = ((struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *)o);
  p->__pyx_base.__pyx_base.__pyx_base.__pyx_vtab = (struct __pyx_vtabstruct_8obitools_5align_8_dynamic_DynamicProgramming*)__pyx_vtabptr_8obitools_5align_12_gprofilenws_GProfileNWS;
  p->alignment1 = Py_None; Py_INCREF(Py_None);
  p->alignment2 = Py_None; Py_INCREF(Py_None);
  return o;
}

static void __pyx_tp_dealloc_8obitools_5align_12_gprofilenws_GProfileNWS(PyObject *o) {
  struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *p = (struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *)o;
  Py_XDECREF(p->alignment1);
  Py_XDECREF(p->alignment2);
  __pyx_ptype_8obitools_5align_11_profilenws_ProfileNWS->tp_dealloc(o);
}

static int __pyx_tp_traverse_8obitools_5align_12_gprofilenws_GProfileNWS(PyObject *o, visitproc v, void *a) {
  int e;
  struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *p = (struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *)o;
  if (__pyx_ptype_8obitools_5align_11_profilenws_ProfileNWS->tp_traverse) {
    e = __pyx_ptype_8obitools_5align_11_profilenws_ProfileNWS->tp_traverse(o, v, a); if (e) return e;
  }
  if (p->alignment1) {
    e = (*v)(p->alignment1, a); if (e) return e;
  }
  if (p->alignment2) {
    e = (*v)(p->alignment2, a); if (e) return e;
  }
  return 0;
}

static int __pyx_tp_clear_8obitools_5align_12_gprofilenws_GProfileNWS(PyObject *o) {
  struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *p = (struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS *)o;
  PyObject* tmp;
  if (__pyx_ptype_8obitools_5align_11_profilenws_ProfileNWS->tp_clear) {
    __pyx_ptype_8obitools_5align_11_profilenws_ProfileNWS->tp_clear(o);
  }
  tmp = ((PyObject*)p->alignment1);
  p->alignment1 = Py_None; Py_INCREF(Py_None);
  Py_XDECREF(tmp);
  tmp = ((PyObject*)p->alignment2);
  p->alignment2 = Py_None; Py_INCREF(Py_None);
  Py_XDECREF(tmp);
  return 0;
}

static PyMethodDef __pyx_methods_8obitools_5align_12_gprofilenws_GProfileNWS[] = {
  {0, 0, 0, 0}
};

static PyNumberMethods __pyx_tp_as_number_GProfileNWS = {
  0, /*nb_add*/
  0, /*nb_subtract*/
  0, /*nb_multiply*/
  #if PY_MAJOR_VERSION < 3
  0, /*nb_divide*/
  #endif
  0, /*nb_remainder*/
  0, /*nb_divmod*/
  0, /*nb_power*/
  0, /*nb_negative*/
  0, /*nb_positive*/
  0, /*nb_absolute*/
  0, /*nb_nonzero*/
  0, /*nb_invert*/
  0, /*nb_lshift*/
  0, /*nb_rshift*/
  0, /*nb_and*/
  0, /*nb_xor*/
  0, /*nb_or*/
  #if PY_MAJOR_VERSION < 3
  0, /*nb_coerce*/
  #endif
  0, /*nb_int*/
  #if PY_MAJOR_VERSION < 3
  0, /*nb_long*/
  #else
  0, /*reserved*/
  #endif
  0, /*nb_float*/
  #if PY_MAJOR_VERSION < 3
  0, /*nb_oct*/
  #endif
  #if PY_MAJOR_VERSION < 3
  0, /*nb_hex*/
  #endif
  0, /*nb_inplace_add*/
  0, /*nb_inplace_subtract*/
  0, /*nb_inplace_multiply*/
  #if PY_MAJOR_VERSION < 3
  0, /*nb_inplace_divide*/
  #endif
  0, /*nb_inplace_remainder*/
  0, /*nb_inplace_power*/
  0, /*nb_inplace_lshift*/
  0, /*nb_inplace_rshift*/
  0, /*nb_inplace_and*/
  0, /*nb_inplace_xor*/
  0, /*nb_inplace_or*/
  0, /*nb_floor_divide*/
  0, /*nb_true_divide*/
  0, /*nb_inplace_floor_divide*/
  0, /*nb_inplace_true_divide*/
  #if PY_VERSION_HEX >= 0x02050000
  0, /*nb_index*/
  #endif
};

static PySequenceMethods __pyx_tp_as_sequence_GProfileNWS = {
  0, /*sq_length*/
  0, /*sq_concat*/
  0, /*sq_repeat*/
  0, /*sq_item*/
  0, /*sq_slice*/
  0, /*sq_ass_item*/
  0, /*sq_ass_slice*/
  0, /*sq_contains*/
  0, /*sq_inplace_concat*/
  0, /*sq_inplace_repeat*/
};

static PyMappingMethods __pyx_tp_as_mapping_GProfileNWS = {
  0, /*mp_length*/
  0, /*mp_subscript*/
  0, /*mp_ass_subscript*/
};

static PyBufferProcs __pyx_tp_as_buffer_GProfileNWS = {
  #if PY_MAJOR_VERSION < 3
  0, /*bf_getreadbuffer*/
  #endif
  #if PY_MAJOR_VERSION < 3
  0, /*bf_getwritebuffer*/
  #endif
  #if PY_MAJOR_VERSION < 3
  0, /*bf_getsegcount*/
  #endif
  #if PY_MAJOR_VERSION < 3
  0, /*bf_getcharbuffer*/
  #endif
  #if PY_VERSION_HEX >= 0x02060000
  0, /*bf_getbuffer*/
  #endif
  #if PY_VERSION_HEX >= 0x02060000
  0, /*bf_releasebuffer*/
  #endif
};

static PyTypeObject __pyx_type_8obitools_5align_12_gprofilenws_GProfileNWS = {
  PyVarObject_HEAD_INIT(0, 0)
  __Pyx_NAMESTR("obitools.align._gprofilenws.GProfileNWS"), /*tp_name*/
  sizeof(struct __pyx_obj_8obitools_5align_12_gprofilenws_GProfileNWS), /*tp_basicsize*/
  0, /*tp_itemsize*/
  __pyx_tp_dealloc_8obitools_5align_12_gprofilenws_GProfileNWS, /*tp_dealloc*/
  0, /*tp_print*/
  0, /*tp_getattr*/
  0, /*tp_setattr*/
  #if PY_MAJOR_VERSION < 3
  0, /*tp_compare*/
  #else
  0, /*reserved*/
  #endif
  0, /*tp_repr*/
  &__pyx_tp_as_number_GProfileNWS, /*tp_as_number*/
  &__pyx_tp_as_sequence_GProfileNWS, /*tp_as_sequence*/
  &__pyx_tp_as_mapping_GProfileNWS, /*tp_as_mapping*/
  0, /*tp_hash*/
  __pyx_pw_8obitools_5align_12_gprofilenws_11GProfileNWS_1__call__, /*tp_call*/
  0, /*tp_str*/
  0, /*tp_getattro*/
  0, /*tp_setattro*/
  &__pyx_tp_as_buffer_GProfileNWS, /*tp_as_buffer*/
  Py_TPFLAGS_DEFAULT|Py_TPFLAGS_CHECKTYPES|Py_TPFLAGS_HAVE_NEWBUFFER|Py_TPFLAGS_BASETYPE|Py_TPFLAGS_HAVE_GC, /*tp_flags*/
  0, /*tp_doc*/
  __pyx_tp_traverse_8obitools_5align_12_gprofilenws_GProfileNWS, /*tp_traverse*/
  __pyx_tp_clear_8obitools_5align_12_gprofilenws_GProfileNWS, /*tp_clear*/
  0, /*tp_richcompare*/
  0, /*tp_weaklistoffset*/
  0, /*tp_iter*/
  0, /*tp_iternext*/
  __pyx_methods_8obitools_5align_12_gprofilenws_GProfileNWS, /*tp_methods*/
  0, /*tp_members*/
  0, /*tp_getset*/
  0, /*tp_base*/
  0, /*tp_dict*/
  0, /*tp_descr_get*/
  0, /*tp_descr_set*/
  0, /*tp_dictoffset*/
  0, /*tp_init*/
  0, /*tp_alloc*/
  __pyx_tp_new_8obitools_5align_12_gprofilenws_GProfileNWS, /*tp_new*/
  0, /*tp_free*/
  0, /*tp_is_gc*/
  0, /*tp_bases*/
  0, /*tp_mro*/
  0, /*tp_cache*/
  0, /*tp_subclasses*/
  0, /*tp_weaklist*/
  0, /*tp_del*/
  #if PY_VERSION_HEX >= 0x02060000
  0, /*tp_version_tag*/
  #endif
};

static PyMethodDef __pyx_methods[] = {
  {0, 0, 0, 0}
};

#if PY_MAJOR_VERSION >= 3
static struct PyModuleDef __pyx_moduledef = {
    PyModuleDef_HEAD_INIT,
    __Pyx_NAMESTR("_gprofilenws"),
    __Pyx_DOCSTR(__pyx_k_1), /* 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__A, __pyx_k__A, sizeof(__pyx_k__A), 0, 0, 1, 1},
  {&__pyx_n_s__C, __pyx_k__C, sizeof(__pyx_k__C), 0, 0, 1, 1},
  {&__pyx_n_s__Eg, __pyx_k__Eg, sizeof(__pyx_k__Eg), 0, 0, 1, 1},
  {&__pyx_n_s__G, __pyx_k__G, sizeof(__pyx_k__G), 0, 0, 1, 1},
  {&__pyx_n_s__Og, __pyx_k__Og, sizeof(__pyx_k__Og), 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__pseudo, __pyx_k__pseudo, sizeof(__pyx_k__pseudo), 0, 0, 1, 1},
  {&__pyx_n_s__pseudocounts, __pyx_k__pseudocounts, sizeof(__pyx_k__pseudocounts), 0, 0, 1, 1},
  {&__pyx_n_s__range, __pyx_k__range, sizeof(__pyx_k__range), 0, 0, 1, 1},
  {&__pyx_n_s__size, __pyx_k__size, sizeof(__pyx_k__size), 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 = 55; __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 = 55; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  #endif
  __pyx_builtin_range = __Pyx_GetName(__pyx_b, __pyx_n_s__range); if (!__pyx_builtin_range) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 73; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  return 0;
  __pyx_L1_error:;
  return -1;
}

static int __Pyx_InitCachedConstants(void) {
  __Pyx_RefNannyDeclarations
  __Pyx_RefNannySetupContext("__Pyx_InitCachedConstants", 0);
  __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;};
  return 0;
  __pyx_L1_error:;
  return -1;
}

#if PY_MAJOR_VERSION < 3
PyMODINIT_FUNC init_gprofilenws(void); /*proto*/
PyMODINIT_FUNC init_gprofilenws(void)
#else
PyMODINIT_FUNC PyInit__gprofilenws(void); /*proto*/
PyMODINIT_FUNC PyInit__gprofilenws(void)
#endif
{
  PyObject *__pyx_t_1 = 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__gprofilenws(void)", 0);
  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_CyFunction_USED
  if (__Pyx_CyFunction_init() < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  #endif
  #ifdef __Pyx_FusedFunction_USED
  if (__pyx_FusedFunction_init() < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  #endif
  #ifdef __Pyx_Generator_USED
  if (__pyx_Generator_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("_gprofilenws"), __pyx_methods, __Pyx_DOCSTR(__pyx_k_1), 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_obitools__align___gprofilenws) {
    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 ---*/
  __pyx_ptype_8obitools_5align_11_profilenws_ProfileNWS = __Pyx_ImportType("obitools.align._profilenws", "ProfileNWS", sizeof(struct __pyx_obj_8obitools_5align_11_profilenws_ProfileNWS), 1); if (unlikely(!__pyx_ptype_8obitools_5align_11_profilenws_ProfileNWS)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_vtabptr_8obitools_5align_11_profilenws_ProfileNWS = (struct __pyx_vtabstruct_8obitools_5align_11_profilenws_ProfileNWS*)__Pyx_GetVtable(__pyx_ptype_8obitools_5align_11_profilenws_ProfileNWS->tp_dict); if (unlikely(!__pyx_vtabptr_8obitools_5align_11_profilenws_ProfileNWS)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 1; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_vtabptr_8obitools_5align_12_gprofilenws_GProfileNWS = &__pyx_vtable_8obitools_5align_12_gprofilenws_GProfileNWS;
  __pyx_vtable_8obitools_5align_12_gprofilenws_GProfileNWS.__pyx_base = *__pyx_vtabptr_8obitools_5align_11_profilenws_ProfileNWS;
  __pyx_vtable_8obitools_5align_12_gprofilenws_GProfileNWS.__pyx_base.__pyx_base.matchScore = (double (*)(struct __pyx_obj_8obitools_5align_4_nws_NWS *, int, int))__pyx_f_8obitools_5align_12_gprofilenws_11GProfileNWS_matchScore;
  __pyx_type_8obitools_5align_12_gprofilenws_GProfileNWS.tp_base = __pyx_ptype_8obitools_5align_11_profilenws_ProfileNWS;
  if (PyType_Ready(&__pyx_type_8obitools_5align_12_gprofilenws_GProfileNWS) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 11; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (__Pyx_SetVtable(__pyx_type_8obitools_5align_12_gprofilenws_GProfileNWS.tp_dict, __pyx_vtabptr_8obitools_5align_12_gprofilenws_GProfileNWS) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 11; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  if (__Pyx_SetAttrString(__pyx_m, "GProfileNWS", (PyObject *)&__pyx_type_8obitools_5align_12_gprofilenws_GProfileNWS) < 0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 11; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_ptype_8obitools_5align_12_gprofilenws_GProfileNWS = &__pyx_type_8obitools_5align_12_gprofilenws_GProfileNWS;
  /*--- Type import code ---*/
  __pyx_ptype_8obitools_5align_8_dynamic_DynamicProgramming = __Pyx_ImportType("obitools.align._dynamic", "DynamicProgramming", sizeof(struct __pyx_obj_8obitools_5align_8_dynamic_DynamicProgramming), 1); if (unlikely(!__pyx_ptype_8obitools_5align_8_dynamic_DynamicProgramming)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 63; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_vtabptr_8obitools_5align_8_dynamic_DynamicProgramming = (struct __pyx_vtabstruct_8obitools_5align_8_dynamic_DynamicProgramming*)__Pyx_GetVtable(__pyx_ptype_8obitools_5align_8_dynamic_DynamicProgramming->tp_dict); if (unlikely(!__pyx_vtabptr_8obitools_5align_8_dynamic_DynamicProgramming)) {__pyx_filename = __pyx_f[1]; __pyx_lineno = 63; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_ptype_8obitools_5align_4_nws_NWS = __Pyx_ImportType("obitools.align._nws", "NWS", sizeof(struct __pyx_obj_8obitools_5align_4_nws_NWS), 1); if (unlikely(!__pyx_ptype_8obitools_5align_4_nws_NWS)) {__pyx_filename = __pyx_f[2]; __pyx_lineno = 3; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_vtabptr_8obitools_5align_4_nws_NWS = (struct __pyx_vtabstruct_8obitools_5align_4_nws_NWS*)__Pyx_GetVtable(__pyx_ptype_8obitools_5align_4_nws_NWS->tp_dict); if (unlikely(!__pyx_vtabptr_8obitools_5align_4_nws_NWS)) {__pyx_filename = __pyx_f[2]; __pyx_lineno = 3; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_ptype_8obitools_7profile_8_profile__MemIntArray = __Pyx_ImportType("obitools.profile._profile", "_MemIntArray", sizeof(struct __pyx_obj_8obitools_7profile_8_profile__MemIntArray), 1); if (unlikely(!__pyx_ptype_8obitools_7profile_8_profile__MemIntArray)) {__pyx_filename = __pyx_f[3]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_vtabptr_8obitools_7profile_8_profile__MemIntArray = (struct __pyx_vtabstruct_8obitools_7profile_8_profile__MemIntArray*)__Pyx_GetVtable(__pyx_ptype_8obitools_7profile_8_profile__MemIntArray->tp_dict); if (unlikely(!__pyx_vtabptr_8obitools_7profile_8_profile__MemIntArray)) {__pyx_filename = __pyx_f[3]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_ptype_8obitools_7profile_8_profile_DNAProfile = __Pyx_ImportType("obitools.profile._profile", "DNAProfile", sizeof(struct __pyx_obj_8obitools_7profile_8_profile_DNAProfile), 1); if (unlikely(!__pyx_ptype_8obitools_7profile_8_profile_DNAProfile)) {__pyx_filename = __pyx_f[3]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  __pyx_vtabptr_8obitools_7profile_8_profile_DNAProfile = (struct __pyx_vtabstruct_8obitools_7profile_8_profile_DNAProfile*)__Pyx_GetVtable(__pyx_ptype_8obitools_7profile_8_profile_DNAProfile->tp_dict); if (unlikely(!__pyx_vtabptr_8obitools_7profile_8_profile_DNAProfile)) {__pyx_filename = __pyx_f[3]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;}
  /*--- Variable import code ---*/
  /*--- Function import code ---*/
  /*--- Execution code ---*/

  /* "obitools/align/_gprofilenws.pyx":1
 * '''             # <<<<<<<<<<<<<<
 * Created on 16 Feb. 2011
 * 
 */
  __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);
  if (__pyx_m) {
    __Pyx_AddTraceback("init obitools.align._gprofilenws", __pyx_clineno, __pyx_lineno, __pyx_filename);
    Py_DECREF(__pyx_m); __pyx_m = 0;
  } else if (!PyErr_Occurred()) {
    PyErr_SetString(PyExc_ImportError, "init obitools.align._gprofilenws");
  }
  __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_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)
{
    PyObject *key = 0, *value = 0;
    Py_ssize_t pos = 0;
    PyObject*** name;
    PyObject*** first_kw_arg = argnames + num_pos_args;
    while (PyDict_Next(kwds, &pos, &key, &value)) {
        name = first_kw_arg;
        while (*name && (**name != key)) name++;
        if (*name) {
            values[name-argnames] = value;
        } else {
            #if PY_MAJOR_VERSION < 3
            if (unlikely(!PyString_CheckExact(key)) && unlikely(!PyString_Check(key))) {
            #else
            if (unlikely(!PyUnicode_Check(key))) {
            #endif
                goto invalid_keyword_type;
            } else {
                for (name = first_kw_arg; *name; name++) {
                    #if PY_MAJOR_VERSION >= 3
                    if (PyUnicode_GET_SIZE(**name) == PyUnicode_GET_SIZE(key) &&
                        PyUnicode_Compare(**name, key) == 0) break;
                    #else
                    if (PyString_GET_SIZE(**name) == PyString_GET_SIZE(key) &&
                        _PyString_Eq(**name, key)) break;
                    #endif
                }
                if (*name) {
                    values[name-argnames] = value;
                } else {
                    for (name=argnames; name != first_kw_arg; name++) {
                        if (**name == key) goto arg_passed_twice;
                        #if PY_MAJOR_VERSION >= 3
                        if (PyUnicode_GET_SIZE(**name) == PyUnicode_GET_SIZE(key) &&
                            PyUnicode_Compare(**name, key) == 0) goto arg_passed_twice;
                        #else
                        if (PyString_GET_SIZE(**name) == PyString_GET_SIZE(key) &&
                            _PyString_Eq(**name, key)) goto arg_passed_twice;
                        #endif
                    }
                    if (kwds2) {
                        if (unlikely(PyDict_SetItem(kwds2, key, value))) goto bad;
                    } else {
                        goto invalid_keyword;
                    }
                }
            }
        }
    }
    return 0;
arg_passed_twice:
    __Pyx_RaiseDoubleKeywordsError(function_name, **name);
    goto bad;
invalid_keyword_type:
    PyErr_Format(PyExc_TypeError,
        "%s() keywords must be strings", function_name);
    goto bad;
invalid_keyword:
    PyErr_Format(PyExc_TypeError,
    #if PY_MAJOR_VERSION < 3
        "%s() got an unexpected keyword argument '%s'",
        function_name, PyString_AsString(key));
    #else
        "%s() got an unexpected keyword argument '%U'",
        function_name, key);
    #endif
bad:
    return -1;
}

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 CYTHON_INLINE int __Pyx_TypeTest(PyObject *obj, PyTypeObject *type) {
    if (unlikely(!type)) {
        PyErr_Format(PyExc_SystemError, "Missing type object");
        return 0;
    }
    if (likely(PyObject_TypeCheck(obj, type)))
        return 1;
    PyErr_Format(PyExc_TypeError, "Cannot convert %.200s to %.200s",
                 Py_TYPE(obj)->tp_name, type->tp_name);
    return 0;
}



static CYTHON_INLINE unsigned char __Pyx_PyInt_AsUnsignedChar(PyObject* x) {
    const unsigned char neg_one = (unsigned char)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
    if (sizeof(unsigned char) < sizeof(long)) {
        long val = __Pyx_PyInt_AsLong(x);
        if (unlikely(val != (long)(unsigned char)val)) {
            if (!unlikely(val == -1 && PyErr_Occurred())) {
                PyErr_SetString(PyExc_OverflowError,
                    (is_unsigned && unlikely(val < 0)) ?
                    "can't convert negative value to unsigned char" :
                    "value too large to convert to unsigned char");
            }
            return (unsigned char)-1;
        }
        return (unsigned char)val;
    }
    return (unsigned char)__Pyx_PyInt_AsUnsignedLong(x);
}

static CYTHON_INLINE unsigned short __Pyx_PyInt_AsUnsignedShort(PyObject* x) {
    const unsigned short neg_one = (unsigned short)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
    if (sizeof(unsigned short) < sizeof(long)) {
        long val = __Pyx_PyInt_AsLong(x);
        if (unlikely(val != (long)(unsigned short)val)) {
            if (!unlikely(val == -1 && PyErr_Occurred())) {
                PyErr_SetString(PyExc_OverflowError,
                    (is_unsigned && unlikely(val < 0)) ?
                    "can't convert negative value to unsigned short" :
                    "value too large to convert to unsigned short");
            }
            return (unsigned short)-1;
        }
        return (unsigned short)val;
    }
    return (unsigned short)__Pyx_PyInt_AsUnsignedLong(x);
}

static CYTHON_INLINE unsigned int __Pyx_PyInt_AsUnsignedInt(PyObject* x) {
    const unsigned int neg_one = (unsigned int)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
    if (sizeof(unsigned int) < sizeof(long)) {
        long val = __Pyx_PyInt_AsLong(x);
        if (unlikely(val != (long)(unsigned int)val)) {
            if (!unlikely(val == -1 && PyErr_Occurred())) {
                PyErr_SetString(PyExc_OverflowError,
                    (is_unsigned && unlikely(val < 0)) ?
                    "can't convert negative value to unsigned int" :
                    "value too large to convert to unsigned int");
            }
            return (unsigned int)-1;
        }
        return (unsigned int)val;
    }
    return (unsigned int)__Pyx_PyInt_AsUnsignedLong(x);
}

static CYTHON_INLINE char __Pyx_PyInt_AsChar(PyObject* x) {
    const char neg_one = (char)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
    if (sizeof(char) < sizeof(long)) {
        long val = __Pyx_PyInt_AsLong(x);
        if (unlikely(val != (long)(char)val)) {
            if (!unlikely(val == -1 && PyErr_Occurred())) {
                PyErr_SetString(PyExc_OverflowError,
                    (is_unsigned && unlikely(val < 0)) ?
                    "can't convert negative value to char" :
                    "value too large to convert to char");
            }
            return (char)-1;
        }
        return (char)val;
    }
    return (char)__Pyx_PyInt_AsLong(x);
}

static CYTHON_INLINE short __Pyx_PyInt_AsShort(PyObject* x) {
    const short neg_one = (short)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
    if (sizeof(short) < sizeof(long)) {
        long val = __Pyx_PyInt_AsLong(x);
        if (unlikely(val != (long)(short)val)) {
            if (!unlikely(val == -1 && PyErr_Occurred())) {
                PyErr_SetString(PyExc_OverflowError,
                    (is_unsigned && unlikely(val < 0)) ?
                    "can't convert negative value to short" :
                    "value too large to convert to short");
            }
            return (short)-1;
        }
        return (short)val;
    }
    return (short)__Pyx_PyInt_AsLong(x);
}

static CYTHON_INLINE int __Pyx_PyInt_AsInt(PyObject* x) {
    const int neg_one = (int)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
    if (sizeof(int) < sizeof(long)) {
        long val = __Pyx_PyInt_AsLong(x);
        if (unlikely(val != (long)(int)val)) {
            if (!unlikely(val == -1 && PyErr_Occurred())) {
                PyErr_SetString(PyExc_OverflowError,
                    (is_unsigned && unlikely(val < 0)) ?
                    "can't convert negative value to int" :
                    "value too large to convert to int");
            }
            return (int)-1;
        }
        return (int)val;
    }
    return (int)__Pyx_PyInt_AsLong(x);
}

static CYTHON_INLINE signed char __Pyx_PyInt_AsSignedChar(PyObject* x) {
    const signed char neg_one = (signed char)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
    if (sizeof(signed char) < sizeof(long)) {
        long val = __Pyx_PyInt_AsLong(x);
        if (unlikely(val != (long)(signed char)val)) {
            if (!unlikely(val == -1 && PyErr_Occurred())) {
                PyErr_SetString(PyExc_OverflowError,
                    (is_unsigned && unlikely(val < 0)) ?
                    "can't convert negative value to signed char" :
                    "value too large to convert to signed char");
            }
            return (signed char)-1;
        }
        return (signed char)val;
    }
    return (signed char)__Pyx_PyInt_AsSignedLong(x);
}

static CYTHON_INLINE signed short __Pyx_PyInt_AsSignedShort(PyObject* x) {
    const signed short neg_one = (signed short)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
    if (sizeof(signed short) < sizeof(long)) {
        long val = __Pyx_PyInt_AsLong(x);
        if (unlikely(val != (long)(signed short)val)) {
            if (!unlikely(val == -1 && PyErr_Occurred())) {
                PyErr_SetString(PyExc_OverflowError,
                    (is_unsigned && unlikely(val < 0)) ?
                    "can't convert negative value to signed short" :
                    "value too large to convert to signed short");
            }
            return (signed short)-1;
        }
        return (signed short)val;
    }
    return (signed short)__Pyx_PyInt_AsSignedLong(x);
}

static CYTHON_INLINE signed int __Pyx_PyInt_AsSignedInt(PyObject* x) {
    const signed int neg_one = (signed int)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
    if (sizeof(signed int) < sizeof(long)) {
        long val = __Pyx_PyInt_AsLong(x);
        if (unlikely(val != (long)(signed int)val)) {
            if (!unlikely(val == -1 && PyErr_Occurred())) {
                PyErr_SetString(PyExc_OverflowError,
                    (is_unsigned && unlikely(val < 0)) ?
                    "can't convert negative value to signed int" :
                    "value too large to convert to signed int");
            }
            return (signed int)-1;
        }
        return (signed int)val;
    }
    return (signed int)__Pyx_PyInt_AsSignedLong(x);
}

static CYTHON_INLINE int __Pyx_PyInt_AsLongDouble(PyObject* x) {
    const int neg_one = (int)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
    if (sizeof(int) < sizeof(long)) {
        long val = __Pyx_PyInt_AsLong(x);
        if (unlikely(val != (long)(int)val)) {
            if (!unlikely(val == -1 && PyErr_Occurred())) {
                PyErr_SetString(PyExc_OverflowError,
                    (is_unsigned && unlikely(val < 0)) ?
                    "can't convert negative value to int" :
                    "value too large to convert to int");
            }
            return (int)-1;
        }
        return (int)val;
    }
    return (int)__Pyx_PyInt_AsLong(x);
}

static CYTHON_INLINE unsigned long __Pyx_PyInt_AsUnsignedLong(PyObject* x) {
    const unsigned long neg_one = (unsigned long)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
#if PY_VERSION_HEX < 0x03000000
    if (likely(PyInt_Check(x))) {
        long val = PyInt_AS_LONG(x);
        if (is_unsigned && unlikely(val < 0)) {
            PyErr_SetString(PyExc_OverflowError,
                            "can't convert negative value to unsigned long");
            return (unsigned long)-1;
        }
        return (unsigned long)val;
    } else
#endif
    if (likely(PyLong_Check(x))) {
        if (is_unsigned) {
            if (unlikely(Py_SIZE(x) < 0)) {
                PyErr_SetString(PyExc_OverflowError,
                                "can't convert negative value to unsigned long");
                return (unsigned long)-1;
            }
            return (unsigned long)PyLong_AsUnsignedLong(x);
        } else {
            return (unsigned long)PyLong_AsLong(x);
        }
    } else {
        unsigned long val;
        PyObject *tmp = __Pyx_PyNumber_Int(x);
        if (!tmp) return (unsigned long)-1;
        val = __Pyx_PyInt_AsUnsignedLong(tmp);
        Py_DECREF(tmp);
        return val;
    }
}

static CYTHON_INLINE unsigned PY_LONG_LONG __Pyx_PyInt_AsUnsignedLongLong(PyObject* x) {
    const unsigned PY_LONG_LONG neg_one = (unsigned PY_LONG_LONG)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
#if PY_VERSION_HEX < 0x03000000
    if (likely(PyInt_Check(x))) {
        long val = PyInt_AS_LONG(x);
        if (is_unsigned && unlikely(val < 0)) {
            PyErr_SetString(PyExc_OverflowError,
                            "can't convert negative value to unsigned PY_LONG_LONG");
            return (unsigned PY_LONG_LONG)-1;
        }
        return (unsigned PY_LONG_LONG)val;
    } else
#endif
    if (likely(PyLong_Check(x))) {
        if (is_unsigned) {
            if (unlikely(Py_SIZE(x) < 0)) {
                PyErr_SetString(PyExc_OverflowError,
                                "can't convert negative value to unsigned PY_LONG_LONG");
                return (unsigned PY_LONG_LONG)-1;
            }
            return (unsigned PY_LONG_LONG)PyLong_AsUnsignedLongLong(x);
        } else {
            return (unsigned PY_LONG_LONG)PyLong_AsLongLong(x);
        }
    } else {
        unsigned PY_LONG_LONG val;
        PyObject *tmp = __Pyx_PyNumber_Int(x);
        if (!tmp) return (unsigned PY_LONG_LONG)-1;
        val = __Pyx_PyInt_AsUnsignedLongLong(tmp);
        Py_DECREF(tmp);
        return val;
    }
}

static CYTHON_INLINE long __Pyx_PyInt_AsLong(PyObject* x) {
    const long neg_one = (long)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
#if PY_VERSION_HEX < 0x03000000
    if (likely(PyInt_Check(x))) {
        long val = PyInt_AS_LONG(x);
        if (is_unsigned && unlikely(val < 0)) {
            PyErr_SetString(PyExc_OverflowError,
                            "can't convert negative value to long");
            return (long)-1;
        }
        return (long)val;
    } else
#endif
    if (likely(PyLong_Check(x))) {
        if (is_unsigned) {
            if (unlikely(Py_SIZE(x) < 0)) {
                PyErr_SetString(PyExc_OverflowError,
                                "can't convert negative value to long");
                return (long)-1;
            }
            return (long)PyLong_AsUnsignedLong(x);
        } else {
            return (long)PyLong_AsLong(x);
        }
    } else {
        long val;
        PyObject *tmp = __Pyx_PyNumber_Int(x);
        if (!tmp) return (long)-1;
        val = __Pyx_PyInt_AsLong(tmp);
        Py_DECREF(tmp);
        return val;
    }
}

static CYTHON_INLINE PY_LONG_LONG __Pyx_PyInt_AsLongLong(PyObject* x) {
    const PY_LONG_LONG neg_one = (PY_LONG_LONG)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
#if PY_VERSION_HEX < 0x03000000
    if (likely(PyInt_Check(x))) {
        long val = PyInt_AS_LONG(x);
        if (is_unsigned && unlikely(val < 0)) {
            PyErr_SetString(PyExc_OverflowError,
                            "can't convert negative value to PY_LONG_LONG");
            return (PY_LONG_LONG)-1;
        }
        return (PY_LONG_LONG)val;
    } else
#endif
    if (likely(PyLong_Check(x))) {
        if (is_unsigned) {
            if (unlikely(Py_SIZE(x) < 0)) {
                PyErr_SetString(PyExc_OverflowError,
                                "can't convert negative value to PY_LONG_LONG");
                return (PY_LONG_LONG)-1;
            }
            return (PY_LONG_LONG)PyLong_AsUnsignedLongLong(x);
        } else {
            return (PY_LONG_LONG)PyLong_AsLongLong(x);
        }
    } else {
        PY_LONG_LONG val;
        PyObject *tmp = __Pyx_PyNumber_Int(x);
        if (!tmp) return (PY_LONG_LONG)-1;
        val = __Pyx_PyInt_AsLongLong(tmp);
        Py_DECREF(tmp);
        return val;
    }
}

static CYTHON_INLINE signed long __Pyx_PyInt_AsSignedLong(PyObject* x) {
    const signed long neg_one = (signed long)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
#if PY_VERSION_HEX < 0x03000000
    if (likely(PyInt_Check(x))) {
        long val = PyInt_AS_LONG(x);
        if (is_unsigned && unlikely(val < 0)) {
            PyErr_SetString(PyExc_OverflowError,
                            "can't convert negative value to signed long");
            return (signed long)-1;
        }
        return (signed long)val;
    } else
#endif
    if (likely(PyLong_Check(x))) {
        if (is_unsigned) {
            if (unlikely(Py_SIZE(x) < 0)) {
                PyErr_SetString(PyExc_OverflowError,
                                "can't convert negative value to signed long");
                return (signed long)-1;
            }
            return (signed long)PyLong_AsUnsignedLong(x);
        } else {
            return (signed long)PyLong_AsLong(x);
        }
    } else {
        signed long val;
        PyObject *tmp = __Pyx_PyNumber_Int(x);
        if (!tmp) return (signed long)-1;
        val = __Pyx_PyInt_AsSignedLong(tmp);
        Py_DECREF(tmp);
        return val;
    }
}

static CYTHON_INLINE signed PY_LONG_LONG __Pyx_PyInt_AsSignedLongLong(PyObject* x) {
    const signed PY_LONG_LONG neg_one = (signed PY_LONG_LONG)-1, const_zero = 0;
    const int is_unsigned = neg_one > const_zero;
#if PY_VERSION_HEX < 0x03000000
    if (likely(PyInt_Check(x))) {
        long val = PyInt_AS_LONG(x);
        if (is_unsigned && unlikely(val < 0)) {
            PyErr_SetString(PyExc_OverflowError,
                            "can't convert negative value to signed PY_LONG_LONG");
            return (signed PY_LONG_LONG)-1;
        }
        return (signed PY_LONG_LONG)val;
    } else
#endif
    if (likely(PyLong_Check(x))) {
        if (is_unsigned) {
            if (unlikely(Py_SIZE(x) < 0)) {
                PyErr_SetString(PyExc_OverflowError,
                                "can't convert negative value to signed PY_LONG_LONG");
                return (signed PY_LONG_LONG)-1;
            }
            return (signed PY_LONG_LONG)PyLong_AsUnsignedLongLong(x);
        } else {
            return (signed PY_LONG_LONG)PyLong_AsLongLong(x);
        }
    } else {
        signed PY_LONG_LONG val;
        PyObject *tmp = __Pyx_PyNumber_Int(x);
        if (!tmp) return (signed PY_LONG_LONG)-1;
        val = __Pyx_PyInt_AsSignedLongLong(tmp);
        Py_DECREF(tmp);
        return val;
    }
}

static int __Pyx_check_binary_version(void) {
    char ctversion[4], rtversion[4];
    PyOS_snprintf(ctversion, 4, "%d.%d", PY_MAJOR_VERSION, PY_MINOR_VERSION);
    PyOS_snprintf(rtversion, 4, "%s", Py_GetVersion());
    if (ctversion[0] != rtversion[0] || ctversion[2] != rtversion[2]) {
        char message[200];
        PyOS_snprintf(message, sizeof(message),
                      "compiletime version %s of module '%.100s' "
                      "does not match runtime version %s",
                      ctversion, __Pyx_MODULE_NAME, rtversion);
        #if PY_VERSION_HEX < 0x02050000
        return PyErr_Warn(NULL, message);
        #else
        return PyErr_WarnEx(NULL, message, 1);
        #endif
    }
    return 0;
}

#ifndef __PYX_HAVE_RT_ImportType
#define __PYX_HAVE_RT_ImportType
static PyTypeObject *__Pyx_ImportType(const char *module_name, const char *class_name,
    size_t size, int strict)
{
    PyObject *py_module = 0;
    PyObject *result = 0;
    PyObject *py_name = 0;
    char warning[200];
    py_module = __Pyx_ImportModule(module_name);
    if (!py_module)
        goto bad;
    py_name = __Pyx_PyIdentifier_FromString(class_name);
    if (!py_name)
        goto bad;
    result = PyObject_GetAttr(py_module, py_name);
    Py_DECREF(py_name);
    py_name = 0;
    Py_DECREF(py_module);
    py_module = 0;
    if (!result)
        goto bad;
    if (!PyType_Check(result)) {
        PyErr_Format(PyExc_TypeError,
            "%s.%s is not a type object",
            module_name, class_name);
        goto bad;
    }
    if (!strict && (size_t)((PyTypeObject *)result)->tp_basicsize > size) {
        PyOS_snprintf(warning, sizeof(warning),
            "%s.%s size changed, may indicate binary incompatibility",
            module_name, class_name);
        #if PY_VERSION_HEX < 0x02050000
        if (PyErr_Warn(NULL, warning) < 0) goto bad;
        #else
        if (PyErr_WarnEx(NULL, warning, 0) < 0) goto bad;
        #endif
    }
    else if ((size_t)((PyTypeObject *)result)->tp_basicsize != size) {
        PyErr_Format(PyExc_ValueError,
            "%s.%s has the wrong size, try recompiling",
            module_name, class_name);
        goto bad;
    }
    return (PyTypeObject *)result;
bad:
    Py_XDECREF(py_module);
    Py_XDECREF(result);
    return NULL;
}
#endif

#ifndef __PYX_HAVE_RT_ImportModule
#define __PYX_HAVE_RT_ImportModule
static PyObject *__Pyx_ImportModule(const char *name) {
    PyObject *py_name = 0;
    PyObject *py_module = 0;
    py_name = __Pyx_PyIdentifier_FromString(name);
    if (!py_name)
        goto bad;
    py_module = PyImport_Import(py_name);
    Py_DECREF(py_name);
    return py_module;
bad:
    Py_XDECREF(py_name);
    return 0;
}
#endif

static void* __Pyx_GetVtable(PyObject *dict) {
    void* ptr;
    PyObject *ob = PyMapping_GetItemString(dict, (char *)"__pyx_vtable__");
    if (!ob)
        goto bad;
#if PY_VERSION_HEX >= 0x02070000 && !(PY_MAJOR_VERSION==3&&PY_MINOR_VERSION==0)
    ptr = PyCapsule_GetPointer(ob, 0);
#else
    ptr = PyCObject_AsVoidPtr(ob);
#endif
    if (!ptr && !PyErr_Occurred())
        PyErr_SetString(PyExc_RuntimeError, "invalid vtable found for imported type");
    Py_DECREF(ob);
    return ptr;
bad:
    Py_XDECREF(ob);
    return NULL;
}

static int __Pyx_SetVtable(PyObject *dict, void *vtable) {
#if PY_VERSION_HEX >= 0x02070000 && !(PY_MAJOR_VERSION==3&&PY_MINOR_VERSION==0)
    PyObject *ob = PyCapsule_New(vtable, 0, 0);
#else
    PyObject *ob = PyCObject_FromVoidPtr(vtable, 0);
#endif
    if (!ob)
        goto bad;
    if (PyDict_SetItemString(dict, "__pyx_vtable__", ob) < 0)
        goto bad;
    Py_DECREF(ob);
    return 0;
bad:
    Py_XDECREF(ob);
    return -1;
}

static int __pyx_bisect_code_objects(__Pyx_CodeObjectCacheEntry* entries, int count, int code_line) {
    int start = 0, mid = 0, end = count - 1;
    if (end >= 0 && code_line > entries[end].code_line) {
        return count;
    }
    while (start < end) {
        mid = (start + end) / 2;
        if (code_line < entries[mid].code_line) {
            end = mid;
        } else if (code_line > entries[mid].code_line) {
             start = mid + 1;
        } else {
            return mid;
        }
    }
    if (code_line <= entries[mid].code_line) {
        return mid;
    } else {
        return mid + 1;
    }
}
static PyCodeObject *__pyx_find_code_object(int code_line) {
    PyCodeObject* code_object;
    int pos;
    if (unlikely(!code_line) || unlikely(!__pyx_code_cache.entries)) {
        return NULL;
    }
    pos = __pyx_bisect_code_objects(__pyx_code_cache.entries, __pyx_code_cache.count, code_line);
    if (unlikely(pos >= __pyx_code_cache.count) || unlikely(__pyx_code_cache.entries[pos].code_line != code_line)) {
        return NULL;
    }
    code_object = __pyx_code_cache.entries[pos].code_object;
    Py_INCREF(code_object);
    return code_object;
}
static void __pyx_insert_code_object(int code_line, PyCodeObject* code_object) {
    int pos, i;
    __Pyx_CodeObjectCacheEntry* entries = __pyx_code_cache.entries;
    if (unlikely(!code_line)) {
        return;
    }
    if (unlikely(!entries)) {
        entries = (__Pyx_CodeObjectCacheEntry*)PyMem_Malloc(64*sizeof(__Pyx_CodeObjectCacheEntry));
        if (likely(entries)) {
            __pyx_code_cache.entries = entries;
            __pyx_code_cache.max_count = 64;
            __pyx_code_cache.count = 1;
            entries[0].code_line = code_line;
            entries[0].code_object = code_object;
            Py_INCREF(code_object);
        }
        return;
    }
    pos = __pyx_bisect_code_objects(__pyx_code_cache.entries, __pyx_code_cache.count, code_line);
    if ((pos < __pyx_code_cache.count) && unlikely(__pyx_code_cache.entries[pos].code_line == code_line)) {
        PyCodeObject* tmp = entries[pos].code_object;
        entries[pos].code_object = code_object;
        Py_DECREF(tmp);
        return;
    }
    if (__pyx_code_cache.count == __pyx_code_cache.max_count) {
        int new_max = __pyx_code_cache.max_count + 64;
        entries = (__Pyx_CodeObjectCacheEntry*)PyMem_Realloc(
            __pyx_code_cache.entries, new_max*sizeof(__Pyx_CodeObjectCacheEntry));
        if (unlikely(!entries)) {
            return;
        }
        __pyx_code_cache.entries = entries;
        __pyx_code_cache.max_count = new_max;
    }
    for (i=__pyx_code_cache.count; i>pos; i--) {
        entries[i] = entries[i-1];
    }
    entries[pos].code_line = code_line;
    entries[pos].code_object = code_object;
    __pyx_code_cache.count++;
    Py_INCREF(code_object);
}

#include "compile.h"
#include "frameobject.h"
#include "traceback.h"
static PyCodeObject* __Pyx_CreateCodeObjectForTraceback(
            const char *funcname, int c_line,
            int py_line, const char *filename) {
    PyCodeObject *py_code = 0;
    PyObject *py_srcfile = 0;
    PyObject *py_funcname = 0;
    #if PY_MAJOR_VERSION < 3
    py_srcfile = PyString_FromString(filename);
    #else
    py_srcfile = PyUnicode_FromString(filename);
    #endif
    if (!py_srcfile) goto bad;
    if (c_line) {
        #if PY_MAJOR_VERSION < 3
        py_funcname = PyString_FromFormat( "%s (%s:%d)", funcname, __pyx_cfilenm, c_line);
        #else
        py_funcname = PyUnicode_FromFormat( "%s (%s:%d)", funcname, __pyx_cfilenm, c_line);
        #endif
    }
    else {
        #if PY_MAJOR_VERSION < 3
        py_funcname = PyString_FromString(funcname);
        #else
        py_funcname = PyUnicode_FromString(funcname);
        #endif
    }
    if (!py_funcname) goto bad;
    py_code = __Pyx_PyCode_New(
        0,            /*int argcount,*/
        0,            /*int kwonlyargcount,*/
        0,            /*int nlocals,*/
        0,            /*int stacksize,*/
        0,            /*int flags,*/
        __pyx_empty_bytes, /*PyObject *code,*/
        __pyx_empty_tuple, /*PyObject *consts,*/
        __pyx_empty_tuple, /*PyObject *names,*/
        __pyx_empty_tuple, /*PyObject *varnames,*/
        __pyx_empty_tuple, /*PyObject *freevars,*/
        __pyx_empty_tuple, /*PyObject *cellvars,*/
        py_srcfile,   /*PyObject *filename,*/
        py_funcname,  /*PyObject *name,*/
        py_line,      /*int firstlineno,*/
        __pyx_empty_bytes  /*PyObject *lnotab*/
    );
    Py_DECREF(py_srcfile);
    Py_DECREF(py_funcname);
    return py_code;
bad:
    Py_XDECREF(py_srcfile);
    Py_XDECREF(py_funcname);
    return NULL;
}
static void __Pyx_AddTraceback(const char *funcname, int c_line,
                               int py_line, const char *filename) {
    PyCodeObject *py_code = 0;
    PyObject *py_globals = 0;
    PyFrameObject *py_frame = 0;
    py_code = __pyx_find_code_object(c_line ? c_line : py_line);
    if (!py_code) {
        py_code = __Pyx_CreateCodeObjectForTraceback(
            funcname, c_line, py_line, filename);
        if (!py_code) goto bad;
        __pyx_insert_code_object(c_line ? c_line : py_line, py_code);
    }
    py_globals = PyModule_GetDict(__pyx_m);
    if (!py_globals) goto bad;
    py_frame = PyFrame_New(
        PyThreadState_GET(), /*PyThreadState *tstate,*/
        py_code,             /*PyCodeObject *code,*/
        py_globals,          /*PyObject *globals,*/
        0                    /*PyObject *locals*/
    );
    if (!py_frame) goto bad;
    py_frame->f_lineno = py_line;
    PyTraceBack_Here(py_frame);
bad:
    Py_XDECREF(py_code);
    Py_XDECREF(py_frame);
}

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


/* Type Conversion Functions */

static CYTHON_INLINE int __Pyx_PyObject_IsTrue(PyObject* x) {
   int is_true = x == Py_True;
   if (is_true | (x == Py_False) | (x == Py_None)) return is_true;
   else return PyObject_IsTrue(x);
}

static CYTHON_INLINE PyObject* __Pyx_PyNumber_Int(PyObject* x) {
  PyNumberMethods *m;
  const char *name = NULL;
  PyObject *res = NULL;
#if PY_VERSION_HEX < 0x03000000
  if (PyInt_Check(x) || PyLong_Check(x))
#else
  if (PyLong_Check(x))
#endif
    return Py_INCREF(x), x;
  m = Py_TYPE(x)->tp_as_number;
#if PY_VERSION_HEX < 0x03000000
  if (m && m->nb_int) {
    name = "int";
    res = PyNumber_Int(x);
  }
  else if (m && m->nb_long) {
    name = "long";
    res = PyNumber_Long(x);
  }
#else
  if (m && m->nb_int) {
    name = "int";
    res = PyNumber_Long(x);
  }
#endif
  if (res) {
#if PY_VERSION_HEX < 0x03000000
    if (!PyInt_Check(res) && !PyLong_Check(res)) {
#else
    if (!PyLong_Check(res)) {
#endif
      PyErr_Format(PyExc_TypeError,
                   "__%s__ returned non-%s (type %.200s)",
                   name, name, Py_TYPE(res)->tp_name);
      Py_DECREF(res);
      return NULL;
    }
  }
  else if (!PyErr_Occurred()) {
    PyErr_SetString(PyExc_TypeError,
                    "an integer is required");
  }
  return res;
}

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

static CYTHON_INLINE PyObject * __Pyx_PyInt_FromSize_t(size_t ival) {
#if PY_VERSION_HEX < 0x02050000
   if (ival <= LONG_MAX)
       return PyInt_FromLong((long)ival);
   else {
       unsigned char *bytes = (unsigned char *) &ival;
       int one = 1; int little = (int)*(unsigned char*)&one;
       return _PyLong_FromByteArray(bytes, sizeof(size_t), little, 0);
   }
#else
   return PyInt_FromSize_t(ival);
#endif
}

static CYTHON_INLINE size_t __Pyx_PyInt_AsSize_t(PyObject* x) {
   unsigned PY_LONG_LONG val = __Pyx_PyInt_AsUnsignedLongLong(x);
   if (unlikely(val == (unsigned PY_LONG_LONG)-1 && PyErr_Occurred())) {
       return (size_t)-1;
   } else if (unlikely(val != (unsigned PY_LONG_LONG)(size_t)val)) {
       PyErr_SetString(PyExc_OverflowError,
                       "value too large to convert to size_t");
       return (size_t)-1;
   }
   return (size_t)val;
}


#endif /* Py_PYTHON_H */
