Chris,
Take a look in the petsc branch barry/wirth-fusion-materials
https://bitbucket.org/petsc/petsc/branch/barry%2Fwirth-fusion-materials
in the file src/ts/examples/tutorials/advection-diffusion-reaction/ex10.c
feel free to ask any questions.
Barry
On Oct 7, 2013, at 3:03 AM, Christophe Ortiz <[email protected]> wrote:
> Hi all
>
> I need some suggestions to design a code with PETSc.
>
> I want to solve a 1D problem composed of several diffusion equations and a
> lot of ODEs (thousands).
>
> - After discretization, the diffusion equations have terms in i-1, i and i+1.
>
> - The diffusion equations have additional terms due to couping with all the
> ODEs. These terms are non-linear.
>
> - The terms of the ODEs are local (only depend on node i) but are non-linear.
>
> Any advice to design the code ?
>
> How should I distribute the terms between IFunction and RHSFunction ?
>
> Any special attention to DMDA ? Should I declare just one DMDA with dof >>1 ?
>
> Many thanks in advance !
> Christophe
>
