47 lines
1.3 KiB
Diff
47 lines
1.3 KiB
Diff
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 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;
|
|
|
|
+ /*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,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)
|
|
|