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?” -gideon
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
Yes, is there and example for this? -gideon
On Jun 23, 2015, at 11:05 PM, Matthew Knepley <[email protected]> wrote:
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
On Tue, Jun 23, 2015 at 10:30 PM, Gideon Simpson <[email protected]> wrote:
Yes, is there and example for this?
src/ts/examples/tutorials/advection-diffusion-reaction/ex5.c Thanks, Matt
-gideon
On Jun 23, 2015, at 11:05 PM, Matthew Knepley <[email protected]> wrote:
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
Would the same work for non-structured problems, for instance for power flow equations? 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
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
participants (3)
-
domenico lahaye -
Gideon Simpson -
Matthew Knepley