Re: [Libmesh-users] question on Trapezoidal integration rule

[Michael P. <mpo...@pu...>]

On 11/07/2018 02:59 PM, Stogner, Roy H wrote:
> On Wed, 7 Nov 2018, Michael Povolotskyi wrote:
>
>> if I use the trapezoidal rule:
>>
>> QTrap qrule(dim)
>>
>> Can I assume the quadrature points are ordered in the same way at the element
>> nodes?
> Doesn't look that way, in general, at first glance.  We build the quad
> and hex QTrap rules via tensor products, but that's not how we order
> nodes.
> ---
> Roy
Thank you,
then how to solve the following problem.

I have an element that has n_nodes;
It has the same number of quadrature points that are at the nodal 
positions (trapezoidal rule is assumed).
How to find the node number that corresponds to a quadrature point?
Should I compare coordinates?

Thank you,
Michael.