1) First of all, can you describe in a bit more detailed way the usage of AppCtx class of matfree.py module to solve ODE systems (determined as in rober.py), without jacobian initialisation, in other words how can change rober.py to solve this issue?
2) Does THETA integration implement time step adaptation? 3) Suppose i have a large ODE system, how can i implement multiprocessor (parallel) integration in a way similar with those (function definition and plotting) in rober.py? -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100926/983112dd/attachment.htm>