ffc team mailing list archive
-
ffc team
-
Mailing list archive
-
Message #01458
Re: defining sudbomain integrals
On Wed, Jan 30, 2008 at 09:23:42PM -0500, 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?
>
> Thanks,
> Jake
Looks like an obvious bug that should be easy to fix. I will take a
look.
--
Anders
Follow ups
References