[petsc-users] jacobians for multicomponent problems

[Matthew Knepley]

On Wed, Jun 24, 2015 at 2:39 AM, domenico lahaye <domenico_lahaye at yahoo.com>
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 <knepley at gmail.com>
> *To:* Gideon Simpson <gideon.simpson at gmail.com>
> *Cc:* petsc-users <petsc-users at 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 at 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
>
>
>


-- 
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150624/7220ab3a/attachment.html>