I forgot to multiply the value of the coefficient at the nodes by the value of the basis function at the quadrature point

[gmarkall]

In the quadrature loop nest, we get:

__shared__ double c_q0[24];
for(int i_g = 0; i_g < 6; i_g++)
{
  for(int i_d_0 = 0; i_d_0 < 2; i_d_0++)
  {
    for(int i_d_1 = 0; i_d_1 < 2; i_d_1++)
    {
      (c_q0[((i_g + (6 * i_d_0)) + (2 * (6 * i_d_1)))] = 0.0);
      for(int i_r_0 = 0; i_r_0 < 3; i_r_0++)
      {
        (c_q0[((i_g + (6 * i_d_0)) + (2 * (6 * i_d_1)))] += c0[((i_d_1 + (2 * i_d_0)) + (2 * (2 * i_r_0)))]);
      };
    };
  };
};

The RHS of the expression should be:

c0[((i_d_1 + (2 * i_d_0)) + (2 * (2 * i_r_0)))]) * CG0[subscript for current basis fn at current quad point]

[gmarkall]

Fixed in [master 5f3b7f6]