Would the same work for non-structured problems, for instance for power flow equations?

Thx. Domenico Lahaye.

From: Matthew Knepley <knepley@gmail.com>
To: Gideon Simpson <gideon.simpson@gmail.com>
Cc: petsc-users <petsc-users@mcs.anl.gov>
Sent: Wednesday, June 24, 2015 5:05 AM
Subject: Re: [petsc-users] jacobians for multicomponent problems

On Tue, Jun 23, 2015 at 7:52 PM, Gideon Simpson <gideon.simpson@gmail.com> wrote:

Suppose I have a problem for which the field variable has multiple components, in this case real and imaginary parts, and I have created a DM with two degrees of freedom. When assembling a Jacobian associated with a nonlinear problem on this data, is there any data management accounting for the different components that makes this “easier?”

DMDA knows about components so that MatStencil has a .c slot for indices. Is that what you mean?

Thanks,

Matt

-gideon

What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
-- Norbert Wiener