[petsc-users] Vector form of TS solvers

[Chris Bradley]

Dear PETSc,
I was wondering if you had any advice for using PETSc TS solvers to
solve a vector of ODEs. PETSc TS solvers solve DAEs of the form

f(t,u,u^dot,a) = g(u,t) = 0 (for my case)

where f is a user provided function, u is a vector of state variables,
u^dot a vector of the derivatives of state variables and a is a vector
of parameters (passed through the context).

I need to solve a large number of these DAEs (one for every node in a
bigger finite element problem). Is there a vector form of the TS solvers
that can solve

F(t,U,U^dot) = G(U,t) = 0

where F is now a vector, U is a matrix of state variables and U^dot is a
matrix of the derivatives of U. Each row of the system corresponds to
the standard PETSc TS problem above and is independent of every other row.

If I understand the TS manual correctly then if I wanted to do this in
PETSc I would have to create a TS object and state and derivative
vectors objects for each row? All these PETSc objects for each row would
seem like a big overhead if there was a very large number of rows
(potentially tens of millions for my problem).

Thanks for your help
Chris
-- 
Chris Bradley
Auckland Bioengineering Institute,
University of Auckland,
Private Bag 92019,
New Zealand.

Tel: +64 9 373 7599 xtn 89924
Fax: +64 9 367 7157