← Back to team overview

ffc team mailing list archive

Re: General / Special quadrature rules in FFC

 

On Mon, Jul 06, 2009 at 04:35:45PM +0200, Anders Logg wrote:
> On Mon, Jul 06, 2009 at 03:34:29PM +0200, Kristian Oelgaard wrote:
> > Quoting Robert Kirby <robert.c.kirby@xxxxxxxxx>:
> > 
> > > Such rules may best be put into FIAT, and then propogated into ffc.
> > 
> > That would be great, the infrastructure is in place on the FFC side, for
> > quadrature representation, it is just a matter of deciding on a naming
> > convention and then have FIAT accepting an additional string argument when
> > generating the rules.
> > 
> > I also know that Anders would like to move a lot of other things (like the
> > finite element classes) from the ffc/fem module to FIAT in order to shrink the
> > code base, but I'll let him comment on the details.
> 
> Almost. I would like to be able to use the FIAT FiniteElement class
> directly from FFC without wrapping it. For that to work, some (small)
> changes will be needed on the FIAT side (and some on the FFC side).
> 
> One thing that we currently handle in FFC which FIAT can't handle is
> the MixedElement class. This needs to tabulate basis functions for a
> list of elements with zeros inserted (taking the tensor product). Rob
> has not wanted to put this into FIAT before, but maybe he can be
> convinced... :-)

And while we're on the subject, should it be

  TH = P2 + P1 (what we have now)

or

  TH = P2 * P1

?

Someone mentioned it to me recently, I think it was Kristian.

-- 
Anders

Attachment: signature.asc
Description: Digital signature


Follow ups

References