ffc team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Wed, Feb 17, 2010 at 11:19:53AM +0000, Garth N. Wells wrote:
> I've run into a tensor contraction bug. The jit log contains:
>
>  In member function ‘virtual void
> form_96980ae3ebe60648ee99493c681869c0377c12ac_exterior_facet_integral_0_0::tabulate_tensor(double*,
> const double* const*, const ufc::cell&, unsigned int) const’:
> form_96980ae3ebe60648ee99493c681869c0377c12ac.h:1952: error: ‘detJ’ was
> not declared in this scope
> error: command 'gcc' failed with exit status 1
>
> It seems to related to linear forms which have only an exterior facet
> integral. In the generated code I see:
>
>     . . . .
>     const double J_11 = x[2][1] - x[0][1];
>
>     // Compute determinant of Jacobian
>
>     // Compute inverse of Jacobian
>
>     // Get vertices on edge
>     static unsigned int edge_vertices[3][2] = {{1, 2}, {0, 2}, {0, 1}};
>     . . . .
>
> Which looks like the detJ has been stripped out, despite being used
> later in the code.
>
> Garth
>
> _______________________________________________
> Mailing list: https://launchpad.net/~ffc
> Post to     : ffc@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~ffc
> More help   : https://help.launchpad.net/ListHelp

ok. Do you have a simple test that fails?

--
Anders

Attachment: signature.asc
Description: Digital signature