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-18 16:10:40 +01:00
parent 161deccb61
commit a1291cc4bd
1 changed files with 17 additions and 1 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/cstd.c.wrk wrk/ode-initval2/cstd.c
--- wrk/ode-initval2/cstd.c.wrk 2013-01-18 16:02:19.566833381 +0100
+++ wrk/ode-initval2/cstd.c 2013-01-18 16:04:15.023570525 +0100
+++ wrk/ode-initval2/cstd.c 2013-01-18 16:10:37.126101180 +0100
@@ -85,6 +85,7 @@ std_control_hadjust (void *vstate, size_
const double y[], const double yerr[], const double yp[],
double *h)
@ -9,6 +9,22 @@ diff -up wrk/ode-initval2/cstd.c.wrk wrk/ode-initval2/cstd.c
std_control_state_t *state = (std_control_state_t *) vstate;
const double eps_abs = state->eps_abs;
@@ -109,6 +110,7 @@ std_control_hadjust (void *vstate, size_
if (rmax > 1.1)
{
+ printf("skracujem\n");
/* decrease step, no more than factor of 5, but a fraction S more
than scaling suggests (for better accuracy) */
double r = S / pow (rmax, 1.0 / ord);
@@ -122,6 +124,7 @@ std_control_hadjust (void *vstate, size_
}
else if (rmax < 0.5)
{
+ printf("predlzujem");
/* increase step, no more than factor of 5 */
double r = S / pow (rmax, 1.0 / (ord + 1.0));
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-17 20:00:18.350164952 +0100