Hi Roy, 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. Best, Viviana. On Thu, Dec 14, 2017 at 9:03 AM, Roy Stogner <royst...@ices.utexas.edu> wrote: > > On Wed, 13 Dec 2017, Viviana Palacio Betancur wrote: > > I am trying to assemble the system's matrix for a mesh composed of TET10 >> (bulk) bounded by TRI6 (surface) elements but I get an error related to an >> unsupported type of elements. The problem solves for a single variable. >> > > This should be fine. We do something similar in fem_system_ex4. > > When I try to assemble the matrix, I get the error: *"ERROR: >> Unsupported type: 4."* >> > > You'll understand that error better if you look at the line number or > the full stack trace to get context. > > I don't understand if I should create separate quadrature rules >> because there are two types of elements. >> > > You'll need separate quadrature rules incidentally, because you'll > need separate FE objects (note that they're templated around element > dimension) and each FE object will need its own quadrature rule. > > But then the elements are attached to the same variable, so how >> would I be able to distinguish between them? >> > > You can test Elem::dim() for each. > --- > 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
