<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Oct 18, 2016 at 7:38 AM, Julian Andrej <span dir="ltr"><<a href="mailto:juan@tf.uni-kiel.de" target="_blank">juan@tf.uni-kiel.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
i have general question about PetscFE. When i want to assemble certain<br>
parts of physics separately, how can i do that? I basically want to<br>
assemble matrices/vectors from the weak forms on the same DM (and<br>
avoid copying the DM) and use them afterwards. Is there a convenient<br>
way for doing that?<br>
<br>
The "workflow" i'm approaching is something like:<br>
<br>
- Setup the DM<br>
- Setup discretization (spaces and quadrature) for each weak form i<br>
want to compute<br>
- Compute just the weak form i want right now for a specific<br>
discretization and field.<br>
<br>
The reason is i need certain parts of the "complete" Jacobian for<br>
computations of eigenproblems and like to avoid computing those more<br>
often than needed.<br></blockquote><div><br></div><div>The way I envision this working is to use DMCreateSubDM(). It should extract</div><div>everything correctly for the subset of fields you select. However, I have not</div><div>extensively tested, so if something is wrong let me know.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Regards<br>
<span class="HOEnZb"><font color="#888888">Julian<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>