Re: [Libmesh-users] Zero quadrature points in QGauss

[Vijay S. M. <vi...@gm...>]

Yes, I use TRI3 with first order Lagrange basis functions. The code is
part of a much bigger system. I am right now trying to create a test
case to see if I can replicate the issue. I will get back to you with
the code soon. Meanwhile, if you have any suggestions that could help
me, that would be fantastic.

Vijay

On Fri, May 1, 2009 at 9:06 AM, John Peterson <jwp...@gm...> wrote:
> On Fri, May 1, 2009 at 9:03 AM, Vijay S. Mahadevan <vi...@gm...> wrote:
>> Hi,
>>
>> This might sound ridiculous but I was wondering if there is a scenario
>> when the quadrature rule you created for say, 2 dimensions and SIXTH
>> order returns zero points ? Is there ever a need for initing
>> quadrature rules explicitly ?
>>
>> To give you a background, I have an ImplicitSystem and an
>> ExplicitSystem in my EquationSystems. I want to update the solution in
>> ExplicitSystem (solution_exp) based on the solution in the
>> ImplicitSystem (solution_imp). Hence, I loop over the local elements
>> in the mesh, perform user code and compute solution_exp. But at the
>> beginning of this loop, I attach the quadrature rule to the FEType of
>> variable in ImplicitSystem which returns zero points (N_Q_Points=0).
>> The order and dimension variables in the rule are correct but just
>> that the points have not been computed.
>>
>> I have tried to debug this code for a while but I cannot get past the
>> zero quadrature points. Any ideas or suggestions to get further on
>> this will be extremely helpful.
>>
>> I am terribly sorry if I was vague or confusing or if this is
>> downright absurd. I spent an entire day trying to get to the root of
>> this problem but got nowhere. Please feel free to let me know if you
>> need more details. Thanks a bunch in advance.
>
> Which element type, triangles?  We may need you to post your code or
> at least email it to one of us to take a closer look...
>
> --
> John
>