rpms/gsl
3
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

temporary diagnostic

Browse Source
This commit is contained in:
Frantisek Kluknavsky 2013-01-23 17:26:53 +01:00
parent ad9e05d57f
commit bcbc6cc329
1 changed files with 13 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
wrk.patch
View File

@ -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:19:08.012618826 +0100
+++ wrk/ode-initval2/evolve.c 2013-01-23 17:26:39.628412797 +0100
@@ -28,6 +28,8 @@
#include "odeiv_util.h"
@ -99,7 +99,15 @@ 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;
@@ -206,6 +211,7 @@ try_step:
@@ -167,6 +172,7 @@ try_step:
step_status =
gsl_odeiv2_step_apply (step, t0, h0, y, e->yerr, e->dydt_in,
e->dydt_out, dydt);
+ printf("po stepapply vonku\n%d\n", step->type->order (step->state));
}
else
{
@@ -206,6 +212,7 @@ try_step:
DBL_MEMCPY (y, e->y0, dydt->dimension);
e->failed_steps++;
@ -107,7 +115,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
goto try_step;
}
else
@@ -228,16 +234,19 @@ try_step:
@@ -228,16 +235,19 @@ try_step:
{
*t = t0 + h0;
}
@ -129,7 +137,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
if (hadjust_status == GSL_ODEIV_HADJ_DEC)
{
@@ -254,6 +263,8 @@ try_step:
@@ -254,6 +264,8 @@ try_step:
DBL_MEMCPY (y, e->y0, dydt->dimension);
e->failed_steps++;
@ -138,7 +146,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
goto try_step;
}
else
@@ -279,6 +290,9 @@ try_step:
@@ -279,6 +291,9 @@ try_step:
*h = h0;
}