 
         double *const qr     = start;          /* meqn   */
         double *const qd     = qr + meqn;      /* meqn   */
-        double *const bpdq   = qd + meqn;    /* meqn   */
+        double *const bpdq   = qd + meqn;      /* meqn   */
         double *const aux1   = bpdq + meqn;    /* 2*maux */
         double *const aux2   = aux1 + 2*maux;  /* 2*maux */
         double *const aux3   = aux2 + 2*maux;  /* 2*maux */


        double *const qr     = start;          /* meqn   */
        double *const qd     = qr + meqn;      /* meqn   */
        double *const bpdq   = qd + meqn;    /* meqn   */
        double *const aux1   = bpdq + meqn;    /* 2*maux */
        double *const aux2   = aux1 + 2*maux;  /* 2*maux */
        double *const aux3   = aux2 + 2*maux;  /* 2*maux */


        double *const qr     = start;          /* meqn   */
        double *const ql     = qr + meqn;      /* meqn   */
        double *const apdq   = ql + meqn;      /* meqn   */
        double *const aux1   = apdq + meqn;    /* 2*maux */
        double *const aux2   = aux1 + 2*maux;  /* 2*maux */
        double *const aux3   = aux2 + 2*maux;  /* 2*maux */
