dolfin team mailing list archive

[Johan Hake <johan.hake@xxxxxxxxx>]

Hello!

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. 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.