[petsc-users] Setting up a PETSc section for field-split

[Matthew Knepley]

On Fri, Jan 24, 2014 at 5:10 PM, Luc Berger-Vergiat <lb2653 at columbia.edu>wrote:

> Hi all,
> I want to use PETSc as a solver for an FEM problem: modelization of shear
> bands using a 4 fields mixed formulation.
>
> 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.
>
> As of now I have issues when I try to assign a point and its associated
> degree of freedom to a field using: PetscSectionSetFieldDof.
> Is this the correct way to associate a dof/point to a field?
>

You have to set the number of fields before the chart. I am updating the
docs.

 Thanks,

    Matt


> I attached an example code and its makefile.
>
> --
> Best,
> Luc
>
>


-- 
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/20140125/2fb62fc7/attachment.html>