message: Accidentally called GenerateIntegrationRule rather than IntegrationRules::Get
diff: @@ -1595,20 +1595,20 @@ IntegrationRule *IntegrationRules::TetrahedronIntegrationRule(int Order)
  // Integration rules for reference prism
  IntegrationRule *IntegrationRules::PrismIntegrationRule(int Order)
  {
  - IntegrationRule * irt = GenerateIntegrationRule(Geometry::TRIANGLE, Order);
  - IntegrationRule * irs = GenerateIntegrationRule(Geometry::SEGMENT, Order);
  - int nt = irt->GetNPoints();
  - int ns = irs->GetNPoints();
  + const IntegrationRule & irt = Get(Geometry::TRIANGLE, Order);
  + const IntegrationRule & irs = Get(Geometry::SEGMENT, Order);
  + int nt = irt.GetNPoints();
  + int ns = irs.GetNPoints();
  AllocIntRule(PrismIntRules, Order);
  PrismIntRules[Order] = new IntegrationRule(nt * ns);
  for (int ks=0; ks<ns; ks++)
  {
  - const IntegrationPoint & ips = irs->IntPoint(ks);
  + const IntegrationPoint & ips = irs.IntPoint(ks);
  for (int kt=0; kt<nt; kt++)
  {
  int kp = ks * nt + kt;
  - const IntegrationPoint & ipt = irt->IntPoint(kt);
  + const IntegrationPoint & ipt = irt.IntPoint(kt);
  IntegrationPoint & ipp = PrismIntRules[Order]->IntPoint(kp);
  ipp.x = ipt.x;
  ipp.y = ipt.y;
  