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:52:12.141920679 +0100
++++ wrk/ode-initval2/msbdf.c	2013-01-23 17:56:07.543365882 +0100
 @@ -43,6 +43,8 @@
     framework.
  */
@@ -244,6 +244,15 @@ diff -up wrk/ode-initval2/msbdf.c.wrk wrk/ode-initval2/msbdf.c
    if (state->ordwait == 0)
      {
        msbdf_eval_order (abscor, tempvec, svec, errcoeff, dim, errlev,
+@@ -1626,7 +1641,7 @@ msbdf_apply (void *vstate, size_t dim, d
+                         state->ordp1coeffprev, ordp2coeff,
+                         hprev, h, z, &(state->ord), &(state->ordwait));
+     }
+-
++printf("ord after eval\n%d\n", state->ord); 
+   /* Undo scaling of abscor for possible order increase on next step */
+   {
+     size_t i;
 @@ -1675,7 +1690,9 @@ msbdf_apply (void *vstate, size_t dim, d
    printf ("-- nJ=%d, nM=%d\n", (int) state->nJ, (int) state->nM);
  #endif