diff options
Diffstat (limited to 'packages/sundials/src')
-rw-r--r-- | packages/sundials/src/Numeric/Sundials/ARKode/ODE.hs | 19 |
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++) { |