[petsc-users] Can I use TS routines for operator split formulation

[Jed Brown]

Shriram Srinivasan <shriram at ualberta.ca> writes:
> Au = (k grad(u)_x,  is a discretization of the flux in the x direction, 
> while B is for flux in the y direction. These matrices are assembled by 
> solving local elliptic problems in each cell following a  multiscale 
> finite volume method.
>   To begin, the operator split is done as follows:
>
> u_t + Au = f ( considering only x dir fluxes) as (u* - u_n)/tau   + A u* 
> = f and solve for u*
>
> then u_t + B u = f (considering only y dir fluxes) as (u - u*)/tau + B u 
> = f and solve for u at current level.

This is ADI "alternating direction implicit" and tends not to be very
accurate, but more importantly, is a pretty terrible parallel algorithm
that isn't used much any more.  Unless you are implementing this for
nostalgia purposes, I recommend using multigrid (algebraic or geometric)
to solve the coupled problem.  This will be fast and eliminate the
splitting error.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140609/f42e1cea/attachment.pgp>