[petsc-dev] [petsc-users] Recursive fieldsplit

Jed Brown jedbrown at mcs.anl.gov
Mon May 20 18:10:12 CDT 2013


Matthew Knepley <knepley at gmail.com> writes:

> Actually I now think I was wrong. Let me try to go through the use case,
> and you tell
> me where I am wrong:
>
>   1) We start with a system with block size B
>
>   2) We split it into k parts, so that \sum_k b_k = B, and b_k is the
> blocksize of each
>       component. Isn't this how FS for structured grids works?

Are you mixing up the size of the split with the size of the pieces?
That code path is only used when all the scalar splits are activated,
otherwise.



More information about the petsc-dev mailing list