[deal.II] Re: DoFTools::make_periodicity_constraints error when working with FEColleciton

Hi Peter, The assert is no longer triggered in serial (independently of the triangulation is globally refined or not) but it is still triggered in parallel (Line 2336 ). I also tested in release mode. In

Re: [deal.II] shape function gradient at arbitrary points in the element

Just to make sure that I chose a reasonable approach: //create FEValues object whenever I need it, that is, 'number_of_cells*number_of_qps_per_cell' times FEValues<2> fe_values(dof_handler.get_fe(), Quadrature<2>(my_point_in_ref_coords), update_gradients);

Re: [deal.II] shape function gradient at arbitrary points in the element

On 7/7/22 02:34, Simon Wiesheier wrote: Just to make sure that I chose a reasonable approach: //create FEValues object whenever I need it, that is, 'number_of_cells*number_of_qps_per_cell' times FEValues<2> fe_values(dof_handler.get_fe(),   

Re: [deal.II] error while running code dependent on step-44

On 7/7/22 09:46, ME20D503 NEWTON wrote: *** Caution: EXTERNAL Sender *** Dear all, i am getting error while running a code based on step-44. i am not able to understand the error output, can anyone please help me out ? ME20D503 NEWTON: As I think others have already pointed out, many error

[deal.II] error while running code dependent on step-44

Dear all, i am getting error while running a code based on step-44. i am not able to understand the error output, can anyone please help me out ? thanks in advance. Exception on processing: An error occurred

Re: [deal.II] 3D Periodic BCs

Jack, step-45 https://www.dealii.org/current/doxygen/deal.II/step_45.html provides an example for dealing with periodic boundary conditions. Even though it's a 2D example, the interface is the same in 3D. Assuming that the arguments for your set_periodic_boundary_condition function are similar

[deal.II] Re: Using AMR with MatrixFree

The problem is that ECL is not working with AMR, yet. A technical issue is that one cannot that easily loop over 2*dim faces and compute the fluxes for each but also has to go over the subfaces when the neighbor is refined. One might be able to hide this in FEFaceEvaluation and load the values

[deal.II] 3D Periodic BCs

I am trying to set periodic boundary conditions for a 3D problem like *velocity->set_periodic_boundary_condition(0, 1, 2); velocity->set_periodic_boundary_condition(2, 3, 1); pressure->set_periodic_boundary_condition(0, 1, 2); pressure->set_periodic_boundary_condition(2, 3, 1);* The