[deal.II] extrude triangulation with n_slices = 1

2023-03-29 Thread Greg Wang
Hi there, We want to construct a 3D triangulation by extruding a 2D triangulation (one that potentially contains hanging nodes) and we only want one slice/layer of mesh on the extrusion direction. Looking around in the GridGenerator namespace led me to the extrude_triangulation function. It’

Re: [deal.II] writing data to a file changes program output although this file is never used

2023-03-29 Thread Simon Wiesheier
Of course, Approach B is the right way to go if I want to produce reliable results. After plenty of hours of debugging, I found a possible source for the differences: The member function called by my assembly routine has the following signature void fun(const Tensor<2,3 & F, const d

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

2023-03-29 Thread jose.a...@gmail.com
Hi dealii community, I think I found where the problem lies. In order to reproduce the assertion that gets thrown, i.e., An error occurred in line <1690> of file in function const dealii::FiniteElement& dealii::DoFAccessor::get_fe(unsigned int) const [with int structdim = 1; int dim =