[gsl/temp_odecrash] temporary diagnostic

[Frantisek Kluknavsky]

commit 973e1d57de55bea6a85bc7b41b08db6664259820
Author: Frantisek Kluknavsky <fkluknav at redhat.com>
Date:   Wed Jan 16 14:02:55 2013 +0100

    temporary diagnostic

 wrk.patch |  117 +------------------------------------------------------------
 1 files changed, 1 insertions(+), 116 deletions(-)
---
diff --git a/wrk.patch b/wrk.patch
index 00a4665..8b13789 100644
--- a/wrk.patch
+++ b/wrk.patch
@@ -1,116 +1 @@
-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
-@@ -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;
- 
-+  /*int i;
-+  printf("temp states\n");
-+  for (i=0; i<15; ++i) {
-+  
-+	  printf("%lu\n", *((unsigned long *)(&(y[i]))));
-+  }
-+
-+	  printf("%lu\n", *((unsigned long *)(&(rg2))));
-+	  printf("%lu\n", *((unsigned long *)(&(ls2))));
-+	  printf("%lu\n", *((unsigned long *)(&(rg3))));
-+	  printf("%lu\n", *((unsigned long *)(&(ls3))));
-+*/
-   return GSL_SUCCESS;
- }
- 
-@@ -1264,14 +1276,21 @@ 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);
- 
--#ifdef DEBUG
--      printf ("%.5e %.5e %.5e %d\n", t, y[0], y[1],
--              gsl_odeiv2_step_order (d->s));
--#endif
-+      printf("krok %d\n\n", ++poc);
-+      //int k;
-+      //for (k=0; k<15; ++k) {
-+	//      printf(" %lu",(*(unsigned long int *)(&(y[k]))));
-+      //}
-+      //printf("\n");
-+//#ifdef DEBUG
-+      //printf ("%.5e %.5e %.5e %d\n", t, y[0], y[1],
-+      printf("%d\n", 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 ());
- }
+