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 19:35:20 +01:00
parent 671930b335
commit 21c386cc36
1 changed files with 14 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
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:05:07.589629814 +0100 +++ wrk/ode-initval2/evolve.c 2013-01-17 19:35:12.739355883 +0100
@@ -28,6 +28,8 @@ @@ -28,6 +28,8 @@
#include "odeiv_util.h" #include "odeiv_util.h"
@ -19,7 +19,17 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
const double t0 = *t; const double t0 = *t;
double h0 = *h; double h0 = *h;
int step_status; int step_status;
@@ -206,6 +210,7 @@ try_step: @@ -151,7 +155,8 @@ gsl_odeiv2_evolve_apply (gsl_odeiv2_evol
}
try_step:
-
+printf("po navasti\n");
+PD(h0);
if ((dt >= 0.0 && h0 > dt) || (dt < 0.0 && h0 < dt))
{
h0 = dt;
@@ -206,6 +211,7 @@ try_step:
DBL_MEMCPY (y, e->y0, dydt->dimension); DBL_MEMCPY (y, e->y0, dydt->dimension);
e->failed_steps++; e->failed_steps++;
@ -27,7 +37,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
goto try_step; goto try_step;
} }
else else
@@ -254,6 +259,7 @@ try_step: @@ -254,6 +260,7 @@ try_step:
DBL_MEMCPY (y, e->y0, dydt->dimension); DBL_MEMCPY (y, e->y0, dydt->dimension);
e->failed_steps++; e->failed_steps++;
@ -35,7 +45,7 @@ diff -up wrk/ode-initval2/evolve.c.wrk wrk/ode-initval2/evolve.c
goto try_step; goto try_step;
} }
else else
@@ -279,6 +285,9 @@ try_step: @@ -279,6 +286,9 @@ try_step:
*h = h0; *h = h0;
} }