[deal.II] tutorial step-36 failed with PETSc error number 76 while calling a SLEPc function

2019-09-10 Thread Pragalv Karki
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

Re: [deal.II] Stabilization for Q1Q0 Stokes by Bochev [2006]

2019-09-10 Thread Bruno Blais
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

Re: [deal.II] Re: Extracting volume data to the boundary

2019-09-10 Thread Alberto Salvadori
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