static GOData *
gog_xy_dropbar_plot_axis_get_bounds (GogPlot *plot, GogAxisType axis,
GogPlotBoundInfo *bounds)
{
GogXYDropBarPlot *model = GOG_XY_DROPBAR_PLOT (plot);
if ((model->horizontal && axis == GOG_AXIS_Y) || (!model->horizontal && axis == GOG_AXIS_X)) {
GSList *ptr;
bounds->val.minima = model->x.minima;
bounds->val.maxima = model->x.maxima;
bounds->is_discrete = model->x.minima > model->x.maxima ||
!go_finite (model->x.minima) ||
!go_finite (model->x.maxima);
if (bounds->fmt == NULL && model->x.fmt != NULL)
bounds->fmt = go_format_ref (model->x.fmt);
if (model->x.date_conv)
bounds->date_conv = model->x.date_conv;
for (ptr = plot->series; ptr != NULL ; ptr = ptr->next)
if (gog_series_is_valid (GOG_SERIES (ptr->data)))
return GOG_SERIES (ptr->data)->values[0].data;
return NULL;
}
if ((model->horizontal && axis == GOG_AXIS_X) || (!model->horizontal && axis == GOG_AXIS_Y)) {
bounds->val.minima = model->y.minima;
bounds->val.maxima = model->y.maxima;
if (bounds->fmt == NULL && model->y.fmt != NULL)
bounds->fmt = go_format_ref (model->y.fmt);
if (model->y.date_conv)
bounds->date_conv = model->y.date_conv;
}
return NULL;
}
static gboolean
find_element (GogView *view, double cx, double cy, double x, double y,
unsigned int *index, GogPieSeries **series)
{
GogPiePlot *pie = GOG_PIE_PLOT (view->model);
GSList *ptr;
double *vals, scale, len = 0, theta;
*series = NULL;
*index = 0;
for (ptr = pie->base.series ; ptr != NULL ; ptr = ptr->next)
if (gog_series_is_valid (GOG_SERIES (*series = ptr->data)))
break;
if (ptr == NULL)
return FALSE;
theta = (atan2 (y - cy, x - cx)
* 180 / M_PI - pie->initial_angle + 90.) / pie->span / 3.6;
if (theta < 0)
theta += 1.;
vals = go_data_get_values ((*series)->base.values[1].data);
scale = 1 / (*series)->total;
for (*index = 0 ; *index < (*series)->base.num_elements; (*index)++) {
len = vals[*index] * scale;
if (len < 0.)
len = pie->show_negatives? -len: 0.;
if (go_finite (len) && len > 1e-3) {
theta -= len;
if (theta < 0)
break;
}
}
return TRUE;
}
static void
gog_xy_dropbar_view_render (GogView *view, GogViewAllocation const *bbox)
{
GogXYDropBarPlot const *model = GOG_XY_DROPBAR_PLOT (view->model);
GogPlot *plot = GOG_PLOT (model);
GogSeries const *series;
GogAxisMap *x_map, *y_map, *val_map, *index_map;
GogViewAllocation work;
double *pos_vals, *start_vals, *end_vals;
unsigned i;
GSList *ptr;
unsigned n, tmp, num_series;
GOStyle *neg_style;
for (num_series = 0, ptr = plot->series ; ptr != NULL ; ptr = ptr->next, num_series++);
if (num_series < 1)
return;
x_map = gog_axis_map_new (GOG_PLOT (model)->axis[0],
view->allocation.x, view->allocation.w);
y_map = gog_axis_map_new (GOG_PLOT (model)->axis[1], view->allocation.y + view->allocation.h,
-view->allocation.h);
if (!(gog_axis_map_is_valid (x_map) &&
gog_axis_map_is_valid (y_map))) {
gog_axis_map_free (x_map);
gog_axis_map_free (y_map);
return;
}
work.w = view->allocation.w * model->width / 100.;
for (ptr = plot->series ; ptr != NULL ; ptr = ptr->next) {
series = ptr->data;
if (!gog_series_is_valid (GOG_SERIES (series)))
continue;
neg_style = go_style_dup ((GOG_STYLED_OBJECT (series))->style);
neg_style->line.color ^= 0xffffff00;
neg_style->fill.pattern.back ^= 0xffffff00;
neg_style->fill.pattern.fore ^= 0xffffff00;
pos_vals = go_data_get_values (series->values[0].data);
n = go_data_get_vector_size (series->values[1].data);
start_vals = go_data_get_values (series->values[1].data);
tmp = go_data_get_vector_size (series->values[1].data);
if (n > tmp)
n = tmp;
end_vals = go_data_get_values (series->values[2].data);
tmp = go_data_get_vector_size (series->values[2].data);
if (n > tmp)
n = tmp;
if (model->horizontal) {
index_map = y_map;
val_map = x_map;
} else {
index_map = x_map;
val_map = y_map;
}
for (i = 0; i < n; i++) {
work.x = pos_vals[i];
work.y = start_vals[i];
work.h = end_vals[i] - work.y;
if (!gog_axis_map_finite (index_map, work.x) ||
!gog_axis_map_finite (val_map, start_vals[i]) ||
!gog_axis_map_finite (val_map, end_vals[i]))
continue;
gog_renderer_push_style (view->renderer, (start_vals[i] <= end_vals[i])?
GOG_STYLED_OBJECT (series)->style: neg_style);
barcol_draw_rect (view->renderer, model->horizontal, x_map, y_map, &work);
gog_renderer_pop_style (view->renderer);
}
g_object_unref (neg_style);
}
gog_axis_map_free (x_map);
gog_axis_map_free (y_map);
}
static void
gog_xy_dropbar_plot_update (GogObject *obj)
{
GogXYDropBarPlot *model = GOG_XY_DROPBAR_PLOT (obj);
GogSeries const *series = NULL;
double x_min, x_max, y_min, y_max, tmp_min, tmp_max;
GSList *ptr;
unsigned xaxis, yaxis;
if (model->horizontal) {
xaxis = 1;
yaxis = 0;
} else {
xaxis = 0;
yaxis = 1;
}
x_min = y_min = DBL_MAX;
x_max = y_max = -DBL_MAX;
gog_xy_dropbar_plot_clear_formats (model);
for (ptr = model->base.series ; ptr != NULL ; ptr = ptr->next) {
series = ptr->data;
if (!gog_series_is_valid (GOG_SERIES (series)))
continue;
go_data_get_bounds (series->values[1].data, &tmp_min, &tmp_max);
if (y_min > tmp_min) y_min = tmp_min;
if (y_max < tmp_max) y_max = tmp_max;
if (model->y.fmt == NULL) {
model->y.fmt = go_data_preferred_fmt (series->values[1].data);
model->y.date_conv = go_data_date_conv (series->values[1].data);
}
go_data_get_bounds (series->values[2].data, &tmp_min, &tmp_max);
if (y_min > tmp_min)
y_min = tmp_min;
if (y_max < tmp_max)
y_max = tmp_max;
if (series->values[0].data != NULL) {
go_data_get_bounds (series->values[0].data, &tmp_min, &tmp_max);
if (!go_finite (tmp_min) || !go_finite (tmp_max) ||
tmp_min > tmp_max) {
tmp_min = 0;
tmp_max = go_data_get_vector_size (series->values[1].data);
} else if (model->x.fmt == NULL) {
model->x.fmt = go_data_preferred_fmt (series->values[0].data);
model->x.date_conv = go_data_date_conv (series->values[0].data);
}
} else {
tmp_min = 0;
tmp_max = go_data_get_vector_size (series->values[1].data);
}
if (x_min > tmp_min)
x_min = tmp_min;
if (x_max < tmp_max)
x_max = tmp_max;
}
/* make room for bar width, this is approximate since real width is larger */
tmp_max = (x_max - x_min) * model->width / 200.;
x_min -= tmp_max;
x_max += tmp_max;
if (model->x.minima != x_min || model->x.maxima != x_max) {
model->x.minima = x_min;
model->x.maxima = x_max;
gog_axis_bound_changed (model->base.axis[xaxis], GOG_OBJECT (model));
}
if (model->y.minima != y_min || model->y.maxima != y_max) {
model->y.minima = y_min;
model->y.maxima = y_max;
gog_axis_bound_changed (model->base.axis[yaxis], GOG_OBJECT (model));
}
gog_object_emit_changed (GOG_OBJECT (obj), FALSE);
if (gog_xy_dropbar_parent_klass->update)
gog_xy_dropbar_parent_klass->update (obj);
}
static void
gog_pie_view_render (GogView *view, GogViewAllocation const *bbox)
{
GogPiePlot const *model = GOG_PIE_PLOT (view->model);
GogPieSeries const *series = NULL;
double separated_cx, separated_cy, cx, cy, r, theta, len, *vals, scale;
double default_sep;
unsigned elem, k;
GOPath *path;
GogTheme *theme = gog_object_get_theme (GOG_OBJECT (model));
GOStyle *style, *white_style = NULL, *elt_style = NULL;
GSList *ptr;
unsigned num_series = 0;
unsigned index;
double center_radius;
double center_size = 0.0;
double r_ext, r_int, r_tot;
double outline_width_max;
gboolean has_hole;
double separation_max, separation;
GogPieSeriesElement *gpse;
GList const *overrides;
GogShowNegsMode mode = model->show_negatives;
gboolean negative;
/* compute number of valid series */
for (ptr = model->base.series ; ptr != NULL ; ptr = ptr->next) {
series = ptr->data;
if (!gog_series_is_valid (GOG_SERIES (ptr->data)))
continue;
num_series++;
}
if (num_series <= 0)
return;
separation_max = .0;
outline_width_max = .0;
if ((style = go_styled_object_get_style (GO_STYLED_OBJECT (series))))
outline_width_max = gog_renderer_line_size (view->renderer, style->line.width);
if (mode == GOG_SHOW_NEGS_WHITE) {
white_style = go_style_dup (style);
white_style->fill.type = GO_STYLE_FILL_PATTERN;
go_pattern_set_solid (&white_style->fill.pattern, GO_COLOR_WHITE);
}
for (overrides = gog_series_get_overrides (GOG_SERIES (series));
overrides != NULL;
overrides = overrides->next) {
separation = GOG_PIE_SERIES_ELEMENT (overrides->data)->separation;
if (separation_max < separation)
separation_max = separation;
style = go_styled_object_get_style (GO_STYLED_OBJECT (overrides->data));
if (outline_width_max < style->line.width)
outline_width_max = style->line.width;
}
if (separation_max < -model->default_separation)
separation_max = -model->default_separation;
if (GOG_IS_RING_PLOT (model))
center_size = GOG_RING_PLOT(model)->center_size;
else if (num_series > 1)
num_series = 1;
/* calculate things when span < 100. */
if (model->span < 100.) {
double xmin, xmax, ymin, ymax, ratio;
double begin, cur, end;
if (num_series == 1) {
begin = (model->initial_angle + series->initial_angle) * M_PI / 180. - M_PI / 2.;
end = begin + model->span * M_PI / 50.;
} else {
/* WARNING: this code has not been checked */
ptr = model->base.series;
while (!gog_series_is_valid (GOG_SERIES (ptr->data)))
ptr = ptr->next;
series = ptr->data;
begin = end = series->initial_angle;
for (ptr = ptr->next; ptr != NULL; ptr = ptr->next) {
if (!gog_series_is_valid (GOG_SERIES (ptr->data)))
continue;
series = ptr->data;
cur = series->initial_angle;
if (cur < begin)
begin = cur;
if (cur > end)
end = cur;
}
begin = (model->initial_angle + begin) * M_PI / 180. - M_PI / 2.;
end = ((model->initial_angle + end) / 180. + model->span / 50.) * M_PI - M_PI / 2.;
}
cur = ceil (begin / M_PI * 2 - 1e-10) * M_PI / 2;
xmin = xmax = cos (begin);
ymin = ymax = sin (begin);
while (cur < end) {
cx = cos (cur);
cy = sin (cur);
if (cx > xmax)
xmax = cx;
if (cx < xmin)
xmin = cx;
if (cy > ymax)
ymax = cy;
if (cy < ymin)
ymin = cy;
cur += M_PI / 2;
}
cx = cos (end);
cy = sin (end);
if (cx > xmax)
xmax = cx;
if (cx < xmin)
xmin = cx;
if (cy > ymax)
ymax = cy;
if (cy < ymin)
ymin = cy;
/* we ensure that the center will be visible */
if (xmin > 0.)
xmin = 0.;
if (xmax < 0.)
xmax = 0.;
if (ymin > 0.)
ymin = 0.;
if (ymax < 0.)
ymax = 0.;
ratio = (ymax - ymin) / (xmax - xmin);
if (view->allocation.h > view->allocation.w * ratio) {
r_tot = view->allocation.w * MAX (xmax, -xmin) / (xmax - xmin);
cx = view->allocation.x - view->allocation.w * xmin / (xmax - xmin);
cy = view->allocation.y +
(view->allocation.h + view->allocation.w * ratio) / 2.
- view->allocation.w * ratio * ymax / (ymax - ymin);
} else {
r_tot = view->allocation.h * MAX (ymax, -ymin) / (ymax - ymin);
cx = view->allocation.x +
(view->allocation.w - view->allocation.h / ratio) / 2.
- view->allocation.h / ratio * xmin / (xmax - xmin);
cy = view->allocation.y - view->allocation.h * ymin / (ymax - ymin);
}
r_tot /= 1. + model->default_separation + separation_max;
} else {
/* centre things */
cx = view->allocation.x + view->allocation.w/2.;
cy = view->allocation.y + view->allocation.h/2.;
r_tot = view->allocation.h;
if (r_tot > view->allocation.w)
r_tot = view->allocation.w;
r_tot /= 2. * (1. + model->default_separation + separation_max);
}
path = go_path_new ();
default_sep = r_tot * model->default_separation;
center_radius = r_tot * center_size;
r = r_tot * (1. - center_size);
elem = model->base.index_num;
index = 1;
for (ptr = model->base.series ; ptr != NULL ; ptr = ptr->next) {
series = ptr->data;
if (!gog_series_is_valid (GOG_SERIES (series)))
continue;
if (index > num_series) /* people snuck extra series into a pie */
break;
if (num_series == index)
r -= outline_width_max / 2.0;
has_hole = center_radius > 0. || index > 1;
r_int = center_radius + r * ((double)index - 1.0) / (double)num_series;
r_ext = center_radius + r * (double)index / (double)num_series;
theta = (model->initial_angle + series->initial_angle) * M_PI / 180. - M_PI / 2.;
scale = 2 * M_PI / 100 * model->span / series->total;
vals = go_data_get_values (series->base.values[1].data);
style = GOG_STYLED_OBJECT (series)->style;
if (model->base.vary_style_by_element || mode == GOG_SHOW_NEGS_WHITE)
style = go_style_dup (style);
gog_renderer_push_style (view->renderer, style);
overrides = gog_series_get_overrides (GOG_SERIES (series));
for (k = 0 ; k < series->base.num_elements; k++) {
len = vals[k] * scale;
negative = len < 0;
if (negative) {
if (mode == GOG_SHOW_NEGS_SKIP) {
if ((overrides != NULL) &&
(GOG_SERIES_ELEMENT (overrides->data)->index == k))
overrides = overrides->next;
continue;
}
len = -len;
}
if (!go_finite (len) || len < 1e-3) {
if ((overrides != NULL) &&
(GOG_SERIES_ELEMENT (overrides->data)->index == k))
overrides = overrides->next;
continue;
}
gpse = NULL;
if ((overrides != NULL) &&
(GOG_SERIES_ELEMENT (overrides->data)->index == k)) {
gpse = GOG_PIE_SERIES_ELEMENT (overrides->data);
overrides = overrides->next;
if (negative && mode == GOG_SHOW_NEGS_WHITE) {
elt_style = go_style_dup (go_styled_object_get_style (
GO_STYLED_OBJECT (gpse)));
elt_style->fill.type = GO_STYLE_FILL_PATTERN;
go_pattern_set_solid (&elt_style->fill.pattern, GO_COLOR_WHITE);
gog_renderer_push_style (view->renderer, elt_style);
} else
gog_renderer_push_style (view->renderer,
go_styled_object_get_style (
GO_STYLED_OBJECT (gpse)));
} else {
if (negative && mode == GOG_SHOW_NEGS_WHITE)
gog_renderer_push_style (view->renderer, white_style);
else if (model->base.vary_style_by_element)
gog_theme_fillin_style (theme, style, GOG_OBJECT (series),
model->base.index_num + k, GO_STYLE_FILL);
}
/* only separate the outer ring */
separated_cx = cx;
separated_cy = cy;
if (num_series == index && (default_sep > 0. || gpse != NULL)) {
separation = default_sep;
if (gpse != NULL)
separation += gpse->separation * r_tot;
separated_cx += separation * cos (theta + len/2.);
separated_cy += separation * sin (theta + len/2.);
}
theta += len;
go_path_ring_wedge (path, separated_cx, separated_cy,
r_ext, r_ext,
has_hole ? r_int : 0.0,
has_hole ? r_int : 0.0,
theta - len,
theta);
gog_renderer_draw_shape (view->renderer, path);
go_path_clear (path);
if (gpse != NULL || (negative && mode == GOG_SHOW_NEGS_WHITE)) {
gog_renderer_pop_style (view->renderer);
if (elt_style) {
g_object_unref (elt_style);
elt_style = NULL;
}
}
}
gog_renderer_pop_style (view->renderer);
if (model->base.vary_style_by_element || mode == GOG_SHOW_NEGS_WHITE)
g_object_unref (style);
if (white_style)
g_object_unref (white_style);
index ++;
}
go_path_free (path);
}
static int
gog_pie_view_get_data_at_point (GogPlotView *view, double x, double y, GogSeries **series)
{
GogPiePlot const *model = GOG_PIE_PLOT (view->base.model);
GogPieSeries const *pseries = NULL;
double r_tot, r_cur, r_int, r_ext, cx, cy, r, th, th0, theta, scale, *vals;
double separated_cx, separated_cy;
unsigned int index, elem, k;
double center_radius;
double center_size = 0.0;
unsigned num_series = 0;
double default_sep, len;
double separation_max = 0., separation;
GSList *ptr;
double outline_width_max = 0.;
gboolean has_hole;
GogPieSeriesElement *gpse;
GList const *overrides;
GogShowNegsMode mode = model->show_negatives;
gboolean negative;
GOStyle *style;
*series = NULL;
/* compute number of valid series */
for (ptr = model->base.series ; ptr != NULL ; ptr = ptr->next) {
if (!gog_series_is_valid (GOG_SERIES (ptr->data)))
continue;
pseries = ptr->data;
if ((style = go_styled_object_get_style (GO_STYLED_OBJECT (pseries))))
outline_width_max = MAX (outline_width_max, gog_renderer_line_size (view->base.renderer, style->line.width));
for (overrides = gog_series_get_overrides (GOG_SERIES (pseries));
overrides != NULL;
overrides = overrides->next) {
separation = GOG_PIE_SERIES_ELEMENT (overrides->data)->separation;
if (separation_max < separation)
separation_max = separation;
style = go_styled_object_get_style (GO_STYLED_OBJECT (overrides->data));
if (outline_width_max < style->line.width)
outline_width_max = style->line.width;
}
num_series++;
}
if (separation_max < -model->default_separation)
separation_max = -model->default_separation;
if (num_series <= 0)
return -1;
if (GOG_IS_RING_PLOT (model))
center_size = GOG_RING_PLOT(model)->center_size;
else if (num_series > 1)
num_series = 1;
/* calculate things when span < 100. */
if (model->span < 100.) {
double xmin, xmax, ymin, ymax, ratio;
double begin, cur, end;
if (num_series == 1) {
pseries = model->base.series->data;
begin = (model->initial_angle + pseries->initial_angle) * M_PI / 180. - M_PI / 2.;
end = begin + model->span * M_PI / 50.;
} else {
/* WARNING: this code has not been checked */
ptr = model->base.series;
while (!gog_series_is_valid (GOG_SERIES (ptr->data)))
ptr = ptr->next;
pseries = ptr->data;
begin = end = pseries->initial_angle;
for (ptr = ptr->next; ptr != NULL; ptr = ptr->next) {
if (!gog_series_is_valid (GOG_SERIES (ptr->data)))
continue;
pseries = ptr->data;
cur = pseries->initial_angle;
if (cur < begin)
begin = cur;
if (cur > end)
end = cur;
}
begin = (model->initial_angle + begin) * M_PI / 180. - M_PI / 2.;
end = ((model->initial_angle + end) / 180. + model->span / 50.) * M_PI - M_PI / 2.;
}
cur = ceil (begin / M_PI * 2 - 1e-10) * M_PI / 2;
xmin = xmax = cos (begin);
ymin = ymax = sin (begin);
while (cur < end) {
cx = cos (cur);
cy = sin (cur);
if (cx > xmax)
xmax = cx;
if (cx < xmin)
xmin = cx;
if (cy > ymax)
ymax = cy;
if (cy < ymin)
ymin = cy;
cur += M_PI / 2;
}
cx = cos (end);
cy = sin (end);
if (cx > xmax)
xmax = cx;
if (cx < xmin)
xmin = cx;
if (cy > ymax)
ymax = cy;
if (cy < ymin)
ymin = cy;
/* we ensure that the center will be visible */
if (xmin > 0.)
xmin = 0.;
if (xmax < 0.)
xmax = 0.;
if (ymin > 0.)
ymin = 0.;
if (ymax < 0.)
ymax = 0.;
ratio = (ymax - ymin) / (xmax - xmin);
if (view->base.allocation.h > view->base.allocation.w * ratio) {
r_tot = view->base.allocation.w * MAX (xmax, -xmin) / (xmax - xmin);
cx = view->base.allocation.x - view->base.allocation.w * xmin / (xmax - xmin);
cy = view->base.allocation.y +
(view->base.allocation.h + view->base.allocation.w * ratio) / 2.
- view->base.allocation.w * ratio * ymax / (ymax - ymin);
} else {
r_tot = view->base.allocation.h * MAX (ymax, -ymin) / (ymax - ymin);
cx = view->base.allocation.x +
(view->base.allocation.w - view->base.allocation.h / ratio) / 2.
- view->base.allocation.h / ratio * xmin / (xmax - xmin);
cy = view->base.allocation.y - view->base.allocation.h * ymin / (ymax - ymin);
}
r_tot /= 1. + model->default_separation + separation_max;
} else {
/* centre things */
cx = view->base.allocation.x + view->base.allocation.w/2.;
cy = view->base.allocation.y + view->base.allocation.h/2.;
r_tot = view->base.allocation.h;
if (r_tot > view->base.allocation.w)
r_tot = view->base.allocation.w;
r_tot /= 2. * (1. + model->default_separation + separation_max);
}
r = hypot (x - cx, y - cy);
default_sep = r_tot * model->default_separation;
center_radius = r_tot * center_size;
r_cur = r_tot * (1. - center_size);
if (r < center_radius)
return -1;
elem = model->base.index_num;
index = 1;
th0 = (model->initial_angle + pseries->initial_angle) * M_PI / 180. - M_PI / 2.;
th = atan2 (y - cy, x - cx);
if (th < th0)
th += 2 * M_PI;
for (ptr = model->base.series ; ptr != NULL ; ptr = ptr->next) {
pseries = ptr->data;
if (!gog_series_is_valid (GOG_SERIES (pseries)))
continue;
if (index > num_series) /* people snuck extra series into a pie */
break;
if (num_series == index)
r_cur -= outline_width_max / 2.0;
has_hole = center_radius > 0. || index > 1;
r_int = (has_hole)? (center_radius + r_cur * ((double) index - 1.0) / (double) num_series): 0.;
r_ext = center_radius + r_cur * (double) index / (double) num_series;
theta = th0;
scale = 2 * M_PI / 100 * model->span / pseries->total;
vals = go_data_get_values (pseries->base.values[1].data);
overrides = gog_series_get_overrides (GOG_SERIES (pseries));
for (k = 0 ; k < pseries->base.num_elements; k++) {
len = vals[k] * scale;
negative = len < 0;
if (negative) {
if (mode == GOG_SHOW_NEGS_SKIP) {
if ((overrides != NULL) &&
(GOG_SERIES_ELEMENT (overrides->data)->index == k))
overrides = overrides->next;
continue;
}
len = -len;
}
if (!go_finite (len) || len < 1e-3) {
if ((overrides != NULL) &&
(GOG_SERIES_ELEMENT (overrides->data)->index == k))
overrides = overrides->next;
continue;
}
gpse = NULL;
if ((overrides != NULL) &&
(GOG_SERIES_ELEMENT (overrides->data)->index == k)) {
gpse = GOG_PIE_SERIES_ELEMENT (overrides->data);
overrides = overrides->next;
}
/* only separate the outer ring */
separated_cx = cx;
separated_cy = cy;
if (num_series == index && (default_sep > 0. || gpse != NULL)) {
separation = default_sep;
if (gpse != NULL)
separation += gpse->separation * r_tot;
separated_cx += separation * cos (theta + len/2.);
separated_cy += separation * sin (theta + len/2.);
r = hypot (x - separated_cx, y - separated_cy);
th = atan2 (y - separated_cy, x - separated_cx);
if (th < th0)
th += 2 * M_PI;
}
theta += len;
if (r > r_int && r <= r_ext && th > theta - len && th <= theta) {
*series = GOG_SERIES (pseries);
return k;
}
if (gpse) {
r = hypot (x - cx, y - cy);
th = atan2 (y - cy, x - cx);
if (th < th0)
th += 2 * M_PI;
}
}
index ++;
}
return -1;
}
