summaryrefslogtreecommitdiff
path: root/packages/gsl
diff options
context:
space:
mode:
authorDominic Steinitz <dominic@steinitz.org>2017-09-12 17:20:55 +0100
committerDominic Steinitz <dominic@steinitz.org>2017-09-12 17:20:55 +0100
commit8c8b8e452974897fe9272fcbcbf7dd9d623feaf0 (patch)
treeda02d92106f17940e8fdcb4cee436c0cabbd502e /packages/gsl
parent86f539fe5edf0377b40146c7c0feaa43c9780720 (diff)
Provide better error message
Diffstat (limited to 'packages/gsl')
-rw-r--r--packages/gsl/src/Numeric/GSL/gsl-ode.c15
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];