diff --git a/wrk.patch b/wrk.patch index ef9bb32..656bf4b 100644 --- a/wrk.patch +++ b/wrk.patch @@ -70,7 +70,7 @@ diff -up wrk/ode-initval2/cstd.c.wrk wrk/ode-initval2/cstd.c } diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c --- wrk/ode-initval2/evolve.c.wrk 2013-01-17 18:48:06.473062671 +0100 -+++ wrk/ode-initval2/evolve.c 2013-01-23 17:19:08.012618826 +0100 ++++ wrk/ode-initval2/evolve.c 2013-01-23 17:26:39.628412797 +0100 @@ -28,6 +28,8 @@ #include "odeiv_util.h" @@ -99,7 +99,15 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c if ((dt >= 0.0 && h0 > dt) || (dt < 0.0 && h0 < dt)) { h0 = dt; -@@ -206,6 +211,7 @@ try_step: +@@ -167,6 +172,7 @@ try_step: + step_status = + gsl_odeiv2_step_apply (step, t0, h0, y, e->yerr, e->dydt_in, + e->dydt_out, dydt); ++ printf("po stepapply vonku\n%d\n", step->type->order (step->state)); + } + else + { +@@ -206,6 +212,7 @@ try_step: DBL_MEMCPY (y, e->y0, dydt->dimension); e->failed_steps++; @@ -107,7 +115,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c goto try_step; } else -@@ -228,16 +234,19 @@ try_step: +@@ -228,16 +235,19 @@ try_step: { *t = t0 + h0; } @@ -129,7 +137,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c if (hadjust_status == GSL_ODEIV_HADJ_DEC) { -@@ -254,6 +263,8 @@ try_step: +@@ -254,6 +264,8 @@ try_step: DBL_MEMCPY (y, e->y0, dydt->dimension); e->failed_steps++; @@ -138,7 +146,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c goto try_step; } else -@@ -279,6 +290,9 @@ try_step: +@@ -279,6 +291,9 @@ try_step: *h = h0; }