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
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
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