Re: [deal.II] Calculate cell center distance from a boundary

2022-02-17 Thread vachan potluri
> > Hello, > Here is the PR https://github.com/dealii/dealii/pull/13394 that adds the > new wrappers for ArborX > Best, > Bruno Thank you very much! Didn't expect it to come so fast :) ! -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

Re: [deal.II] Re: MPI, synchronize processes

2022-02-17 Thread Wolfgang Bangerth
On 2/17/22 09:22, Uclus Heis wrote: I still had problems as I first copy the array and then I store it in a matrix for diffeerent frequencies.The result I got was differet whene using few process compared to using 1 single process. I added the following code and now works, is it right? It

[deal.II] Re: MPI, synchronize processes

2022-02-17 Thread Uclus Heis
Dear Bruno, I still had problems as I first copy the array and then I store it in a matrix for diffeerent frequencies.The result I got was differet whene using few process compared to using 1 single process. I added the following code and now works, is it right? { static

Re: [deal.II] Impose values inside a material/ not on the boundary

2022-02-17 Thread Wolfgang Bangerth
On 2/15/22 07:32, 'Markus Mehnert' via deal.II User Group wrote: Thank you for your response. I identified the correct dofs by trial and error which worked (9 dofs for a scalar quantity with quadratic polyorder). However, when I use /"fe_face.system_to_base_index().first.first" /to

Re: [deal.II] Triangulation from lower-dimensional object

2022-02-17 Thread Wolfgang Bangerth
On 2/17/22 02:48, Konrad Simon wrote: Can I somehow extract a Triangulation form a Triangulation by providing, for example, the boundary id on the Triangulation level? Yes! Take a look at GridGenerator::extract_boundary_mesh(). Best W. --

Re: [deal.II] Coupled poroelastic-elastic problem

2022-02-17 Thread Wolfgang Bangerth
On 2/17/22 01:50, Мария Бронзова wrote: I took the step-46 as an example for matching two physical domains, everything works fine, apart from the interface boundary condition. In my case I need to strongly enforce the continuity of the solid displacement vector. In other words I need to set

Re: [deal.II] Calculate cell center distance from a boundary

2022-02-17 Thread Bruno Turcksin
Hello, Here is the PR https://github.com/dealii/dealii/pull/13394 that adds the new wrappers for ArborX Best, Bruno On Saturday, February 12, 2022 at 11:58:35 PM UTC-5 vachanpo...@gmail.com wrote: > Drs. Bruno and Wolfgang, > *Actually what Wolfgang is suggesting is very doable. We have

[deal.II] Re: MPI, synchronize processes

2022-02-17 Thread Bruno Turcksin
Hello, You don't need to synchronize anything the function has an MPI_Barrier. This is the reason for the comment. If one processor call the function but the others don't, the code will hang. Best, Bruno On Thursday, February 17, 2022 at 4:01:32 AM UTC-5 uclu...@gmail.com wrote: > Dear

Re: [deal.II] Way to use the erf function with MuParser?

2022-02-17 Thread blais...@gmail.com
Awesome, I'll make a PR for it. On Wednesday, February 16, 2022 at 11:33:58 p.m. UTC-5 Wolfgang Bangerth wrote: > On 2/16/22 20:16, blais...@gmail.com wrote: > > > > It seems MuParser does not, by default, allow for the use of the error > > function (erf). We are solving a Stefan problem

[deal.II] Triangulation from lower-dimensional object

2022-02-17 Thread Konrad Simon
Dear all, I was recently looking for a solution to the following problem: I have a triangulation, e.g., of a unit cube. I need to solve a lower-dimensional (FEM-like) problem on all faces with BCs on edges and for each face I need to solve a problem on all edges with BCs at vertices. Can I

[deal.II] MPI, synchronize processes

2022-02-17 Thread Uclus Heis
Dear all, I am copying a distributed vector using the following function: Vector & Vector < Number >::operator=(const PETScWrappers::VectorBase

[deal.II] Coupled poroelastic-elastic problem

2022-02-17 Thread Мария Бронзова
Dear All, I have been working with acoustic porous media and would like to extend my script to deal with two-layered acoustic absorbers, in this case I am interested in adding an elastic layer on top of my porous media. I took the step-46 as an example for matching two physical domains,