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:10:25 +01:00
parent e5bf834570
commit c7b509e575
1 changed files with 6 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
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 16:42:55.405445512 +0100
+++ wrk/ode-initval2/evolve.c 2013-01-23 17:10:19.668348408 +0100
@@ -28,6 +28,8 @@
#include "odeiv_util.h"
@ -107,7 +107,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
goto try_step;
}
else
@@ -228,7 +234,8 @@ try_step:
@@ -228,16 +234,19 @@ try_step:
{
*t = t0 + h0;
}
@ -117,7 +117,10 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
if (con != NULL)
{
/* Check error and attempt to adjust the step. */
@@ -238,6 +245,8 @@ try_step:
double h_old = h0;
-
+ printf("pred hadjust vonku\n%d\n", s->type->order (s->state));
const int hadjust_status
=
gsl_odeiv2_control_hadjust (con, step, y, e->yerr, e->dydt_out, &h0);