From 7b294c2aef79e8cb5619ad1f17253d65464d7ca2 Mon Sep 17 00:00:00 2001 From: Frantisek Kluknavsky Date: Mon, 21 Jan 2013 11:16:56 +0100 Subject: [PATCH] temporary diagnostic --- wrk.patch | 16 +++++++++++----- 1 file changed, 11 insertions(+), 5 deletions(-) diff --git a/wrk.patch b/wrk.patch index 52b2daa..1845e0c 100644 --- a/wrk.patch +++ b/wrk.patch @@ -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-21 11:01:43.136808297 +0100 ++++ wrk/ode-initval2/cstd.c 2013-01-21 11:16:35.168162631 +0100 @@ -85,6 +85,7 @@ std_control_hadjust (void *vstate, size_ const double y[], const double yerr[], const double yp[], double *h) @@ -9,7 +9,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,6 +110,7 @@ std_control_hadjust (void *vstate, size_ +@@ -109,10 +110,12 @@ std_control_hadjust (void *vstate, size_ if (rmax > 1.1) { @@ -17,15 +17,21 @@ diff -up wrk/ode-initval2/cstd.c.wrk wrk/ode-initval2/cstd.c /* 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_ +- ++printf ("r\n"); ++PD(r); + if (r < 0.2) + r = 0.2; + +@@ -122,6 +125,7 @@ std_control_hadjust (void *vstate, size_ } else if (rmax < 0.5) { -+ printf("predlzujem"); ++ printf("predlzujem\n"); /* increase step, no more than factor of 5 */ double r = S / pow (rmax, 1.0 / (ord + 1.0)); -@@ -137,6 +140,7 @@ std_control_hadjust (void *vstate, size_ +@@ -137,6 +141,7 @@ std_control_hadjust (void *vstate, size_ } else {