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
