ffc team mailing list archive
-
ffc team
-
Mailing list archive
-
Message #02716
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