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?
--
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.
> >
>
signature.asc
Description: Digital signature
_______________________________________________ Mailing list: https://launchpad.net/~ffc Post to : ffc@lists.launchpad.net Unsubscribe : https://launchpad.net/~ffc More help : https://help.launchpad.net/ListHelp
