ffc team mailing list archive

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

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