[gsl/temp_odecrash] temporary diagnostic

Frantisek Kluknavsky fkluknav at fedoraproject.org
Wed Jan 23 16:30:40 UTC 2013


commit df114a92224aafee5f107880e86e5d734d3774a2
Author: Frantisek Kluknavsky <fkluknav at redhat.com>
Date:   Wed Jan 23 17:30:32 2013 +0100

    temporary diagnostic

 wrk.patch |   16 ++++++++++------
 1 files changed, 10 insertions(+), 6 deletions(-)
---
diff --git a/wrk.patch b/wrk.patch
index 656bf4b..459ad8c 100644
--- a/wrk.patch
+++ b/wrk.patch
@@ -70,7 +70,7 @@ diff -up wrk/ode-initval2/cstd.c.wrk wrk/ode-initval2/cstd.c
      }
 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-23 17:26:39.628412797 +0100
++++ wrk/ode-initval2/evolve.c	2013-01-23 17:30:24.141816833 +0100
 @@ -28,6 +28,8 @@
  
  #include "odeiv_util.h"
@@ -99,7 +99,11 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
    if ((dt >= 0.0 && h0 > dt) || (dt < 0.0 && h0 < dt))
      {
        h0 = dt;
-@@ -167,6 +172,7 @@ try_step:
+@@ -164,9 +169,11 @@ try_step:
+ 
+   if (step->type->can_use_dydt_in)
+     {
++	    printf("pred stepapply vonku\n%d\n",  step->type->order (step->state));
        step_status =
          gsl_odeiv2_step_apply (step, t0, h0, y, e->yerr, e->dydt_in,
                                 e->dydt_out, dydt);
@@ -107,7 +111,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
      }
    else
      {
-@@ -206,6 +212,7 @@ try_step:
+@@ -206,6 +213,7 @@ try_step:
              
              DBL_MEMCPY (y, e->y0, dydt->dimension);
              e->failed_steps++;
@@ -115,7 +119,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
              goto try_step;
            }
          else
-@@ -228,16 +235,19 @@ try_step:
+@@ -228,16 +236,19 @@ try_step:
      {
        *t = t0 + h0;
      }
@@ -137,7 +141,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
  
        if (hadjust_status == GSL_ODEIV_HADJ_DEC)
          {
-@@ -254,6 +264,8 @@ try_step:
+@@ -254,6 +265,8 @@ try_step:
  
                DBL_MEMCPY (y, e->y0, dydt->dimension);
                e->failed_steps++;
@@ -146,7 +150,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
                goto try_step;
              }
            else
-@@ -279,6 +291,9 @@ try_step:
+@@ -279,6 +292,9 @@ try_step:
        *h = h0;
      }
  


More information about the scm-commits mailing list