alloc_rule_pts()

static int alloc_rule_pts ( rule *  r, unsigned  num_regions  )

static

{
     if (num_regions > r->num_regions) {
          free(r->pts);
          r->pts = r->vals = NULL;
          r->num_regions = 0;
          num_regions *= 2; /* allocate extra so that
                               repeatedly calling alloc_rule_pts with
                               growing num_regions only needs
                               a logarithmic number of allocations */
          r->pts = (double *) malloc(sizeof(double) * 
                                     (num_regions
                                      * r->num_points * (r->dim + r->fdim)));
          if (r->fdim + r->dim > 0 && !r->pts) return FAILURE;
          r->vals = r->pts + num_regions * r->num_points * r->dim;
          r->num_regions = num_regions;
     }
     return SUCCESS;
}