dolfin team mailing list archive

[Johan Hake <johan.hake@xxxxxxxxx>]

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.

Johan