> On 22 Dec 2016, at 17:10, cmhame...@gmail.com wrote: > > The thing is though that I don't actually need to integrate at these points. > The method still employ a 4node quad element in 2d or an 8 node hex element > in 3d but only uses the element centroid to calculate the deformation > gradient, no quadrature is actually done at that point so I feel like I need > to interpolate otherwise the way deal makes the FESystem objects I'll now > have quad points at all centroids that should be integrated, or am I > misunderstanding something here?
Nothing stops you from using different quadrature rules with different FEValues objects simultaneously. One pair with centroid will be auxiliary and not participate in integration directly (so you won’t use the weight associated to your single point quadrature rule anywhere). Whereas for integration on your linear elements you would most likely be using QGauss<dim> quarature(degree+1); Regards, Denis. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
