1226
1226
// Compute circumradius, assuming triangle is embedded in 2D.
1229
1231
// Array of quadrature weights.
1230
1232
static const double W16[16] = {0.023568368, 0.035388068, 0.022584049, 0.0054232259, 0.044185089, 0.066344216, 0.042339725, 0.01016726, 0.044185089, 0.066344216, 0.042339725, 0.01016726, 0.023568368, 0.035388068, 0.022584049, 0.0054232259};
1231
1233
// Quadrature points on the UFC reference element: (0.065466995, 0.057104196), (0.050210123, 0.27684301), (0.028912084, 0.58359043), (0.0097037851, 0.86024014), (0.31116455, 0.057104196), (0.23864866, 0.27684301), (0.1374191, 0.58359043), (0.04612208, 0.86024014), (0.63173125, 0.057104196), (0.48450833, 0.27684301), (0.27899046, 0.58359043), (0.093637784, 0.86024014), (0.87742881, 0.057104196), (0.67294686, 0.27684301), (0.38749748, 0.58359043), (0.13005608, 0.86024014)