Re: [Libmesh-users] Question about hybrid mesh

[Viviana P. B. <vpa...@uc...>]

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 <roy...@ic...>
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
>