[deal.II] Re: neighboring faces on the boundary

Erik, you might want to have a look at https://groups.google.com/d/msg/dealii/wG6aJdXHweA/IcsUjxj5dFsJ Basically, I would loop through all the neighbors recursively only considering cells that have a vertex with the face you are looking for in common. You can access the DoFs on faces just as fo

Re: [deal.II] Re: How to shift a ConstraintMatrix that initialized with IndexSet?

Ce Qin, This finally made its way into the developer version via PR #3722 . The approach we finally came up with is to also shift the local indices when shifting the ConstraintMatrix objects.[1] Merging two ConstraintMatrix objects also merges the lo

Re: [deal.II] Re: neighboring faces on the boundary

Also this page might be useful https://www.dealii.org/8.4.1/doxygen/deal.II/structGeometryInfo.html It shows how the faces within a cell are topologically related to one another. Best praveen -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https

[deal.II] Re: Synchronization issue for parallel computation of material forces

Maybe this helps: How do you compute "reaction forces" within deal.II? Did you ever compute reaction forces and write them as output within a parallel framework? Best, Seyed Ali -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.g

[deal.II] Re: Synchronization issue for parallel computation of material forces

Seyed, your results look like the ghost values are wrong. My first guess would be that you are still using a completely distributed vector (one initialized with locally_owned_dofs only) instead of a ghosted one (initialized with locally_owned_dofs and locally_relevant_dofs). Can you confirm tha

[deal.II] Re: neighboring faces on the boundary

Dear Erik, We actually already have a function that can do 1/2 the job for you, namely GridTools:: compute_active_cell_halo_layer . In this instance your predicate would be whether the inpu

Re: [deal.II] Synchronization issue for parallel computation of material forces

Dear Seyed, Daniel’s been nice enough to furnish you with more information, but at this point one of your deleted comments is quite telling to me: > Any ideas what could cause such a behavior? I am really not an expert yet to > dive deeply inside the deal.II parallel structure. As we explain

[deal.II] Re: Synchronization issue for parallel computation of material forces

@Daniel: Thanks. I just tried Jean-Paul's approach again, but with a new procedure. It worked. If I could, I would mark both of your comments as the best answer. Thank you everyone, again :) I will not forget your kind assistance and hopefully help you also soon, if possible.

[deal.II] Re: neighboring faces on the boundary

> > We actually already have a function that can do 1/2 the job for you, > namely GridTools:: compute_active_cell_halo_layer > . > > In this instance your predicate would be whether the inpu

[deal.II] Re: neighboring faces on the boundary

> > Since Erik explicitly asks for a way to "get the degrees of freedom living > on the set of faces", I would still just modify the referenced code to > return face_iterators and access the corresponding DoFs using them. > You're 100% correct, I misread the question (well, rather I focussed on