summaryrefslogtreecommitdiff
path: root/packages/sundials/src
diff options
context:
space:
mode:
Diffstat (limited to 'packages/sundials/src')
-rw-r--r--packages/sundials/src/Numeric/Sundials/ARKode/ODE.hs19
1 files changed, 0 insertions, 19 deletions
diff --git a/packages/sundials/src/Numeric/Sundials/ARKode/ODE.hs b/packages/sundials/src/Numeric/Sundials/ARKode/ODE.hs
index 2ede918..44b724e 100644
--- a/packages/sundials/src/Numeric/Sundials/ARKode/ODE.hs
+++ b/packages/sundials/src/Numeric/Sundials/ARKode/ODE.hs
@@ -263,25 +263,6 @@ solveOdeC method relTol absTol fun f0 ts = unsafePerformIO $ do
263 flag = ARKodeSetIRKTableNum(arkode_mem, $(int method)); 263 flag = ARKodeSetIRKTableNum(arkode_mem, $(int method));
264 if (check_flag(&flag, "ARKode", 1)) return 1; 264 if (check_flag(&flag, "ARKode", 1)) return 1;
265 265
266 int s, q, p;
267 realtype *ai = (realtype *)malloc(ARK_S_MAX * ARK_S_MAX * sizeof(realtype));
268 realtype *ae = (realtype *)malloc(ARK_S_MAX * ARK_S_MAX * sizeof(realtype));
269 realtype *ci = (realtype *)malloc(ARK_S_MAX * sizeof(realtype));
270 realtype *ce = (realtype *)malloc(ARK_S_MAX * sizeof(realtype));
271 realtype *bi = (realtype *)malloc(ARK_S_MAX * sizeof(realtype));
272 realtype *be = (realtype *)malloc(ARK_S_MAX * sizeof(realtype));
273 realtype *b2i = (realtype *)malloc(ARK_S_MAX * sizeof(realtype));
274 realtype *b2e = (realtype *)malloc(ARK_S_MAX * sizeof(realtype));
275 flag = ARKodeGetCurrentButcherTables(arkode_mem, &s, &q, &p, ai, ae, ci, ce, bi, be, b2i, b2e);
276 if (check_flag(&flag, "ARKode", 1)) return 1;
277 fprintf(stderr, "s = %d, q = %d, p = %d\n", s, q, p);
278 for (i = 0; i < s; i++) {
279 for (j = 0; j < s; j++) {
280 fprintf(stderr, "ai[%d,%d] = %f", i, j, ai[i * ARK_S_MAX + j]);
281 }
282 fprintf(stderr, "\n");
283 }
284
285 /* Main time-stepping loop: calls ARKode to perform the integration */ 266 /* Main time-stepping loop: calls ARKode to perform the integration */
286 /* Stops when the final time has been reached */ 267 /* Stops when the final time has been reached */
287 for (i = 1; i < $(int nTs); i++) { 268 for (i = 1; i < $(int nTs); i++) {