temporary diagnostic

wrk.patch

@@ -162,7 +162,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
  
 diff -up wrk/ode-initval2/msbdf.c.wrk wrk/ode-initval2/msbdf.c
 --- wrk/ode-initval2/msbdf.c.wrk	2013-01-09 10:35:45.259960403 +0100
-+++ wrk/ode-initval2/msbdf.c	2013-01-23 17:56:07.543365882 +0100
++++ wrk/ode-initval2/msbdf.c	2013-01-23 18:03:53.982234648 +0100
 @@ -43,6 +43,8 @@
     framework.
  */
@@ -196,6 +196,18 @@ diff -up wrk/ode-initval2/msbdf.c.wrk wrk/ode-initval2/msbdf.c
    return GSL_SUCCESS;
  }
  
+@@ -1007,9 +1017,9 @@ msbdf_eval_order (gsl_vector * abscor, g
+   if (ordm1est > ordest && ordm1est > ordp1est && ordm1est > min_incr)
+     {
+       *ord -= 1;
+-#ifdef DEBUG
++//#ifdef DEBUG
+       printf ("-- eval_order order DECREASED to %d\n", (int) *ord);
+-#endif
++//#endif
+     }
+ 
+   else if (ordp1est > ordest && ordp1est > ordm1est && ordp1est > min_incr)
 @@ -1082,7 +1092,8 @@ msbdf_apply (void *vstate, size_t dim, d
               const gsl_odeiv2_system * sys)
  {