[gsl/temp_odecrash] temporary diagnostic

[Frantisek Kluknavsky]

commit a1291cc4bd726e9bb5a22296a9563efdf7b89853
Author: Frantisek Kluknavsky <fkluknav at redhat.com>
Date:   Fri Jan 18 16:10:40 2013 +0100

    temporary diagnostic

 wrk.patch |   18 +++++++++++++++++-
 1 files changed, 17 insertions(+), 1 deletions(-)
---
diff --git a/wrk.patch b/wrk.patch
index 6557fb1..43c24fb 100644
--- a/wrk.patch
+++ b/wrk.patch
@@ -1,6 +1,6 @@
 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-18 16:04:15.023570525 +0100
++++ wrk/ode-initval2/cstd.c	2013-01-18 16:10:37.126101180 +0100
 @@ -85,6 +85,7 @@ std_control_hadjust (void *vstate, size_
                       const double y[], const double yerr[], const double yp[],
                       double *h)
@@ -9,6 +9,22 @@ 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,6 +110,7 @@ std_control_hadjust (void *vstate, size_
+ 
+   if (rmax > 1.1)
+     {
++	    printf("skracujem\n");
+       /* decrease step, no more than factor of 5, but a fraction S more
+          than scaling suggests (for better accuracy) */
+       double r = S / pow (rmax, 1.0 / ord);
+@@ -122,6 +124,7 @@ std_control_hadjust (void *vstate, size_
+     }
+   else if (rmax < 0.5)
+     {
++	    printf("predlzujem");
+       /* increase step, no more than factor of 5 */
+       double r = S / pow (rmax, 1.0 / (ord + 1.0));
+ 
 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-17 20:00:18.350164952 +0100