[aspect-devel] periodic domains
Wolfgang Bangerth
bangerth at math.tamu.edu
Wed Sep 11 11:25:10 PDT 2013
> I am somewhat confused about the different versions for
> make_periodicity_constraints (and I think we need to clean those up!).
> There is no version of make_periodicity_constraints that accepts a
>
> std::vector<std_cxx1x::tuple< typename
> parallel::distributed::Triangulation<dim>::cell_iterator, unsigned int
> /*face index*/,
> typename
> parallel::distributed::Triangulation<dim>::cell_iterator, unsigned int
> /*face index*/> >
>
> coming out of identify_periodic_face_pairs() right now. The closest one is
>
> (const FaceIterator &face_1, const typename identity< FaceIterator
>> ::type &face_2,::ConstraintMatrix &constraint_matrix, const
> ComponentMask &component_mask=ComponentMask(), const bool
> face_orientation=true, const bool face_flip=false, const bool
> face_rotation=false)
>
> but I don't know how to construct the bool flags (or if this is even
> possible in general, because identify_periodic_face_pairs() does not
> allow for flipping the direction for the periodic boundary to say make
> a moebius strip). This might be necessary for a half shell (not sure
> though).
Move this to the deal.II lists.
>
> So, I think what we need to return in the interface is something like
> std::vector of
> std::tuple of
> cell_iterator1, face_index1, cell_iterator2, face_index2, bool flip
>
> Or am I wrong?
Yes, I think that's the right interface for GeometryModel::Interface.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bangerth at math.tamu.edu
www: http://www.math.tamu.edu/~bangerth/
More information about the Aspect-devel
mailing list