message: Adding comment on pyramid integration rule implementation
diff: @@ -1659,6 +1659,10 @@ IntegrationRule *IntegrationRules::TetrahedronIntegrationRule(int Order)
  // Integration rules for reference pyramid
  IntegrationRule *IntegrationRules::PyramidIntegrationRule(int Order)
  {
  + // This is a simple integration rule adapted from an integration
  + // rule for a cube which seems to be adequate for now. When we
  + // implement high order finite elements for pyramids we should
  + // revisit this and see if we can improve upon it.
  IntegrationRule * irc = GenerateIntegrationRule(Geometry::CUBE, Order);
  int npts = irc->GetNPoints();
  AllocIntRule(PyramidIntRules, Order);
  @@ -1671,7 +1675,6 @@ IntegrationRule *IntegrationRules::PyramidIntegrationRule(int Order)
  ipp.x = ipc.x * (1.0 - ipc.z);
  ipp.y = ipc.y * (1.0 - ipc.z);
  ipp.z = ipc.z;
  - // ipp.weight = ipc.weight * pow(1.0 - ipc.z, 2);
  ipp.weight = ipc.weight / 3.0;
  }
  return PyramidIntRules[Order];
  