ffc team mailing list archive

[Jake Ostien <tostien@xxxxxxxxx>]

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?

Thanks,
Jake