[gsl/temp_odecrash] temporary diagnostic
Frantisek Kluknavsky
fkluknav at fedoraproject.org
Wed Jan 16 14:26:52 UTC 2013
commit 829d7f4a32f242c2dd6aec7af2730010bc1f55f1
Author: Frantisek Kluknavsky <fkluknav at redhat.com>
Date: Wed Jan 16 15:26:32 2013 +0100
temporary diagnostic
wrk.patch | 92 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++-
1 files changed, 91 insertions(+), 1 deletions(-)
---
diff --git a/wrk.patch b/wrk.patch
index 7124513..f0fb758 100644
--- a/wrk.patch
+++ b/wrk.patch
@@ -1,7 +1,7 @@
diff -up wrk/ode-initval2/msbdf.c.wrk wrk/ode-initval2/msbdf.c
diff -up wrk/ode-initval2/test.c.wrk wrk/ode-initval2/test.c
--- wrk/ode-initval2/test.c.wrk 2013-01-09 10:48:22.051415928 +0100
-+++ wrk/ode-initval2/test.c 2013-01-16 13:45:32.047667419 +0100
++++ wrk/ode-initval2/test.c 2013-01-16 15:26:17.505976018 +0100
@@ -923,6 +923,18 @@ rhs_ringmod (double t, const double y[],
f[13] = (-y[0] + uin1 - (ri + rg1) * y[13]) / ls1;
f[14] = (-y[1] - (rc + rg1) * y[14]) / ls1;
@@ -20,3 +20,93 @@ diff -up wrk/ode-initval2/test.c.wrk wrk/ode-initval2/test.c
+*/
return GSL_SUCCESS;
}
+
+@@ -1264,13 +1276,20 @@ sys_driver (const gsl_odeiv2_step_type *
+ nfe = 0;
+ nje = 0;
+
++ int poc=0;
+ while (t < t1)
+ {
+ s = gsl_odeiv2_evolve_apply (d->e, d->c, d->s, sys, &t, t1, &h, y);
+
++ printf("\nkrok %d\n", ++poc);
++ int k;
++ for (k=0; k<15; ++k) {
++ printf(" %lu\n",(*(unsigned long int *)(&(y[k]))));
++ }
++ //printf("\n");
+ #ifdef DEBUG
+ printf ("%.5e %.5e %.5e %d\n", t, y[0], y[1],
+- gsl_odeiv2_step_order (d->s));
++ gsl_odeiv2_step_order (d->s));
+ #endif
+
+ if (s != GSL_SUCCESS)
+@@ -1959,7 +1978,7 @@ test_extreme_problems (void)
+
+ /* Loop over problems */
+
+- for (p = 0; p < CONST_EXTREME_NPROB; p++)
++ for (p = 2; p < CONST_EXTREME_NPROB; p++)
+ {
+ /* Initialize */
+
+@@ -2001,12 +2020,13 @@ test_extreme_problems (void)
+
+ /* Call each solver for the problem */
+
+- for (i = 0; steppers[i] != 0; i++)
++ for (i = 1; steppers[i] != 0; i++)
+ {
++ printf("spustam driver p=%ld i=%ld\n", p, i);
+ int s = sys_driver (steppers[i], prob[p], start[p], end[p],
+ initstepsize[p], &y[sd[p] * i],
+ epsabs[p], epsrel[p], probname[p]);
+-
++ printf("koniec drivera\n");
+ if (s != GSL_SUCCESS)
+ {
+ printf ("start=%.5e, initstepsize=%.5e\n", start[p],
+@@ -2450,7 +2470,7 @@ main (void)
+
+ /* Basic tests for all steppers */
+
+- for (i = 0; p[i].type != 0; i++)
++ /*for (i = 0; p[i].type != 0; i++)
+ {
+ test_stepper (p[i].type);
+ }
+@@ -2469,9 +2489,9 @@ main (void)
+ test_stepsize_fail (p[i].type, p[i].h);
+ test_user_break (p[i].type, p[i].h);
+ }
+-
++*/
+ /* Derivative test for explicit steppers */
+-
++/*
+ explicit_stepper[0].type = gsl_odeiv2_step_rk4;
+ explicit_stepper[0].h = 1.0e-3;
+ explicit_stepper[1].type = gsl_odeiv2_step_rk2;
+@@ -2491,16 +2511,16 @@ main (void)
+ test_stepfn (explicit_stepper[i].type);
+ test_stepfn2 (explicit_stepper[i].type);
+ }
+-
++*/
+ /* Special tests */
+
+- test_nonstiff_problems ();
++ // test_nonstiff_problems ();
+
+- test_stiff_problems ();
++ // test_stiff_problems ();
+
+ test_extreme_problems ();
+
+- test_driver ();
++ //test_driver ();
+
+ exit (gsl_test_summary ());
+ }
More information about the scm-commits
mailing list