[deal.II] problem with CylindricalManifold<3> after removing cells from a triangulation

2021-04-23 Thread Simon
Dear all, I created a mesh using GridGenerator::hyper_cube_with_cylindrical_hole and removed after that half of the total number of the cells in order to make use of symmetry. Then I set a SphericalManifold<2> (2D) respectively CylindricalManifold<3> (3D) in order to get a true circle (2D) re

Re: [deal.II] interpolate FE_Nelelec

2021-04-23 Thread Konrad Simon
Hi John, Maybe I can give some hints of how I would approach this problem (these are just some quick thoughts): Write this problem as a vector Laplace problem curl(nu*curl(A)) - grad(div(A)) = J which you then have to write as a system of PDEs. Note, this can only be the same as curl(nu*curl(

[deal.II] Re: integrating over boundaries

2021-04-23 Thread Sylvain Mathonnière
By the way, I am aware than I forgot to take the square root of heat_flux_cell_norm. I am also aware that using *input_data.solution_gradients[face_index][d]* is weird in this context since the first indices should refer to a vertex indices of the cell and not the face index and it only works

[deal.II] integrating over boundaries

2021-04-23 Thread Sylvain Mathonnière
Hello everyone ! I have a very basic question but somehow I did not manage to find the answer in tutorials, neither looking at the documentation or in the this user group. Or more probably I did not understand it when I saw it since I am fairly new to deal.II and to C++. *My main goal :* I am