xref: /petsc/src/sys/classes/draw/impls/x/xops.c (revision 387ef4f7680df1f160fdf549863d1c7ac10fb432)
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                                       PetscDrawSave_X,
600                                       PetscDrawSetSave_X,
601                                       0,
602                                       PetscDrawArrow_X,
603                                       PetscDrawCoordinateToPixel_X,
604                                       PetscDrawPixelToCoordinate_X,
605                                       PetscDrawPointPixel_X,
606                                       0};
607 
608 
609 #undef __FUNCT__
610 #define __FUNCT__ "PetscDrawGetSingleton_X"
611 static PetscErrorCode PetscDrawGetSingleton_X(PetscDraw draw,PetscDraw *sdraw)
612 {
613   PetscDraw_X    *Xwin = (PetscDraw_X*)draw->data,*sXwin;
614   PetscErrorCode ierr;
615 
616   PetscFunctionBegin;
617   ierr = PetscDrawCreate(PETSC_COMM_SELF,draw->display,draw->title,draw->x,draw->y,draw->w,draw->h,sdraw);CHKERRQ(ierr);
618   ierr = PetscObjectChangeTypeName((PetscObject)*sdraw,PETSC_DRAW_X);CHKERRQ(ierr);
619   ierr = PetscMemcpy((*sdraw)->ops,&DvOps,sizeof(DvOps));CHKERRQ(ierr);
620 
621   if (draw->popup) {
622     PetscBool isnull;
623     ierr = PetscDrawIsNull(draw->popup,&isnull);CHKERRQ(ierr);
624     if (isnull) {ierr = PetscDrawOpenNull(PETSC_COMM_SELF,&(*sdraw)->popup);CHKERRQ(ierr);}
625     else        {ierr = PetscDrawGetSingleton(draw->popup,&(*sdraw)->popup);CHKERRQ(ierr);}
626   }
627 
628   (*sdraw)->pause   = draw->pause;
629   (*sdraw)->coor_xl = draw->coor_xl;
630   (*sdraw)->coor_xr = draw->coor_xr;
631   (*sdraw)->coor_yl = draw->coor_yl;
632   (*sdraw)->coor_yr = draw->coor_yr;
633   (*sdraw)->port_xl = draw->port_xl;
634   (*sdraw)->port_xr = draw->port_xr;
635   (*sdraw)->port_yl = draw->port_yl;
636   (*sdraw)->port_yr = draw->port_yr;
637 
638   /* share drawables (windows and/or pixmap) from the parent draw */
639   ierr = PetscNewLog(*sdraw,&sXwin);CHKERRQ(ierr);
640   (*sdraw)->data = (void*)sXwin;
641   ierr = PetscDrawXiInit(sXwin,draw->display);CHKERRQ(ierr);
642   if (Xwin->win) {
643     ierr = PetscDrawXiQuickWindowFromWindow(sXwin,Xwin->win);CHKERRQ(ierr);
644     sXwin->drw = Xwin->drw; /* XXX If the window is ever resized, this is wrong! */
645   } else if (Xwin->drw) {
646     ierr = PetscDrawXiColormap(sXwin);CHKERRQ(ierr);
647     sXwin->drw = Xwin->drw;
648   }
649   (*sdraw)->x = sXwin->x; (*sdraw)->y = sXwin->y;
650   (*sdraw)->w = sXwin->w; (*sdraw)->h = sXwin->h;
651   PetscFunctionReturn(0);
652 }
653 
654 #undef __FUNCT__
655 #define __FUNCT__ "PetscDrawRestoreSingleton_X"
656 static PetscErrorCode PetscDrawRestoreSingleton_X(PetscDraw draw,PetscDraw *sdraw)
657 {
658   PetscErrorCode ierr;
659 
660   PetscFunctionBegin;
661   if (draw->popup && (*sdraw)->popup) {
662     PetscBool isdrawx;
663     PetscDraw_X *pXwin = (PetscDraw_X*)draw->popup->data;
664     PetscDraw_X *sXwin = (PetscDraw_X*)(*sdraw)->popup->data;
665     ierr = PetscObjectTypeCompare((PetscObject)draw->popup,PETSC_DRAW_X,&isdrawx);CHKERRQ(ierr);
666     if (!isdrawx) goto finally;
667     ierr = PetscObjectTypeCompare((PetscObject)(*sdraw)->popup,PETSC_DRAW_X,&isdrawx);CHKERRQ(ierr);
668     if (!isdrawx) goto finally;
669     if (sXwin->win == pXwin->win) {
670       ierr = PetscDrawRestoreSingleton(draw->popup,&(*sdraw)->popup);CHKERRQ(ierr);
671     }
672   }
673 finally:
674   ierr = PetscDrawDestroy(sdraw);CHKERRQ(ierr);
675   PetscFunctionReturn(0);
676 }
677 
678 #undef __FUNCT__
679 #define __FUNCT__ "PetscDrawXGetDisplaySize_Private"
680 static PetscErrorCode PetscDrawXGetDisplaySize_Private(const char name[],int *width,int *height)
681 {
682   Display *display;
683 
684   PetscFunctionBegin;
685   display = XOpenDisplay(name);
686   if (!display) {
687     *width  = *height = 0;
688     SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_LIB,"Unable to open display on %s\n\
689     Make sure your COMPUTE NODES are authorized to connect\n\
690     to this X server and either your DISPLAY variable\n\
691     is set or you use the -display name option\n",name);
692   }
693   *width  = (int)DisplayWidth(display,DefaultScreen(display));
694   *height = (int)DisplayHeight(display,DefaultScreen(display));
695   XCloseDisplay(display);
696   PetscFunctionReturn(0);
697 }
698 
699 /*MC
700      PETSC_DRAW_X  - PETSc graphics device that uses either X windows or its virtual version Xvfb
701 
702    Options Database Keys:
703 +  -display <display> - sets the display to use
704 .  -x_virtual - forces use of a X virtual display Xvfb that will not display anything but -draw_save will still work.
705                 Xvfb is automatically started up in PetscSetDisplay() with this option
706 .  -draw_size w,h - percentage of screeen (either 1, .5, .3, .25), or size in pixels
707 .  -geometry x,y,w,h - set location and size in pixels
708 .  -draw_virtual - do not open a window (draw on a pixmap), -draw_save will still work
709 -  -draw_double_buffer - avoid window flickering (draw on pixmap and flush to window)
710 
711    Level: beginner
712 
713 .seealso:  PetscDrawOpenX(), PetscDrawSetDisplay(), PetscDrawSetFromOptions()
714 
715 M*/
716 
717 #undef __FUNCT__
718 #define __FUNCT__ "PetscDrawCreate_X"
719 PETSC_EXTERN PetscErrorCode PetscDrawCreate_X(PetscDraw draw)
720 {
721   PetscDraw_X    *Xwin;
722   PetscErrorCode ierr;
723   PetscMPIInt    rank;
724   int            x = draw->x,y = draw->y,w = draw->w,h = draw->h;
725   static int     xavailable = 0,yavailable = 0,ybottom = 0,xmax = 0,ymax = 0;
726   PetscBool      set,dvirtual = PETSC_FALSE,doublebuffer = PETSC_TRUE;
727   PetscInt       xywh[4],osize = 4,nsizes=2;
728   PetscReal      sizes[2] = {.3,.3};
729 
730   PetscFunctionBegin;
731   /* get the display variable */
732   if (!draw->display) {
733     ierr = PetscMalloc1(256,&draw->display);CHKERRQ(ierr);
734     ierr = PetscGetDisplay(draw->display,256);CHKERRQ(ierr);
735   }
736 
737   /* initialize the display size */
738   if (!xmax) {
739     ierr = PetscDrawXGetDisplaySize_Private(draw->display,&xmax,&ymax);
740     /* if some processors fail on this and others succed then this is a problem ! */
741     if (ierr) {
742       (*PetscErrorPrintf)("PETSc unable to use X windows\nproceeding without graphics\n");
743       ierr = PetscDrawSetType(draw,PETSC_DRAW_NULL);CHKERRQ(ierr);
744       PetscFunctionReturn(0);
745     }
746   }
747 
748   /* allow user to set size of drawable */
749   ierr = PetscOptionsGetRealArray(((PetscObject)draw)->options,((PetscObject)draw)->prefix,"-draw_size",sizes,&nsizes,&set);CHKERRQ(ierr);
750   if (set && nsizes == 1 && sizes[0] > 1.0) sizes[1] = sizes[0];
751   if (set) {
752     if (sizes[0] > 1.0)       w = (int)sizes[0];
753     else if (sizes[0] == 1.0) w = PETSC_DRAW_FULL_SIZE;
754     else if (sizes[0] == .5)  w = PETSC_DRAW_HALF_SIZE;
755     else if (sizes[0] == .3)  w = PETSC_DRAW_THIRD_SIZE;
756     else if (sizes[0] == .25) w = PETSC_DRAW_QUARTER_SIZE;
757     if (sizes[1] > 1.0)       h = (int)sizes[1];
758     else if (sizes[1] == 1.0) h = PETSC_DRAW_FULL_SIZE;
759     else if (sizes[1] == .5)  h = PETSC_DRAW_HALF_SIZE;
760     else if (sizes[1] == .3)  h = PETSC_DRAW_THIRD_SIZE;
761     else if (sizes[1] == .25) h = PETSC_DRAW_QUARTER_SIZE;
762   }
763   if (w == PETSC_DECIDE || w == PETSC_DEFAULT) w = draw->w = 300;
764   if (h == PETSC_DECIDE || h == PETSC_DEFAULT) h = draw->h = 300;
765   switch (w) {
766   case PETSC_DRAW_FULL_SIZE:    w = draw->w = (xmax - 10);   break;
767   case PETSC_DRAW_HALF_SIZE:    w = draw->w = (xmax - 20)/2; break;
768   case PETSC_DRAW_THIRD_SIZE:   w = draw->w = (xmax - 30)/3; break;
769   case PETSC_DRAW_QUARTER_SIZE: w = draw->w = (xmax - 40)/4; break;
770   }
771   switch (h) {
772   case PETSC_DRAW_FULL_SIZE:    h = draw->h = (ymax - 10);   break;
773   case PETSC_DRAW_HALF_SIZE:    h = draw->h = (ymax - 20)/2; break;
774   case PETSC_DRAW_THIRD_SIZE:   h = draw->h = (ymax - 30)/3; break;
775   case PETSC_DRAW_QUARTER_SIZE: h = draw->h = (ymax - 40)/4; break;
776   }
777 
778   ierr = PetscOptionsGetBool(((PetscObject)draw)->options,((PetscObject)draw)->prefix,"-draw_virtual",&dvirtual,NULL);CHKERRQ(ierr);
779 
780   if (!dvirtual) {
781 
782     /* allow user to set location and size of window */
783     xywh[0] = x; xywh[1] = y; xywh[2] = w; xywh[3] = h;
784     ierr = PetscOptionsGetIntArray(((PetscObject)draw)->options,((PetscObject)draw)->prefix,"-geometry",xywh,&osize,NULL);CHKERRQ(ierr);
785     x = (int)xywh[0]; y = (int)xywh[1]; w = (int)xywh[2]; h = (int)xywh[3];
786     if (w == PETSC_DECIDE || w == PETSC_DEFAULT) w = 300;
787     if (h == PETSC_DECIDE || h == PETSC_DEFAULT) h = 300;
788     draw->x = x; draw->y = y; draw->w = w; draw->h = h;
789 
790     if (draw->x == PETSC_DECIDE || draw->y == PETSC_DECIDE) {
791       /*
792        PETSc tries to place windows starting in the upper left corner
793         and moving across to the right.
794 
795        0,0+----------------------------------------------+
796        |  Region used so far  +xavailable,yavailable  |
797        |                      |                       |
798        |                      |                       |
799        +--------------------- +ybottom                |
800        |                                              |
801        |                                              |
802        +----------------------------------------------+xmax,ymax
803 
804        */
805       /*  First: can we add it to the right? */
806       if (xavailable+w+10 <= xmax) {
807         x       = xavailable;
808         y       = yavailable;
809         ybottom = PetscMax(ybottom,y + h + 30);
810       } else {
811         /* No, so add it below on the left */
812         x          = 0;
813         xavailable = 0;
814         y          = ybottom;
815         yavailable = ybottom;
816         ybottom    = ybottom + h + 30;
817       }
818     }
819     /* update available region */
820     xavailable = PetscMax(xavailable,x + w + 10);
821     if (xavailable >= xmax) {
822       xavailable = 0;
823       yavailable = yavailable + h + 30;
824       ybottom    = yavailable;
825     }
826     if (yavailable >= ymax) {
827       y          = 0;
828       yavailable = 0;
829       ybottom    = 0;
830     }
831 
832   } /* endif(!dvirtual) */
833 
834   ierr = MPI_Comm_rank(PetscObjectComm((PetscObject)draw),&rank);CHKERRQ(ierr);
835   if (!rank && (w <= 0 || h <= 0)) SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Negative window width or height");
836 
837   ierr = PetscNewLog(draw,&Xwin);CHKERRQ(ierr);
838   ierr = PetscMemcpy(draw->ops,&DvOps,sizeof(DvOps));CHKERRQ(ierr);
839   draw->data = (void*)Xwin;
840 
841   ierr = PetscDrawXiInit(Xwin,draw->display);CHKERRQ(ierr);
842   if (!dvirtual) {
843     Xwin->x = x; Xwin->y = y;
844     Xwin->w = w; Xwin->h = h;
845     if (!rank) {ierr = PetscDrawXiQuickWindow(Xwin,draw->title,x,y,w,h);CHKERRQ(ierr);}
846     ierr = MPI_Bcast(&Xwin->win,1,MPI_UNSIGNED_LONG,0,PetscObjectComm((PetscObject)draw));CHKERRQ(ierr);
847     if (rank) {ierr = PetscDrawXiQuickWindowFromWindow(Xwin,Xwin->win);CHKERRQ(ierr);}
848   } else {
849     Xwin->x = 0; Xwin->y = 0;
850     Xwin->w = w; Xwin->h = h;
851     ierr = PetscDrawXiColormap(Xwin);CHKERRQ(ierr);
852     if (!rank) {ierr = PetscDrawXiQuickPixmap(Xwin);CHKERRQ(ierr);}
853     ierr = MPI_Bcast(&Xwin->drw,1,MPI_UNSIGNED_LONG,0,PetscObjectComm((PetscObject)draw));CHKERRQ(ierr);
854   }
855   draw->x = Xwin->x; draw->y = Xwin->y;
856   draw->w = Xwin->w; draw->h = Xwin->h;
857 
858   ierr = PetscOptionsGetBool(((PetscObject)draw)->options,((PetscObject)draw)->prefix,"-draw_double_buffer",&doublebuffer,NULL);CHKERRQ(ierr);
859   if (doublebuffer) {ierr = PetscDrawSetDoubleBuffer(draw);CHKERRQ(ierr);}
860   PetscFunctionReturn(0);
861 }
862