[gsl/temp_odecrash] temporary diagnostic

Frantisek Kluknavsky fkluknav at fedoraproject.org
Mon Jan 21 10:19:06 UTC 2013


commit 7b294c2aef79e8cb5619ad1f17253d65464d7ca2
Author: Frantisek Kluknavsky <fkluknav at redhat.com>
Date:   Mon Jan 21 11:16:56 2013 +0100

    temporary diagnostic

 wrk.patch |   16 +++++++++++-----
 1 files changed, 11 insertions(+), 5 deletions(-)
---
diff --git a/wrk.patch b/wrk.patch
index 52b2daa..1845e0c 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:01:43.136808297 +0100
++++ wrk/ode-initval2/cstd.c	2013-01-21 11:16:35.168162631 +0100
 @@ -85,6 +85,7 @@ std_control_hadjust (void *vstate, size_
                       const double y[], const double yerr[], const double yp[],
                       double *h)
@@ -9,7 +9,7 @@ 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_
+@@ -109,10 +110,12 @@ std_control_hadjust (void *vstate, size_
  
    if (rmax > 1.1)
      {
@@ -17,15 +17,21 @@ diff -up wrk/ode-initval2/cstd.c.wrk wrk/ode-initval2/cstd.c
        /* 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_
+-
++printf ("r\n");
++PD(r);
+       if (r < 0.2)
+         r = 0.2;
+ 
+@@ -122,6 +125,7 @@ std_control_hadjust (void *vstate, size_
      }
    else if (rmax < 0.5)
      {
-+	    printf("predlzujem");
++	    printf("predlzujem\n");
        /* increase step, no more than factor of 5 */
        double r = S / pow (rmax, 1.0 / (ord + 1.0));
  
-@@ -137,6 +140,7 @@ std_control_hadjust (void *vstate, size_
+@@ -137,6 +141,7 @@ std_control_hadjust (void *vstate, size_
      }
    else
      {


More information about the scm-commits mailing list