[gsl/temp_odecrash] temporary diagnostic

[Frantisek Kluknavsky]

commit 3d85b9fe791adfe9ba159b2f5a112f6a77ced5bd
Author: Frantisek Kluknavsky <fkluknav at redhat.com>
Date:   Thu Jan 17 14:20:57 2013 +0100

    temporary diagnostic

 gsl.spec  |    3 +-
 wrk.patch |   87 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++-
 2 files changed, 87 insertions(+), 3 deletions(-)
---
diff --git a/gsl.spec b/gsl.spec
index e22f9f2..0b385fa 100644
--- a/gsl.spec
+++ b/gsl.spec
@@ -34,7 +34,6 @@ The gsl-devel package contains the header files necessary for
 developing programs using the GSL (GNU Scientific Library).
 
 %prep
-#uname -a
 %setup -q
 %patch0 -p1 -b .lib64
 %patch1 -p1 -b .libs
@@ -46,7 +45,7 @@ mv THANKS.aux THANKS
 
 %build
 %ifarch %{ix86}
-%configure CFLAGS="$CFLAGS -fgnu89-inline -ffp-contract=off -msse2"
+%configure CFLAGS="$CFLAGS -fgnu89-inline -ffp-contract=off -msse2 -mfpmath=sse"
 %else
 %configure CFLAGS="$CFLAGS -fgnu89-inline -ffp-contract=off"
 %endif
diff --git a/wrk.patch b/wrk.patch
index 03677e2..8bdd0a1 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-17 13:21:54.742361043 +0100
++++ wrk/ode-initval2/test.c	2013-01-17 14:17:53.342044766 +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;
@@ -21,3 +21,88 @@ diff -up wrk/ode-initval2/test.c.wrk wrk/ode-initval2/test.c
    return GSL_SUCCESS;
  }
  
+@@ -1264,10 +1276,17 @@ 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));
+@@ -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 ());
+ }