[gsl/temp_odecrash] temporary diagnostic

[Frantisek Kluknavsky]

commit ab84e8557b74cbab761d5234004a28123c2b7d83
Author: Frantisek Kluknavsky <fkluknav at redhat.com>
Date:   Mon Jan 21 11:28:33 2013 +0100

    temporary diagnostic

 wrk.patch |   23 ++++++++++++++---------
 1 files changed, 14 insertions(+), 9 deletions(-)
---
diff --git a/wrk.patch b/wrk.patch
index d21866b..b778fa8 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-21 11:18:36.213897243 +0100
++++ wrk/ode-initval2/cstd.c	2013-01-21 11:28:19.857402239 +0100
 @@ -80,11 +80,14 @@ std_control_init (void *vstate,
    return GSL_SUCCESS;
  }
@@ -16,21 +16,26 @@ 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 +112,12 @@ std_control_hadjust (void *vstate, size_
+@@ -109,9 +112,17 @@ 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);
--
-+printf ("r\n");
-+PD(r);
+-      double r = S / pow (rmax, 1.0 / ord);
++	 double pom = pow (rmax, 1.0 / ord);
++      double r = S / pom;
++printf ("S\n");
++PD(S);
++printf("rmax\n");
++PD(rmax);
++printf("pom\n");
++PD(pom);
+ 
        if (r < 0.2)
          r = 0.2;
- 
-@@ -122,6 +127,7 @@ std_control_hadjust (void *vstate, size_
+@@ -122,6 +133,7 @@ std_control_hadjust (void *vstate, size_
      }
    else if (rmax < 0.5)
      {
@@ -38,7 +43,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 +143,7 @@ std_control_hadjust (void *vstate, size_
+@@ -137,6 +149,7 @@ std_control_hadjust (void *vstate, size_
      }
    else
      {