autodiffutils_matlabMATLAB_wrap.cxx

/* ----------------------------------------------------------------------------
 * This file was automatically generated by SWIG (http://www.swig.org).
 * Version 3.0.9
 *
 * This file is not intended to be easily readable and contains a number of
 * coding conventions designed to improve portability and efficiency. Do not make
 * changes to this file unless you know what you are doing--modify the SWIG
 * interface file instead.
 * ----------------------------------------------------------------------------- */

#define SWIGMATLAB
#define SWIG_name_d      "autodiffutils"
#define SWIG_name        autodiffutils

#define SWIG_op_prefix        "op_"
#define SWIG_pkg_name        "autodiffutils"
#define SwigVar_mxArray mxArray*


#ifdef __cplusplus
/* SwigValueWrapper is described in swig.swg */
template<typename T> class SwigValueWrapper {
  struct SwigMovePointer {
    T *ptr;
    SwigMovePointer(T *p) : ptr(p) { }
    ~SwigMovePointer() { delete ptr; }
    SwigMovePointer& operator=(SwigMovePointer& rhs) { T* oldptr = ptr; ptr = 0; delete oldptr; ptr = rhs.ptr; rhs.ptr = 0; return *this; }
  } pointer;
  SwigValueWrapper& operator=(const SwigValueWrapper<T>& rhs);
  SwigValueWrapper(const SwigValueWrapper<T>& rhs);
public:
  SwigValueWrapper() : pointer(0) { }
  SwigValueWrapper& operator=(const T& t) { SwigMovePointer tmp(new T(t)); pointer = tmp; return *this; }
  operator T&() const { return *pointer.ptr; }
  T *operator&() { return pointer.ptr; }
};

template <typename T> T SwigValueInit() {
  return T();
}
#endif

/* -----------------------------------------------------------------------------
 *  This section contains generic SWIG labels for method/variable
 *  declarations/attributes, and other compiler dependent labels.
 * ----------------------------------------------------------------------------- */

/* template workaround for compilers that cannot correctly implement the C++ standard */
#ifndef SWIGTEMPLATEDISAMBIGUATOR
# if defined(__SUNPRO_CC) && (__SUNPRO_CC <= 0x560)
#  define SWIGTEMPLATEDISAMBIGUATOR template
# elif defined(__HP_aCC)
/* Needed even with `aCC -AA' when `aCC -V' reports HP ANSI C++ B3910B A.03.55 */
/* If we find a maximum version that requires this, the test would be __HP_aCC <= 35500 for A.03.55 */
#  define SWIGTEMPLATEDISAMBIGUATOR template
# else
#  define SWIGTEMPLATEDISAMBIGUATOR
# endif
#endif

/* inline attribute */
#ifndef SWIGINLINE
# if defined(__cplusplus) || (defined(__GNUC__) && !defined(__STRICT_ANSI__))
#   define SWIGINLINE inline
# else
#   define SWIGINLINE
# endif
#endif

/* attribute recognised by some compilers to avoid 'unused' warnings */
#ifndef SWIGUNUSED
# if defined(__GNUC__)
#   if !(defined(__cplusplus)) || (__GNUC__ > 3 || (__GNUC__ == 3 && __GNUC_MINOR__ >= 4))
#     define SWIGUNUSED __attribute__ ((__unused__))
#   else
#     define SWIGUNUSED
#   endif
# elif defined(__ICC)
#   define SWIGUNUSED __attribute__ ((__unused__))
# else
#   define SWIGUNUSED
# endif
#endif

#ifndef SWIG_MSC_UNSUPPRESS_4505
# if defined(_MSC_VER)
#   pragma warning(disable : 4505) /* unreferenced local function has been removed */
# endif
#endif

#ifndef SWIGUNUSEDPARM
# ifdef __cplusplus
#   define SWIGUNUSEDPARM(p)
# else
#   define SWIGUNUSEDPARM(p) p SWIGUNUSED
# endif
#endif

/* internal SWIG method */
#ifndef SWIGINTERN
# define SWIGINTERN static SWIGUNUSED
#endif

/* internal inline SWIG method */
#ifndef SWIGINTERNINLINE
# define SWIGINTERNINLINE SWIGINTERN SWIGINLINE
#endif

/* exporting methods */
#if defined(__GNUC__)
#  if (__GNUC__ >= 4) || (__GNUC__ == 3 && __GNUC_MINOR__ >= 4)
#    ifndef GCC_HASCLASSVISIBILITY
#      define GCC_HASCLASSVISIBILITY
#    endif
#  endif
#endif

#ifndef SWIGEXPORT
# if defined(_WIN32) || defined(__WIN32__) || defined(__CYGWIN__)
#   if defined(STATIC_LINKED)
#     define SWIGEXPORT
#   else
#     define SWIGEXPORT __declspec(dllexport)
#   endif
# else
#   if defined(__GNUC__) && defined(GCC_HASCLASSVISIBILITY)
#     define SWIGEXPORT __attribute__ ((visibility("default")))
#   else
#     define SWIGEXPORT
#   endif
# endif
#endif

/* calling conventions for Windows */
#ifndef SWIGSTDCALL
# if defined(_WIN32) || defined(__WIN32__) || defined(__CYGWIN__)
#   define SWIGSTDCALL __stdcall
# else
#   define SWIGSTDCALL
# endif
#endif

/* Deal with Microsoft's attempt at deprecating C standard runtime functions */
#if !defined(SWIG_NO_CRT_SECURE_NO_DEPRECATE) && defined(_MSC_VER) && !defined(_CRT_SECURE_NO_DEPRECATE)
# define _CRT_SECURE_NO_DEPRECATE
#endif

/* Deal with Microsoft's attempt at deprecating methods in the standard C++ library */
#if !defined(SWIG_NO_SCL_SECURE_NO_DEPRECATE) && defined(_MSC_VER) && !defined(_SCL_SECURE_NO_DEPRECATE)
# define _SCL_SECURE_NO_DEPRECATE
#endif

/* Deal with Apple's deprecated 'AssertMacros.h' from Carbon-framework */
#if defined(__APPLE__) && !defined(__ASSERT_MACROS_DEFINE_VERSIONS_WITHOUT_UNDERSCORES)
# define __ASSERT_MACROS_DEFINE_VERSIONS_WITHOUT_UNDERSCORES 0
#endif

/* Intel's compiler complains if a variable which was never initialised is
 * cast to void, which is a common idiom which we use to indicate that we
 * are aware a variable isn't used.  So we just silence that warning.
 * See: https://github.com/swig/swig/issues/192 for more discussion.
 */
#ifdef __INTEL_COMPILER
# pragma warning disable 592
#endif


#include <mex.h>
#include <math.h> /* need this for floor() */
#include <string.h> /* need this for strcmp() */
#include <stdarg.h> /* need this for SWIG_Error below */

/* -----------------------------------------------------------------------------
 * swigrun.swg
 *
 * This file contains generic C API SWIG runtime support for pointer
 * type checking.
 * ----------------------------------------------------------------------------- */

/* This should only be incremented when either the layout of swig_type_info changes,
   or for whatever reason, the runtime changes incompatibly */
#define SWIG_RUNTIME_VERSION "4"

/* define SWIG_TYPE_TABLE_NAME as "SWIG_TYPE_TABLE" */
#ifdef SWIG_TYPE_TABLE
# define SWIG_QUOTE_STRING(x) #x
# define SWIG_EXPAND_AND_QUOTE_STRING(x) SWIG_QUOTE_STRING(x)
# define SWIG_TYPE_TABLE_NAME SWIG_EXPAND_AND_QUOTE_STRING(SWIG_TYPE_TABLE)
#else
# define SWIG_TYPE_TABLE_NAME
#endif

/*
  You can use the SWIGRUNTIME and SWIGRUNTIMEINLINE macros for
  creating a static or dynamic library from the SWIG runtime code.
  In 99.9% of the cases, SWIG just needs to declare them as 'static'.

  But only do this if strictly necessary, ie, if you have problems
  with your compiler or suchlike.
*/

#ifndef SWIGRUNTIME
# define SWIGRUNTIME SWIGINTERN
#endif

#ifndef SWIGRUNTIMEINLINE
# define SWIGRUNTIMEINLINE SWIGRUNTIME SWIGINLINE
#endif

/*  Generic buffer size */
#ifndef SWIG_BUFFER_SIZE
# define SWIG_BUFFER_SIZE 1024
#endif

/* Flags for pointer conversions */
#define SWIG_POINTER_DISOWN        0x1
#define SWIG_CAST_NEW_MEMORY       0x2

/* Flags for new pointer objects */
#define SWIG_POINTER_OWN           0x1


/*
   Flags/methods for returning states.

   The SWIG conversion methods, as ConvertPtr, return an integer
   that tells if the conversion was successful or not. And if not,
   an error code can be returned (see swigerrors.swg for the codes).

   Use the following macros/flags to set or process the returning
   states.

   In old versions of SWIG, code such as the following was usually written:

     if (SWIG_ConvertPtr(obj,vptr,ty.flags) != -1) {
       // success code
     } else {
       //fail code
     }

   Now you can be more explicit:

    int res = SWIG_ConvertPtr(obj,vptr,ty.flags);
    if (SWIG_IsOK(res)) {
      // success code
    } else {
      // fail code
    }

   which is the same really, but now you can also do

    Type *ptr;
    int res = SWIG_ConvertPtr(obj,(void **)(&ptr),ty.flags);
    if (SWIG_IsOK(res)) {
      // success code
      if (SWIG_IsNewObj(res) {
        ...
	delete *ptr;
      } else {
        ...
      }
    } else {
      // fail code
    }

   I.e., now SWIG_ConvertPtr can return new objects and you can
   identify the case and take care of the deallocation. Of course that
   also requires SWIG_ConvertPtr to return new result values, such as

      int SWIG_ConvertPtr(obj, ptr,...) {
        if (<obj is ok>) {
          if (<need new object>) {
            *ptr = <ptr to new allocated object>;
            return SWIG_NEWOBJ;
          } else {
            *ptr = <ptr to old object>;
            return SWIG_OLDOBJ;
          }
        } else {
          return SWIG_BADOBJ;
        }
      }

   Of course, returning the plain '0(success)/-1(fail)' still works, but you can be
   more explicit by returning SWIG_BADOBJ, SWIG_ERROR or any of the
   SWIG errors code.

   Finally, if the SWIG_CASTRANK_MODE is enabled, the result code
   allows to return the 'cast rank', for example, if you have this

       int food(double)
       int fooi(int);

   and you call

      food(1)   // cast rank '1'  (1 -> 1.0)
      fooi(1)   // cast rank '0'

   just use the SWIG_AddCast()/SWIG_CheckState()
*/

#define SWIG_OK                    (0)
#define SWIG_ERROR                 (-1)
#define SWIG_IsOK(r)               (r >= 0)
#define SWIG_ArgError(r)           ((r != SWIG_ERROR) ? r : SWIG_TypeError)

/* The CastRankLimit says how many bits are used for the cast rank */
#define SWIG_CASTRANKLIMIT         (1 << 8)
/* The NewMask denotes the object was created (using new/malloc) */
#define SWIG_NEWOBJMASK            (SWIG_CASTRANKLIMIT  << 1)
/* The TmpMask is for in/out typemaps that use temporal objects */
#define SWIG_TMPOBJMASK            (SWIG_NEWOBJMASK << 1)
/* Simple returning values */
#define SWIG_BADOBJ                (SWIG_ERROR)
#define SWIG_OLDOBJ                (SWIG_OK)
#define SWIG_NEWOBJ                (SWIG_OK | SWIG_NEWOBJMASK)
#define SWIG_TMPOBJ                (SWIG_OK | SWIG_TMPOBJMASK)
/* Check, add and del mask methods */
#define SWIG_AddNewMask(r)         (SWIG_IsOK(r) ? (r | SWIG_NEWOBJMASK) : r)
#define SWIG_DelNewMask(r)         (SWIG_IsOK(r) ? (r & ~SWIG_NEWOBJMASK) : r)
#define SWIG_IsNewObj(r)           (SWIG_IsOK(r) && (r & SWIG_NEWOBJMASK))
#define SWIG_AddTmpMask(r)         (SWIG_IsOK(r) ? (r | SWIG_TMPOBJMASK) : r)
#define SWIG_DelTmpMask(r)         (SWIG_IsOK(r) ? (r & ~SWIG_TMPOBJMASK) : r)
#define SWIG_IsTmpObj(r)           (SWIG_IsOK(r) && (r & SWIG_TMPOBJMASK))

/* Cast-Rank Mode */
#if defined(SWIG_CASTRANK_MODE)
#  ifndef SWIG_TypeRank
#    define SWIG_TypeRank             unsigned long
#  endif
#  ifndef SWIG_MAXCASTRANK            /* Default cast allowed */
#    define SWIG_MAXCASTRANK          (2)
#  endif
#  define SWIG_CASTRANKMASK          ((SWIG_CASTRANKLIMIT) -1)
#  define SWIG_CastRank(r)           (r & SWIG_CASTRANKMASK)
SWIGINTERNINLINE int SWIG_AddCast(int r) {
  return SWIG_IsOK(r) ? ((SWIG_CastRank(r) < SWIG_MAXCASTRANK) ? (r + 1) : SWIG_ERROR) : r;
}
SWIGINTERNINLINE int SWIG_CheckState(int r) {
  return SWIG_IsOK(r) ? SWIG_CastRank(r) + 1 : 0;
}
#else /* no cast-rank mode */
#  define SWIG_AddCast(r) (r)
#  define SWIG_CheckState(r) (SWIG_IsOK(r) ? 1 : 0)
#endif


#include <string.h>

#ifdef __cplusplus
extern "C" {
#endif

typedef void *(*swig_converter_func)(void *, int *);
typedef struct swig_type_info *(*swig_dycast_func)(void **);

/* Structure to store information on one type */
typedef struct swig_type_info {
  const char             *name;			/* mangled name of this type */
  const char             *str;			/* human readable name of this type */
  swig_dycast_func        dcast;		/* dynamic cast function down a hierarchy */
  struct swig_cast_info  *cast;			/* linked list of types that can cast into this type */
  void                   *clientdata;		/* language specific type data */
  int                    owndata;		/* flag if the structure owns the clientdata */
} swig_type_info;

/* Structure to store a type and conversion function used for casting */
typedef struct swig_cast_info {
  swig_type_info         *type;			/* pointer to type that is equivalent to this type */
  swig_converter_func     converter;		/* function to cast the void pointers */
  struct swig_cast_info  *next;			/* pointer to next cast in linked list */
  struct swig_cast_info  *prev;			/* pointer to the previous cast */
} swig_cast_info;

/* Structure used to store module information
 * Each module generates one structure like this, and the runtime collects
 * all of these structures and stores them in a circularly linked list.*/
typedef struct swig_module_info {
  swig_type_info         **types;		/* Array of pointers to swig_type_info structures that are in this module */
  size_t                 size;		        /* Number of types in this module */
  struct swig_module_info *next;		/* Pointer to next element in circularly linked list */
  swig_type_info         **type_initial;	/* Array of initially generated type structures */
  swig_cast_info         **cast_initial;	/* Array of initially generated casting structures */
  void                    *clientdata;		/* Language specific module data */
} swig_module_info;

/*
  Compare two type names skipping the space characters, therefore
  "char*" == "char *" and "Class<int>" == "Class<int >", etc.

  Return 0 when the two name types are equivalent, as in
  strncmp, but skipping ' '.
*/
SWIGRUNTIME int
SWIG_TypeNameComp(const char *f1, const char *l1,
		  const char *f2, const char *l2) {
  for (;(f1 != l1) && (f2 != l2); ++f1, ++f2) {
    while ((*f1 == ' ') && (f1 != l1)) ++f1;
    while ((*f2 == ' ') && (f2 != l2)) ++f2;
    if (*f1 != *f2) return (*f1 > *f2) ? 1 : -1;
  }
  return (int)((l1 - f1) - (l2 - f2));
}

/*
  Check type equivalence in a name list like <name1>|<name2>|...
  Return 0 if equal, -1 if nb < tb, 1 if nb > tb
*/
SWIGRUNTIME int
SWIG_TypeCmp(const char *nb, const char *tb) {
  int equiv = 1;
  const char* te = tb + strlen(tb);
  const char* ne = nb;
  while (equiv != 0 && *ne) {
    for (nb = ne; *ne; ++ne) {
      if (*ne == '|') break;
    }
    equiv = SWIG_TypeNameComp(nb, ne, tb, te);
    if (*ne) ++ne;
  }
  return equiv;
}

/*
  Check type equivalence in a name list like <name1>|<name2>|...
  Return 0 if not equal, 1 if equal
*/
SWIGRUNTIME int
SWIG_TypeEquiv(const char *nb, const char *tb) {
  return SWIG_TypeCmp(nb, tb) == 0 ? 1 : 0;
}

/*
  Check the typename
*/
SWIGRUNTIME swig_cast_info *
SWIG_TypeCheck(const char *c, swig_type_info *ty) {
  if (ty) {
    swig_cast_info *iter = ty->cast;
    while (iter) {
      if (strcmp(iter->type->name, c) == 0) {
        if (iter == ty->cast)
          return iter;
        /* Move iter to the top of the linked list */
        iter->prev->next = iter->next;
        if (iter->next)
          iter->next->prev = iter->prev;
        iter->next = ty->cast;
        iter->prev = 0;
        if (ty->cast) ty->cast->prev = iter;
        ty->cast = iter;
        return iter;
      }
      iter = iter->next;
    }
  }
  return 0;
}

/*
  Identical to SWIG_TypeCheck, except strcmp is replaced with a pointer comparison
*/
SWIGRUNTIME swig_cast_info *
SWIG_TypeCheckStruct(swig_type_info *from, swig_type_info *ty) {
  if (ty) {
    swig_cast_info *iter = ty->cast;
    while (iter) {
      if (iter->type == from) {
        if (iter == ty->cast)
          return iter;
        /* Move iter to the top of the linked list */
        iter->prev->next = iter->next;
        if (iter->next)
          iter->next->prev = iter->prev;
        iter->next = ty->cast;
        iter->prev = 0;
        if (ty->cast) ty->cast->prev = iter;
        ty->cast = iter;
        return iter;
      }
      iter = iter->next;
    }
  }
  return 0;
}

/*
  Cast a pointer up an inheritance hierarchy
*/
SWIGRUNTIMEINLINE void *
SWIG_TypeCast(swig_cast_info *ty, void *ptr, int *newmemory) {
  return ((!ty) || (!ty->converter)) ? ptr : (*ty->converter)(ptr, newmemory);
}

/*
   Dynamic pointer casting. Down an inheritance hierarchy
*/
SWIGRUNTIME swig_type_info *
SWIG_TypeDynamicCast(swig_type_info *ty, void **ptr) {
  swig_type_info *lastty = ty;
  if (!ty || !ty->dcast) return ty;
  while (ty && (ty->dcast)) {
    ty = (*ty->dcast)(ptr);
    if (ty) lastty = ty;
  }
  return lastty;
}

/*
  Return the name associated with this type
*/
SWIGRUNTIMEINLINE const char *
SWIG_TypeName(const swig_type_info *ty) {
  return ty->name;
}

/*
  Return the pretty name associated with this type,
  that is an unmangled type name in a form presentable to the user.
*/
SWIGRUNTIME const char *
SWIG_TypePrettyName(const swig_type_info *type) {
  /* The "str" field contains the equivalent pretty names of the
     type, separated by vertical-bar characters.  We choose
     to print the last name, as it is often (?) the most
     specific. */
  if (!type) return NULL;
  if (type->str != NULL) {
    const char *last_name = type->str;
    const char *s;
    for (s = type->str; *s; s++)
      if (*s == '|') last_name = s+1;
    return last_name;
  }
  else
    return type->name;
}

/*
   Set the clientdata field for a type
*/
SWIGRUNTIME void
SWIG_TypeClientData(swig_type_info *ti, void *clientdata) {
  swig_cast_info *cast = ti->cast;
  /* if (ti->clientdata == clientdata) return; */
  ti->clientdata = clientdata;

  while (cast) {
    if (!cast->converter) {
      swig_type_info *tc = cast->type;
      if (!tc->clientdata) {
	SWIG_TypeClientData(tc, clientdata);
      }
    }
    cast = cast->next;
  }
}
SWIGRUNTIME void
SWIG_TypeNewClientData(swig_type_info *ti, void *clientdata) {
  SWIG_TypeClientData(ti, clientdata);
  ti->owndata = 1;
}

/*
  Search for a swig_type_info structure only by mangled name
  Search is a O(log #types)

  We start searching at module start, and finish searching when start == end.
  Note: if start == end at the beginning of the function, we go all the way around
  the circular list.
*/
SWIGRUNTIME swig_type_info *
SWIG_MangledTypeQueryModule(swig_module_info *start,
                            swig_module_info *end,
		            const char *name) {
  swig_module_info *iter = start;
  do {
    if (iter->size) {
      size_t l = 0;
      size_t r = iter->size - 1;
      do {
	/* since l+r >= 0, we can (>> 1) instead (/ 2) */
	size_t i = (l + r) >> 1;
	const char *iname = iter->types[i]->name;
	if (iname) {
	  int compare = strcmp(name, iname);
	  if (compare == 0) {
	    return iter->types[i];
	  } else if (compare < 0) {
	    if (i) {
	      r = i - 1;
	    } else {
	      break;
	    }
	  } else if (compare > 0) {
	    l = i + 1;
	  }
	} else {
	  break; /* should never happen */
	}
      } while (l <= r);
    }
    iter = iter->next;
  } while (iter != end);
  return 0;
}

/*
  Search for a swig_type_info structure for either a mangled name or a human readable name.
  It first searches the mangled names of the types, which is a O(log #types)
  If a type is not found it then searches the human readable names, which is O(#types).

  We start searching at module start, and finish searching when start == end.
  Note: if start == end at the beginning of the function, we go all the way around
  the circular list.
*/
SWIGRUNTIME swig_type_info *
SWIG_TypeQueryModule(swig_module_info *start,
                     swig_module_info *end,
		     const char *name) {
  /* STEP 1: Search the name field using binary search */
  swig_type_info *ret = SWIG_MangledTypeQueryModule(start, end, name);
  if (ret) {
    return ret;
  } else {
    /* STEP 2: If the type hasn't been found, do a complete search
       of the str field (the human readable name) */
    swig_module_info *iter = start;
    do {
      size_t i = 0;
      for (; i < iter->size; ++i) {
	if (iter->types[i]->str && (SWIG_TypeEquiv(iter->types[i]->str, name)))
	  return iter->types[i];
      }
      iter = iter->next;
    } while (iter != end);
  }

  /* neither found a match */
  return 0;
}

/*
   Pack binary data into a string
*/
SWIGRUNTIME char *
SWIG_PackData(char *c, void *ptr, size_t sz) {
  static const char hex[17] = "0123456789abcdef";
  const unsigned char *u = (unsigned char *) ptr;
  const unsigned char *eu =  u + sz;
  for (; u != eu; ++u) {
    unsigned char uu = *u;
    *(c++) = hex[(uu & 0xf0) >> 4];
    *(c++) = hex[uu & 0xf];
  }
  return c;
}

/*
   Unpack binary data from a string
*/
SWIGRUNTIME const char *
SWIG_UnpackData(const char *c, void *ptr, size_t sz) {
  unsigned char *u = (unsigned char *) ptr;
  const unsigned char *eu = u + sz;
  for (; u != eu; ++u) {
    char d = *(c++);
    unsigned char uu;
    if ((d >= '0') && (d <= '9'))
      uu = ((d - '0') << 4);
    else if ((d >= 'a') && (d <= 'f'))
      uu = ((d - ('a'-10)) << 4);
    else
      return (char *) 0;
    d = *(c++);
    if ((d >= '0') && (d <= '9'))
      uu |= (d - '0');
    else if ((d >= 'a') && (d <= 'f'))
      uu |= (d - ('a'-10));
    else
      return (char *) 0;
    *u = uu;
  }
  return c;
}

/*
   Pack 'void *' into a string buffer.
*/
SWIGRUNTIME char *
SWIG_PackVoidPtr(char *buff, void *ptr, const char *name, size_t bsz) {
  char *r = buff;
  if ((2*sizeof(void *) + 2) > bsz) return 0;
  *(r++) = '_';
  r = SWIG_PackData(r,&ptr,sizeof(void *));
  if (strlen(name) + 1 > (bsz - (r - buff))) return 0;
  strcpy(r,name);
  return buff;
}

SWIGRUNTIME const char *
SWIG_UnpackVoidPtr(const char *c, void **ptr, const char *name) {
  if (*c != '_') {
    if (strcmp(c,"NULL") == 0) {
      *ptr = (void *) 0;
      return name;
    } else {
      return 0;
    }
  }
  return SWIG_UnpackData(++c,ptr,sizeof(void *));
}

SWIGRUNTIME char *
SWIG_PackDataName(char *buff, void *ptr, size_t sz, const char *name, size_t bsz) {
  char *r = buff;
  size_t lname = (name ? strlen(name) : 0);
  if ((2*sz + 2 + lname) > bsz) return 0;
  *(r++) = '_';
  r = SWIG_PackData(r,ptr,sz);
  if (lname) {
    strncpy(r,name,lname+1);
  } else {
    *r = 0;
  }
  return buff;
}

SWIGRUNTIME const char *
SWIG_UnpackDataName(const char *c, void *ptr, size_t sz, const char *name) {
  if (*c != '_') {
    if (strcmp(c,"NULL") == 0) {
      memset(ptr,0,sz);
      return name;
    } else {
      return 0;
    }
  }
  return SWIG_UnpackData(++c,ptr,sz);
}

#ifdef __cplusplus
}
#endif

/*  Errors in SWIG */
#define  SWIG_UnknownError    	   -1
#define  SWIG_IOError        	   -2
#define  SWIG_RuntimeError   	   -3
#define  SWIG_IndexError     	   -4
#define  SWIG_TypeError      	   -5
#define  SWIG_DivisionByZero 	   -6
#define  SWIG_OverflowError  	   -7
#define  SWIG_SyntaxError    	   -8
#define  SWIG_ValueError     	   -9
#define  SWIG_SystemError    	   -10
#define  SWIG_AttributeError 	   -11
#define  SWIG_MemoryError    	   -12
#define  SWIG_NullReferenceError   -13



/* Last error */
static int SWIG_lasterror_code = 0;
static char SWIG_lasterror_msg[1024];
SWIGRUNTIME void SWIG_Error(int code, const char *format, ...) {
  SWIG_lasterror_code = code;
  va_list arg;
  va_start(arg, format);
  vsnprintf(SWIG_lasterror_msg, sizeof(SWIG_lasterror_msg), format, arg);
  va_end(arg);
}

SWIGRUNTIME const char* SWIG_ErrorType(int code) {
  switch (code) {
  case SWIG_MemoryError:
    return "SWIG:MemoryError";
  case SWIG_IOError:
    return "SWIG:IOError";
  case SWIG_RuntimeError:
    return "SWIG:RuntimeError";
  case SWIG_IndexError:
    return "SWIG:IndexError";
  case SWIG_TypeError:
    return "SWIG:TypeError";
  case SWIG_DivisionByZero:
    return "SWIG:DivisionByZero";
  case SWIG_OverflowError:
    return "SWIG:OverflowError";
  case SWIG_SyntaxError:
    return "SWIG:SyntaxError";
  case SWIG_ValueError:
    return "SWIG:ValueError";
  case SWIG_SystemError:
    return "SWIG:SystemError";
  case SWIG_AttributeError:
    return "SWIG:AttributeError";
  }
  return "SWIG:UnknownError";
}

SWIGRUNTIME bool SWIG_check_num_args(const char *func_name, int num_args, int max_args, int min_args, int varargs) {
  if (num_args > max_args && !varargs) {
    SWIG_Error(SWIG_RuntimeError, "function %s takes at most %i arguments", func_name, max_args);
  } else if (num_args < min_args) {
    SWIG_Error(SWIG_RuntimeError, "function %s requires at least %i arguments", func_name, min_args);
  } else {
    return true;
  }
  return false;
}

#define SWIG_fail                                       goto fail

#define SWIG_Matlab_ConvertPtr(obj, pptr, type, flags)  SWIG_Matlab_ConvertPtrAndOwn(obj, pptr, type, flags, 0)
#define SWIG_ConvertPtr(obj, pptr, type, flags)         SWIG_Matlab_ConvertPtr(obj, pptr, type, flags)
#define SWIG_ConvertPtrAndOwn(obj,pptr,type,flags,own)  SWIG_Matlab_ConvertPtrAndOwn(obj, pptr, type, flags, own)
#define SWIG_ConvertPtr(obj, pptr, type, flags)         SWIG_Matlab_ConvertPtr(obj, pptr, type, flags)
#define SWIG_NewPointerObj(ptr, type, flags)            SWIG_Matlab_NewPointerObj(ptr, type, flags)
#define swig_owntype                                    int

#define SWIG_ConvertPacked(obj, ptr, sz, ty)            SWIG_Matlab_ConvertPacked(obj, ptr, sz, ty)
#define SWIG_NewPackedObj(ptr, sz, type)                SWIG_Matlab_NewPackedObj(ptr, sz, type)

#define SWIG_ConvertFunctionPtr(obj, pptr, type)        SWIG_ConvertPtr(obj, pptr, type, 0)
#define SWIG_NewFunctionPtrObj(ptr, type)               SWIG_NewPointerObj(ptr, type, 0)

#define SWIG_ConvertMember(obj, ptr, sz, ty)            SWIG_Matlab_ConvertPacked(obj, ptr, sz, ty)
#define SWIG_NewMemberObj(ptr, sz, type)                SWIG_Matlab_NewPackedObj(ptr, sz, type)

#define SWIG_GetModule(clientdata)                      SWIG_Matlab_GetModule(clientdata)
#define SWIG_SetModule(clientdata, pointer)             SWIG_Matlab_SetModule(pointer);
#define SWIG_MODULE_CLIENTDATA_TYPE void*

#define Matlab_Error_Occurred() 0
#define SWIG_Matlab_AddErrorMsg(msg) {;}

SWIGRUNTIME swig_module_info *SWIG_Matlab_GetModule(void *clientdata);
SWIGRUNTIME void SWIG_Matlab_SetModule(swig_module_info *pointer);

/* For backward compatibility only */
#define SWIG_POINTER_EXCEPTION  0
#define SWIG_arg_fail(arg)      0

/* Pointer to C/C++ object, including type information and reference counters */
typedef struct SwigPtr {
  void *ptr;
  swig_type_info *type;
  int own;
  void *self;
  struct SwigPtr* next;
} SwigPtr;
static swig_module_info *saved_swig_module = 0;

/* Add new pointer object to the memory allocator
   Returns the location of the new element by reference and a return flag by value.
 */
int SWIG_Matlab_NewPointer(SwigPtr** swig_ptr, void *ptr, swig_type_info *type, int own) {
  /* Check input arguments */
  if (!swig_ptr || !type) return 1;

  /* Allocate memory */
  SwigPtr *obj=(SwigPtr *)malloc(sizeof(SwigPtr));
  if (!obj) return 1;

  /* Initialize object */
  obj->ptr = ptr;
  obj->type = type;
  obj->own = own;
  obj->self = obj;
  obj->next = 0;
  
  /* Return memory reference and successful return */
  *swig_ptr = obj;
  return 0;
}

/* Free pointer object
   Returns a void* with the object to be freed (if any)
 */
void* SWIG_Matlab_DeletePointer(SwigPtr* p) {
  /* Quick return if null pointer? */
  if (p==0) return 0;

  /* Clear and free */
  void* ret = p->own ? p->ptr : 0;
  p->ptr=0;
  p->type=0;
  p->own=0;
  p->self=0;
  p->next=0;
  free(p);
  return ret;
}

/* This should probably be moved somewhere else */
#ifdef __cplusplus
extern "C"
#endif /* cplusplus */
void SWIG_Matlab_ExitFcn(void) {
  /* Unused at the moment */
}

/* Runtime API implementation */
SWIGRUNTIME mxArray* SWIG_Matlab_NewPointerObj(void *ptr, swig_type_info *type, int flags) {
  /* Ownership? */
  int own = (flags & SWIG_POINTER_OWN) ? SWIG_POINTER_OWN : 0;

  /* Allocate a pointer object */  
  SwigPtr* swig_ptr;
  if (SWIG_Matlab_NewPointer(&swig_ptr, ptr, type, own)) {
    mexErrMsgIdAndTxt("SWIG:NewPointerObj","Cannot allocate pointer");
  }

  /* Create a SwigRef instance */
  mxArray* pm = 0;
  mxArray* ex = mexCallMATLABWithTrap(1, &pm, 0, 0, "SwigRef");
  if(ex) mexErrMsgIdAndTxt("SWIG:NewPointerObj", "Cannot create SwigRef instance");

  /* Set pointer property */
  mxArray *pm_ptr = mxCreateNumericMatrix(1, 1, mxUINT64_CLASS, mxREAL);
  *(uint64_T *)mxGetData(pm_ptr) = (uint64_T)swig_ptr;
  if(!pm_ptr) mexErrMsgIdAndTxt("SWIG:NewPointerObj", "mxCreateNumericMatrix failed");
  mxSetProperty(pm, 0, "swigPtr", pm_ptr);
  mxDestroyArray(pm_ptr);

  /* Create a proxy class of the desired type */
  if (type->clientdata) {
    const char* matlab_classname = (const char *)type->clientdata;
    mxArray* pm_test = 0;
    ex = mexCallMATLABWithTrap(1, &pm_test, 1, &pm, matlab_classname);
    if (ex) {
      mexWarnMsgIdAndTxt("SWIG:NewPointerObj",
                         "NewPointerObj problem. creating SwigRef as opposed to %s", matlab_classname);
      mxDestroyArray(ex);
    } else {
      mxDestroyArray(pm);
      pm = pm_test;
    }
  }
  
  return pm;
}

SWIGRUNTIME int SWIG_Matlab_isOwned(mxArray* pm) {
/* Get index, represented as a uint64_T scalar */
  mxArray *pm_ptr = mxGetProperty(pm, 0, "swigPtr");
  if(!pm_ptr || mxGetNumberOfElements(pm_ptr) != 1 ||
     mxGetClassID(pm_ptr) != mxUINT64_CLASS || mxIsComplex(pm_ptr)) {
    if(pm_ptr) mxDestroyArray(pm_ptr);
    return SWIG_ERROR;
  }
  SwigPtr* p = (SwigPtr*)(*(uint64_T*)mxGetData(pm_ptr));
  mxDestroyArray(pm_ptr);

  return p ? p->own : 0;
}

SWIGRUNTIME SwigPtr* SWIG_Matlab_getSwigPtr(mxArray* pm) {
/* Get index, represented as a uint64_T scalar */
  mxArray *pm_ptr = mxGetProperty(pm, 0, "swigPtr");
  if(!pm_ptr || mxGetNumberOfElements(pm_ptr) != 1 ||
     mxGetClassID(pm_ptr) != mxUINT64_CLASS || mxIsComplex(pm_ptr)) {
    if(pm_ptr) mxDestroyArray(pm_ptr);
    return 0;
  }
  SwigPtr* p = (SwigPtr*)(*(uint64_T*)mxGetData(pm_ptr));
  mxDestroyArray(pm_ptr);

  /* Consistency check */
  if (p->self != p) return 0;

  /* Get pointer object */
  return p;
}

SWIGRUNTIME int SWIG_Matlab_ConvertPtrAndOwn(mxArray* pm, void **ptr, swig_type_info *type, int flags, int *own) {
  /* Get index, represented as a uint64_T scalar */
  mxArray *pm_ptr = mxGetProperty(pm, 0, "swigPtr");
  if(!pm_ptr || mxGetNumberOfElements(pm_ptr) != 1 ||
     mxGetClassID(pm_ptr) != mxUINT64_CLASS || mxIsComplex(pm_ptr)){
    if(pm_ptr) mxDestroyArray(pm_ptr);
    return SWIG_ERROR;
  }
  SwigPtr* swig_ptr = (SwigPtr*)(*(uint64_T*)mxGetData(pm_ptr));
  mxDestroyArray(pm_ptr);

  /* Get pointer object */
  if (!swig_ptr) return SWIG_ERROR;

  swig_cast_info *tc;
  bool same_type = false;
  SwigPtr* p=swig_ptr;
  while (p) {
    /* Loop over types */
    swig_type_info *p_type = p->type;

    /* No type cast needed */
    if (p_type==type || !type) {
      same_type = true;
      break;
    }

    /* Check if type conversion is possible */
    tc = SWIG_TypeCheck(p_type->name, type);
    if (tc) break;

    /* Check next type */
    p = p->next;
  }

  /* Conversion not possible */
  if (!same_type && !tc) return SWIG_ERROR;

  /* Get data */
  if (ptr) *ptr = p->ptr;
  if (own) *own = p->own;
    
  /* Set ownership marker to false? */
  if(flags & SWIG_POINTER_DISOWN){
    SWIG_Matlab_DeletePointer(swig_ptr);
  }

  /* No type cast needed */
  if (same_type || !ptr) {
    return SWIG_OK;
  }

  /* Do type-cast */
  int newmemory = 0;
  *ptr = SWIG_TypeCast(tc, *ptr, &newmemory);
  if (newmemory == SWIG_CAST_NEW_MEMORY) {
    if (own) {
      *own = *own | SWIG_CAST_NEW_MEMORY;
    } else {
      /* badly formed typemap which will lead to a memory leak - it must set and use own to delete *ptr */
      mexErrMsgIdAndTxt("SWIG:ConvertPointerObj","SWIG error in converting pointer for type %s. Bad typemap?",type->name);
      return SWIG_ERROR;
    }
  }
  return SWIG_OK;
}

SWIGRUNTIME swig_module_info *SWIG_Matlab_GetModule(void *SWIGUNUSEDPARM(clientdata)) {
  if (!saved_swig_module) {
    /* First call for this module, fetch pointer from MATLAB */
    mxArray* mem = 0;
    int flag = mexCallMATLAB(1, &mem, 0, 0, "SwigMem");
    if (flag || !mem) mexErrMsgIdAndTxt("SWIG:GetModule", "Cannot call SwigMem");
    /* On first call, mem is an empty matrix */
    if (mxGetNumberOfElements(mem)!= 0) {
      if (!mxIsStruct(mem)) {
        mexErrMsgIdAndTxt("SWIG:GetModule", "Corrupted memory");
      }
      mxArray *module = mxGetField(mem, 0,
        "swig_runtime_data_type_pointer" SWIG_RUNTIME_VERSION SWIG_TYPE_TABLE_NAME);
      if (module && mxGetNumberOfElements(module) == 1 && mxGetClassID(module) == mxUINT64_CLASS && !mxIsComplex(module)) {
        /* Cast memory as a pointer */
        saved_swig_module = (swig_module_info*)(*(uint64_T*)mxGetData(module));
      }
    }
    mxDestroyArray(mem);
  }
  return saved_swig_module;
}

SWIGRUNTIME void SWIG_Matlab_SetModule(swig_module_info *pointer) {
  mxArray* mem = 0;
  const char *fields[1] = {"swig_runtime_data_type_pointer" SWIG_RUNTIME_VERSION SWIG_TYPE_TABLE_NAME};
  int flag = mexCallMATLAB(1, &mem, 0, 0, "SwigMem");
  if (flag || !mem) mexErrMsgIdAndTxt("SWIG:SetModule", "Cannot call SwigMem");
  /* Under normal circumstances, swigMem is an empty matrix */
  /* but is not required to be for a valid call to SWIG_Matlab_SetModule() */
  if (mxGetNumberOfElements(mem) == 0) {
    mxDestroyArray(mem);
    mem = mxCreateStructMatrix(1,1,1,fields);
    if(!mem) mexErrMsgIdAndTxt("SWIG:SetModule","mxCreateStructMatrix failed");
    mxArray *module = mxCreateNumericMatrix(1, 1, mxUINT64_CLASS, mxREAL);
    if(!module) mexErrMsgIdAndTxt("SWIG:SetModule","mxCreateNumericMatrix failed");
    *(uint64_T *)mxGetData(module) = (uint64_T)pointer;
    mxSetField(mem, 0, fields[0], module);
  } else {
    int nfields = 0;
    int fieldNum = 0;
    /* Add to existing mem if not the same */
    if (!mxIsStruct(mem)) {
      mexErrMsgIdAndTxt("SWIG:SetModule", "Corrupted memory");
    }
    nfields = mxGetNumberOfFields(mem);
    /* Check if same type table */
    for (int ii=0; ii<nfields; ++ii){
      if(strcmp(mxGetFieldNameByNumber(mem, ii), fields[0]) == 0) {
        mxDestroyArray(mem);
        return;
      }
    }
    fieldNum = mxAddField(mem, fields[0]);
    if (fieldNum < 0) mexErrMsgIdAndTxt("SWIG:SetModule", "Error adding field to SwigMem");
    mxArray *module = mxCreateNumericMatrix(1, 1, mxUINT64_CLASS, mxREAL);
    if(!module) mexErrMsgIdAndTxt("SWIG:SetModule","mxCreateNumericMatrix failed");
    *(uint64_T *)mxGetData(module) = (uint64_T)pointer;
    mxSetFieldByNumber(mem, 0, fieldNum, module);
  }
  /* Set variable in SwigMem */
  flag = mexCallMATLAB(0, 0, 1, &mem, "SwigMem");
  mxDestroyArray(mem);
  if (flag) {
    mexErrMsgIdAndTxt("SWIG:SetModule", "Could not set variable in SwigMem.\n");
  }
  saved_swig_module = pointer;
}

/* HACK HACK */
#define SWIG_Matlab_SetConstant(dummy1,dummy2,pm) (pm)

#ifdef __cplusplus
#include <streambuf>
#include <iostream>
namespace swig {
  // Stream buffer to allow redirecting output to MATLAB
  class SWIG_Matlab_streambuf : public std::streambuf {
  public:
    SWIG_Matlab_streambuf() {}
  protected:
    virtual int_type overflow(int_type ch) {
      if(ch != traits_type::eof()) {
        mexPrintf("%c", static_cast<char>(ch));
      }
      return ch;
    }
    virtual std::streamsize xsputn(const char* s, std::streamsize num) {
      // Pass straight to mexPrintf without buffering
      mexPrintf("%.*s", static_cast<int>(num), s);
      return num;
    }
  };

  // Instantiation
  static SWIG_Matlab_streambuf SWIG_Matlab_buf;
} // namespace swig
#endif /* cplusplus */



#define SWIG_exception_fail(code, msg) do { SWIG_Error(code, msg); SWIG_fail; } while(0) 

#define SWIG_contract_assert(expr, msg) if (!(expr)) { SWIG_Error(SWIG_RuntimeError, msg); SWIG_fail; } else 



  #define SWIG_exception(code, msg) do { SWIG_Error(code, msg); SWIG_fail;; } while(0) 


/* -------- TYPES TABLE (BEGIN) -------- */

#define SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t swig_types[0]
#define SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t swig_types[1]
#define SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t swig_types[2]
#define SWIGTYPE_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_t_Eigen__Dynamic_1_t swig_types[3]
#define SWIGTYPE_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_1_t swig_types[4]
#define SWIGTYPE_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_Eigen__Dynamic_t swig_types[5]
#define SWIGTYPE_p_char swig_types[6]
static swig_type_info *swig_types[8];
static swig_module_info swig_module = {swig_types, 7, 0, 0, 0, 0};
#define SWIG_TypeQuery(name) SWIG_TypeQueryModule(&swig_module, &swig_module, name)
#define SWIG_MangledTypeQuery(name) SWIG_MangledTypeQueryModule(&swig_module, &swig_module, name)

/* -------- TYPES TABLE (END) -------- */


#define SWIGVERSION 0x030009 
#define SWIG_VERSION SWIGVERSION


#define SWIG_as_voidptr(a) const_cast< void * >(static_cast< const void * >(a)) 
#define SWIG_as_voidptrptr(a) ((void)SWIG_as_voidptr(*a),reinterpret_cast< void** >(a)) 


#include <stdexcept>


  #include <Eigen/Core>
  #include <vector>


#ifdef SWIGPYTHON
  #define SWIG_FILE_WITH_INIT
  #include <Python.h>
#endif
#include <unsupported/Eigen/AutoDiff>
#include <iostream>



template <typename DerType, int RowsAtCompileTime, int ColsAtCompileTime>
class AutoDiffWrapper: public Eigen::Matrix<Eigen::AutoDiffScalar<DerType>, RowsAtCompileTime, ColsAtCompileTime> {

  typedef Eigen::Matrix<Eigen::AutoDiffScalar<DerType>, RowsAtCompileTime, ColsAtCompileTime> BaseMatrix;

public:
  AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime>() {}

  AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime>(const Eigen::Matrix<Eigen::AutoDiffScalar<DerType>, RowsAtCompileTime, ColsAtCompileTime>& x): Eigen::Matrix<Eigen::AutoDiffScalar<DerType>, RowsAtCompileTime, ColsAtCompileTime>(x) {};

  AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime>(const Eigen::Matrix<double, RowsAtCompileTime, ColsAtCompileTime> &value, const Eigen::MatrixXd &derivatives) {
    this->resize(value.rows(), value.cols());
    if (derivatives.rows() != value.size()) {
      throw std::runtime_error("derivatives must have one row for every element in value");
    }
    for (size_t i=0; i < value.rows(); i++) {
      for (size_t j=0; j < value.cols(); j++) {
        this->coeffRef(i, j) = Eigen::AutoDiffScalar<Eigen::VectorXd>(value(i,j), derivatives.row(i + value.rows() * j));
      }
    }
  }

  Eigen::Matrix<double, RowsAtCompileTime, ColsAtCompileTime> value() {
    Eigen::Matrix<double, RowsAtCompileTime, ColsAtCompileTime> result(this->rows(), this->cols());
    for (size_t i=0; i < this->size(); i++) {
      result(i) = this->coeffRef(i).value();
    }
    return result;
  }

  Eigen::MatrixXd derivatives() {
    Eigen::MatrixXd result(this->size(), this->size() > 0 ? this->coeffRef(0).derivatives().size() : 0);
    for (int i=0; i < this->size(); i++) {
      if (this->coeffRef(i).derivatives().size() > result.cols()) {
        result.conservativeResize(result.rows(), this->coeffRef(i).derivatives().size());
      }
      for (int j=0; j < this->coeffRef(i).derivatives().size(); j++) {
        result(i, j) = this->coeffRef(i).derivatives()(j);
      }
    }
    return result;
  }

  AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime> operator+ (const AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime>& other) {
    return BaseMatrix::operator+(other).eval();
  }
  AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime> operator- (const AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime>& other) {
    return BaseMatrix::operator-(other).eval();
  }
  AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime> arrayMultiply (const AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime>& other) {
    return this->array().operator*(other.array()).matrix().eval();
  }
  AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime> arrayDivide (const AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime>& other) {
    return this->array().operator/(other.array()).matrix().eval();
  }

  AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime> operator+ (double other) {
    return this->array().operator+(other).matrix().eval();
  }
  AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime> operator- (double other) {
    return this->array().operator-(other).matrix().eval();
  }
  AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime> operator* (double other) {
    return BaseMatrix::operator*(other).eval();
  }
  AutoDiffWrapper<DerType, RowsAtCompileTime, ColsAtCompileTime> operator/ (double other) {
    return BaseMatrix::operator/(other).eval();
  }

  int rows() const {
    return BaseMatrix::rows();
  }

  int cols() const {
    return BaseMatrix::cols();
  }

  int size() const {
    return BaseMatrix::size();
  }
};



  template <typename T> mxClassID MxType() {return mxUNKNOWN_CLASS;};

  template <class Derived>
  bool ConvertFromMatlabToEigenMatrix(Eigen::MatrixBase<Derived>* out, const mxArray* in)
  {
    size_t rows = 0;
    size_t cols = 0;

    // Do a bunch of validation and type checking
    if (!mxIsNumeric(in)) {
      mexPrintf("input is not a matlab numeric array");
      return false;
    } // check data type
    else if (mxGetClassID(in) != MxType<typename Derived::Scalar>()) {
      mexPrintf("Type mismatch between Matlab and Eigen objects.");
      return false;
    }
    // check dimensions
    else if (mxGetNumberOfDimensions(in) > 2) {
      mexPrintf("Eigen only supports 1D or 2D arrays, but I got a Matlab array with more than 2 dimensions");
      return false;
    }
    else if (mxGetNumberOfDimensions(in) == 1) {
      rows = mxGetNumberOfElements(in);
      cols = 1;
      if ((Derived::RowsAtCompileTime != Eigen::Dynamic) && (Derived::RowsAtCompileTime != rows)) {
        mexPrintf("Row dimension mismatch between Matlab and Eigen objects (1D).");
        return false;
      }
      else if ((Derived::ColsAtCompileTime != Eigen::Dynamic) && (Derived::ColsAtCompileTime != 1)) {
        mexPrintf("Column dimension mismatch between Matlab and Eigen objects (1D).");
        return false;
      }
    }
    else if (mxGetNumberOfDimensions(in) == 2) {
      rows = mxGetM(in);
      cols = mxGetN(in);
      if ((Derived::RowsAtCompileTime != Eigen::Dynamic) && (Derived::RowsAtCompileTime != rows)) {
        mexPrintf("Row dimension mismatch between Matlab and Eigen objects (2D).");
        return false;
      }
      else if ((Derived::ColsAtCompileTime != Eigen::Dynamic) && (Derived::ColsAtCompileTime != cols)) {
        mexPrintf("Column dimension mismatch between Matlab and Eigen objects (2D).");
        return false;
      }
    }

    // Extract the data
    out->derived().setZero(rows, cols);
    typename Derived::Scalar* data = static_cast<typename Derived::Scalar*>(mxGetData(in));
    for (size_t i=0; i < rows; i++) {
      for (size_t j=0; j < cols; j++) {
        out->coeffRef(i,j) = data[i + rows*j];
      }
    }

    return true;
  };

  // Copy values from an Eigen matrix into an *existing* Matlab matrix
  template<class Derived>
  bool CopyFromEigenToMatlabMatrix(mxArray* out, Eigen::MatrixBase<Derived>* in) {
    size_t rows;
    size_t cols;

    // check object type
    if (!mxIsNumeric(out)) {
      mexPrintf("Argout matrix must be a numeric matlab matrix");
      return false;
    }
    // check scalar type
    else if (mxGetClassID(out) != MxType<typename Derived::Scalar>()) {
      mexPrintf("Argout matrix data type does not match the scalar type of the Eigen matrix");
      return false;
    }
    else if (mxGetNumberOfDimensions(out) > 2) {
      mexPrintf("I can't handle Matlab matrices with more than 2 dimensions");
      return false;
    }
    rows = mxGetM(out);
    cols = mxGetN(out);

    if ((Derived::RowsAtCompileTime != Eigen::Dynamic) && (Derived::RowsAtCompileTime != rows))
    {
      mexPrintf("Row dimension mismatch between Matlab and Eigen objects");
      return false;
    }
    else if ((Derived::ColsAtCompileTime != Eigen::Dynamic) && (Derived::ColsAtCompileTime != cols))
    {
      mexPrintf("Column dimension mismatch between Matlab and Eigen objects");
      return false;
    }

    typename Derived::Scalar* data = static_cast<typename Derived::Scalar*>(mxGetData(out));

    for (size_t i=0; i < rows; i++) {
      for (size_t j= 0; j < cols; j++) {
        data[i + j * rows] = in->coeff(i,j);
      }
    }

    return true;
  };

  template<class Derived>
  bool ConvertFromEigenToMatlabMatrix(mxArray** out, Eigen::MatrixBase<Derived>* in) {
    size_t rows = in->rows();
    size_t cols = in->cols();

    if (MxType<typename Derived::Scalar>() == mxUNKNOWN_CLASS) {
      mexPrintf("I don't know how to convert this scalar type to a Matlab matrix");
      return false;
    }

    *out = mxCreateNumericMatrix(rows, cols, MxType<typename Derived::Scalar>(), mxREAL);

    typename Derived::Scalar* data = static_cast<typename Derived::Scalar*>(mxGetData(*out));

    for (size_t i=0; i < rows; i++) {
      for (size_t j=0; j < cols; j++) {
        data[i + j * rows] = in->coeff(i,j);
      }
    }

    return true;
  }

  template<> mxClassID MxType<double>() {return mxDOUBLE_CLASS;};
  template<> mxClassID MxType<int>() {return mxINT32_CLASS;};


  SWIGINTERN int SWIG_AsVal_double (mxArray* pm, double* val)
    {
      if(!mxIsDouble(pm) || mxGetNumberOfElements(pm)!=1) return SWIG_TypeError;
      if (val) *val = mxGetScalar(pm);
      return SWIG_OK;
    }


  SWIGINTERNINLINE mxArray* SWIG_From_long    (long value)
    {    
      mxArray* out = mxCreateNumericMatrix(1,1,mxINT64_CLASS,mxREAL);
      *((int64_T*)mxGetData(out)) = (int64_T)value;
      return out;
    }


SWIGINTERNINLINE mxArray*
SWIG_From_int  (int value)
{    
  return SWIG_From_long  (value);
}

int _wrap_new_AutoDiffVectorDynamic__SWIG_0 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *result = 0 ;
  
  if (!SWIG_check_num_args("new_AutoDiffVectorDynamic",argc,0,0,0)) {
    SWIG_fail;
  }
  {
    try {
      result = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *)new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 1 |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_new_AutoDiffVectorDynamic__SWIG_1 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,1 > *arg1 = 0 ;
  void *argp1 ;
  int res1 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *result = 0 ;
  
  if (!SWIG_check_num_args("new_AutoDiffVectorDynamic",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1, SWIGTYPE_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_1_t,  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "new_AutoDiffVectorDynamic" "', argument " "1"" of type '" "Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,1 > const &""'"); 
  }
  if (!argp1) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "new_AutoDiffVectorDynamic" "', argument " "1"" of type '" "Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,1 > const &""'"); 
  }
  arg1 = reinterpret_cast< Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *)new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >((Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,1 > const &)*arg1);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 1 |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_new_AutoDiffVectorDynamic__SWIG_2 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  Eigen::Matrix< double,Eigen::Dynamic,1 > *arg1 = 0 ;
  Eigen::MatrixXd *arg2 = 0 ;
  Eigen::Matrix< double,Eigen::Dynamic,1 > temp1 ;
  Eigen::MatrixXd temp2 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *result = 0 ;
  
  if (!SWIG_check_num_args("new_AutoDiffVectorDynamic",argc,2,2,0)) {
    SWIG_fail;
  }
  {
    // In: const&
    if (!ConvertFromMatlabToEigenMatrix<Eigen::Matrix<double,Eigen::Dynamic,1> >(&temp1, argv[0]))
    SWIG_fail;
    arg1 = &temp1;
  }
  {
    // In: const&
    if (!ConvertFromMatlabToEigenMatrix<Eigen::MatrixXd >(&temp2, argv[1]))
    SWIG_fail;
    arg2 = &temp2;
  }
  {
    try {
      result = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *)new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >((Eigen::Matrix< double,Eigen::Dynamic,1 > const &)*arg1,(Eigen::MatrixXd const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 1 |  0 );
  if (_out) --resc, *resv++ = _out;
  
  
  return 0;
fail:
  return 1;
}


int _wrap_new_AutoDiffVectorDynamic (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc == 0) {
    return _wrap_new_AutoDiffVectorDynamic__SWIG_0(resc,resv,argc,argv);
  }
  if (argc == 1) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_1_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      return _wrap_new_AutoDiffVectorDynamic__SWIG_1(resc,resv,argc,argv);
    }
  }
  if (argc == 2) {
    int _v;
    {
      _v = mxIsNumeric(argv[0]);
    }
    if (_v) {
      {
        _v = mxIsNumeric(argv[1]);
      }
      if (_v) {
        return _wrap_new_AutoDiffVectorDynamic__SWIG_2(resc,resv,argc,argv);
      }
    }
  }
  
  SWIG_Error(SWIG_RuntimeError, "No matching function for overload function 'new_AutoDiffVectorDynamic'."
    "  Possible C/C++ prototypes are:\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >::AutoDiffWrapper()\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >::AutoDiffWrapper(Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,1 > const &)\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >::AutoDiffWrapper(Eigen::Matrix< double,Eigen::Dynamic,1 > const &,Eigen::MatrixXd const &)\n");
  return 1;
}


int _wrap_AutoDiffVectorDynamic_value (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  Eigen::Matrix< double,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_value",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_value" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (arg1)->value();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  {
    if (!ConvertFromEigenToMatlabMatrix<Eigen::Matrix<double,Eigen::Dynamic,1> >(&_out, &result))
    SWIG_fail;
  }
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_derivatives (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  Eigen::MatrixXd result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_derivatives",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_derivatives" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (arg1)->derivatives();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  {
    if (!ConvertFromEigenToMatlabMatrix<Eigen::MatrixXd >(&_out, &result))
    SWIG_fail;
  }
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_plus__SWIG_0 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_plus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_plus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffVectorDynamic_plus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffVectorDynamic_plus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp2);
  {
    try {
      result = (arg1)->operator +((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_minus__SWIG_0 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_minus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_minus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffVectorDynamic_minus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffVectorDynamic_minus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp2);
  {
    try {
      result = (arg1)->operator -((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_arrayMultiply (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_arrayMultiply",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_arrayMultiply" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffVectorDynamic_arrayMultiply" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffVectorDynamic_arrayMultiply" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp2);
  {
    try {
      result = (arg1)->arrayMultiply((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_arrayDivide (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_arrayDivide",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_arrayDivide" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffVectorDynamic_arrayDivide" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffVectorDynamic_arrayDivide" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp2);
  {
    try {
      result = (arg1)->arrayDivide((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_plus__SWIG_1 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_plus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_plus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffVectorDynamic_plus" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator +(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_plus (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      void *vptr = 0;
      int res = SWIG_ConvertPtr(argv[1], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0);
      _v = SWIG_CheckState(res);
      if (_v) {
        return _wrap_AutoDiffVectorDynamic_plus__SWIG_0(resc,resv,argc,argv);
      }
    }
  }
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      {
        int res = SWIG_AsVal_double(argv[1], NULL);
        _v = SWIG_CheckState(res);
      }
      if (_v) {
        return _wrap_AutoDiffVectorDynamic_plus__SWIG_1(resc,resv,argc,argv);
      }
    }
  }
  
  SWIG_Error(SWIG_RuntimeError, "No matching function for overload function 'AutoDiffVectorDynamic_plus'."
    "  Possible C/C++ prototypes are:\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >::operator +(AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &)\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >::operator +(double)\n");
  return 1;
}


int _wrap_AutoDiffVectorDynamic_minus__SWIG_1 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_minus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_minus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffVectorDynamic_minus" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator -(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_minus (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      void *vptr = 0;
      int res = SWIG_ConvertPtr(argv[1], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0);
      _v = SWIG_CheckState(res);
      if (_v) {
        return _wrap_AutoDiffVectorDynamic_minus__SWIG_0(resc,resv,argc,argv);
      }
    }
  }
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      {
        int res = SWIG_AsVal_double(argv[1], NULL);
        _v = SWIG_CheckState(res);
      }
      if (_v) {
        return _wrap_AutoDiffVectorDynamic_minus__SWIG_1(resc,resv,argc,argv);
      }
    }
  }
  
  SWIG_Error(SWIG_RuntimeError, "No matching function for overload function 'AutoDiffVectorDynamic_minus'."
    "  Possible C/C++ prototypes are:\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >::operator -(AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const &)\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >::operator -(double)\n");
  return 1;
}


int _wrap_AutoDiffVectorDynamic_mtimes (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_mtimes",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_mtimes" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffVectorDynamic_mtimes" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator *(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_mrdivide (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_mrdivide",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_mrdivide" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffVectorDynamic_mrdivide" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator /(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_rows (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  int result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_rows",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_rows" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (int)((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const *)arg1)->rows();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_From_int(static_cast< int >(result));
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_cols (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  int result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_cols",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_cols" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (int)((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const *)arg1)->cols();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_From_int(static_cast< int >(result));
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorDynamic_size (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  int result;
  
  if (!SWIG_check_num_args("AutoDiffVectorDynamic_size",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorDynamic_size" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (int)((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > const *)arg1)->size();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_From_int(static_cast< int >(result));
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_delete_AutoDiffVectorDynamic (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  
  int is_owned;
  if (!SWIG_check_num_args("delete_AutoDiffVectorDynamic",argc,1,1,0)) {
    SWIG_fail;
  }
  is_owned = SWIG_Matlab_isOwned(argv[0]);
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, SWIG_POINTER_DISOWN |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "delete_AutoDiffVectorDynamic" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > * >(argp1);
  if (is_owned) {
    {
      try {
        delete arg1;
      } catch (const std::exception& e) {
        SWIG_exception(SWIG_RuntimeError, e.what());
      } catch (...) {
        SWIG_exception(SWIG_RuntimeError, "Unknown error");
      }
    }
  }
  _out = (mxArray*)0;
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_new_AutoDiffVectorMax73__SWIG_0 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  mxArray * _out;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *result = 0 ;
  
  if (!SWIG_check_num_args("new_AutoDiffVectorMax73",argc,0,0,0)) {
    SWIG_fail;
  }
  {
    try {
      result = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *)new AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 1 |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_new_AutoDiffVectorMax73__SWIG_1 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 > >,Eigen::Dynamic,1 > *arg1 = 0 ;
  void *argp1 ;
  int res1 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *result = 0 ;
  
  if (!SWIG_check_num_args("new_AutoDiffVectorMax73",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1, SWIGTYPE_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_t_Eigen__Dynamic_1_t,  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "new_AutoDiffVectorMax73" "', argument " "1"" of type '" "Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 > >,Eigen::Dynamic,1 > const &""'"); 
  }
  if (!argp1) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "new_AutoDiffVectorMax73" "', argument " "1"" of type '" "Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 > >,Eigen::Dynamic,1 > const &""'"); 
  }
  arg1 = reinterpret_cast< Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 > >,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *)new AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >((Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 > >,Eigen::Dynamic,1 > const &)*arg1);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 1 |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_new_AutoDiffVectorMax73__SWIG_2 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  Eigen::Matrix< double,Eigen::Dynamic,1 > *arg1 = 0 ;
  Eigen::MatrixXd *arg2 = 0 ;
  Eigen::Matrix< double,Eigen::Dynamic,1 > temp1 ;
  Eigen::MatrixXd temp2 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *result = 0 ;
  
  if (!SWIG_check_num_args("new_AutoDiffVectorMax73",argc,2,2,0)) {
    SWIG_fail;
  }
  {
    // In: const&
    if (!ConvertFromMatlabToEigenMatrix<Eigen::Matrix<double,Eigen::Dynamic,1> >(&temp1, argv[0]))
    SWIG_fail;
    arg1 = &temp1;
  }
  {
    // In: const&
    if (!ConvertFromMatlabToEigenMatrix<Eigen::MatrixXd >(&temp2, argv[1]))
    SWIG_fail;
    arg2 = &temp2;
  }
  {
    try {
      result = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *)new AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >((Eigen::Matrix< double,Eigen::Dynamic,1 > const &)*arg1,(Eigen::MatrixXd const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 1 |  0 );
  if (_out) --resc, *resv++ = _out;
  
  
  return 0;
fail:
  return 1;
}


int _wrap_new_AutoDiffVectorMax73 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc == 0) {
    return _wrap_new_AutoDiffVectorMax73__SWIG_0(resc,resv,argc,argv);
  }
  if (argc == 1) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_t_Eigen__Dynamic_1_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      return _wrap_new_AutoDiffVectorMax73__SWIG_1(resc,resv,argc,argv);
    }
  }
  if (argc == 2) {
    int _v;
    {
      _v = mxIsNumeric(argv[0]);
    }
    if (_v) {
      {
        _v = mxIsNumeric(argv[1]);
      }
      if (_v) {
        return _wrap_new_AutoDiffVectorMax73__SWIG_2(resc,resv,argc,argv);
      }
    }
  }
  
  SWIG_Error(SWIG_RuntimeError, "No matching function for overload function 'new_AutoDiffVectorMax73'."
    "  Possible C/C++ prototypes are:\n"
    "    AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >::AutoDiffWrapper()\n"
    "    AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >::AutoDiffWrapper(Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 > >,Eigen::Dynamic,1 > const &)\n"
    "    AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >::AutoDiffWrapper(Eigen::Matrix< double,Eigen::Dynamic,1 > const &,Eigen::MatrixXd const &)\n");
  return 1;
}


int _wrap_AutoDiffVectorMax73_value (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  Eigen::Matrix< double,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_value",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_value" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (arg1)->value();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  {
    if (!ConvertFromEigenToMatlabMatrix<Eigen::Matrix<double,Eigen::Dynamic,1> >(&_out, &result))
    SWIG_fail;
  }
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_derivatives (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  Eigen::MatrixXd result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_derivatives",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_derivatives" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (arg1)->derivatives();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  {
    if (!ConvertFromEigenToMatlabMatrix<Eigen::MatrixXd >(&_out, &result))
    SWIG_fail;
  }
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_plus__SWIG_0 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_plus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_plus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffVectorMax73_plus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffVectorMax73_plus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp2);
  {
    try {
      result = (arg1)->operator +((AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_minus__SWIG_0 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_minus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_minus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffVectorMax73_minus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffVectorMax73_minus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp2);
  {
    try {
      result = (arg1)->operator -((AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_arrayMultiply (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_arrayMultiply",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_arrayMultiply" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffVectorMax73_arrayMultiply" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffVectorMax73_arrayMultiply" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp2);
  {
    try {
      result = (arg1)->arrayMultiply((AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_arrayDivide (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_arrayDivide",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_arrayDivide" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffVectorMax73_arrayDivide" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffVectorMax73_arrayDivide" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp2);
  {
    try {
      result = (arg1)->arrayDivide((AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_plus__SWIG_1 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_plus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_plus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffVectorMax73_plus" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator +(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_plus (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      void *vptr = 0;
      int res = SWIG_ConvertPtr(argv[1], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0);
      _v = SWIG_CheckState(res);
      if (_v) {
        return _wrap_AutoDiffVectorMax73_plus__SWIG_0(resc,resv,argc,argv);
      }
    }
  }
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      {
        int res = SWIG_AsVal_double(argv[1], NULL);
        _v = SWIG_CheckState(res);
      }
      if (_v) {
        return _wrap_AutoDiffVectorMax73_plus__SWIG_1(resc,resv,argc,argv);
      }
    }
  }
  
  SWIG_Error(SWIG_RuntimeError, "No matching function for overload function 'AutoDiffVectorMax73_plus'."
    "  Possible C/C++ prototypes are:\n"
    "    AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >::operator +(AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &)\n"
    "    AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >::operator +(double)\n");
  return 1;
}


int _wrap_AutoDiffVectorMax73_minus__SWIG_1 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_minus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_minus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffVectorMax73_minus" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator -(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_minus (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      void *vptr = 0;
      int res = SWIG_ConvertPtr(argv[1], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0);
      _v = SWIG_CheckState(res);
      if (_v) {
        return _wrap_AutoDiffVectorMax73_minus__SWIG_0(resc,resv,argc,argv);
      }
    }
  }
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      {
        int res = SWIG_AsVal_double(argv[1], NULL);
        _v = SWIG_CheckState(res);
      }
      if (_v) {
        return _wrap_AutoDiffVectorMax73_minus__SWIG_1(resc,resv,argc,argv);
      }
    }
  }
  
  SWIG_Error(SWIG_RuntimeError, "No matching function for overload function 'AutoDiffVectorMax73_minus'."
    "  Possible C/C++ prototypes are:\n"
    "    AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >::operator -(AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const &)\n"
    "    AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >::operator -(double)\n");
  return 1;
}


int _wrap_AutoDiffVectorMax73_mtimes (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_mtimes",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_mtimes" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffVectorMax73_mtimes" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator *(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_mrdivide (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_mrdivide",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_mrdivide" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffVectorMax73_mrdivide" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator /(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >(static_cast< const AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_rows (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  int result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_rows",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_rows" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (int)((AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const *)arg1)->rows();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_From_int(static_cast< int >(result));
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_cols (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  int result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_cols",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_cols" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (int)((AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const *)arg1)->cols();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_From_int(static_cast< int >(result));
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffVectorMax73_size (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  int result;
  
  if (!SWIG_check_num_args("AutoDiffVectorMax73_size",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffVectorMax73_size" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  {
    try {
      result = (int)((AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > const *)arg1)->size();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_From_int(static_cast< int >(result));
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_delete_AutoDiffVectorMax73 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *arg1 = (AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  
  int is_owned;
  if (!SWIG_check_num_args("delete_AutoDiffVectorMax73",argc,1,1,0)) {
    SWIG_fail;
  }
  is_owned = SWIG_Matlab_isOwned(argv[0]);
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, SWIG_POINTER_DISOWN |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "delete_AutoDiffVectorMax73" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > * >(argp1);
  if (is_owned) {
    {
      try {
        delete arg1;
      } catch (const std::exception& e) {
        SWIG_exception(SWIG_RuntimeError, e.what());
      } catch (...) {
        SWIG_exception(SWIG_RuntimeError, "Unknown error");
      }
    }
  }
  _out = (mxArray*)0;
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_new_AutoDiffMatrixDynamic__SWIG_0 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *result = 0 ;
  
  if (!SWIG_check_num_args("new_AutoDiffMatrixDynamic",argc,0,0,0)) {
    SWIG_fail;
  }
  {
    try {
      result = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *)new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 1 |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_new_AutoDiffMatrixDynamic__SWIG_1 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,Eigen::Dynamic > *arg1 = 0 ;
  void *argp1 ;
  int res1 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *result = 0 ;
  
  if (!SWIG_check_num_args("new_AutoDiffMatrixDynamic",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1, SWIGTYPE_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_Eigen__Dynamic_t,  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "new_AutoDiffMatrixDynamic" "', argument " "1"" of type '" "Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,Eigen::Dynamic > const &""'"); 
  }
  if (!argp1) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "new_AutoDiffMatrixDynamic" "', argument " "1"" of type '" "Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,Eigen::Dynamic > const &""'"); 
  }
  arg1 = reinterpret_cast< Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  {
    try {
      result = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *)new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >((Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,Eigen::Dynamic > const &)*arg1);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 1 |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_new_AutoDiffMatrixDynamic__SWIG_2 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  Eigen::Matrix< double,Eigen::Dynamic,Eigen::Dynamic > *arg1 = 0 ;
  Eigen::MatrixXd *arg2 = 0 ;
  Eigen::Matrix< double,Eigen::Dynamic,Eigen::Dynamic > temp1 ;
  Eigen::MatrixXd temp2 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *result = 0 ;
  
  if (!SWIG_check_num_args("new_AutoDiffMatrixDynamic",argc,2,2,0)) {
    SWIG_fail;
  }
  {
    // In: const&
    if (!ConvertFromMatlabToEigenMatrix<Eigen::Matrix<double,Eigen::Dynamic,Eigen::Dynamic> >(&temp1, argv[0]))
    SWIG_fail;
    arg1 = &temp1;
  }
  {
    // In: const&
    if (!ConvertFromMatlabToEigenMatrix<Eigen::MatrixXd >(&temp2, argv[1]))
    SWIG_fail;
    arg2 = &temp2;
  }
  {
    try {
      result = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *)new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >((Eigen::Matrix< double,Eigen::Dynamic,Eigen::Dynamic > const &)*arg1,(Eigen::MatrixXd const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 1 |  0 );
  if (_out) --resc, *resv++ = _out;
  
  
  return 0;
fail:
  return 1;
}


int _wrap_new_AutoDiffMatrixDynamic (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc == 0) {
    return _wrap_new_AutoDiffMatrixDynamic__SWIG_0(resc,resv,argc,argv);
  }
  if (argc == 1) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_Eigen__Dynamic_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      return _wrap_new_AutoDiffMatrixDynamic__SWIG_1(resc,resv,argc,argv);
    }
  }
  if (argc == 2) {
    int _v;
    {
      _v = mxIsNumeric(argv[0]);
    }
    if (_v) {
      {
        _v = mxIsNumeric(argv[1]);
      }
      if (_v) {
        return _wrap_new_AutoDiffMatrixDynamic__SWIG_2(resc,resv,argc,argv);
      }
    }
  }
  
  SWIG_Error(SWIG_RuntimeError, "No matching function for overload function 'new_AutoDiffMatrixDynamic'."
    "  Possible C/C++ prototypes are:\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >::AutoDiffWrapper()\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >::AutoDiffWrapper(Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,Eigen::Dynamic > const &)\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >::AutoDiffWrapper(Eigen::Matrix< double,Eigen::Dynamic,Eigen::Dynamic > const &,Eigen::MatrixXd const &)\n");
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_value (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  Eigen::Matrix< double,Eigen::Dynamic,Eigen::Dynamic > result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_value",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_value" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  {
    try {
      result = (arg1)->value();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  {
    if (!ConvertFromEigenToMatlabMatrix<Eigen::Matrix<double,Eigen::Dynamic,Eigen::Dynamic> >(&_out, &result))
    SWIG_fail;
  }
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_derivatives (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  Eigen::MatrixXd result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_derivatives",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_derivatives" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  {
    try {
      result = (arg1)->derivatives();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  {
    if (!ConvertFromEigenToMatlabMatrix<Eigen::MatrixXd >(&_out, &result))
    SWIG_fail;
  }
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_plus__SWIG_0 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_plus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_plus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffMatrixDynamic_plus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffMatrixDynamic_plus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp2);
  {
    try {
      result = (arg1)->operator +((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_minus__SWIG_0 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_minus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_minus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffMatrixDynamic_minus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffMatrixDynamic_minus" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp2);
  {
    try {
      result = (arg1)->operator -((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_arrayMultiply (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_arrayMultiply",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_arrayMultiply" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffMatrixDynamic_arrayMultiply" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffMatrixDynamic_arrayMultiply" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp2);
  {
    try {
      result = (arg1)->arrayMultiply((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_arrayDivide (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg2 = 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  void *argp2 ;
  int res2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_arrayDivide",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_arrayDivide" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  res2 = SWIG_ConvertPtr(argv[1], &argp2, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t,  0 );
  if (!SWIG_IsOK(res2)) {
    SWIG_exception_fail(SWIG_ArgError(res2), "in method '" "AutoDiffMatrixDynamic_arrayDivide" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &""'"); 
  }
  if (!argp2) {
    SWIG_exception_fail(SWIG_ValueError, "invalid null reference " "in method '" "AutoDiffMatrixDynamic_arrayDivide" "', argument " "2"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &""'"); 
  }
  arg2 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp2);
  {
    try {
      result = (arg1)->arrayDivide((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &)*arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_plus__SWIG_1 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_plus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_plus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffMatrixDynamic_plus" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator +(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_plus (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      void *vptr = 0;
      int res = SWIG_ConvertPtr(argv[1], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0);
      _v = SWIG_CheckState(res);
      if (_v) {
        return _wrap_AutoDiffMatrixDynamic_plus__SWIG_0(resc,resv,argc,argv);
      }
    }
  }
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      {
        int res = SWIG_AsVal_double(argv[1], NULL);
        _v = SWIG_CheckState(res);
      }
      if (_v) {
        return _wrap_AutoDiffMatrixDynamic_plus__SWIG_1(resc,resv,argc,argv);
      }
    }
  }
  
  SWIG_Error(SWIG_RuntimeError, "No matching function for overload function 'AutoDiffMatrixDynamic_plus'."
    "  Possible C/C++ prototypes are:\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >::operator +(AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &)\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >::operator +(double)\n");
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_minus__SWIG_1 (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_minus",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_minus" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffMatrixDynamic_minus" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator -(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_minus (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      void *vptr = 0;
      int res = SWIG_ConvertPtr(argv[1], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0);
      _v = SWIG_CheckState(res);
      if (_v) {
        return _wrap_AutoDiffMatrixDynamic_minus__SWIG_0(resc,resv,argc,argv);
      }
    }
  }
  if (argc == 2) {
    int _v;
    void *vptr = 0;
    int res = SWIG_ConvertPtr(argv[0], &vptr, SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0);
    _v = SWIG_CheckState(res);
    if (_v) {
      {
        int res = SWIG_AsVal_double(argv[1], NULL);
        _v = SWIG_CheckState(res);
      }
      if (_v) {
        return _wrap_AutoDiffMatrixDynamic_minus__SWIG_1(resc,resv,argc,argv);
      }
    }
  }
  
  SWIG_Error(SWIG_RuntimeError, "No matching function for overload function 'AutoDiffMatrixDynamic_minus'."
    "  Possible C/C++ prototypes are:\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >::operator -(AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const &)\n"
    "    AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >::operator -(double)\n");
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_mtimes (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_mtimes",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_mtimes" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffMatrixDynamic_mtimes" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator *(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_mrdivide (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  double arg2 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  double val2 ;
  int ecode2 = 0 ;
  mxArray * _out;
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_mrdivide",argc,2,2,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_mrdivide" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  ecode2 = SWIG_AsVal_double(argv[1], &val2);
  if (!SWIG_IsOK(ecode2)) {
    SWIG_exception_fail(SWIG_ArgError(ecode2), "in method '" "AutoDiffMatrixDynamic_mrdivide" "', argument " "2"" of type '" "double""'");
  } 
  arg2 = static_cast< double >(val2);
  {
    try {
      result = (arg1)->operator /(arg2);
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_NewPointerObj((new AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >(static_cast< const AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic >& >(result))), SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, SWIG_POINTER_OWN |  0 );
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_rows (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  int result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_rows",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_rows" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  {
    try {
      result = (int)((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const *)arg1)->rows();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_From_int(static_cast< int >(result));
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_cols (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  int result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_cols",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_cols" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  {
    try {
      result = (int)((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const *)arg1)->cols();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_From_int(static_cast< int >(result));
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_AutoDiffMatrixDynamic_size (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  int result;
  
  if (!SWIG_check_num_args("AutoDiffMatrixDynamic_size",argc,1,1,0)) {
    SWIG_fail;
  }
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0 |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "AutoDiffMatrixDynamic_size" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  {
    try {
      result = (int)((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > const *)arg1)->size();
    } catch (const std::exception& e) {
      SWIG_exception(SWIG_RuntimeError, e.what());
    } catch (...) {
      SWIG_exception(SWIG_RuntimeError, "Unknown error");
    }
  }
  _out = SWIG_From_int(static_cast< int >(result));
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}


int _wrap_delete_AutoDiffMatrixDynamic (int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *arg1 = (AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) 0 ;
  void *argp1 = 0 ;
  int res1 = 0 ;
  mxArray * _out;
  
  int is_owned;
  if (!SWIG_check_num_args("delete_AutoDiffMatrixDynamic",argc,1,1,0)) {
    SWIG_fail;
  }
  is_owned = SWIG_Matlab_isOwned(argv[0]);
  res1 = SWIG_ConvertPtr(argv[0], &argp1,SWIGTYPE_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, SWIG_POINTER_DISOWN |  0 );
  if (!SWIG_IsOK(res1)) {
    SWIG_exception_fail(SWIG_ArgError(res1), "in method '" "delete_AutoDiffMatrixDynamic" "', argument " "1"" of type '" "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *""'"); 
  }
  arg1 = reinterpret_cast< AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > * >(argp1);
  if (is_owned) {
    {
      try {
        delete arg1;
      } catch (const std::exception& e) {
        SWIG_exception(SWIG_RuntimeError, e.what());
      } catch (...) {
        SWIG_exception(SWIG_RuntimeError, "Unknown error");
      }
    }
  }
  _out = (mxArray*)0;
  if (_out) --resc, *resv++ = _out;
  return 0;
fail:
  return 1;
}



/* -------- TYPE CONVERSION AND EQUIVALENCE RULES (BEGIN) -------- */

static void *_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_tTo_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_t_Eigen__Dynamic_1_t(void *x, int *SWIGUNUSEDPARM(newmemory)) {
    return (void *)((Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 > >,Eigen::Dynamic,1 > *)  ((AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *) x));
}
static void *_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_tTo_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_1_t(void *x, int *SWIGUNUSEDPARM(newmemory)) {
    return (void *)((Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,1 > *)  ((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *) x));
}
static void *_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_tTo_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_Eigen__Dynamic_t(void *x, int *SWIGUNUSEDPARM(newmemory)) {
    return (void *)((Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,Eigen::Dynamic > *)  ((AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *) x));
}
static swig_type_info _swigt__p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t = {"_p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t", "AutoDiffWrapper< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 >,Eigen::Dynamic,1 > *", 0, 0, (void*)"autodiffutils.AutoDiffVectorMax73", 0};
static swig_type_info _swigt__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t = {"_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t", "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,1 > *", 0, 0, (void*)"autodiffutils.AutoDiffVectorDynamic", 0};
static swig_type_info _swigt__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t = {"_p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t", "AutoDiffWrapper< Eigen::VectorXd,Eigen::Dynamic,Eigen::Dynamic > *", 0, 0, (void*)"autodiffutils.AutoDiffMatrixDynamic", 0};
static swig_type_info _swigt__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_t_Eigen__Dynamic_1_t = {"_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_t_Eigen__Dynamic_1_t", "Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::Matrix< double,Eigen::Dynamic,1,0,73 > >,Eigen::Dynamic,1 > *", 0, 0, (void*)0, 0};
static swig_type_info _swigt__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_1_t = {"_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_1_t", "Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,1 > *", 0, 0, (void*)0, 0};
static swig_type_info _swigt__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_Eigen__Dynamic_t = {"_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_Eigen__Dynamic_t", "Eigen::Matrix< Eigen::AutoDiffScalar< Eigen::VectorXd >,Eigen::Dynamic,Eigen::Dynamic > *", 0, 0, (void*)0, 0};
static swig_type_info _swigt__p_char = {"_p_char", "char *", 0, 0, (void*)0, 0};

static swig_type_info *swig_type_initial[] = {
  &_swigt__p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t,
  &_swigt__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t,
  &_swigt__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t,
  &_swigt__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_t_Eigen__Dynamic_1_t,
  &_swigt__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_1_t,
  &_swigt__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_Eigen__Dynamic_t,
  &_swigt__p_char,
};

static swig_cast_info _swigc__p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t[] = {  {&_swigt__p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, 0, 0, 0},{0, 0, 0, 0}};
static swig_cast_info _swigc__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t[] = {  {&_swigt__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, 0, 0, 0},{0, 0, 0, 0}};
static swig_cast_info _swigc__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t[] = {  {&_swigt__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, 0, 0, 0},{0, 0, 0, 0}};
static swig_cast_info _swigc__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_t_Eigen__Dynamic_1_t[] = {  {&_swigt__p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t, _p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_tTo_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_t_Eigen__Dynamic_1_t, 0, 0},  {&_swigt__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_t_Eigen__Dynamic_1_t, 0, 0, 0},{0, 0, 0, 0}};
static swig_cast_info _swigc__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_1_t[] = {  {&_swigt__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t, _p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_tTo_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_1_t, 0, 0},  {&_swigt__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_1_t, 0, 0, 0},{0, 0, 0, 0}};
static swig_cast_info _swigc__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_Eigen__Dynamic_t[] = {  {&_swigt__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t, _p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_tTo_p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_Eigen__Dynamic_t, 0, 0},  {&_swigt__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_Eigen__Dynamic_t, 0, 0, 0},{0, 0, 0, 0}};
static swig_cast_info _swigc__p_char[] = {  {&_swigt__p_char, 0, 0, 0},{0, 0, 0, 0}};

static swig_cast_info *swig_cast_initial[] = {
  _swigc__p_AutoDiffWrapperT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_Eigen__Dynamic_1_t,
  _swigc__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_1_t,
  _swigc__p_AutoDiffWrapperT_Eigen__VectorXd_Eigen__Dynamic_Eigen__Dynamic_t,
  _swigc__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__MatrixT_double_Eigen__Dynamic_1_0_73_t_t_Eigen__Dynamic_1_t,
  _swigc__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_1_t,
  _swigc__p_Eigen__MatrixT_Eigen__AutoDiffScalarT_Eigen__VectorXd_t_Eigen__Dynamic_Eigen__Dynamic_t,
  _swigc__p_char,
};


/* -------- TYPE CONVERSION AND EQUIVALENCE RULES (END) -------- */

/* -----------------------------------------------------------------------------
 * Type initialization:
 * This problem is tough by the requirement that no dynamic
 * memory is used. Also, since swig_type_info structures store pointers to
 * swig_cast_info structures and swig_cast_info structures store pointers back
 * to swig_type_info structures, we need some lookup code at initialization.
 * The idea is that swig generates all the structures that are needed.
 * The runtime then collects these partially filled structures.
 * The SWIG_InitializeModule function takes these initial arrays out of
 * swig_module, and does all the lookup, filling in the swig_module.types
 * array with the correct data and linking the correct swig_cast_info
 * structures together.
 *
 * The generated swig_type_info structures are assigned statically to an initial
 * array. We just loop through that array, and handle each type individually.
 * First we lookup if this type has been already loaded, and if so, use the
 * loaded structure instead of the generated one. Then we have to fill in the
 * cast linked list. The cast data is initially stored in something like a
 * two-dimensional array. Each row corresponds to a type (there are the same
 * number of rows as there are in the swig_type_initial array). Each entry in
 * a column is one of the swig_cast_info structures for that type.
 * The cast_initial array is actually an array of arrays, because each row has
 * a variable number of columns. So to actually build the cast linked list,
 * we find the array of casts associated with the type, and loop through it
 * adding the casts to the list. The one last trick we need to do is making
 * sure the type pointer in the swig_cast_info struct is correct.
 *
 * First off, we lookup the cast->type name to see if it is already loaded.
 * There are three cases to handle:
 *  1) If the cast->type has already been loaded AND the type we are adding
 *     casting info to has not been loaded (it is in this module), THEN we
 *     replace the cast->type pointer with the type pointer that has already
 *     been loaded.
 *  2) If BOTH types (the one we are adding casting info to, and the
 *     cast->type) are loaded, THEN the cast info has already been loaded by
 *     the previous module so we just ignore it.
 *  3) Finally, if cast->type has not already been loaded, then we add that
 *     swig_cast_info to the linked list (because the cast->type) pointer will
 *     be correct.
 * ----------------------------------------------------------------------------- */

#ifdef __cplusplus
extern "C" {
#if 0
} /* c-mode */
#endif
#endif

#if 0
#define SWIGRUNTIME_DEBUG
#endif


SWIGRUNTIME void
SWIG_InitializeModule(void *clientdata) {
  size_t i;
  swig_module_info *module_head, *iter;
  int init;

  /* check to see if the circular list has been setup, if not, set it up */
  if (swig_module.next==0) {
    /* Initialize the swig_module */
    swig_module.type_initial = swig_type_initial;
    swig_module.cast_initial = swig_cast_initial;
    swig_module.next = &swig_module;
    init = 1;
  } else {
    init = 0;
  }

  /* Try and load any already created modules */
  module_head = SWIG_GetModule(clientdata);
  if (!module_head) {
    /* This is the first module loaded for this interpreter */
    /* so set the swig module into the interpreter */
    SWIG_SetModule(clientdata, &swig_module);
  } else {
    /* the interpreter has loaded a SWIG module, but has it loaded this one? */
    iter=module_head;
    do {
      if (iter==&swig_module) {
        /* Our module is already in the list, so there's nothing more to do. */
        return;
      }
      iter=iter->next;
    } while (iter!= module_head);

    /* otherwise we must add our module into the list */
    swig_module.next = module_head->next;
    module_head->next = &swig_module;
  }

  /* When multiple interpreters are used, a module could have already been initialized in
     a different interpreter, but not yet have a pointer in this interpreter.
     In this case, we do not want to continue adding types... everything should be
     set up already */
  if (init == 0) return;

  /* Now work on filling in swig_module.types */
#ifdef SWIGRUNTIME_DEBUG
  printf("SWIG_InitializeModule: size %d\n", swig_module.size);
#endif
  for (i = 0; i < swig_module.size; ++i) {
    swig_type_info *type = 0;
    swig_type_info *ret;
    swig_cast_info *cast;

#ifdef SWIGRUNTIME_DEBUG
    printf("SWIG_InitializeModule: type %d %s\n", i, swig_module.type_initial[i]->name);
#endif

    /* if there is another module already loaded */
    if (swig_module.next != &swig_module) {
      type = SWIG_MangledTypeQueryModule(swig_module.next, &swig_module, swig_module.type_initial[i]->name);
    }
    if (type) {
      /* Overwrite clientdata field */
#ifdef SWIGRUNTIME_DEBUG
      printf("SWIG_InitializeModule: found type %s\n", type->name);
#endif
      if (swig_module.type_initial[i]->clientdata) {
	type->clientdata = swig_module.type_initial[i]->clientdata;
#ifdef SWIGRUNTIME_DEBUG
      printf("SWIG_InitializeModule: found and overwrite type %s \n", type->name);
#endif
      }
    } else {
      type = swig_module.type_initial[i];
    }

    /* Insert casting types */
    cast = swig_module.cast_initial[i];
    while (cast->type) {

      /* Don't need to add information already in the list */
      ret = 0;
#ifdef SWIGRUNTIME_DEBUG
      printf("SWIG_InitializeModule: look cast %s\n", cast->type->name);
#endif
      if (swig_module.next != &swig_module) {
        ret = SWIG_MangledTypeQueryModule(swig_module.next, &swig_module, cast->type->name);
#ifdef SWIGRUNTIME_DEBUG
	if (ret) printf("SWIG_InitializeModule: found cast %s\n", ret->name);
#endif
      }
      if (ret) {
	if (type == swig_module.type_initial[i]) {
#ifdef SWIGRUNTIME_DEBUG
	  printf("SWIG_InitializeModule: skip old type %s\n", ret->name);
#endif
	  cast->type = ret;
	  ret = 0;
	} else {
	  /* Check for casting already in the list */
	  swig_cast_info *ocast = SWIG_TypeCheck(ret->name, type);
#ifdef SWIGRUNTIME_DEBUG
	  if (ocast) printf("SWIG_InitializeModule: skip old cast %s\n", ret->name);
#endif
	  if (!ocast) ret = 0;
	}
      }

      if (!ret) {
#ifdef SWIGRUNTIME_DEBUG
	printf("SWIG_InitializeModule: adding cast %s\n", cast->type->name);
#endif
        if (type->cast) {
          type->cast->prev = cast;
          cast->next = type->cast;
        }
        type->cast = cast;
      }
      cast++;
    }
    /* Set entry in modules->types array equal to the type */
    swig_module.types[i] = type;
  }
  swig_module.types[i] = 0;

#ifdef SWIGRUNTIME_DEBUG
  printf("**** SWIG_InitializeModule: Cast List ******\n");
  for (i = 0; i < swig_module.size; ++i) {
    int j = 0;
    swig_cast_info *cast = swig_module.cast_initial[i];
    printf("SWIG_InitializeModule: type %d %s\n", i, swig_module.type_initial[i]->name);
    while (cast->type) {
      printf("SWIG_InitializeModule: cast type %s\n", cast->type->name);
      cast++;
      ++j;
    }
  printf("---- Total casts: %d\n",j);
  }
  printf("**** SWIG_InitializeModule: Cast List ******\n");
#endif
}

/* This function will propagate the clientdata field of type to
* any new swig_type_info structures that have been added into the list
* of equivalent types.  It is like calling
* SWIG_TypeClientData(type, clientdata) a second time.
*/
SWIGRUNTIME void
SWIG_PropagateClientData(void) {
  size_t i;
  swig_cast_info *equiv;
  static int init_run = 0;

  if (init_run) return;
  init_run = 1;

  for (i = 0; i < swig_module.size; i++) {
    if (swig_module.types[i]->clientdata) {
      equiv = swig_module.types[i]->cast;
      while (equiv) {
        if (!equiv->converter) {
          if (equiv->type && !equiv->type->clientdata)
            SWIG_TypeClientData(equiv->type, swig_module.types[i]->clientdata);
        }
        equiv = equiv->next;
      }
    }
  }
}

#ifdef __cplusplus
#if 0
{ /* c-mode */
#endif
}
#endif


void SWIG_Matlab_LoadModule() {
  /* Only call once */
  static int moduleIsLoaded = 0;
  if (moduleIsLoaded) return;
  moduleIsLoaded = 1;

  /* Prevent unloading this file until MATLAB exits */
  mexLock();

  /* Exit function (called when unloaded) */
  mexAtExit(SWIG_Matlab_ExitFcn);

  /* Load dependent modules and initialize */

SWIG_InitializeModule(0);

}


const char* swigConstantName_(int con_id) {
  switch (con_id) {
  default: return 0;
  }
}

int swigConstantName(int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc!=1 || !mxIsDouble(*argv) || mxGetNumberOfElements(*argv)!=1) {
    SWIG_Error(SWIG_RuntimeError, "This mex file should only be called from inside the .m files generated by SWIG. The input should be the constant ID.");
    return 1;
  }
  if (resc!=1) {
    SWIG_Error(SWIG_RuntimeError, "The function should have one output.");
    return 1;
  }
  int con_id = (int)mxGetScalar(*argv);
  const char* s = swigConstantName_(con_id);
  if (s==0) {
    SWIG_Error(SWIG_RuntimeError, "No such constant ID.");
    return 1;
  }
  *resv = mxCreateString(s);
  return 0;
}

extern "C"
int swigConstant(int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (--argc < 0 || !mxIsDouble(*argv) || mxGetNumberOfElements(*argv)!=1) {
    SWIG_Error(SWIG_RuntimeError, "This function should only be called from inside the .m files generated by SWIG. First input should be the constant ID .");
    return 1;
  }
  int con_id = (int)mxGetScalar(*argv++);
  switch (con_id) {
  default:
    SWIG_Error(SWIG_RuntimeError, "No such constant.");
    return 1;
  }
  return 0;
}

const char* swigFunctionName_(int fcn_id) {
  switch(fcn_id) {
  case 0: return "swigConstant";
  case 1: return "swigFunctionName";
  case 2: return "swigConstantName";
  case 3: return "swigThis";
  case 4: return "swigTouch";
  case 5: return "new_AutoDiffVectorDynamic";
  case 6: return "AutoDiffVectorDynamic_value";
  case 7: return "AutoDiffVectorDynamic_derivatives";
  case 8: return "AutoDiffVectorDynamic_arrayMultiply";
  case 9: return "AutoDiffVectorDynamic_arrayDivide";
  case 10: return "AutoDiffVectorDynamic_plus";
  case 11: return "AutoDiffVectorDynamic_minus";
  case 12: return "AutoDiffVectorDynamic_mtimes";
  case 13: return "AutoDiffVectorDynamic_mrdivide";
  case 14: return "AutoDiffVectorDynamic_rows";
  case 15: return "AutoDiffVectorDynamic_cols";
  case 16: return "AutoDiffVectorDynamic_size";
  case 17: return "delete_AutoDiffVectorDynamic";
  case 18: return "new_AutoDiffVectorMax73";
  case 19: return "AutoDiffVectorMax73_value";
  case 20: return "AutoDiffVectorMax73_derivatives";
  case 21: return "AutoDiffVectorMax73_arrayMultiply";
  case 22: return "AutoDiffVectorMax73_arrayDivide";
  case 23: return "AutoDiffVectorMax73_plus";
  case 24: return "AutoDiffVectorMax73_minus";
  case 25: return "AutoDiffVectorMax73_mtimes";
  case 26: return "AutoDiffVectorMax73_mrdivide";
  case 27: return "AutoDiffVectorMax73_rows";
  case 28: return "AutoDiffVectorMax73_cols";
  case 29: return "AutoDiffVectorMax73_size";
  case 30: return "delete_AutoDiffVectorMax73";
  case 31: return "new_AutoDiffMatrixDynamic";
  case 32: return "AutoDiffMatrixDynamic_value";
  case 33: return "AutoDiffMatrixDynamic_derivatives";
  case 34: return "AutoDiffMatrixDynamic_arrayMultiply";
  case 35: return "AutoDiffMatrixDynamic_arrayDivide";
  case 36: return "AutoDiffMatrixDynamic_plus";
  case 37: return "AutoDiffMatrixDynamic_minus";
  case 38: return "AutoDiffMatrixDynamic_mtimes";
  case 39: return "AutoDiffMatrixDynamic_mrdivide";
  case 40: return "AutoDiffMatrixDynamic_rows";
  case 41: return "AutoDiffMatrixDynamic_cols";
  case 42: return "AutoDiffMatrixDynamic_size";
  case 43: return "delete_AutoDiffMatrixDynamic";
  default: return 0;
  }
}

int swigFunctionName(int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc!=1 || !mxIsDouble(*argv) || mxGetNumberOfElements(*argv)!=1) {
    SWIG_Error(SWIG_RuntimeError, "This mex file should only be called from inside the .m files generated by SWIG. The input should be the function ID.");
    return 1;
  }
  if (resc!=1) {
    SWIG_Error(SWIG_RuntimeError, "The function should have one output.");
    return 1;
  }
  int fcn_id = (int)mxGetScalar(*argv);
  const char* s = swigFunctionName_(fcn_id);
  if (s==0) {
    SWIG_Error(SWIG_RuntimeError, "No such function ID.");
    return 1;
  }
  *resv = mxCreateString(s);
  return 0;
}

int swigThis(int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc!=1 || resc!=1) {
    SWIG_Error(SWIG_RuntimeError, "swigThis should have one input and one output.");
    return 1;
  }
  SwigPtr* swig_ptr = SWIG_Matlab_getSwigPtr(argv[0]);
  if (!swig_ptr) {
    SWIG_Error(SWIG_RuntimeError, "The argument should be a SWIG class.");
    return 1;
  }
  resv[0] = mxCreateNumericMatrix(1, 1, mxUINT64_CLASS, mxREAL);
  *(uint64_T *)mxGetData(resv[0]) = (uint64_T)swig_ptr->ptr;
  return 0;
}

int swigTouch(int resc, mxArray *resv[], int argc, mxArray *argv[]) {
  if (argc!=0 || resc!=0) {
    SWIG_Error(SWIG_RuntimeError, "swigTouch should have no inputs or outputs.");
    return 1;
  }
  return 0;
}

extern "C"
void mexFunction(int resc, mxArray *resv[], int argc, const mxArray *argv[]) {
  /* Initialize module if first call */
  SWIG_Matlab_LoadModule();

  if (--argc < 0 || !mxIsDouble(*argv) || mxGetNumberOfElements(*argv)!=1)
    mexErrMsgTxt("This mex file should only be called from inside the .m files generated by SWIG. First input should be the function ID .");
  int fcn_id = (int)mxGetScalar(*argv++);
  int flag=0;
  switch (fcn_id) {
  case 0: flag=swigConstant(resc,resv,argc,(mxArray**)(argv)); break;
  case 1: flag=swigFunctionName(resc,resv,argc,(mxArray**)(argv)); break;
  case 2: flag=swigConstantName(resc,resv,argc,(mxArray**)(argv)); break;
  case 3: flag=swigThis(resc,resv,argc,(mxArray**)(argv)); break;
  case 4: flag=swigTouch(resc,resv,argc,(mxArray**)(argv)); break;
  case 5: flag=_wrap_new_AutoDiffVectorDynamic(resc,resv,argc,(mxArray**)(argv)); break;
  case 6: flag=_wrap_AutoDiffVectorDynamic_value(resc,resv,argc,(mxArray**)(argv)); break;
  case 7: flag=_wrap_AutoDiffVectorDynamic_derivatives(resc,resv,argc,(mxArray**)(argv)); break;
  case 8: flag=_wrap_AutoDiffVectorDynamic_arrayMultiply(resc,resv,argc,(mxArray**)(argv)); break;
  case 9: flag=_wrap_AutoDiffVectorDynamic_arrayDivide(resc,resv,argc,(mxArray**)(argv)); break;
  case 10: flag=_wrap_AutoDiffVectorDynamic_plus(resc,resv,argc,(mxArray**)(argv)); break;
  case 11: flag=_wrap_AutoDiffVectorDynamic_minus(resc,resv,argc,(mxArray**)(argv)); break;
  case 12: flag=_wrap_AutoDiffVectorDynamic_mtimes(resc,resv,argc,(mxArray**)(argv)); break;
  case 13: flag=_wrap_AutoDiffVectorDynamic_mrdivide(resc,resv,argc,(mxArray**)(argv)); break;
  case 14: flag=_wrap_AutoDiffVectorDynamic_rows(resc,resv,argc,(mxArray**)(argv)); break;
  case 15: flag=_wrap_AutoDiffVectorDynamic_cols(resc,resv,argc,(mxArray**)(argv)); break;
  case 16: flag=_wrap_AutoDiffVectorDynamic_size(resc,resv,argc,(mxArray**)(argv)); break;
  case 17: flag=_wrap_delete_AutoDiffVectorDynamic(resc,resv,argc,(mxArray**)(argv)); break;
  case 18: flag=_wrap_new_AutoDiffVectorMax73(resc,resv,argc,(mxArray**)(argv)); break;
  case 19: flag=_wrap_AutoDiffVectorMax73_value(resc,resv,argc,(mxArray**)(argv)); break;
  case 20: flag=_wrap_AutoDiffVectorMax73_derivatives(resc,resv,argc,(mxArray**)(argv)); break;
  case 21: flag=_wrap_AutoDiffVectorMax73_arrayMultiply(resc,resv,argc,(mxArray**)(argv)); break;
  case 22: flag=_wrap_AutoDiffVectorMax73_arrayDivide(resc,resv,argc,(mxArray**)(argv)); break;
  case 23: flag=_wrap_AutoDiffVectorMax73_plus(resc,resv,argc,(mxArray**)(argv)); break;
  case 24: flag=_wrap_AutoDiffVectorMax73_minus(resc,resv,argc,(mxArray**)(argv)); break;
  case 25: flag=_wrap_AutoDiffVectorMax73_mtimes(resc,resv,argc,(mxArray**)(argv)); break;
  case 26: flag=_wrap_AutoDiffVectorMax73_mrdivide(resc,resv,argc,(mxArray**)(argv)); break;
  case 27: flag=_wrap_AutoDiffVectorMax73_rows(resc,resv,argc,(mxArray**)(argv)); break;
  case 28: flag=_wrap_AutoDiffVectorMax73_cols(resc,resv,argc,(mxArray**)(argv)); break;
  case 29: flag=_wrap_AutoDiffVectorMax73_size(resc,resv,argc,(mxArray**)(argv)); break;
  case 30: flag=_wrap_delete_AutoDiffVectorMax73(resc,resv,argc,(mxArray**)(argv)); break;
  case 31: flag=_wrap_new_AutoDiffMatrixDynamic(resc,resv,argc,(mxArray**)(argv)); break;
  case 32: flag=_wrap_AutoDiffMatrixDynamic_value(resc,resv,argc,(mxArray**)(argv)); break;
  case 33: flag=_wrap_AutoDiffMatrixDynamic_derivatives(resc,resv,argc,(mxArray**)(argv)); break;
  case 34: flag=_wrap_AutoDiffMatrixDynamic_arrayMultiply(resc,resv,argc,(mxArray**)(argv)); break;
  case 35: flag=_wrap_AutoDiffMatrixDynamic_arrayDivide(resc,resv,argc,(mxArray**)(argv)); break;
  case 36: flag=_wrap_AutoDiffMatrixDynamic_plus(resc,resv,argc,(mxArray**)(argv)); break;
  case 37: flag=_wrap_AutoDiffMatrixDynamic_minus(resc,resv,argc,(mxArray**)(argv)); break;
  case 38: flag=_wrap_AutoDiffMatrixDynamic_mtimes(resc,resv,argc,(mxArray**)(argv)); break;
  case 39: flag=_wrap_AutoDiffMatrixDynamic_mrdivide(resc,resv,argc,(mxArray**)(argv)); break;
  case 40: flag=_wrap_AutoDiffMatrixDynamic_rows(resc,resv,argc,(mxArray**)(argv)); break;
  case 41: flag=_wrap_AutoDiffMatrixDynamic_cols(resc,resv,argc,(mxArray**)(argv)); break;
  case 42: flag=_wrap_AutoDiffMatrixDynamic_size(resc,resv,argc,(mxArray**)(argv)); break;
  case 43: flag=_wrap_delete_AutoDiffMatrixDynamic(resc,resv,argc,(mxArray**)(argv)); break;
  default: flag=1, SWIG_Error(SWIG_RuntimeError, "No function id %d.", fcn_id);
  }
  if (flag) {
    mexErrMsgIdAndTxt(SWIG_ErrorType(SWIG_lasterror_code), SWIG_lasterror_msg);
  }
}