[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