Hi Wolfgang > > 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.
Thanks a lot. This is exactly how I ended up solving the problem and all works fine Andrew > > 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
