static void
update_data (GtkPlotData *data, gboolean new_range)
{
if(new_range && data->is_function)
GTK_PLOT_SURFACE(data)->recalc_dt = TRUE;
GTK_PLOT_DATA_CLASS(parent_class)->update(data, new_range);
}
GtkWidget*
gtk_plot_surface_new_function (GtkPlotFunc3D function)
{
GtkWidget *dataset;
dataset = gtk_type_new (gtk_plot_surface_get_type ());
gtk_plot_surface_construct_function(GTK_PLOT_SURFACE(dataset), function);
return dataset;
}
static gboolean
gtk_plot_surface_add_to_plot (GtkPlotData *data, GtkPlot *plot)
{
if(!data->is_function && data->num_points > 0){
gtk_plot_surface_build_mesh(GTK_PLOT_SURFACE(data));
}
data->redraw_pending = TRUE;
return TRUE;
}
void
gtk_plot_surface_destroy(GtkObject *object)
{
GtkPlotSurface *surface;
surface = GTK_PLOT_SURFACE(object);
if(surface->dt)
gtk_object_destroy(GTK_OBJECT(surface->dt));
if ( GTK_OBJECT_CLASS (parent_class)->destroy )
(* GTK_OBJECT_CLASS (parent_class)->destroy) (object);
}
static void
gtk_plot_surface_draw_private (GtkPlotData *data)
{
GtkPlot *plot;
GtkPlotSurface *surface;
g_return_if_fail(GTK_IS_PLOT_SURFACE(data));
surface = GTK_PLOT_SURFACE(data);
data = GTK_PLOT_DATA(surface);
g_return_if_fail(GTK_PLOT_DATA(data)->plot != NULL);
g_return_if_fail(GTK_IS_PLOT(GTK_PLOT_DATA(data)->plot));
if(data->redraw_pending) gtk_plot_surface_update_range(data);
plot = GTK_PLOT(data->plot);
gtk_plot_surface_real_draw(surface);
}
int main(int argc, char *argv[]){
GtkWidget *window1;
GtkWidget *vbox1;
GtkWidget *scrollw1;
GtkWidget *canvas;
GtkWidget *button;
GtkWidget *surface;
GtkPlotCanvasChild *child;
gint page_width, page_height;
gfloat scale = 1.;
gint n;
static gdouble *x, *y, *z;
page_width = GTK_PLOT_LETTER_W * scale;
page_height = GTK_PLOT_LETTER_H * scale;
gtk_init(&argc,&argv);
window1=gtk_window_new(GTK_WINDOW_TOPLEVEL);
gtk_window_set_title(GTK_WINDOW(window1), "GtkPlot3D Demo");
gtk_widget_set_usize(window1,550,650);
gtk_container_border_width(GTK_CONTAINER(window1),0);
gtk_signal_connect (GTK_OBJECT (window1), "destroy",
GTK_SIGNAL_FUNC (quit), NULL);
vbox1=gtk_vbox_new(FALSE,0);
gtk_container_add(GTK_CONTAINER(window1),vbox1);
gtk_widget_show(vbox1);
scrollw1=gtk_scrolled_window_new(NULL, NULL);
gtk_container_border_width(GTK_CONTAINER(scrollw1),0);
gtk_scrolled_window_set_policy(GTK_SCROLLED_WINDOW(scrollw1),
GTK_POLICY_ALWAYS,GTK_POLICY_ALWAYS);
gtk_box_pack_start(GTK_BOX(vbox1),scrollw1, TRUE, TRUE,0);
gtk_widget_show(scrollw1);
canvas = gtk_plot_canvas_new(page_width, page_height, 1.);
GTK_PLOT_CANVAS_SET_FLAGS(GTK_PLOT_CANVAS(canvas), GTK_PLOT_CANVAS_DND_FLAGS);
gtk_scrolled_window_add_with_viewport(GTK_SCROLLED_WINDOW(scrollw1), canvas);
/*
gdk_color_parse("light blue", &color);
gdk_color_alloc(gtk_widget_get_colormap(canvas), &color);
gtk_plot_canvas_set_background(GTK_PLOT_CANVAS(canvas), &color);
*/
gtk_widget_show(canvas);
plot3d = GTK_PLOT3D(new_layer(canvas, GTK_TYPE_PLOT3D));
child = gtk_plot_canvas_plot_new(GTK_PLOT(plot3d));
gtk_plot_canvas_put_child(GTK_PLOT_CANVAS(canvas), child, .16, .05, .65, .45);
gtk_widget_show(GTK_WIDGET(plot3d));
x = (gdouble *)g_malloc(1600 * sizeof(gdouble));
y = (gdouble *)g_malloc(1600 * sizeof(gdouble));
z = (gdouble *)g_malloc(1600 * sizeof(gdouble));
n = 0;
/*
for(nx = 0; nx < 40; nx++)
for(ny = 0; ny < 40; ny++)
{
x[n]=0.025*(gdouble)nx;
y[n]=0.025*(gdouble)ny;
z[n]=cos(((x[n]-0.5)*(x[n]-0.5) + (y[n]-0.5)*(y[n]-0.5))*24) / 4. + .5;
n++;
}
surface = gtk_plot_csurface_new();
gtk_plot_surface_set_points(GTK_PLOT_SURFACE(surface),
x, y, z, NULL, NULL, NULL, 40, 40);
gtk_plot_data_set_legend(GTK_PLOT_DATA(surface), "cos ((r-r\\s0\\N)\\S2\\N)");
*/
/*
for(nx = 0; nx < 10; nx++)
for(ny = 0; ny < 10; ny++)
{
x[n]=(gdouble)0.1*(gdouble)nx;
y[n]=(gdouble)0.1*(gdouble)ny;
z[n]=(x[n]-0.5)*(x[n]-0.5) + (y[n]-0.5)*(y[n]-0.5);
z[n]=x[n]*x[n] + y[n]*y[n];
n++;
}
surface = gtk_plot_csurface_new();
gtk_plot_surface_set_points(GTK_PLOT_SURFACE(surface),
x, y, z, NULL, NULL, NULL, 10, 10);
*/
gtk_plot3d_set_minor_ticks(GTK_PLOT3D(plot3d), GTK_PLOT_AXIS_X, 3);
gtk_plot3d_set_minor_ticks(GTK_PLOT3D(plot3d), GTK_PLOT_AXIS_Y, 3);
gtk_plot3d_minor_grids_set_visible(GTK_PLOT3D(plot3d), FALSE, FALSE, FALSE);
surface = gtk_plot_csurface_new_function((GtkPlotFunc3D) function);
gtk_plot_surface_set_xstep(GTK_PLOT_SURFACE(surface), .05);
gtk_plot_surface_set_ystep(GTK_PLOT_SURFACE(surface), .05);
gtk_plot_data_set_legend(GTK_PLOT_DATA(surface), "cos ((r-r\\s0\\N)\\S2\\N)");
gtk_plot_add_data(GTK_PLOT(plot3d), GTK_PLOT_DATA(surface));
gtk_widget_show(surface);
gtk_plot_data_set_gradient(GTK_PLOT_DATA(surface), .2, .8, 6, 0);
/*
gtk_plot_surface_set_transparent(GTK_PLOT_SURFACE(surface), TRUE);
*/
gtk_plot_surface_set_transparent(GTK_PLOT_SURFACE(surface), FALSE);
gtk_plot_csurface_set_projection(GTK_PLOT_CSURFACE(surface), GTK_PLOT_PROJECT_EMPTY);
GTK_PLOT_CSURFACE(surface)->levels_line.line_style = GTK_PLOT_LINE_NONE;
GTK_PLOT_CSURFACE(surface)->sublevels_line.line_style = GTK_PLOT_LINE_NONE;
/*******************/
active_plot = new_layer(canvas, GTK_TYPE_PLOT);
child = gtk_plot_canvas_plot_new(GTK_PLOT(active_plot));
gtk_plot_canvas_put_child(GTK_PLOT_CANVAS(canvas), child, .26, .54, .65, .90);
gtk_plot_axis_set_minor_ticks(gtk_plot_get_axis(GTK_PLOT(active_plot), GTK_PLOT_AXIS_LEFT), 1);
gtk_plot_axis_set_minor_ticks(gtk_plot_get_axis(GTK_PLOT(active_plot), GTK_PLOT_AXIS_BOTTOM), 1);
surface = gtk_plot_csurface_new_function((GtkPlotFunc3D) function);
gtk_plot_surface_set_xstep(GTK_PLOT_SURFACE(surface), .050);
gtk_plot_surface_set_ystep(GTK_PLOT_SURFACE(surface), .050);
gtk_plot_data_set_legend(GTK_PLOT_DATA(surface), "cos ((r-r\\s0\\N)\\S2\\N)");
gtk_plot_add_data(GTK_PLOT(active_plot), GTK_PLOT_DATA(surface));
gtk_widget_show(surface);
gtk_plot_data_set_gradient(GTK_PLOT_DATA(surface), .2, .8, 3, 0);
gtk_plot_surface_set_grid_visible(GTK_PLOT_SURFACE(surface), FALSE);
gtk_plot_surface_set_transparent(GTK_PLOT_SURFACE(surface), TRUE);
gtk_plot_csurface_set_projection(GTK_PLOT_CSURFACE(surface), GTK_PLOT_PROJECT_FULL);
/*
gtk_plot_csurface_set_lines_visible(GTK_PLOT_CSURFACE(surface), FALSE);
*/
button = gtk_button_new_with_label("Rotate X");
gtk_fixed_put(GTK_FIXED(canvas), button, 80, 0);
gtk_widget_show(button);
gtk_signal_connect(GTK_OBJECT(button), "clicked",
GTK_SIGNAL_FUNC(rotate_x), canvas);
button = gtk_button_new_with_label("Rotate Y");
gtk_fixed_put(GTK_FIXED(canvas), button, 160, 0);
gtk_widget_show(button);
gtk_signal_connect(GTK_OBJECT(button), "clicked",
GTK_SIGNAL_FUNC(rotate_y), canvas);
button = gtk_button_new_with_label("Rotate Z");
gtk_fixed_put(GTK_FIXED(canvas), button, 240, 0);
gtk_widget_show(button);
gtk_signal_connect(GTK_OBJECT(button), "clicked",
GTK_SIGNAL_FUNC(rotate_z), canvas);
gtk_signal_connect(GTK_OBJECT(canvas), "select_item",
(GtkSignalFunc) select_item, NULL);
gtk_widget_show(window1);
gtk_plot_canvas_export_ps(GTK_PLOT_CANVAS(canvas), "democsurface.ps", 0, 0,
GTK_PLOT_LETTER);
gtk_main();
return(0);
}
int main(int argc, char *argv[]){
GtkWidget *window1;
GtkWidget *vbox1;
GtkWidget *scrollw1;
GtkWidget *canvas;
GtkWidget *button;
GtkWidget *surface;
gdouble *x,*y,*z,*a;
gint cnt= 0;
gint page_width, page_height;
gfloat scale = 1.;
char buffer[1000];
FILE *f;
gdouble xmin=1e99;
gdouble xmax=-1e99;
gdouble ymin=1e99;
gdouble ymax=-1e99;
gdouble zmin=1e99;
gdouble zmax=-1e99;
gdouble dx,dy,dz;
gint n, nx, ny;
GtkPlotCanvasChild *child;
page_width = GTK_PLOT_LETTER_W * scale;
page_height = GTK_PLOT_LETTER_H * scale;
gtk_init(&argc,&argv);
window1=gtk_window_new(GTK_WINDOW_TOPLEVEL);
gtk_window_set_title(GTK_WINDOW(window1), "GtkPlot3D Demo");
gtk_widget_set_usize(window1,550,650);
gtk_container_border_width(GTK_CONTAINER(window1),0);
gtk_signal_connect (GTK_OBJECT (window1), "destroy",
GTK_SIGNAL_FUNC (quit), NULL);
vbox1=gtk_vbox_new(FALSE,0);
gtk_container_add(GTK_CONTAINER(window1),vbox1);
scrollw1=gtk_scrolled_window_new(NULL, NULL);
gtk_container_border_width(GTK_CONTAINER(scrollw1),0);
gtk_scrolled_window_set_policy(GTK_SCROLLED_WINDOW(scrollw1),
GTK_POLICY_ALWAYS,GTK_POLICY_ALWAYS);
gtk_box_pack_start(GTK_BOX(vbox1),scrollw1, TRUE, TRUE,0);
canvas = gtk_plot_canvas_new(page_width, page_height, 1.);
GTK_PLOT_CANVAS_SET_FLAGS(GTK_PLOT_CANVAS(canvas), GTK_PLOT_CANVAS_DND_FLAGS);
gtk_scrolled_window_add_with_viewport(GTK_SCROLLED_WINDOW(scrollw1), canvas);
/*
gdk_color_parse("light blue", &color);
gdk_color_alloc(gtk_widget_get_colormap(canvas), &color);
gtk_plot_canvas_set_background(GTK_PLOT_CANVAS(canvas), &color);
*/
active_plot = new_layer(canvas);
child = gtk_plot_canvas_plot_new(GTK_PLOT(active_plot));
gtk_plot_canvas_put_child(GTK_PLOT_CANVAS(canvas), child, .10, .06, .85, .85);
x = (gdouble *)g_malloc(1600 * sizeof(gdouble));
y = (gdouble *)g_malloc(1600 * sizeof(gdouble));
z = (gdouble *)g_malloc(1600 * sizeof(gdouble));
a = (gdouble *)g_malloc(1600 * sizeof(gdouble));
n = 0;
for(nx = 0; nx < 40; nx++)
for(ny = 0; ny < 40; ny++)
{
x[n]=0.025*(gdouble)nx;
y[n]=0.025*(gdouble)ny;
z[n]=cos(((x[n]-0.5)*(x[n]-0.5) + (y[n]-0.5)*(y[n]-0.5))*24) / 4. + .5;
a[n]=(x[n]-0.5)*(x[n]-0.5) + (y[n]-0.5)*(y[n]-0.5);
n++;
}
surface = gtk_plot_surface_new();
gtk_plot_surface_set_points(GTK_PLOT_SURFACE(surface),
x, y, z, NULL, NULL, NULL, 40, 40);
gtk_plot_data_set_legend(GTK_PLOT_DATA(surface), "cos ((r-r\\s0\\N)\\S2\\N)");
gtk_plot_surface_use_amplitud(GTK_PLOT_SURFACE(surface), TRUE);
gtk_plot_data_set_a(GTK_PLOT_DATA(surface), a);
gtk_plot_data_gradient_set_visible(GTK_PLOT_DATA(surface), TRUE);
/*
for(nx = 0; nx < 10; nx++)
for(ny = 0; ny < 10; ny++)
{
x[n]=(gdouble)0.1*(gdouble)nx;
y[n]=(gdouble)0.1*(gdouble)ny;
z[n]=(x[n]-0.5)*(x[n]-0.5) + (y[n]-0.5)*(y[n]-0.5);
n++;
}
surface = gtk_plot_surface_new();
gtk_plot_surface_set_points(GTK_PLOT_SURFACE(surface),
x, y, z, NULL, NULL, NULL, 10, 10);
*/
gtk_plot3d_set_minor_ticks(GTK_PLOT3D(active_plot), GTK_PLOT_AXIS_X, 1);
gtk_plot3d_set_minor_ticks(GTK_PLOT3D(active_plot), GTK_PLOT_AXIS_Y, 1);
gtk_plot3d_show_ticks(GTK_PLOT3D(active_plot),
GTK_PLOT_SIDE_XY,
GTK_PLOT_TICKS_OUT, GTK_PLOT_TICKS_OUT);
gtk_plot3d_show_ticks(GTK_PLOT3D(active_plot),
GTK_PLOT_SIDE_XZ,
GTK_PLOT_TICKS_OUT, GTK_PLOT_TICKS_OUT);
gtk_plot3d_show_ticks(GTK_PLOT3D(active_plot),
GTK_PLOT_SIDE_YX,
GTK_PLOT_TICKS_OUT, GTK_PLOT_TICKS_OUT);
gtk_plot3d_show_ticks(GTK_PLOT3D(active_plot),
GTK_PLOT_SIDE_YZ,
GTK_PLOT_TICKS_OUT, GTK_PLOT_TICKS_OUT);
gtk_plot3d_show_ticks(GTK_PLOT3D(active_plot),
GTK_PLOT_SIDE_ZX,
GTK_PLOT_TICKS_OUT, GTK_PLOT_TICKS_OUT);
gtk_plot3d_show_ticks(GTK_PLOT3D(active_plot),
GTK_PLOT_SIDE_ZY,
GTK_PLOT_TICKS_OUT, GTK_PLOT_TICKS_OUT);
/* if a file was specified on the commandline, try to open and display it: */
if (argc!=2 || !(f=fopen(argv[1],"r"))) {
/*
x= y= z= NULL;
surface = gtk_plot_surface_new_function((GtkPlotFunc3D) function);
gtk_plot_surface_set_xstep(GTK_PLOT_SURFACE(surface), .025);
gtk_plot_surface_set_ystep(GTK_PLOT_SURFACE(surface), .025);
gtk_plot_data_set_legend(GTK_PLOT_DATA(surface), "cos ((r-r\\s0\\N)\\S2\\N)");
*/
gtk_plot_add_data(GTK_PLOT(active_plot), GTK_PLOT_DATA(surface));
} else {
x= g_new(gdouble,MAXNODES);
y= g_new(gdouble,MAXNODES);
z= g_new(gdouble,MAXNODES);
while (fgets(buffer,1000,f) && cnt<MAXNODES)
if (sscanf(buffer,"%lf %lf %lf", &x[cnt], &y[cnt], &z[cnt])==3) {
z[cnt]*=10; /* REMOVE THIS BEFORE FLIGHT */
if (xmin>x[cnt]) xmin= x[cnt];
if (xmax<x[cnt]) xmax= x[cnt];
if (ymin>y[cnt]) ymin= y[cnt];
if (ymax<y[cnt]) ymax= y[cnt];
if (zmin>z[cnt]) zmin= z[cnt];
if (zmax<z[cnt]) zmax= z[cnt];
cnt++;
}
fclose(f);
dx= (xmax-xmin)*.02;
dy= (ymax-ymin)*.02;
dz= (zmax-zmin)*.02;
/* start the triangulation */
fprintf(stderr,"data ranges from (%g,%g) to (%g,%g)\n",
xmin,ymin,xmax,ymax);
/*
gtk_plot3d_set_xrange(GTK_PLOT3D(active_plot), xmin-dx, xmax+dx);
gtk_plot3d_set_yrange(GTK_PLOT3D(active_plot), ymin-dy, ymax+dy);
gtk_plot3d_set_zrange(GTK_PLOT3D(active_plot), zmin-dz, zmax+dz);
gtk_plot3d_set_ticks(GTK_PLOT3D(active_plot), GTK_PLOT_AXIS_X, (xmax-xmin)/10., 1);
gtk_plot3d_set_ticks(GTK_PLOT3D(active_plot), GTK_PLOT_AXIS_Y, (ymax-ymin)/10., 1);
gtk_plot3d_set_ticks(GTK_PLOT3D(active_plot), GTK_PLOT_AXIS_Z, (zmax-zmin)/10., 1);
*/
surface = gtk_plot_csurface_new();
gtk_plot_add_data(GTK_PLOT(active_plot), GTK_PLOT_DATA(surface));
gtk_plot_surface_set_points(GTK_PLOT_SURFACE(surface),
x, y, z, NULL, NULL, NULL, cnt, 1);
gtk_plot_data_set_legend(GTK_PLOT_DATA(surface), argv[1]);
gtk_plot_data_gradient_autoscale_z(GTK_PLOT_DATA(surface));
}
/*
gtk_plot_surface_use_height_gradient(GTK_PLOT_SURFACE(surface), TRUE);
gtk_plot_data_set_gradient(GTK_PLOT_DATA(surface), .2, .8, 6, 0);
*/
/*
gtk_plot_data_gradient_set_scale(GTK_PLOT_DATA(surface),GTK_PLOT_SCALE_LOG10);
*/
gtk_plot3d_autoscale(GTK_PLOT3D(active_plot));
// gtk_plot3d_set_xrange(GTK_PLOT3D(active_plot), -0.5,1);
// gtk_plot3d_set_yrange(GTK_PLOT3D(active_plot), -0.2,1);
// gtk_plot3d_set_zrange(GTK_PLOT3D(active_plot), 0.,1);
gtk_widget_show(surface);
button = gtk_button_new_with_label("Rotate X");
gtk_fixed_put(GTK_FIXED(canvas), button, 80, 0);
gtk_signal_connect(GTK_OBJECT(button), "clicked",
GTK_SIGNAL_FUNC(rotate_x), canvas);
button = gtk_button_new_with_label("Rotate Y");
gtk_fixed_put(GTK_FIXED(canvas), button, 160, 0);
gtk_signal_connect(GTK_OBJECT(button), "clicked",
GTK_SIGNAL_FUNC(rotate_y), canvas);
button = gtk_button_new_with_label("Rotate Z");
gtk_fixed_put(GTK_FIXED(canvas), button, 240, 0);
/*
gtk_plot3d_set_xfactor(GTK_PLOT3D(active_plot), 0.5);
*/
/*
gtk_plot3d_set_xrange(GTK_PLOT3D(active_plot), 0., 2.0);
*/
gtk_signal_connect(GTK_OBJECT(button), "clicked",
GTK_SIGNAL_FUNC(rotate_z), canvas);
gtk_widget_show_all(window1);
/*
gtk_plot_canvas_export_ps(GTK_PLOT_CANVAS(canvas), "demo3d.ps", 0, 0,
GTK_PLOT_LETTER);
*/
gtk_plot_canvas_export_ps_with_size(GTK_PLOT_CANVAS(canvas), "demo3d.ps",
0, 0, GTK_PLOT_PSPOINTS,
GTK_PLOT_LETTER_W, GTK_PLOT_LETTER_H);
gtk_main();
return(0);
}
static void
gtk_plot_surface_draw_legend(GtkPlotData *data, gint x, gint y)
{
GtkPlotSurface *surface;
GtkPlot *plot = NULL;
GtkPlotText legend;
GdkRectangle area;
gint lascent, ldescent, lheight, lwidth;
gdouble m;
surface = GTK_PLOT_SURFACE(data);
g_return_if_fail(data->plot != NULL);
g_return_if_fail(GTK_IS_PLOT(data->plot));
plot = data->plot;
area.x = GTK_WIDGET(plot)->allocation.x;
area.y = GTK_WIDGET(plot)->allocation.y;
area.width = GTK_WIDGET(plot)->allocation.width;
area.height = GTK_WIDGET(plot)->allocation.height;
m = plot->magnification;
legend = plot->legends_attr;
if(data->legend)
legend.text = data->legend;
else
legend.text = "";
gtk_plot_text_get_size(legend.text, legend.angle, legend.font,
roundint(legend.height * m),
&lwidth, &lheight,
&lascent, &ldescent);
if(!surface->height_gradient){
gtk_plot_pc_set_color(plot->pc, &surface->color);
gtk_plot_pc_draw_rectangle(plot->pc, TRUE,
x, y,
roundint(plot->legends_line_width * m),
lascent + ldescent);
}else{
gdouble level, step;
gint lx = x, lstep;
step = (data->gradient.end - data->gradient.begin) / 10;
lstep = roundint(plot->legends_line_width * m / 10.);
for(level = data->gradient.begin; level < data->gradient.end; level += step){
GdkColor color;
gtk_plot_data_get_gradient_level(data, level, &color);
gtk_plot_pc_set_color(plot->pc, &color);
gtk_plot_pc_draw_rectangle(plot->pc, TRUE,
lx, y,
lstep, lascent + ldescent);
lx += lstep;
}
}
legend.x = (gdouble)(area.x + x + roundint((plot->legends_line_width + 4) * m))
/ (gdouble)area.width;
legend.y = (gdouble)(area.y + y + lascent) / (gdouble)area.height;
gtk_plot_draw_text(plot, legend);
y += 2*lheight;
gtk_plot_data_draw_gradient(data, x, y);
}
static void
gtk_plot_surface_update_range (GtkPlotData *data)
{
GtkPlotSurface *surface;
GtkPlot *plot;
if(!data->redraw_pending) return;
surface = GTK_PLOT_SURFACE(data);
plot = data->plot;
if(!surface->recalc_dt && data->num_points > 0){
gtk_plot_surface_recalc_nodes(surface);
data->redraw_pending = FALSE;
surface->recalc_dt = FALSE;
return;
}
if(!data->is_function && data->num_points > 0){
if(surface->recalc_dt)
gtk_plot_surface_build_mesh(surface);
else
gtk_plot_surface_recalc_nodes(surface);
data->redraw_pending = FALSE;
surface->recalc_dt = FALSE;
return;
}
if(data->is_function){
gdouble xstep, ystep;
gdouble xmin, xmax, ymin, ymax;
gdouble x, y;
gdouble *fx = NULL, *fy = NULL, *fz = NULL;
gdouble *old_x, *old_y, *old_z;
gint nx, ny;
gint npoints;
xmin = GTK_PLOT(plot)->xmin;
xmax = GTK_PLOT(plot)->xmax;
ymin = GTK_PLOT(plot)->ymin;
ymax = GTK_PLOT(plot)->ymax;
xstep = surface->xstep;
surface->nx = roundint((xmax - xmin) / xstep) + 1;
ystep = surface->ystep;
surface->ny = roundint((ymax - ymin) / ystep) + 1;
npoints = surface->nx * surface->ny;
fx = (gdouble *)g_malloc((npoints + 1) * sizeof(gdouble));
fy = (gdouble *)g_malloc((npoints + 1) * sizeof(gdouble));
fz = (gdouble *)g_malloc((npoints + 1) * sizeof(gdouble));
npoints = 0;
y = ymin;
for(ny = 0; ny < surface->ny; ny++)
{
x = xmin;
for(nx = 0; nx < surface->nx; nx++)
{
gboolean error;
fx[npoints] = x;
fy[npoints] = y;
fz[npoints] = data->function3d(plot, data, x, y, &error);
x += xstep;
npoints++;
}
y += ystep;
}
old_x = data->x; old_y = data->y; old_z = data->z;
data->x = fx;
data->y = fy;
data->z = fz;
data->num_points = npoints;
gtk_plot_surface_build_mesh(surface);
g_free(fx);
g_free(fy);
g_free(fz);
data->x = old_x; data->y = old_y; data->z = old_z;
}
data->redraw_pending = FALSE;
surface->recalc_dt = FALSE;
}
