[petsc-users] question about small matrices

[Jed Brown]

"Povolotskyi, Mykhailo via petsc-users" <petsc-users at mcs.anl.gov> writes:

> Hi Matthew,
>
> is it possible to do in principle what I would like to do?

SNES isn't meant to solve tiny independent systems.  (It's just high
overhead for that purpose.)  You can solve many such instances together
by creating a residual function that evaluates them all.  The linear
algebra will be more efficient with that granularity, though all
sub-problems will take the same number of Newton iterations.