ffc team mailing list archive

[Kristian Oelgaard <k.b.oelgaard@xxxxxxxxxx>]

Quoting Jake Ostien <tostien@xxxxxxxxx>:

> And for anyone paying attention, here are the missing files.
> 
> 
> 
> Jake Ostien wrote:
> > Hi,
> >
> > Something doesn't seem right when I try and define multiple integrals 
> > over subdomains in a form file.
> >
> > See the attached test.form and test.txt.  Basically I copied the element 
> > tensor section of the generated header file from two different 
> > compilations of test.form, one where
> >
> > a = v*u*dx  (for a reference)
> >
> > and one with
> >
> > dx0 = Integral("cell", 0)
> > dx1 = Integral("cell", 1)
> >
> > a = v*u*dx0 + 10*v*u*dx1
> >
> > Now I would expect that the cell_integral_0 from the second compilation 
> > to be identical to the reference compilation, and that cell_integral_1 
> > would be just scaled by 10, but that is not what is happening.  It 
> > appears that in the second compilation the two subdomain integrals are 
> > just added together (see test.txt).
> >
> > I didn't get too much out of running ffc with the -d1 flag, but I really 
> > don't know what I am looking for.  Any help on where to look?

I'm not quite sure about what is going on, or what should be going on. It does,
however, look a bit strange that cell_integral_0 also contains the dx1 part of
the form.

The values look OK to me so it looks like a bug in determining which monomial
terms go to which integral? I don't really know where this takes place.

Kristian 

> > Thanks,
> > Jake
> > _______________________________________________
> > FFC-dev mailing list
> > FFC-dev@xxxxxxxxxx
> > http://www.fenics.org/mailman/listinfo/ffc-dev
> >
> >
> >   
> 
>