dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Tue, Apr 26, 2011 at 10:23:35AM -0700, Johan Hake wrote:
> Hello!
>
> (Sorry for the spamming, the previous email was prematurely sent...)
>
> I am about to use DOLFIN to solve a set of distrinct ODEs. The ODEs are a
> result of an operator splitting of a PDE, avoiding an expensive reassemble
> each newton itteration.
>
> Each ODE is potentially small but it needs to be solved for each node on a
> mesh. Each one of the ODEs are decoupled from eachother.
>
> What do you think would be fastest:
>
>   1) To solve each small ODE using a DenseMatrix and a direct solver
>   2) Put all the small dense matrices into a big sparse one and go for an
>      iterative solver.
>
> For the last one I need to construct my own assembler, iterating over the
> degrees of freedom.

I think Option (1) since it can be easily and efficiently
parallelized.

You could also consider checking out the ODECollection class.

--
Anders