ffc team mailing list archive

[Mehdi Nikbakht <m.nikbakht@xxxxxxxxxx>]

On Tue, 2010-02-02 at 10:30 +0100, Anders Logg wrote:
> I tried looking at this but I'm unsure how it should be
> handled. Should a cell_integral class be generated or should a
> surface_integral class be generated?
> 

We handle terms related to surface integral inside a class derived from
ufc::cell_integral. I have started working on updating ffcpum module
which is built against standard ffc. 

Mehdi


> --
> Anders
> 
> 
> On Tue, Jan 26, 2010 at 10:21:51AM +0100, Medhi Nikbakht wrote:
> > Hi Anders,
> > This is the simplest form that I trying to compile.
> >
> > elem_cont = FiniteElement("CG", triangle, 1)
> > elem_discont = ElementRestriction(elem_cont, dc) # or ec[dc]
> >
> > element = elem_cont + elem_discont
> >
> > (vc, vd) = TestFunctions(element)
> > (uc, ud) = TrialFunctions(element)
> >
> > v = vc + vd
> > u = uc + ud
> >
> >
> > k  = Constant(triangle)
> > f  = Coefficient(elem_cont)
> >
> > a = inner(grad(v), grad(u))*dx + k*vd*ud*dc
> > L = v*f*dx
> >
> > What about adding this form as a test for surface integral?
> >
> > Mehdi
> >
> > On Tue, 2010-01-26 at 08:37 +0000, Anders Logg wrote:
> > > The unit test should be a small piece of the form that you are
> > > compiling and it should contain the surface integral.
> > >
> > > Put in in a UFL file.
> > >
> >