message: Fixed GeometryRefiner::RefineInterior to not append existing integration rules.
diff: @@ -1344,8 +1344,8 @@ const IntegrationRule *GeometryRefiner::RefineInterior(Geometry::Type Geom,
  return NULL;
  }
  ir = FindInIntPts(Geom, Times-1);
  - if (ir == NULL)
  - {
  + if (ir) { return ir; }
  +
  ir = new IntegrationRule(Times-1);
  for (int i = 1; i < Times; i++)
  {
  @@ -1354,7 +1354,6 @@ const IntegrationRule *GeometryRefiner::RefineInterior(Geometry::Type Geom,
  ip.y = ip.z = 0.0;
  }
  }
  - }
  break;
  case Geometry::TRIANGLE:
  @@ -1364,8 +1363,8 @@ const IntegrationRule *GeometryRefiner::RefineInterior(Geometry::Type Geom,
  return NULL;
  }
  ir = FindInIntPts(Geom, ((Times-1)*(Times-2))/2);
  - if (ir == NULL)
  - {
  + if (ir) { return ir; }
  +
  ir = new IntegrationRule(((Times-1)*(Times-2))/2);
  for (int k = 0, j = 1; j < Times-1; j++)
  for (int i = 1; i < Times-j; i++, k++)
  @@ -1376,7 +1375,6 @@ const IntegrationRule *GeometryRefiner::RefineInterior(Geometry::Type Geom,
  ip.z = 0.0;
  }
  }
  - }
  break;
  case Geometry::SQUARE:
  @@ -1386,8 +1384,8 @@ const IntegrationRule *GeometryRefiner::RefineInterior(Geometry::Type Geom,
  return NULL;
  }
  ir = FindInIntPts(Geom, (Times-1)*(Times-1));
  - if (ir == NULL)
  - {
  + if (ir) { return ir; }
  +
  ir = new IntegrationRule((Times-1)*(Times-1));
  for (int k = 0, j = 1; j < Times; j++)
  for (int i = 1; i < Times; i++, k++)
  @@ -1398,7 +1396,6 @@ const IntegrationRule *GeometryRefiner::RefineInterior(Geometry::Type Geom,
  ip.z = 0.0;
  }
  }
  - }
  break;
  default:
  