On Wed, Jun 24, 2015 at 2:39 AM, domenico lahaye <[email protected]> wrote:
Would the same work for non-structured problems, for instance for power flow equations?
Yes, PetscSection knows about components for fields, and this carries through to the DMNetwork setup. It will, for instance, use block matrix formats if the block size is constant. Thanks, Matt
Thx. Domenico Lahaye.
------------------------------ *From:* Matthew Knepley <[email protected]> *To:* Gideon Simpson <[email protected]> *Cc:* petsc-users <[email protected]> *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 <[email protected]> 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
-- 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