Andrew, > I've adapted a 2d periodic boundary conditions algorithm. I've sorted the > faces of the element on the opposite sides based on the position of the > centroid of the elements face. The problem is now to match degrees of > freedom on opposite faces. I can determine if they are associated with the > displacements and to which local support point they correspond but this > does not give me a geometric measure to order the dof.
Have you taken a look at step-45 already? The algorithm used in the program itself may not be of much use in your case, but take a look at the "Possibilities for extensions" section at the end of the "Results" section: it has a more generic algorithm that should be what you are looking for. In particular, if you only want to periodicize some of your vector components, you would have to restrict the final look over all dofs_per_face in make_periodicity_constraints_recursively to those degrees of freedom for which fe.face_system_to_component_index(i) is a vector component that you are interested in. Does that help? Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
