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-17 20:00:38 +01:00
parent 9dbe497d9c
commit b6ceb62e81
1 changed files with 12 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
wrk.patch
View File

@ -1,6 +1,6 @@
diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.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.wrk 2013-01-17 18:48:06.473062671 +0100
+++ wrk/ode-initval2/evolve.c 2013-01-17 19:53:40.837541111 +0100 +++ wrk/ode-initval2/evolve.c 2013-01-17 20:00:18.350164952 +0100
@@ -28,6 +28,8 @@ @@ -28,6 +28,8 @@
#include "odeiv_util.h" #include "odeiv_util.h"
@ -47,7 +47,16 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
if (con != NULL) if (con != NULL)
{ {
/* Check error and attempt to adjust the step. */ /* Check error and attempt to adjust the step. */
@@ -254,6 +261,8 @@ try_step: @@ -238,6 +245,8 @@ try_step:
const int hadjust_status
=
gsl_odeiv2_control_hadjust (con, step, y, e->yerr, e->dydt_out, &h0);
+printf("po hadjust\n");
+PD(h0);
if (hadjust_status == GSL_ODEIV_HADJ_DEC)
{
@@ -254,6 +263,8 @@ try_step:
DBL_MEMCPY (y, e->y0, dydt->dimension); DBL_MEMCPY (y, e->y0, dydt->dimension);
e->failed_steps++; e->failed_steps++;
@ -56,7 +65,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
goto try_step; goto try_step;
} }
else else
@@ -279,6 +288,9 @@ try_step: @@ -279,6 +290,9 @@ try_step:
*h = h0; *h = h0;
} }