Re: [Libmesh-devel] New Quadrature Rule

2010-12-02 Thread John Peterson
On Thu, Dec 2, 2010 at 9:20 AM, Boyce Griffith wrote: > > > On 11/16/10 2:24 PM, Derek Gaston wrote: >> >> On Nov 16, 2010, at 10:27 AM, Boyce Griffith wrote: >> >>> Is adding qrule->need_reinit() the only change that you are making to >>> QBase?  I am probably missing something, but it seems like

Re: [Libmesh-devel] New Quadrature Rule

2010-12-02 Thread Boyce Griffith
On 11/16/10 2:24 PM, Derek Gaston wrote: > On Nov 16, 2010, at 10:27 AM, Boyce Griffith wrote: > >> Is adding qrule->need_reinit() the only change that you are making to QBase? >> I am probably missing something, but it seems like allowing the quadrature >> rule to vary from element to element

Re: [Libmesh-devel] New Quadrature Rule

2010-11-16 Thread Derek Gaston
On Nov 16, 2010, at 10:27 AM, Boyce Griffith wrote: > Is adding qrule->need_reinit() the only change that you are making to QBase? > I am probably missing something, but it seems like allowing the quadrature > rule to vary from element to element would also require something like > changing th

Re: [Libmesh-devel] New Quadrature Rule

2010-11-16 Thread Boyce Griffith
On 11/15/10 3:40 PM, Derek Gaston wrote: > On Nov 15, 2010, at 3:38 PM, John Peterson wrote: > >> Are we married to the "has_constant_positions()" name? Maybe >> something like qrule->need_reinit() in the vein of >> FE::shapes_need_reinit(). Just thinking there might potentially be >> other rea

Re: [Libmesh-devel] New Quadrature Rule

2010-11-15 Thread Roy Stogner
On Mon, 15 Nov 2010, John Peterson wrote: > On Mon, Nov 15, 2010 at 2:49 PM, Roy Stogner wrote: >> >> >> On Mon, 15 Nov 2010, Derek Gaston wrote: >> >>> On Nov 15, 2010, at 1:21 PM, Roy Stogner wrote: >>> Does the number of quadrature points change, or just their locations? If it's ju

Re: [Libmesh-devel] New Quadrature Rule

2010-11-15 Thread Derek Gaston
On Nov 15, 2010, at 3:38 PM, John Peterson wrote: > Are we married to the "has_constant_positions()" name? Maybe > something like qrule->need_reinit() in the vein of > FE::shapes_need_reinit(). Just thinking there might potentially be > other reasons quadrature rules need to be recomputed for ea

Re: [Libmesh-devel] New Quadrature Rule

2010-11-15 Thread John Peterson
On Mon, Nov 15, 2010 at 2:49 PM, Roy Stogner wrote: > > > On Mon, 15 Nov 2010, Derek Gaston wrote: > >> On Nov 15, 2010, at 1:21 PM, Roy Stogner wrote: >> >>> Does the number of quadrature points change, or just their locations? >>> If it's just their locations then I think you've probably hit on

Re: [Libmesh-devel] New Quadrature Rule

2010-11-15 Thread Roy Stogner
On Mon, 15 Nov 2010, Derek Gaston wrote: > On Nov 15, 2010, at 1:21 PM, Roy Stogner wrote: > >> Does the number of quadrature points change, or just their locations? >> If it's just their locations then I think you've probably hit on the >> best solution for now. > > Both the number and position

Re: [Libmesh-devel] New Quadrature Rule

2010-11-15 Thread Derek Gaston
On Nov 15, 2010, at 1:21 PM, Roy Stogner wrote: > Does the number of quadrature points change, or just their locations? > If it's just their locations then I think you've probably hit on the > best solution for now. Both the number and positions. I'll go ahead and add something like this for now

Re: [Libmesh-devel] New Quadrature Rule

2010-11-15 Thread Roy Stogner
On Mon, 15 Nov 2010, Derek Gaston wrote: > Would anyone be opposed to me adding a > qrule->has_constant_positions() (that would return true for all of > our current rules and false for the new one I'm making)... or > something like it (I'm open to suggestions!) Does the number of quadrature poi

[Libmesh-devel] New Quadrature Rule

2010-11-15 Thread Derek Gaston
I'm trying to make a new quadrature rule... that is different in that the points it integrates at are different on every element in the mesh. This would seem to be an easy task but I am wondering about FE::reinit(). Reading through it, it seems as if it is missing the case where the quadrature