xref: /petsc/src/sys/classes/draw/impls/x/xops.c (revision 09440f2536211689395dd83abbfd05fbf503dd2c)
1 /*
2     Defines the operations for the X PetscDraw implementation.
3 */
4 
5 #include <../src/sys/classes/draw/impls/x/ximpl.h>         /*I  "petscsys.h" I*/
6 
7 /*
8      These macros transform from the users coordinates to the  X-window pixel coordinates.
9 */
10 #define XTRANS(draw,xwin,x)  (int)(((xwin)->w)*((draw)->port_xl + (((x - (draw)->coor_xl)*((draw)->port_xr - (draw)->port_xl))/((draw)->coor_xr - (draw)->coor_xl))))
11 #define YTRANS(draw,xwin,y)  (int)(((xwin)->h)*(1.0-(draw)->port_yl - (((y - (draw)->coor_yl)*((draw)->port_yr - (draw)->port_yl))/((draw)->coor_yr - (draw)->coor_yl))))
12 
13 #define ITRANS(draw,xwin,i)  (draw)->coor_xl + (i*((draw)->coor_xr - (draw)->coor_xl)/((xwin)->w) - (draw)->port_xl)/((draw)->port_xr - (draw)->port_xl)
14 #define JTRANS(draw,xwin,j)  draw->coor_yl + (((double)j)/xwin->h + draw->port_yl - 1.0)*(draw->coor_yr - draw->coor_yl)/(draw->port_yl - draw->port_yr)
15 
16 #undef __FUNCT__
17 #define __FUNCT__ "PetscDrawSetViewport_X"
18 static PetscErrorCode PetscDrawSetViewport_X(PetscDraw draw,PetscReal xl,PetscReal yl,PetscReal xr,PetscReal yr)
19 {
20   PetscDraw_X *XiWin = (PetscDraw_X*)draw->data;
21   XRectangle  box;
22 
23   PetscFunctionBegin;
24   box.x     = (int)(xl*XiWin->w);     box.y      = (int)((1.0-yr)*XiWin->h);
25   box.width = (int)((xr-xl)*XiWin->w);box.height = (int)((yr-yl)*XiWin->h);
26   XSetClipRectangles(XiWin->disp,XiWin->gc.set,0,0,&box,1,Unsorted);
27   PetscFunctionReturn(0);
28 }
29 
30 #undef __FUNCT__
31 #define __FUNCT__ "PetscDrawCoordinateToPixel_X"
32 static PetscErrorCode PetscDrawCoordinateToPixel_X(PetscDraw draw,PetscReal x,PetscReal y,PetscInt *i,PetscInt *j)
33 {
34   PetscDraw_X *XiWin = (PetscDraw_X*)draw->data;
35 
36   PetscFunctionBegin;
37   *i = XTRANS(draw,XiWin,x);
38   *j = YTRANS(draw,XiWin,y);
39   PetscFunctionReturn(0);
40 }
41 
42 #undef __FUNCT__
43 #define __FUNCT__ "PetscDrawPixelToCoordinate_X"
44 static PetscErrorCode PetscDrawPixelToCoordinate_X(PetscDraw draw,PetscInt i,PetscInt j,PetscReal *x,PetscReal *y)
45 {
46   PetscDraw_X *XiWin = (PetscDraw_X*)draw->data;
47 
48   PetscFunctionBegin;
49   *x = ITRANS(draw,XiWin,i);
50   *y = JTRANS(draw,XiWin,j);
51   PetscFunctionReturn(0);
52 }
53 
54 #undef __FUNCT__
55 #define __FUNCT__ "PetscDrawPoint_X"
56 static PetscErrorCode PetscDrawPoint_X(PetscDraw draw,PetscReal x,PetscReal y,int c)
57 {
58   PetscDraw_X *XiWin = (PetscDraw_X*)draw->data;
59   int         xx,yy,i,j;
60 
61   PetscFunctionBegin;
62   xx = XTRANS(draw,XiWin,x);  yy = YTRANS(draw,XiWin,y);
63   PetscDrawXiSetColor(XiWin,c);
64   for (i=-1; i<2; i++) {
65     for (j=-1; j<2; j++) {
66       XDrawPoint(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,xx+i,yy+j);
67     }
68   }
69   PetscFunctionReturn(0);
70 }
71 
72 #undef __FUNCT__
73 #define __FUNCT__ "PetscDrawPointPixel_X"
74 static PetscErrorCode PetscDrawPointPixel_X(PetscDraw draw,PetscInt x,PetscInt y,int c)
75 {
76   PetscDraw_X *XiWin = (PetscDraw_X*)draw->data;
77 
78   PetscFunctionBegin;
79   PetscDrawXiSetColor(XiWin,c);
80   XDrawPoint(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,x,y);
81   PetscFunctionReturn(0);
82 }
83 
84 #undef __FUNCT__
85 #define __FUNCT__ "PetscDrawLine_X"
86 static PetscErrorCode PetscDrawLine_X(PetscDraw draw,PetscReal xl,PetscReal yl,PetscReal xr,PetscReal yr,int cl)
87 {
88   PetscDraw_X *XiWin = (PetscDraw_X*)draw->data;
89   int         x1,y_1,x2,y2;
90 
91   PetscFunctionBegin;
92   PetscDrawXiSetColor(XiWin,cl);
93   x1  = XTRANS(draw,XiWin,xl);   x2  = XTRANS(draw,XiWin,xr);
94   y_1 = YTRANS(draw,XiWin,yl);   y2  = YTRANS(draw,XiWin,yr);
95   if (x1 == x2 && y_1 == y2) PetscFunctionReturn(0);
96   XDrawLine(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,x1,y_1,x2,y2);
97   PetscFunctionReturn(0);
98 }
99 
100 #undef __FUNCT__
101 #define __FUNCT__ "PetscDrawArrow_X"
102 static PetscErrorCode PetscDrawArrow_X(PetscDraw draw,PetscReal xl,PetscReal yl,PetscReal xr,PetscReal yr,int cl)
103 {
104   PetscDraw_X *XiWin = (PetscDraw_X*)draw->data;
105   int         x1,y_1,x2,y2;
106 
107   PetscFunctionBegin;
108   PetscDrawXiSetColor(XiWin,cl);
109   x1  = XTRANS(draw,XiWin,xl);   x2  = XTRANS(draw,XiWin,xr);
110   y_1 = YTRANS(draw,XiWin,yl);   y2  = YTRANS(draw,XiWin,yr);
111   if (x1 == x2 && y_1 == y2) PetscFunctionReturn(0);
112   XDrawLine(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,x1,y_1,x2,y2);
113   if (x1 == x2 && PetscAbs(y_1 - y2) > 7) {
114     if (y2 > y_1) {
115       XDrawLine(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,x2,y2,x2-3,y2-3);
116       XDrawLine(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,x2,y2,x2+3,y2-3);
117     } else {
118       XDrawLine(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,x2,y2,x2-3,y2+3);
119       XDrawLine(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,x2,y2,x2+3,y2+3);
120     }
121   }
122   PetscFunctionReturn(0);
123 }
124 
125 #undef __FUNCT__
126 #define __FUNCT__ "PetscDrawRectangle_X"
127 static PetscErrorCode PetscDrawRectangle_X(PetscDraw draw,PetscReal xl,PetscReal yl,PetscReal xr,PetscReal yr,int c1,int c2,int c3,int c4)
128 {
129   PetscDraw_X *XiWin = (PetscDraw_X*)draw->data;
130   int         x1,y_1,w,h,c = (c1 + c2 + c3 + c4)/4;
131 
132   PetscFunctionBegin;
133   PetscDrawXiSetColor(XiWin,c);
134   x1  = XTRANS(draw,XiWin,xl);   w  = XTRANS(draw,XiWin,xr) - x1;
135   y_1 = YTRANS(draw,XiWin,yr);   h  = YTRANS(draw,XiWin,yl) - y_1;
136   if (w <= 0) w = 1;
137   if (h <= 0) h = 1;
138   XFillRectangle(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,x1,y_1,w,h);
139   PetscFunctionReturn(0);
140 }
141 
142 #undef __FUNCT__
143 #define __FUNCT__ "PetscDrawEllipse_X"
144 static PetscErrorCode PetscDrawEllipse_X(PetscDraw Win,PetscReal x,PetscReal y,PetscReal a,PetscReal b,int c)
145 {
146   PetscDraw_X *XiWin = (PetscDraw_X*) Win->data;
147   int         xA,yA,w,h;
148 
149   PetscFunctionBegin;
150   PetscDrawXiSetColor(XiWin, c);
151   xA = XTRANS(Win, XiWin, x - a/2.0); w = XTRANS(Win, XiWin, x + a/2.0) - xA;
152   yA = YTRANS(Win, XiWin, y + b/2.0); h = PetscAbs(YTRANS(Win, XiWin, y - b/2.0) - yA);
153   XFillArc(XiWin->disp, PetscDrawXiDrawable(XiWin), XiWin->gc.set, xA, yA, w, h, 0, 23040);
154   PetscFunctionReturn(0);
155 }
156 
157 PETSC_INTERN PetscErrorCode PetscDrawInterpolatedTriangle_X(PetscDraw_X*,int,int,int,int,int,int,int,int,int);
158 
159 #undef __FUNCT__
160 #define __FUNCT__ "PetscDrawTriangle_X"
161 static PetscErrorCode PetscDrawTriangle_X(PetscDraw draw,PetscReal X1,PetscReal Y_1,PetscReal X2,PetscReal Y2,PetscReal X3,PetscReal Y3,int c1,int c2,int c3)
162 {
163   PetscDraw_X    *XiWin = (PetscDraw_X*)draw->data;
164   PetscErrorCode ierr;
165 
166   PetscFunctionBegin;
167   if (c1 == c2 && c2 == c3) {
168     XPoint pt[3];
169     PetscDrawXiSetColor(XiWin,c1);
170     pt[0].x = XTRANS(draw,XiWin,X1);
171     pt[0].y = YTRANS(draw,XiWin,Y_1);
172     pt[1].x = XTRANS(draw,XiWin,X2);
173     pt[1].y = YTRANS(draw,XiWin,Y2);
174     pt[2].x = XTRANS(draw,XiWin,X3);
175     pt[2].y = YTRANS(draw,XiWin,Y3);
176     XFillPolygon(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,pt,3,Convex,CoordModeOrigin);
177   } else {
178     int x1,y_1,x2,y2,x3,y3;
179     x1   = XTRANS(draw,XiWin,X1);
180     y_1  = YTRANS(draw,XiWin,Y_1);
181     x2   = XTRANS(draw,XiWin,X2);
182     y2   = YTRANS(draw,XiWin,Y2);
183     x3   = XTRANS(draw,XiWin,X3);
184     y3   = YTRANS(draw,XiWin,Y3);
185     ierr = PetscDrawInterpolatedTriangle_X(XiWin,x1,y_1,c1,x2,y2,c2,x3,y3,c3);CHKERRQ(ierr);
186   }
187   PetscFunctionReturn(0);
188 }
189 
190 #undef __FUNCT__
191 #define __FUNCT__ "PetscDrawStringSetSize_X"
192 static PetscErrorCode PetscDrawStringSetSize_X(PetscDraw draw,PetscReal x,PetscReal y)
193 {
194   PetscDraw_X    *XiWin = (PetscDraw_X*)draw->data;
195   int            w,h;
196   PetscErrorCode ierr;
197 
198   PetscFunctionBegin;
199   w    = (int)((XiWin->w)*x*(draw->port_xr - draw->port_xl)/(draw->coor_xr - draw->coor_xl));
200   h    = (int)((XiWin->h)*y*(draw->port_yr - draw->port_yl)/(draw->coor_yr - draw->coor_yl));
201   ierr = PetscFree(XiWin->font);CHKERRQ(ierr);
202   ierr = PetscDrawXiFontFixed(XiWin,w,h,&XiWin->font);CHKERRQ(ierr);
203   PetscFunctionReturn(0);
204 }
205 
206 #undef __FUNCT__
207 #define __FUNCT__ "PetscDrawStringGetSize_X"
208 static PetscErrorCode PetscDrawStringGetSize_X(PetscDraw draw,PetscReal *x,PetscReal  *y)
209 {
210   PetscDraw_X *XiWin = (PetscDraw_X*)draw->data;
211   PetscReal   w,h;
212 
213   PetscFunctionBegin;
214   w = XiWin->font->font_w; h = XiWin->font->font_h;
215   if (x) *x = w*(draw->coor_xr - draw->coor_xl)/((XiWin->w)*(draw->port_xr - draw->port_xl));
216   if (y) *y = h*(draw->coor_yr - draw->coor_yl)/((XiWin->h)*(draw->port_yr - draw->port_yl));
217   PetscFunctionReturn(0);
218 }
219 
220 #undef __FUNCT__
221 #define __FUNCT__ "PetscDrawString_X"
222 static PetscErrorCode PetscDrawString_X(PetscDraw draw,PetscReal x,PetscReal y,int c,const char chrs[])
223 {
224   PetscDraw_X    *XiWin = (PetscDraw_X*)draw->data;
225   int            xx,yy;
226   size_t         len;
227   char           *substr;
228   PetscToken     token;
229   PetscErrorCode ierr;
230 
231   PetscFunctionBegin;
232   xx = XTRANS(draw,XiWin,x);
233   yy = YTRANS(draw,XiWin,y);
234   PetscDrawXiSetColor(XiWin,c);
235 
236   ierr = PetscTokenCreate(chrs,'\n',&token);CHKERRQ(ierr);
237   ierr = PetscTokenFind(token,&substr);CHKERRQ(ierr);
238   ierr = PetscStrlen(substr,&len);CHKERRQ(ierr);
239   XDrawString(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,xx,yy - XiWin->font->font_descent,substr,len);
240   ierr = PetscTokenFind(token,&substr);CHKERRQ(ierr);
241   while (substr) {
242     yy  += 4*XiWin->font->font_descent;
243     ierr = PetscStrlen(substr,&len);CHKERRQ(ierr);
244     XDrawString(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,xx,yy - XiWin->font->font_descent,substr,len);
245     ierr = PetscTokenFind(token,&substr);CHKERRQ(ierr);
246   }
247   ierr = PetscTokenDestroy(&token);CHKERRQ(ierr);
248   PetscFunctionReturn(0);
249 }
250 
251 #undef __FUNCT__
252 #define __FUNCT__ "PetscDrawStringVertical_X"
253 static PetscErrorCode PetscDrawStringVertical_X(PetscDraw draw,PetscReal x,PetscReal y,int c,const char chrs[])
254 {
255   PetscDraw_X    *XiWin = (PetscDraw_X*)draw->data;
256   int            xx,yy;
257   char           tmp[2];
258   PetscReal      tw,th;
259   size_t         i,n;
260   PetscErrorCode ierr;
261 
262   PetscFunctionBegin;
263   ierr   = PetscStrlen(chrs,&n);CHKERRQ(ierr);
264   tmp[1] = 0;
265   PetscDrawXiSetColor(XiWin,c);
266   ierr = PetscDrawStringGetSize_X(draw,&tw,&th);CHKERRQ(ierr);
267   xx   = XTRANS(draw,XiWin,x);
268   for (i=0; i<n; i++) {
269     tmp[0] = chrs[i];
270     yy     = YTRANS(draw,XiWin,y-th*i);
271     XDrawString(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set, xx,yy - XiWin->font->font_descent,tmp,1);
272   }
273   PetscFunctionReturn(0);
274 }
275 
276 #undef __FUNCT__
277 #define __FUNCT__ "PetscDrawFlush_X"
278 static PetscErrorCode PetscDrawFlush_X(PetscDraw draw)
279 {
280   PetscDraw_X    *XiWin = (PetscDraw_X*)draw->data;
281   PetscMPIInt    rank;
282   PetscErrorCode ierr;
283 
284   PetscFunctionBegin;
285   ierr = MPI_Comm_rank(PetscObjectComm((PetscObject)draw),&rank);CHKERRQ(ierr);
286 
287   ierr = PetscDrawCollectiveBegin(draw);CHKERRQ(ierr);
288   XSync(XiWin->disp,False);
289   ierr = PetscDrawCollectiveEnd(draw);CHKERRQ(ierr);
290   ierr = MPI_Barrier(PetscObjectComm((PetscObject)draw));CHKERRQ(ierr);
291 
292   ierr = PetscDrawCollectiveBegin(draw);CHKERRQ(ierr);
293   if (!rank && XiWin->drw && XiWin->win) {
294     XCopyArea(XiWin->disp,XiWin->drw,XiWin->win,XiWin->gc.set,0,0,XiWin->w,XiWin->h,0,0);
295     XSync(XiWin->disp,False);
296   }
297   ierr = PetscDrawCollectiveEnd(draw);CHKERRQ(ierr);
298   ierr = MPI_Barrier(PetscObjectComm((PetscObject)draw));CHKERRQ(ierr);
299 
300   ierr = PetscDrawSave(draw);CHKERRQ(ierr);
301   PetscFunctionReturn(0);
302 }
303 
304 #undef __FUNCT__
305 #define __FUNCT__ "PetscDrawClear_X"
306 static PetscErrorCode PetscDrawClear_X(PetscDraw draw)
307 {
308   PetscDraw_X    *XiWin = (PetscDraw_X*)draw->data;
309   PetscMPIInt    rank;
310   PetscErrorCode ierr;
311 
312   PetscFunctionBegin;
313   ierr = MPI_Comm_rank(PetscObjectComm((PetscObject)draw),&rank);CHKERRQ(ierr);
314 
315   ierr = PetscDrawCollectiveBegin(draw);CHKERRQ(ierr);
316   XSync(XiWin->disp,False);
317   ierr = PetscDrawCollectiveEnd(draw);CHKERRQ(ierr);
318   ierr = MPI_Barrier(PetscObjectComm((PetscObject)draw));CHKERRQ(ierr);
319 
320   ierr = PetscDrawCollectiveBegin(draw);CHKERRQ(ierr);
321   if (!rank) {
322     int x = (int)(draw->port_xl*XiWin->w);
323     int y = (int)((1.0-draw->port_yr)*XiWin->h);
324     int w = (int)((draw->port_xr - draw->port_xl)*XiWin->w);
325     int h = (int)((draw->port_yr - draw->port_yl)*XiWin->h);
326     PetscDrawXiSetPixVal(XiWin,XiWin->background);
327     XFillRectangle(XiWin->disp,PetscDrawXiDrawable(XiWin),XiWin->gc.set,x,y,w,h);
328     XSync(XiWin->disp,False);
329   }
330   ierr = PetscDrawCollectiveEnd(draw);CHKERRQ(ierr);
331   ierr = MPI_Barrier(PetscObjectComm((PetscObject)draw));CHKERRQ(ierr);
332   PetscFunctionReturn(0);
333 }
334 
335 #undef __FUNCT__
336 #define __FUNCT__ "PetscDrawSetDoubleBuffer_X"
337 static PetscErrorCode PetscDrawSetDoubleBuffer_X(PetscDraw draw)
338 {
339   PetscDraw_X    *win = (PetscDraw_X*)draw->data;
340   PetscMPIInt    rank;
341   PetscErrorCode ierr;
342 
343   PetscFunctionBegin;
344   if (win->drw) PetscFunctionReturn(0);
345   ierr = MPI_Comm_rank(PetscObjectComm((PetscObject)draw),&rank);CHKERRQ(ierr);
346 
347   ierr = PetscDrawCollectiveBegin(draw);CHKERRQ(ierr);
348   if (!rank) {ierr = PetscDrawXiQuickPixmap(win);CHKERRQ(ierr);}
349   ierr = PetscDrawCollectiveEnd(draw);CHKERRQ(ierr);
350   ierr = MPI_Bcast(&win->drw,1,MPI_UNSIGNED_LONG,0,PetscObjectComm((PetscObject)draw));CHKERRQ(ierr);
351   PetscFunctionReturn(0);
352 }
353 
354 #include <X11/cursorfont.h>
355 
356 #undef __FUNCT__
357 #define __FUNCT__ "PetscDrawGetMouseButton_X"
358 static PetscErrorCode PetscDrawGetMouseButton_X(PetscDraw draw,PetscDrawButton *button,PetscReal *x_user,PetscReal *y_user,PetscReal *x_phys,PetscReal *y_phys)
359 {
360   PetscDraw_X    *win = (PetscDraw_X*)draw->data;
361   Cursor         cursor;
362   XEvent         report;
363   Window         root,child;
364   int            root_x,root_y,px=0,py=0;
365   unsigned int   w,h,border,depth;
366   unsigned int   keys_button;
367   PetscMPIInt    rank;
368   PetscErrorCode ierr;
369 
370   PetscFunctionBegin;
371   *button = PETSC_BUTTON_NONE;
372   if (!win->win) PetscFunctionReturn(0);
373   ierr = MPI_Comm_rank(PetscObjectComm((PetscObject)draw),&rank);CHKERRQ(ierr);
374 
375   ierr = PetscDrawCollectiveBegin(draw);CHKERRQ(ierr);
376   if (rank) goto finally;
377 
378   /* change cursor to indicate input */
379   cursor = XCreateFontCursor(win->disp,XC_hand2); if (!cursor) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_LIB,"Unable to create X cursor");
380   XDefineCursor(win->disp,win->win,cursor);
381   /* wait for mouse button events */
382   XSelectInput(win->disp,win->win,ButtonPressMask|ButtonReleaseMask);
383   while (XCheckTypedEvent(win->disp,ButtonPress,&report));
384   XMaskEvent(win->disp,ButtonReleaseMask,&report);
385   /* get mouse pointer coordinates */
386   XQueryPointer(win->disp,report.xmotion.window,&root,&child,&root_x,&root_y,&px,&py,&keys_button);
387   /* the user may resize the window before pressing the mouse button */
388   XGetGeometry(win->disp,win->win,&root,&root_x,&root_y,&w,&h,&border,&depth);
389   /* cleanup input event handler and cursor  */
390   XSelectInput(win->disp,win->win,NoEventMask);
391   XUndefineCursor(win->disp,win->win);
392   XFreeCursor(win->disp, cursor);
393   XSync(win->disp,False);
394 
395   switch (report.xbutton.button) {
396   case Button1: *button = PETSC_BUTTON_LEFT; break;
397   case Button2: *button = PETSC_BUTTON_CENTER; break;
398   case Button3: *button = PETSC_BUTTON_RIGHT; break;
399   case Button4: *button = PETSC_BUTTON_WHEEL_UP; break;
400   case Button5: *button = PETSC_BUTTON_WHEEL_DOWN; break;
401   }
402   if (report.xbutton.state & ShiftMask) {
403     switch (report.xbutton.button) {
404     case Button1: *button = PETSC_BUTTON_LEFT_SHIFT; break;
405     case Button2: *button = PETSC_BUTTON_CENTER_SHIFT; break;
406     case Button3: *button = PETSC_BUTTON_RIGHT_SHIFT; break;
407     }
408   }
409   if (x_phys) *x_phys = ((double)px)/((double)w);
410   if (y_phys) *y_phys = 1.0 - ((double)py)/((double)h);
411   if (x_user) *x_user = draw->coor_xl + ((((double)px)/((double)w)-draw->port_xl))*(draw->coor_xr - draw->coor_xl)/(draw->port_xr - draw->port_xl);
412   if (y_user) *y_user = draw->coor_yl + ((1.0 - ((double)py)/((double)h)-draw->port_yl))*(draw->coor_yr - draw->coor_yl)/(draw->port_yr - draw->port_yl);
413 
414 finally:
415   ierr = PetscDrawCollectiveEnd(draw);CHKERRQ(ierr);
416   PetscFunctionReturn(0);
417 }
418 
419 #undef __FUNCT__
420 #define __FUNCT__ "PetscDrawPause_X"
421 static PetscErrorCode PetscDrawPause_X(PetscDraw draw)
422 {
423   PetscDraw_X    *win = (PetscDraw_X*)draw->data;
424   PetscErrorCode ierr;
425 
426   PetscFunctionBegin;
427   if (!win->win) PetscFunctionReturn(0);
428   if (draw->pause > 0) PetscSleep(draw->pause);
429   else if (draw->pause == -1) {
430     PetscDrawButton button = PETSC_BUTTON_NONE;
431     ierr = PetscDrawGetMouseButton(draw,&button,NULL,NULL,NULL,NULL);CHKERRQ(ierr);
432     if (button == PETSC_BUTTON_CENTER) draw->pause = 0;
433   }
434   PetscFunctionReturn(0);
435 }
436 
437 #undef __FUNCT__
438 #define __FUNCT__ "PetscDrawGetPopup_X"
439 static PetscErrorCode PetscDrawGetPopup_X(PetscDraw draw,PetscDraw *popup)
440 {
441   PetscDraw_X    *win = (PetscDraw_X*)draw->data;
442   PetscBool      flg  = PETSC_TRUE;
443   PetscErrorCode ierr;
444 
445   PetscFunctionBegin;
446   ierr = PetscOptionsGetBool(((PetscObject)draw)->options,((PetscObject)draw)->prefix,"-draw_popup",&flg,NULL);CHKERRQ(ierr);
447   if (!flg || !win->win) {*popup = NULL; PetscFunctionReturn(0);}
448 
449   ierr = PetscDrawCreate(PetscObjectComm((PetscObject)draw),draw->display,NULL,win->x,win->y+win->h+36,220,220,popup);CHKERRQ(ierr);
450   ierr = PetscObjectSetOptionsPrefix((PetscObject)*popup,"popup_");CHKERRQ(ierr);
451   ierr = PetscObjectAppendOptionsPrefix((PetscObject)*popup,((PetscObject)draw)->prefix);CHKERRQ(ierr);
452   ierr = PetscDrawSetType(*popup,PETSC_DRAW_X);CHKERRQ(ierr);
453   draw->popup = *popup;
454   PetscFunctionReturn(0);
455 }
456 
457 #undef __FUNCT__
458 #define __FUNCT__ "PetscDrawSetTitle_X"
459 static PetscErrorCode PetscDrawSetTitle_X(PetscDraw draw,const char title[])
460 {
461   PetscDraw_X    *win = (PetscDraw_X*)draw->data;
462   PetscMPIInt    rank;
463   PetscErrorCode ierr;
464 
465   PetscFunctionBegin;
466   if (!win->win) PetscFunctionReturn(0);
467   ierr = MPI_Comm_rank(PetscObjectComm((PetscObject)draw),&rank);CHKERRQ(ierr);
468   ierr = PetscDrawCollectiveBegin(draw);CHKERRQ(ierr);
469   if (!rank) {
470     size_t        len;
471     XTextProperty prop;
472     ierr = PetscStrlen(title,&len);CHKERRQ(ierr);
473     XGetWMName(win->disp,win->win,&prop);
474     XFree((void*)prop.value);
475     prop.value  = (unsigned char*)title;
476     prop.nitems = (long)len;
477     XSetWMName(win->disp,win->win,&prop);
478   }
479   ierr = PetscDrawCollectiveEnd(draw);CHKERRQ(ierr);
480   PetscFunctionReturn(0);
481 }
482 
483 #undef __FUNCT__
484 #define __FUNCT__ "PetscDrawResizeWindow_X"
485 static PetscErrorCode PetscDrawResizeWindow_X(PetscDraw draw,int w,int h)
486 {
487   PetscDraw_X    *win = (PetscDraw_X*)draw->data;
488   PetscMPIInt    rank;
489   PetscErrorCode ierr;
490 
491   PetscFunctionBegin;
492   if (!win->win) PetscFunctionReturn(0);
493   ierr = MPI_Comm_rank(PetscObjectComm((PetscObject)draw),&rank);CHKERRQ(ierr);
494 
495   ierr = PetscDrawCollectiveBegin(draw);CHKERRQ(ierr);
496   if (!rank) XResizeWindow(win->disp,win->win,w,h);
497   ierr = PetscDrawCollectiveEnd(draw);CHKERRQ(ierr);
498 
499   ierr = PetscDrawCheckResizedWindow(draw);CHKERRQ(ierr);
500   PetscFunctionReturn(0);
501 }
502 
503 #undef __FUNCT__
504 #define __FUNCT__ "PetscDrawCheckResizedWindow_X"
505 static PetscErrorCode PetscDrawCheckResizedWindow_X(PetscDraw draw)
506 {
507   PetscDraw_X    *win = (PetscDraw_X*)draw->data;
508   int            xywh[4];
509   PetscMPIInt    rank;
510   PetscReal      xl,xr,yl,yr;
511   PetscErrorCode ierr;
512 
513   PetscFunctionBegin;
514   if (!win->win) PetscFunctionReturn(0);
515   ierr = MPI_Comm_rank(PetscObjectComm((PetscObject)draw),&rank);CHKERRQ(ierr);
516 
517   ierr = PetscDrawCollectiveBegin(draw);CHKERRQ(ierr);
518   if (!rank) {
519     Window       root;
520     int          x,y;
521     unsigned int w,h,border,depth;
522     XGetGeometry(win->disp,win->win,&root,&x,&y,&w,&h,&border,&depth);
523     xywh[0] = x; xywh[1] = y; xywh[2] = (int)w; xywh[3] = (int)h;
524   }
525   ierr = PetscDrawCollectiveEnd(draw);CHKERRQ(ierr);
526 
527   ierr = MPI_Bcast(xywh,4,MPI_INT,0,PetscObjectComm((PetscObject)draw));CHKERRQ(ierr);
528   /* record new window position */
529   draw->x = win->x = xywh[0];
530   draw->y = win->y = xywh[1];
531   if (xywh[2] == win->w && xywh[3] == win->h) PetscFunctionReturn(0);
532   /* record new window sizes */
533   draw->w = win->w = xywh[2];
534   draw->h = win->h = xywh[3];
535 
536   ierr = PetscDrawCollectiveBegin(draw);CHKERRQ(ierr);
537   /* free buffer space and create new version (only first processor does this) */
538   if (!rank && win->drw) {ierr = PetscDrawXiQuickPixmap(win);CHKERRQ(ierr);}
539   /* reset the clipping */
540   xl = draw->port_xl; yl = draw->port_yl;
541   xr = draw->port_xr; yr = draw->port_yr;
542   ierr = PetscDrawSetViewport_X(draw,xl,xr,yl,yr);CHKERRQ(ierr);
543   /* try to make sure it is actually done before passing info to all */
544   XSync(win->disp,False);
545   ierr = PetscDrawCollectiveEnd(draw);CHKERRQ(ierr);
546 
547   ierr = MPI_Bcast(&win->drw,1,MPI_UNSIGNED_LONG,0,PetscObjectComm((PetscObject)draw));CHKERRQ(ierr);
548   PetscFunctionReturn(0);
549 }
550 
551 static PetscErrorCode PetscDrawGetSingleton_X(PetscDraw,PetscDraw*);
552 static PetscErrorCode PetscDrawRestoreSingleton_X(PetscDraw,PetscDraw*);
553 
554 #undef __FUNCT__
555 #define __FUNCT__ "PetscDrawDestroy_X"
556 static PetscErrorCode PetscDrawDestroy_X(PetscDraw draw)
557 {
558   PetscDraw_X    *win = (PetscDraw_X*)draw->data;
559   PetscErrorCode ierr;
560 
561   PetscFunctionBegin;
562   ierr = PetscDrawDestroy(&draw->popup);CHKERRQ(ierr);
563   ierr = PetscDrawXiClose(win);CHKERRQ(ierr);
564   ierr = PetscFree(draw->data);CHKERRQ(ierr);
565   PetscFunctionReturn(0);
566 }
567 
568 PETSC_INTERN PetscErrorCode PetscDrawSave_X(PetscDraw);
569 PETSC_INTERN PetscErrorCode PetscDrawSetSave_X(PetscDraw,const char[]);
570 
571 static struct _PetscDrawOps DvOps = { PetscDrawSetDoubleBuffer_X,
572                                       PetscDrawFlush_X,
573                                       PetscDrawLine_X,
574                                       0,
575                                       0,
576                                       PetscDrawPoint_X,
577                                       0,
578                                       PetscDrawString_X,
579                                       PetscDrawStringVertical_X,
580                                       PetscDrawStringSetSize_X,
581                                       PetscDrawStringGetSize_X,
582                                       PetscDrawSetViewport_X,
583                                       PetscDrawClear_X,
584                                       PetscDrawRectangle_X,
585                                       PetscDrawTriangle_X,
586                                       PetscDrawEllipse_X,
587                                       PetscDrawGetMouseButton_X,
588                                       PetscDrawPause_X,
589                                       0,
590                                       0,
591                                       PetscDrawGetPopup_X,
592                                       PetscDrawSetTitle_X,
593                                       PetscDrawCheckResizedWindow_X,
594                                       PetscDrawResizeWindow_X,
595                                       PetscDrawDestroy_X,
596                                       0,
597                                       PetscDrawGetSingleton_X,
598                                       PetscDrawRestoreSingleton_X,
599 #if defined(PETSC_HAVE_AFTERIMAGE)
600                                       PetscDrawSave_X,
601 #else
602                                       0,
603 #endif
604                                       PetscDrawSetSave_X,
605                                       0,
606                                       PetscDrawArrow_X,
607                                       PetscDrawCoordinateToPixel_X,
608                                       PetscDrawPixelToCoordinate_X,
609                                       PetscDrawPointPixel_X,
610                                       0};
611 
612 
613 #undef __FUNCT__
614 #define __FUNCT__ "PetscDrawGetSingleton_X"
615 static PetscErrorCode PetscDrawGetSingleton_X(PetscDraw draw,PetscDraw *sdraw)
616 {
617   PetscDraw_X    *Xwin = (PetscDraw_X*)draw->data,*sXwin;
618   PetscErrorCode ierr;
619 
620   PetscFunctionBegin;
621   ierr = PetscDrawCreate(PETSC_COMM_SELF,draw->display,draw->title,draw->x,draw->y,draw->w,draw->h,sdraw);CHKERRQ(ierr);
622   ierr = PetscObjectChangeTypeName((PetscObject)*sdraw,PETSC_DRAW_X);CHKERRQ(ierr);
623   ierr = PetscMemcpy((*sdraw)->ops,&DvOps,sizeof(DvOps));CHKERRQ(ierr);
624   (*sdraw)->ops->destroy  = NULL;
625   (*sdraw)->ops->getpopup = NULL;
626 
627   (*sdraw)->pause   = draw->pause;
628   (*sdraw)->coor_xl = draw->coor_xl;
629   (*sdraw)->coor_xr = draw->coor_xr;
630   (*sdraw)->coor_yl = draw->coor_yl;
631   (*sdraw)->coor_yr = draw->coor_yr;
632   (*sdraw)->port_xl = draw->port_xl;
633   (*sdraw)->port_xr = draw->port_xr;
634   (*sdraw)->port_yl = draw->port_yl;
635   (*sdraw)->port_yr = draw->port_yr;
636 
637   /* share drawables (windows and/or pixmap) from the parent draw */
638   ierr = PetscNewLog(*sdraw,&sXwin);CHKERRQ(ierr);
639   (*sdraw)->data = (void*)sXwin;
640   ierr = PetscDrawXiInit(sXwin,draw->display);CHKERRQ(ierr);
641   if (Xwin->win) {
642     ierr = PetscDrawXiQuickWindowFromWindow(sXwin,Xwin->win);CHKERRQ(ierr);
643     sXwin->drw = Xwin->drw; /* XXX If the window is ever resized, this is wrong! */
644   } else if (Xwin->drw) {
645     ierr = PetscDrawXiColormap(sXwin);CHKERRQ(ierr);
646     sXwin->drw = Xwin->drw;
647   }
648   (*sdraw)->x = sXwin->x; (*sdraw)->y = sXwin->y;
649   (*sdraw)->w = sXwin->w; (*sdraw)->h = sXwin->h;
650   PetscFunctionReturn(0);
651 }
652 
653 #undef __FUNCT__
654 #define __FUNCT__ "PetscDrawRestoreSingleton_X"
655 static PetscErrorCode PetscDrawRestoreSingleton_X(PetscDraw draw,PetscDraw *sdraw)
656 {
657   PetscDraw_X    *sXwin = (PetscDraw_X*)(*sdraw)->data;
658   PetscErrorCode ierr;
659 
660   PetscFunctionBegin;
661   ierr = PetscFree((*sdraw)->display);CHKERRQ(ierr);
662   ierr = PetscFree((*sdraw)->title);CHKERRQ(ierr);
663   ierr = PetscDrawXiClose(sXwin);CHKERRQ(ierr);
664   ierr = PetscFree(sXwin);CHKERRQ(ierr);
665   ierr = PetscHeaderDestroy(sdraw);CHKERRQ(ierr);
666   PetscFunctionReturn(0);
667 }
668 
669 #undef __FUNCT__
670 #define __FUNCT__ "PetscDrawXGetDisplaySize_Private"
671 static PetscErrorCode PetscDrawXGetDisplaySize_Private(const char name[],int *width,int *height)
672 {
673   Display *display;
674 
675   PetscFunctionBegin;
676   display = XOpenDisplay(name);
677   if (!display) {
678     *width  = *height = 0;
679     SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"Unable to open display on %s\n\
680     Make sure your COMPUTE NODES are authorized to connect\n\
681     to this X server and either your DISPLAY variable\n\
682     is set or you use the -display name option\n",name);
683   }
684   *width  = (int)DisplayWidth(display,DefaultScreen(display));
685   *height = (int)DisplayHeight(display,DefaultScreen(display));
686   XCloseDisplay(display);
687   PetscFunctionReturn(0);
688 }
689 
690 /*MC
691      PETSC_DRAW_X  - PETSc graphics device that uses either X windows or its virtual version Xvfb
692 
693    Options Database Keys:
694 +  -display <display> - sets the display to use
695 .  -x_virtual - forces use of a X virtual display Xvfb that will not display anything but -draw_save will still work.
696                 Xvfb is automatically started up in PetscSetDisplay() with this option
697 .  -draw_size w,h - percentage of screeen (either 1, .5, .3, .25), or size in pixels
698 .  -geometry x,y,w,h - set location and size in pixels
699 .  -draw_virtual - do not open a window (draw on a pixmap), -draw_save will still work
700 -  -draw_double_buffer - avoid window flickering (draw on pixmap and flush to window)
701 
702    Level: beginner
703 
704 .seealso:  PetscDrawOpenX(), PetscDrawSetDisplay(), PetscDrawSetFromOptions()
705 
706 M*/
707 
708 #undef __FUNCT__
709 #define __FUNCT__ "PetscDrawCreate_X"
710 PETSC_EXTERN PetscErrorCode PetscDrawCreate_X(PetscDraw draw)
711 {
712   PetscDraw_X    *Xwin;
713   PetscErrorCode ierr;
714   PetscMPIInt    rank;
715   int            x = draw->x,y = draw->y,w = draw->w,h = draw->h;
716   static int     xavailable = 0,yavailable = 0,ybottom = 0,xmax = 0,ymax = 0;
717   PetscBool      set,dvirtual = PETSC_FALSE,doublebuffer = PETSC_TRUE;
718   PetscInt       xywh[4],osize = 4,nsizes=2;
719   PetscReal      sizes[2] = {.3,.3};
720 
721   PetscFunctionBegin;
722   /* get the display variable */
723   if (!draw->display) {
724     ierr = PetscMalloc1(256,&draw->display);CHKERRQ(ierr);
725     ierr = PetscGetDisplay(draw->display,256);CHKERRQ(ierr);
726   }
727 
728   /* initialize the display size */
729   if (!xmax) {
730     ierr = PetscDrawXGetDisplaySize_Private(draw->display,&xmax,&ymax);
731     /* if some processors fail on this and others succed then this is a problem ! */
732     if (ierr) {
733       (*PetscErrorPrintf)("PETSc unable to use X windows\nproceeding without graphics\n");
734       ierr = PetscDrawSetType(draw,PETSC_DRAW_NULL);CHKERRQ(ierr);
735       PetscFunctionReturn(0);
736     }
737   }
738 
739   /* allow user to set size of drawable */
740   ierr = PetscOptionsGetRealArray(((PetscObject)draw)->options,((PetscObject)draw)->prefix,"-draw_size",sizes,&nsizes,&set);CHKERRQ(ierr);
741   if (set && nsizes == 1 && sizes[0] > 1.0) sizes[1] = sizes[0];
742   if (set) {
743     if (sizes[0] > 1.0)       w = (int)sizes[0];
744     else if (sizes[0] == 1.0) w = PETSC_DRAW_FULL_SIZE;
745     else if (sizes[0] == .5)  w = PETSC_DRAW_HALF_SIZE;
746     else if (sizes[0] == .3)  w = PETSC_DRAW_THIRD_SIZE;
747     else if (sizes[0] == .25) w = PETSC_DRAW_QUARTER_SIZE;
748     if (sizes[1] > 1.0)       h = (int)sizes[1];
749     else if (sizes[1] == 1.0) h = PETSC_DRAW_FULL_SIZE;
750     else if (sizes[1] == .5)  h = PETSC_DRAW_HALF_SIZE;
751     else if (sizes[1] == .3)  h = PETSC_DRAW_THIRD_SIZE;
752     else if (sizes[1] == .25) h = PETSC_DRAW_QUARTER_SIZE;
753   }
754   if (w == PETSC_DECIDE || w == PETSC_DEFAULT) w = draw->w = 300;
755   if (h == PETSC_DECIDE || h == PETSC_DEFAULT) h = draw->h = 300;
756   switch (w) {
757   case PETSC_DRAW_FULL_SIZE:    w = draw->w = (xmax - 10);   break;
758   case PETSC_DRAW_HALF_SIZE:    w = draw->w = (xmax - 20)/2; break;
759   case PETSC_DRAW_THIRD_SIZE:   w = draw->w = (xmax - 30)/3; break;
760   case PETSC_DRAW_QUARTER_SIZE: w = draw->w = (xmax - 40)/4; break;
761   }
762   switch (h) {
763   case PETSC_DRAW_FULL_SIZE:    h = draw->h = (ymax - 10);   break;
764   case PETSC_DRAW_HALF_SIZE:    h = draw->h = (ymax - 20)/2; break;
765   case PETSC_DRAW_THIRD_SIZE:   h = draw->h = (ymax - 30)/3; break;
766   case PETSC_DRAW_QUARTER_SIZE: h = draw->h = (ymax - 40)/4; break;
767   }
768 
769   ierr = PetscOptionsGetBool(((PetscObject)draw)->options,((PetscObject)draw)->prefix,"-draw_virtual",&dvirtual,NULL);CHKERRQ(ierr);
770 
771   if (!dvirtual) {
772 
773     /* allow user to set location and size of window */
774     xywh[0] = x; xywh[1] = y; xywh[2] = w; xywh[3] = h;
775     ierr = PetscOptionsGetIntArray(((PetscObject)draw)->options,((PetscObject)draw)->prefix,"-geometry",xywh,&osize,NULL);CHKERRQ(ierr);
776     x = (int)xywh[0]; y = (int)xywh[1]; w = (int)xywh[2]; h = (int)xywh[3];
777     if (w == PETSC_DECIDE || w == PETSC_DEFAULT) w = 300;
778     if (h == PETSC_DECIDE || h == PETSC_DEFAULT) h = 300;
779     draw->x = x; draw->y = y; draw->w = w; draw->h = h;
780 
781     if (draw->x == PETSC_DECIDE || draw->y == PETSC_DECIDE) {
782       /*
783        PETSc tries to place windows starting in the upper left corner
784         and moving across to the right.
785 
786        0,0+----------------------------------------------+
787        |  Region used so far  +xavailable,yavailable  |
788        |                      |                       |
789        |                      |                       |
790        +--------------------- +ybottom                |
791        |                                              |
792        |                                              |
793        +----------------------------------------------+xmax,ymax
794 
795        */
796       /*  First: can we add it to the right? */
797       if (xavailable+w+10 <= xmax) {
798         x       = xavailable;
799         y       = yavailable;
800         ybottom = PetscMax(ybottom,y + h + 30);
801       } else {
802         /* No, so add it below on the left */
803         x          = 0;
804         xavailable = 0;
805         y          = ybottom;
806         yavailable = ybottom;
807         ybottom    = ybottom + h + 30;
808       }
809     }
810     /* update available region */
811     xavailable = PetscMax(xavailable,x + w + 10);
812     if (xavailable >= xmax) {
813       xavailable = 0;
814       yavailable = yavailable + h + 30;
815       ybottom    = yavailable;
816     }
817     if (yavailable >= ymax) {
818       y          = 0;
819       yavailable = 0;
820       ybottom    = 0;
821     }
822 
823   } /* endif(!dvirtual) */
824 
825   ierr = MPI_Comm_rank(PetscObjectComm((PetscObject)draw),&rank);CHKERRQ(ierr);
826   if (!rank && (w <= 0 || h <= 0)) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Negative window width or height");
827 
828   ierr = PetscNewLog(draw,&Xwin);CHKERRQ(ierr);
829   ierr = PetscMemcpy(draw->ops,&DvOps,sizeof(DvOps));CHKERRQ(ierr);
830   draw->data = (void*)Xwin;
831 
832   ierr = PetscDrawXiInit(Xwin,draw->display);CHKERRQ(ierr);
833   if (!dvirtual) {
834     Xwin->x = x; Xwin->y = y;
835     Xwin->w = w; Xwin->h = h;
836     if (!rank) {ierr = PetscDrawXiQuickWindow(Xwin,draw->title,x,y,w,h);CHKERRQ(ierr);}
837     ierr = MPI_Bcast(&Xwin->win,1,MPI_UNSIGNED_LONG,0,PetscObjectComm((PetscObject)draw));CHKERRQ(ierr);
838     if (rank) {ierr = PetscDrawXiQuickWindowFromWindow(Xwin,Xwin->win);CHKERRQ(ierr);}
839   } else {
840     Xwin->x = 0; Xwin->y = 0;
841     Xwin->w = w; Xwin->h = h;
842     ierr = PetscDrawXiColormap(Xwin);CHKERRQ(ierr);
843     if (!rank) {ierr = PetscDrawXiQuickPixmap(Xwin);CHKERRQ(ierr);}
844     ierr = MPI_Bcast(&Xwin->drw,1,MPI_UNSIGNED_LONG,0,PetscObjectComm((PetscObject)draw));CHKERRQ(ierr);
845   }
846   draw->x = Xwin->x; draw->y = Xwin->y;
847   draw->w = Xwin->w; draw->h = Xwin->h;
848 
849   ierr = PetscOptionsGetBool(((PetscObject)draw)->options,((PetscObject)draw)->prefix,"-draw_double_buffer",&doublebuffer,NULL);CHKERRQ(ierr);
850   if (doublebuffer) {ierr = PetscDrawSetDoubleBuffer(draw);CHKERRQ(ierr);}
851   PetscFunctionReturn(0);
852 }
853