On 11/11/22 07:06, Alex Quinlan wrote:

Thanks very much for your reply.  I'm happy to hear that you think my problems can be resolved without too much difficulty. I'd like to try writing a patch for the mixed 3D meshes.  I have done some development work before, but I have never sent it out of our group; your advice on writing the steps to write a patch would be most welcome!

Alex:
At the most basic level, the steps to contribute patches is outlined here:
  https://dealii.org/participate.html
Let me know if you have questions and I'll be happy to walk you through whatever is missing in these documents.


I will re-review the relevant code and see if I can put together a general plan on how to approach this.

And here too: Feel free to discuss it with me or here on this list.

For the specific code where you ended up with an exception: This is the place where a finite element describes how the degrees of freedom on one of its faces should interact with the degrees of freedom of the finite element on a neighboring cell. For example, if you have a P1 and a P2 element on neighboring cells (let's assume triangles for the moment), then the mid-edge degree of freedom of the P2 element needs to have a specific value to ensure that the resulting finite element field is continuous. In your situation, the function has simply not been taught what these constraints should be if you have a FE_SimplexP(1) element on a tetrahedron's triangular face meet the degrees of freedom of a FE_WedgeP(1) or FE_PyramidP(1) on that common triangular face. The answer is that the degrees of freedom should all be equal to each other, as in some of the other cases handled by the function.

Best
 Wolfgang

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           www: http://www.math.colostate.edu/~bangerth/


--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/a3a8ce3c-30bc-09e6-0541-248d99c5e691%40colostate.edu.

Reply via email to