#5 new
Jed Brown

Implement the vector-valued version of dFS

Reported by Jed Brown | November 21st, 2008 @ 08:38 PM

Starting with a scalar-valued dFS, this should

  • Allow an arbitrary number of boundary conditions
  • Put boundary faces into the global space, determine normal vectors at nodes.

Comments and changes to this ticket

  • Jed Brown

    Jed Brown April 17th, 2009 @ 04:55 PM

    An update on this very old issue. See ticket #10 for implementation of boundary conditions. The only things specific to vector-valued dFS is handling of mixed conditions such as slip. For this, the coordinate system of the global vector should be rotated so that the normal component appears as an explicit degree of freedom. Assembly will drop this value, the residual in this component will always be 0 and the rows and columns of the matrix corresponding to this node will always be the identity.

  • Jed Brown

    Jed Brown April 17th, 2009 @ 08:30 PM

    (from [69efc49fcd9c59a7eff1850a61b24f4b7eb7247a]) Work on handling boundary conditions. [#5] [#10]

    • Provide Global, Dirichlet, and Closure vectors. Global and Dirichlet vectors are VecGhost so they have a local form. Scatters are provided from Closure to the others.

    • Use MAIJ format for element assembly matrices (allows vector-valued spaces)

    • Keep track of boundaries by using ITAPS/MOAB tag conventions instead of private data structure (this is useful since the user needs to see the set structure if they need to relate to geometric entities).

    • Add dFSHomogeneousMode to handle which boundary values should be mapped into expanded space. http://github.com/jedbrown/dohp/...

Please Sign in or create a free account to add a new ticket.

With your very own profile, you can contribute to projects, track your activity, watch tickets, receive and update tickets through your email and much more.

New-ticket Create new ticket

Create your profile

Help contribute to this project by taking a few moments to create your personal profile. Create your profile »

An implementation of the ``dual order hp'' version of the finite element method. This project targets parallel domain-decomposition methods for strongly coupled nonlinear problems with PDE constraints.

People watching this ticket

Referenced by

Pages