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-21 11:18:49 +01:00
parent 7b294c2aef
commit 899ed3c564
1 changed files with 12 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
wrk.patch
View File

@ -1,7 +1,14 @@
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-21 11:16:35.168162631 +0100
@@ -85,6 +85,7 @@ std_control_hadjust (void *vstate, size_
+++ wrk/ode-initval2/cstd.c 2013-01-21 11:18:36.213897243 +0100
@@ -80,11 +80,14 @@ std_control_init (void *vstate,
return GSL_SUCCESS;
}
+#define PD(fl) printf("%lu\n", *((unsigned long *)&(fl)))
+
static int
std_control_hadjust (void *vstate, size_t dim, unsigned int ord,
const double y[], const double yerr[], const double yp[],
double *h)
{
@ -9,7 +16,7 @@ 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,10 +110,12 @@ std_control_hadjust (void *vstate, size_
@@ -109,10 +112,12 @@ std_control_hadjust (void *vstate, size_
if (rmax > 1.1)
{
@ -23,7 +30,7 @@ diff -up wrk/ode-initval2/cstd.c.wrk wrk/ode-initval2/cstd.c
if (r < 0.2)
r = 0.2;
@@ -122,6 +125,7 @@ std_control_hadjust (void *vstate, size_
@@ -122,6 +127,7 @@ std_control_hadjust (void *vstate, size_
}
else if (rmax < 0.5)
{
@ -31,7 +38,7 @@ diff -up wrk/ode-initval2/cstd.c.wrk wrk/ode-initval2/cstd.c
/* increase step, no more than factor of 5 */
double r = S / pow (rmax, 1.0 / (ord + 1.0));
@@ -137,6 +141,7 @@ std_control_hadjust (void *vstate, size_
@@ -137,6 +143,7 @@ std_control_hadjust (void *vstate, size_
}
else
{