From: Roy S. <roy...@ic...> - 2017-12-14 15:04:07
|
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 |