Christophe Ortiz <[email protected]> writes: > 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.
How do you want to extend your code? Is it acceptable to have a centralized struct where you place all variables, or should separate modules be able to "register" additional variables?
pgpvYBGTqz1jY.pgp
Description: PGP signature
