ffc team mailing list archive

["Garth N. Wells" <g.n.wells@xxxxxxxxxx>]

Inside the function evaluate_dof, FFC is generating the code

  // Compute affine mapping x = F(X)
  coordinates[0] = w0*x[0][0] + w1*x[1][0] + w2*x[2][0] + w3*x[3][0];
  coordinates[0] = w0*x[0][1] + w1*x[1][1] + w2*x[2][1] + w3*x[3][1];
  coordinates[0] = w0*x[0][2] + w1*x[1][2] + w2*x[2][2] + w3*x[3][2];

which should be

  // Compute affine mapping x = F(X)
  coordinates[0] = w0*x[0][0] + w1*x[1][0] + w2*x[2][0] + w3*x[3][0];
  coordinates[1] = w0*x[0][1] + w1*x[1][1] + w2*x[2][1] + w3*x[3][1];
  coordinates[2] = w0*x[0][2] + w1*x[1][2] + w2*x[2][2] + w3*x[3][2];

(the index for the array coordinates should be 0, 1 ,2).

Garth