diff --git a/wrk.patch b/wrk.patch
index 1845e0c..d21866b 100644
--- a/wrk.patch
+++ b/wrk.patch
@@ -1,7 +1,14 @@
 diff -up wrk/ode-initval2/cstd.c.wrk wrk/ode-initval2/cstd.c
 --- wrk/ode-initval2/cstd.c.wrk	2013-01-18 16:02:19.566833381 +0100
-+++ wrk/ode-initval2/cstd.c	2013-01-21 11:16:35.168162631 +0100
-@@ -85,6 +85,7 @@ std_control_hadjust (void *vstate, size_
++++ wrk/ode-initval2/cstd.c	2013-01-21 11:18:36.213897243 +0100
+@@ -80,11 +80,14 @@ std_control_init (void *vstate,
+   return GSL_SUCCESS;
+ }
+ 
++#define PD(fl) printf("%lu\n", *((unsigned long *)&(fl)))
++
+ static int
+ std_control_hadjust (void *vstate, size_t dim, unsigned int ord,
                       const double y[], const double yerr[], const double yp[],
                       double *h)
  {
@@ -9,7 +16,7 @@ diff -up wrk/ode-initval2/cstd.c.wrk wrk/ode-initval2/cstd.c
    std_control_state_t *state = (std_control_state_t *) vstate;
  
    const double eps_abs = state->eps_abs;
-@@ -109,10 +110,12 @@ std_control_hadjust (void *vstate, size_
+@@ -109,10 +112,12 @@ std_control_hadjust (void *vstate, size_
  
    if (rmax > 1.1)
      {
@@ -23,7 +30,7 @@ diff -up wrk/ode-initval2/cstd.c.wrk wrk/ode-initval2/cstd.c
        if (r < 0.2)
          r = 0.2;
  
-@@ -122,6 +125,7 @@ std_control_hadjust (void *vstate, size_
+@@ -122,6 +127,7 @@ std_control_hadjust (void *vstate, size_
      }
    else if (rmax < 0.5)
      {
@@ -31,7 +38,7 @@ diff -up wrk/ode-initval2/cstd.c.wrk wrk/ode-initval2/cstd.c
        /* increase step, no more than factor of 5 */
        double r = S / pow (rmax, 1.0 / (ord + 1.0));
  
-@@ -137,6 +141,7 @@ std_control_hadjust (void *vstate, size_
+@@ -137,6 +143,7 @@ std_control_hadjust (void *vstate, size_
      }
    else
      {