Thanks Denis. I guess I still don't quite understand the intimate relationship between all the key data structures in deal.ii quite yet. I thought of this previously but I was hesitant because this seems memory expensive to me but then again I don't know how it is implemented so maybe it is not as expensive as I thought.
On Thursday, December 22, 2016 at 11:14:57 AM UTC-5, Denis Davydov wrote: > > > > On 22 Dec 2016, at 17:10, cmha...@gmail.com <javascript:> 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.
