Work at SourceForge, help us to make it a better place! We have an immediate need for a Support Technician in our San Francisco or Denver office.

## [Libmesh-devel] various types of libMesh::Elem volume

 [Libmesh-devel] various types of libMesh::Elem volume From: Vasilis Vavourakis - 2011-05-27 08:02:52 Attachments: Message as HTML ```hi all :) it happened for me to need to evaluate the "volume" of Face elements in my code through these functions: Real Tri3::volume () const Real Quad4::volume () const ...however, i noticed that the same function is not implemented for Tri6, Quad8, Quad9 and so on... due to the fact that these elements might have curved shaped, therefore none can use an analytic expression to evaluate their "volume". my question: do you think that it would fine in the next version of libMesh to calculate the "volume" through standard Gauss integration: Vol = sum_{i=0}^{nqp}{ JxW[i] } i think that this might work for a decent number of quadrature points. obviously the same rule can be applied for 3D quadratic elements and for the high-order 1D elements too. cheers, Vas ```

 [Libmesh-devel] various types of libMesh::Elem volume From: Vasilis Vavourakis - 2011-05-27 08:02:52 Attachments: Message as HTML ```hi all :) it happened for me to need to evaluate the "volume" of Face elements in my code through these functions: Real Tri3::volume () const Real Quad4::volume () const ...however, i noticed that the same function is not implemented for Tri6, Quad8, Quad9 and so on... due to the fact that these elements might have curved shaped, therefore none can use an analytic expression to evaluate their "volume". my question: do you think that it would fine in the next version of libMesh to calculate the "volume" through standard Gauss integration: Vol = sum_{i=0}^{nqp}{ JxW[i] } i think that this might work for a decent number of quadrature points. obviously the same rule can be applied for 3D quadratic elements and for the high-order 1D elements too. cheers, Vas ```
 Re: [Libmesh-devel] [Libmesh-users] various types of libMesh::Elem volume From: John Peterson - 2011-05-27 13:33:14 ```On Fri, May 27, 2011 at 2:02 AM, Vasilis Vavourakis wrote: > hi all :) > > my question: do you think that it would fine in the next version of libMesh > to calculate the "volume" through standard Gauss integration: > > Vol = sum_{i=0}^{nqp}{ JxW[i] } Yes, this technique is already implemented exactly as you suggest in the base class' Elem::volume(). -- John ```