← Back to team overview

ffc team mailing list archive

Re: Tensors

 

On Wed, Aug 31, 2005 at 09:20:03PM +0200, Johan Jansson wrote:
> On Wed, Aug 31, 2005 at 11:08:42AM -0500, Anders Logg wrote:
> 
> ...
> 
> > Karin Kraft and Johan Jansson have already implemented a prototype
> > version of boundary integrals in DOLFIN as far as I understand, so I
> > only need to implement the remaining portions in FFC. I will probably
> > also need to request some new features from FIAT so I can tabulate
> > basis functions on faces/edges of the reference cells.
> > 
> 
> Yes, we have a DOLFIN implementation in 2D with a manually generated
> FFC-style header file for a reference case. I can send you the
> relevant files with some explanatory comments if/when you start
> working on this. Or if someone else wants to work on it perhaps, maybe
> it's a good idea to try to have several developers working in every
> component so that a single person doesn't get swamped?
> 
> I think Rob has already implemented tabulation over boundary entities,
> check out this message:
> 
> http://www.fenics.org/pipermail/fiat-dev/2005-August/000058.html

Great, I must have missed that.

Everything is prepared for boundary integrals in FFC. You can already
now write *ds, boundary terms are separated from interior terms, and
the DOLFIN format generates code for assembly of both terms.

The missing part is to call FIAT from FFC to tabulate on boundaries
and do the integration properly. I figure it will take me about a week
(from the starting point, which is not today) to figure out all the
details and get it running.

You can post the files on the list and then I (or anyone else) know
where to look.

/Anders



References