diff options
Diffstat (limited to 'packages/gsl/src')
-rw-r--r-- | packages/gsl/src/Numeric/GSL/gsl-ode.c | 15 |
1 files changed, 10 insertions, 5 deletions
diff --git a/packages/gsl/src/Numeric/GSL/gsl-ode.c b/packages/gsl/src/Numeric/GSL/gsl-ode.c index a6bdb55..72c8617 100644 --- a/packages/gsl/src/Numeric/GSL/gsl-ode.c +++ b/packages/gsl/src/Numeric/GSL/gsl-ode.c | |||
@@ -178,11 +178,16 @@ int ode(int method, int control, double h, | |||
178 | status = gsl_odeiv2_driver_apply (d, &t, ti, y); | 178 | status = gsl_odeiv2_driver_apply (d, &t, ti, y); |
179 | 179 | ||
180 | if (status != GSL_SUCCESS) { | 180 | if (status != GSL_SUCCESS) { |
181 | printf ("error in ode, return value=%d\n", status); | 181 | int k; |
182 | break; | 182 | printf ("error in ode, return value=%d\n", status); |
183 | } | 183 | printf("last successful values are:\n"); |
184 | 184 | printf("t = %.5e\n", t); | |
185 | // printf ("%.5e %.5e %.5e\n", t, y[0], y[1]); | 185 | for (k=0; k < xin; k++) |
186 | { | ||
187 | printf("y[%d] = %.5e\n", k, y[k]); | ||
188 | } | ||
189 | break; | ||
190 | } | ||
186 | 191 | ||
187 | for(j=0; j<xin; j++) { | 192 | for(j=0; j<xin; j++) { |
188 | solp[i*xin + j] = y[j]; | 193 | solp[i*xin + j] = y[j]; |