On Thu, 14 Dec 2017, Viviana Palacio Betancur wrote:
Thank you! I was able to fix my problem. I was creating the quadrature rule before looping over the elements. Thus, "ERROR: Unsupported type: 4" appeared when the TRI6 elements appeared. I fixed this by creating the quadrature rule and building the FEBase within the element loop.
I wouldn't recommend that in the long run, for performance reasons. The best thing to do currently is to create an FE for each dimension, or to use FEMContext which does that under the hood. For small problems or prototype code you should be fine, though. --- Roy ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
