Alberto,
be aware of the fact that, when you extract the boundary mesh, the orientation
of the cells may be different w.r.t. to the corresponding face on the
volumetric mesh. In other words, using
FEValues codim_one(…);
FEFaceValuescodim_zero(…);
codim_one.reinit(codim_one_cell);
co
I was trying to run the tutorial 36 but it shows up with following errors:
*[0]PETSC ERROR: - Error Message
--*
[0]PETSC ERROR: Error in external library
[0]PETSC ERROR: Error in LAPACK subroutine steqr: info=0
[0]
I second Wolfgang comment on the fact that Q1Q1 is not difficult to
implement. You can also scale it to arbitrary Qn-Qn elements if you are
interested in higher order.
We have implemented such an approach in our code based on dealii (
https://github.com/lethe-cfd/lethe)
Q1Q1 is already very diffu
Thank you very much, Wolfgang, crystal clear as always.
It is still unclear to me, though, how to reinit() the FEFaceValues
for the cell making use of the
map returned by the extract_boundary_mesh() function. I learned to reinit
FEFaceValues as:
fe_face_values.reinit (cell, face_number);
Is ther