<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Jan 24, 2014 at 5:10 PM, Luc Berger-Vergiat <span dir="ltr"><<a href="mailto:lb2653@columbia.edu" target="_blank">lb2653@columbia.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi all,<br>
I want to use PETSc as a solver for an FEM problem: modelization of shear bands using a 4 fields mixed formulation.<br>
<br>
So far I am just trying to set up a two fields in a Petsc section, assign dof too each field and then set up the section to pass it to a DM. I am taking this approach because in general I want for fieldsplit to work purely on the algebraic level without knowledge of boundary conditions or geometry.<br>

<br>
As of now I have issues when I try to assign a point and its associated degree of freedom to a field using: PetscSectionSetFieldDof.<br>
Is this the correct way to associate a dof/point to a field?<br></blockquote><div><br></div><div>You have to set the number of fields before the chart. I am updating the docs.</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">
I attached an example code and its makefile.<span class="HOEnZb"><font color="#888888"><br>
<br>
-- <br>
Best,<br>
Luc<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>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>