ffc team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

Anders Logg wrote:
> 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?
>

Nope. I'll send you the solver off-line.

Garth

> --
> Anders