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-16 13:45:46 +01:00
parent e165263e72
commit 4ed85fc02a
1 changed files with 6 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
wrk.patch
View File

@ -1,7 +1,7 @@
diff -up wrk/ode-initval2/msbdf.c.wrk wrk/ode-initval2/msbdf.c 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 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.wrk 2013-01-09 10:48:22.051415928 +0100
+++ wrk/ode-initval2/test.c 2013-01-16 13:43:15.929791072 +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[], @@ -923,6 +923,18 @@ rhs_ringmod (double t, const double y[],
f[13] = (-y[0] + uin1 - (ri + rg1) * y[13]) / ls1; f[13] = (-y[0] + uin1 - (ri + rg1) * y[13]) / ls1;
f[14] = (-y[1] - (rc + rg1) * y[14]) / ls1; f[14] = (-y[1] - (rc + rg1) * y[14]) / ls1;
@ -31,6 +31,9 @@ diff -up wrk/ode-initval2/test.c.wrk wrk/ode-initval2/test.c
s = gsl_odeiv2_evolve_apply (d->e, d->c, d->s, sys, &t, t1, &h, y); s = gsl_odeiv2_evolve_apply (d->e, d->c, d->s, sys, &t, t1, &h, y);
-#ifdef DEBUG -#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); + printf("krok %d\n\n", ++poc);
+ //int k; + //int k;
+ //for (k=0; k<15; ++k) { + //for (k=0; k<15; ++k) {
@ -38,9 +41,8 @@ diff -up wrk/ode-initval2/test.c.wrk wrk/ode-initval2/test.c
+ //} + //}
+ //printf("\n"); + //printf("\n");
+//#ifdef DEBUG +//#ifdef DEBUG
printf ("%.5e %.5e %.5e %d\n", t, y[0], y[1], + //printf ("%.5e %.5e %.5e %d\n", t, y[0], y[1],
gsl_odeiv2_step_order (d->s)); + printf("%d\n", gsl_odeiv2_step_order (d->s));
-#endif
+//#endif +//#endif
if (s != GSL_SUCCESS) if (s != GSL_SUCCESS)