Hi all, Cutting to the chase: my problem involves using traceless, symmetric 3x3 tensor fields, and I represent these as 5-component vector finite element fields, where each vector-component is a degree of freedom of the tensor. My domain is a cylinder, and I'm trying to set up periodic boundary conditions on the ends, with one end rotated pi/2 with respect to the other. However, the `make_periodicity_constraints` function assumes an actual vector finite element field (in the sense of having the same dimension as space-dim) when applying the rotation matrix, and so I am unable to use periodic constraints in a straightforward way for my application.
Two questions on this: 1) how difficult would it be to accommodate this particular issue from outside the library? I've been looking through the code and seeing "this is a hack" and "we have to be careful" so I'm a bit wary to try to re-implement `make_periodicity_constraints` on my own. 2) Is it possible (in the future) to be able to apply arbitrary linear constraints between periodically-joined fields? Thanks for any help! - Lucas -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/cb8f7124-2ccc-4db9-9a60-02e5e2f7b75fn%40googlegroups.com.
