ROOT  6.07/01
Reference Guide
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
TPainter3dAlgorithms.h
Go to the documentation of this file.
1 // @(#)root/histpainter:$Id$
2 // Author: Rene Brun, Evgueni Tcherniaev, Olivier Couet 12/12/94
3 
4 /*************************************************************************
5  * Copyright (C) 1995-2000, Rene Brun and Fons Rademakers. *
6  * All rights reserved. *
7  * *
8  * For the licensing terms see $ROOTSYS/LICENSE. *
9  * For the list of contributors see $ROOTSYS/README/CREDITS. *
10  *************************************************************************/
11 
12 #ifndef ROOT_TPainter3dAlgorithms
13 #define ROOT_TPainter3dAlgorithms
14 
15 
16 //////////////////////////////////////////////////////////////////////////
17 // //
18 // TPainter3dAlgorithms //
19 // //
20 // 3D graphics representations package. //
21 // //
22 //////////////////////////////////////////////////////////////////////////
23 
24 #ifndef ROOT_TObject
25 #include "TObject.h"
26 #endif
27 
28 #ifndef ROOT_TAttLine
29 #include "TAttLine.h"
30 #endif
31 
32 #ifndef ROOT_TAttFill
33 #include "TAttFill.h"
34 #endif
35 
36 const Int_t kCARTESIAN = 1;
37 const Int_t kPOLAR = 2;
38 const Int_t kCYLINDRICAL = 3;
39 const Int_t kSPHERICAL = 4;
40 const Int_t kRAPIDITY = 5;
41 
42 class TF3;
43 
44 class TPainter3dAlgorithms : public TObject, public TAttLine, public TAttFill {
45 
46 private:
49  Double_t fRmin[3]; //Lower limits of lego
50  Double_t fRmax[3]; //Upper limits of lego
51  Double_t fU[2000]; //
52  Double_t fD[2000]; //
53  Double_t fT[200]; //
54  Double_t fFunLevel[257]; //Function levels corresponding to colors
55  Double_t fPlines[1200]; //
56  Double_t fAphi[183]; //
58  Double_t fYls[4]; //
59  Double_t fVls[12]; //
67  Int_t fSystem; //Coordinate system
68  Int_t fNT; //
69  Int_t fNlevel; //Number of color levels
70  Int_t fColorLevel[258]; //Color levels corresponding to functions
79  Int_t fMesh; //(=1 if mesh to draw, o otherwise)
81  Int_t fLevelLine[200]; //
83  Int_t fNqs; //
84  Int_t fNStack; //Number of histograms in the stack to be painted
88  Int_t *fRaster; //pointer to raster buffer
89  Int_t fJmask[30]; //
90  Int_t fMask[465]; //
91  Double_t fP8[8][3]; //
92  Double_t fF8[8]; //
93  Double_t fG8[8][3]; //
94  Double_t fFmin; // IsoSurface minimum function value
95  Double_t fFmax; // IsoSurface maximum function value
96  Int_t fNcolor; // Number of colours per Iso surface
97  Int_t fIc1; // Base colour for the 1st Iso Surface
98  Int_t fIc2; // Base colour for the 2nd Iso Surface
99  Int_t fIc3; // Base colour for the 3rd Iso Surface
100 
101  static Int_t fgF3Clipping; // Clipping box is off (0) or on (1)
102  static Double_t fgF3XClip; // Clipping plne along X
103  static Double_t fgF3YClip; // Clipping plne along Y
104  static Double_t fgF3ZClip; // Clipping plne along Y
105  static TF3 *fgCurrentF3; // Pointer to the 3D function to be paint.
106 
107 
108 public:
112 
113 private:
114  DrawFaceFunc_t fDrawFace; //pointer to face drawing function
115  LegoFunc_t fLegoFunction; //pointer to lego function
116  SurfaceFunc_t fSurfaceFunction; //pointer to surface function
117 
118 public:
120  TPainter3dAlgorithms(Double_t *rmin, Double_t *rmax, Int_t system=1);
121  virtual ~TPainter3dAlgorithms();
122  void BackBox(Double_t ang);
123  void ClearRaster();
124  void ColorFunction(Int_t nl, Double_t *fl, Int_t *icl, Int_t &irep);
125  void DefineGridLevels(Int_t ndivz);
126  void DrawFaceGouraudShaded(Int_t *icodes, Double_t xyz[][3], Int_t np, Int_t *iface, Double_t *t);
127  void DrawFaceMode1(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *t);
128  void DrawFaceMode2(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *t);
129  void DrawFaceMode3(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *t);
130  void DrawFaceMove1(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *tt);
131  void DrawFaceMove2(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *tt);
132  void DrawFaceMove3(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *tt);
133  void DrawFaceRaster1(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *tt);
134  void DrawFaceRaster2(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *tt);
135  void FillPolygon(Int_t n, Double_t *p, Double_t *f);
136  void FillPolygonBorder(Int_t nn, Double_t *xy);
138  void FindPartEdge(Double_t *p1, Double_t *p2, Double_t f1, Double_t f2, Double_t fmin, Double_t fmax, Int_t &kpp, Double_t *pp);
139  void FindVisibleLine(Double_t *p1, Double_t *p2, Int_t ntmax, Int_t &nt, Double_t *t);
141  void FrontBox(Double_t ang);
142  void GouraudFunction(Int_t ia, Int_t ib, Double_t *f, Double_t *t);
143  void ImplicitFunction(Double_t *rmin, Double_t *rmax, Int_t nx, Int_t ny, Int_t nz, const char *chopt);
144  void IsoSurface (Int_t ns, Double_t *s, Int_t nx, Int_t ny, Int_t nz, Double_t *x, Double_t *y, Double_t *z, const char *chopt);
147  void LegoCartesian(Double_t ang, Int_t nx, Int_t ny, const char *chopt);
148  void LegoFunction(Int_t ia, Int_t ib, Int_t &nv, Double_t *ab, Double_t *vv, Double_t *t);
149  void LegoPolar(Int_t iordr, Int_t na, Int_t nb, const char *chopt);
150  void LegoCylindrical(Int_t iordr, Int_t na, Int_t nb, const char *chopt);
151  void LegoSpherical(Int_t ipsdr, Int_t iordr, Int_t na, Int_t nb, const char *chopt);
152  void LightSource(Int_t nl, Double_t yl, Double_t xscr, Double_t yscr, Double_t zscr, Int_t &irep);
153  void Luminosity(Double_t *anorm, Double_t &flum);
154  void MarchingCube(Double_t fiso, Double_t p[8][3], Double_t f[8], Double_t g[8][3], Int_t &nnod, Int_t &ntria, Double_t xyz[][3], Double_t grad[][3], Int_t itria[][3]);
155  void MarchingCubeCase00(Int_t k1, Int_t k2, Int_t k3, Int_t k4, Int_t k5, Int_t k6, Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3]);
156  void MarchingCubeCase03(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3]);
157  void MarchingCubeCase04(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3]);
158  void MarchingCubeCase06(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3]);
159  void MarchingCubeCase07(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3]);
160  void MarchingCubeCase10(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3]);
161  void MarchingCubeCase12(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3]);
162  void MarchingCubeCase13(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3]);
163  void MarchingCubeSetTriangles(Int_t ntria, Int_t it[][3], Int_t itria[48][3]);
164  void MarchingCubeMiddlePoint(Int_t nnod, Double_t xyz[52][3], Double_t grad[52][3], Int_t it[][3], Double_t *pxyz, Double_t *pgrad);
166  void MarchingCubeFindNodes(Int_t nnod, Int_t *ie, Double_t xyz[52][3], Double_t grad[52][3]);
168  void SetDrawFace(DrawFaceFunc_t pointer);
169  void SetIsoSurfaceParameters(Double_t fmin, Double_t fmax, Int_t ncolor, Int_t ic1, Int_t ic2, Int_t ic3){fFmin=fmin; fFmax=fmax; fNcolor=ncolor; fIc1=ic1; fIc2=ic2; fIc3=ic3;}
170  void SetLegoFunction(LegoFunc_t pointer);
171  void SetMesh(Int_t mesh=1) {fMesh=mesh;}
172  void SetSurfaceFunction(SurfaceFunc_t pointer);
173  void SetColorDark(Color_t color, Int_t n=0);
174  void SetColorMain(Color_t color, Int_t n=0);
175  void SetEdgeAtt(Color_t color=1, Style_t style=1, Width_t width=1, Int_t n=0);
176  void SideVisibilityDecode(Double_t val, Int_t &iv1, Int_t &iv2, Int_t &iv3, Int_t &iv4, Int_t &iv5, Int_t &iv6, Int_t &ir);
177  void SideVisibilityEncode(Int_t iopt, Double_t phi1, Double_t phi2, Double_t &val);
178  void Spectrum(Int_t nl, Double_t fmin, Double_t fmax, Int_t ic, Int_t idc, Int_t &irep);
179  void SurfaceCartesian(Double_t ang, Int_t nx, Int_t ny, const char *chopt);
180  void SurfacePolar(Int_t iordr, Int_t na, Int_t nb, const char *chopt);
181  void SurfaceCylindrical(Int_t iordr, Int_t na, Int_t nb, const char *chopt);
182  void SurfaceFunction(Int_t ia, Int_t ib, Double_t *f, Double_t *t);
183  void SurfaceSpherical(Int_t ipsdr, Int_t iordr, Int_t na, Int_t nb, const char *chopt);
184  void SurfaceProperty(Double_t qqa, Double_t qqd, Double_t qqs, Int_t nnqs, Int_t &irep);
185  void TestEdge(Double_t del, Double_t xyz[52][3], Int_t i1, Int_t i2, Int_t iface[3], Double_t abcd[4], Int_t &irep);
186  void ZDepth(Double_t xyz[52][3], Int_t &nface, Int_t iface[48][3], Double_t dface[48][6], Double_t abcd[48][4], Int_t *iorder);
187 
188  static void SetF3(TF3 *f3);
189  static void SetF3ClippingBoxOff();
190  static void SetF3ClippingBoxOn(Double_t xclip, Double_t yclip, Double_t zclip);
191 
192  ClassDef(TPainter3dAlgorithms,0) //Hidden line removal package
193 };
194 
195 #endif
196 
197 
198 
const int nx
Definition: kalman.C:16
void SideVisibilityEncode(Int_t iopt, Double_t phi1, Double_t phi2, Double_t &val)
Encode side visibilities and order along R for sector.
void(TPainter3dAlgorithms::* LegoFunc_t)(Int_t, Int_t, Int_t &, Double_t *, Double_t *, Double_t *)
void SurfaceProperty(Double_t qqa, Double_t qqd, Double_t qqs, Int_t nnqs, Int_t &irep)
Set surface property coefficients.
float xmin
Definition: THbookFile.cxx:93
void grad()
Definition: grad.C:17
void MarchingCubeCase12(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3])
Consider case No 12.
void MarchingCubeCase00(Int_t k1, Int_t k2, Int_t k3, Int_t k4, Int_t k5, Int_t k6, Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3])
Consideration of trivial cases: 1,2,5,8,9,11,14.
void LegoCartesian(Double_t ang, Int_t nx, Int_t ny, const char *chopt)
Draw stack of lego-plots in cartesian coordinates.
void MarchingCubeCase06(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3])
Consider case No 6.
void MarchingCubeSetTriangles(Int_t ntria, Int_t it[][3], Int_t itria[48][3])
Set triangles (if parameter IALL=1, all edges will be visible)
short Style_t
Definition: RtypesCore.h:76
void SetColorDark(Color_t color, Int_t n=0)
Store dark color for stack number n.
void SideVisibilityDecode(Double_t val, Int_t &iv1, Int_t &iv2, Int_t &iv3, Int_t &iv4, Int_t &iv5, Int_t &iv6, Int_t &ir)
Decode side visibilities and order along R for sector.
const double a11[11]
Definition: RooMath.cxx:491
void SetIsoSurfaceParameters(Double_t fmin, Double_t fmax, Int_t ncolor, Int_t ic1, Int_t ic2, Int_t ic3)
float ymin
Definition: THbookFile.cxx:93
void MarchingCubeFindNodes(Int_t nnod, Int_t *ie, Double_t xyz[52][3], Double_t grad[52][3])
Find nodes and normales.
tuple f2
Definition: surfaces.py:24
void f3()
Definition: na49.C:50
virtual ~TPainter3dAlgorithms()
Lego default destructor.
void InitRaster(Double_t xmin, Double_t ymin, Double_t xmax, Double_t ymax, Int_t nx, Int_t ny)
Initialize hidden lines removal algorithm (RASTER SCREEN)
void DrawFaceMove3(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *tt)
Draw face - 3rd variant for "MOVING SCREEN" algorithm (draw level lines only)
void SurfaceFunction(Int_t ia, Int_t ib, Double_t *f, Double_t *t)
Service function for Surfaces.
void ColorFunction(Int_t nl, Double_t *fl, Int_t *icl, Int_t &irep)
Set correspondance between function and color levels.
void MarchingCubeCase03(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3])
Consider case No 3.
int Int_t
Definition: RtypesCore.h:41
void FindVisibleDraw(Double_t *r1, Double_t *r2)
Find visible parts of line (draw line)
void MarchingCubeCase10(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3])
Consider case No 10.
void InitMoveScreen(Double_t xmin, Double_t xmax)
Initialize "MOVING SCREEN" method.
TFile * f
static void SetF3(TF3 *f3)
Static function Store pointer to current implicit function.
void DefineGridLevels(Int_t ndivz)
Define the grid levels drawn in the background of surface and lego plots.
void LegoPolar(Int_t iordr, Int_t na, Int_t nb, const char *chopt)
Draw stack of lego-plots in polar coordinates.
void TestEdge(Double_t del, Double_t xyz[52][3], Int_t i1, Int_t i2, Int_t iface[3], Double_t abcd[4], Int_t &irep)
Test edge against face (triangle)
Fill Area Attributes class.
Definition: TAttFill.h:32
Double_t x[n]
Definition: legend1.C:17
void DrawFaceMode3(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *t)
Draw face - 3rd option (draw face for stacked lego plot)
void DrawFaceRaster1(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *tt)
Draw face - 1st variant for "RASTER SCREEN" algorithm (draw face with level lines) ...
TText * tt
Definition: textangle.C:16
void LightSource(Int_t nl, Double_t yl, Double_t xscr, Double_t yscr, Double_t zscr, Int_t &irep)
Set light source.
const int ny
Definition: kalman.C:17
void SetMesh(Int_t mesh=1)
static double p2(double t, double a, double b, double c)
void DrawFaceMode2(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *t)
Draw face - 2nd option (fill in correspondance with function levels)
tuple np
Definition: multifit.py:30
ClassDef(TAttLine, 2)
void ImplicitFunction(Double_t *rmin, Double_t *rmax, Int_t nx, Int_t ny, Int_t nz, const char *chopt)
Draw implicit function FUN(X,Y,Z) = 0 in cartesian coordinates using hidden surface removal algorithm...
void DrawFaceMode1(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *t)
Draw face - 1st variant (2 colors: 1st for external surface, 2nd for internal)
void LegoCylindrical(Int_t iordr, Int_t na, Int_t nb, const char *chopt)
Draw stack of lego-plots in cylindrical coordinates.
Float_t z[5]
Definition: Ifit.C:16
void DrawFaceMove1(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *tt)
Draw face - 1st variant for "MOVING SCREEN" algorithm (draw face with level lines) ...
void ClearRaster()
Clear screen.
const Int_t kSPHERICAL
void DrawFaceMove2(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *tt)
Draw face - 2nd variant for "MOVING SCREEN" algorithm (draw face for stacked lego plot) ...
void MarchingCubeCase13(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3])
Consider case No 13.
short Color_t
Definition: RtypesCore.h:79
void FillPolygon(Int_t n, Double_t *p, Double_t *f)
Fill polygon with function values at vertexes.
void SetDrawFace(DrawFaceFunc_t pointer)
Store pointer to current algorithm to draw faces.
void(TPainter3dAlgorithms::* SurfaceFunc_t)(Int_t, Int_t, Double_t *, Double_t *)
TThread * t[5]
Definition: threadsh1.C:13
float ymax
Definition: THbookFile.cxx:93
void SetColorMain(Color_t color, Int_t n=0)
Store color for stack number n.
void IsoSurface(Int_t ns, Double_t *s, Int_t nx, Int_t ny, Int_t nz, Double_t *x, Double_t *y, Double_t *z, const char *chopt)
Draw set of isosurfaces for a scalar function defined on a grid.
void FrontBox(Double_t ang)
Draw forward faces of surrounding box & axes.
XPoint xy[kMAXMK]
Definition: TGX11.cxx:122
void FindLevelLines(Int_t np, Double_t *f, Double_t *t)
Find level lines for face.
const Int_t kPOLAR
A 3-Dim function with parameters.
Definition: TF3.h:30
void LegoSpherical(Int_t ipsdr, Int_t iordr, Int_t na, Int_t nb, const char *chopt)
Draw stack of lego-plots spheric coordinates.
void LegoFunction(Int_t ia, Int_t ib, Int_t &nv, Double_t *ab, Double_t *vv, Double_t *t)
Service function for Legos.
unsigned int r1[N_CITIES]
Definition: simanTSP.cxx:321
void SurfaceCylindrical(Int_t iordr, Int_t na, Int_t nb, const char *chopt)
Draw surface in cylindrical coordinates.
void MarchingCubeCase07(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3])
Consider case No 7.
void Luminosity(Double_t *anorm, Double_t &flum)
Find surface luminosity at given point.
float xmax
Definition: THbookFile.cxx:93
static double p1(double t, double a, double b)
void SetLegoFunction(LegoFunc_t pointer)
Store pointer to current lego function.
short Width_t
Definition: RtypesCore.h:78
void SurfaceSpherical(Int_t ipsdr, Int_t iordr, Int_t na, Int_t nb, const char *chopt)
Draw surface in spheric coordinates.
void MarchingCubeCase04(Int_t &nnod, Int_t &ntria, Double_t xyz[52][3], Double_t grad[52][3], Int_t itria[48][3])
Consider case No 4.
void MarchingCubeMiddlePoint(Int_t nnod, Double_t xyz[52][3], Double_t grad[52][3], Int_t it[][3], Double_t *pxyz, Double_t *pgrad)
Find middle point of a polygon.
const Int_t kCYLINDRICAL
void(TPainter3dAlgorithms::* DrawFaceFunc_t)(Int_t *, Double_t *, Int_t, Int_t *, Double_t *)
void DrawFaceGouraudShaded(Int_t *icodes, Double_t xyz[][3], Int_t np, Int_t *iface, Double_t *t)
Draw the faces for the Gouraud Shaded Iso surfaces.
double Double_t
Definition: RtypesCore.h:55
The Legos and Surfaces painter class.
void ZDepth(Double_t xyz[52][3], Int_t &nface, Int_t iface[48][3], Double_t dface[48][6], Double_t abcd[48][4], Int_t *iorder)
Z-depth algorithm for set of triangles.
TCanvas * style()
Definition: style.C:1
Double_t y[n]
Definition: legend1.C:17
static void SetF3ClippingBoxOn(Double_t xclip, Double_t yclip, Double_t zclip)
Static function Set the implicit function clipping box "on" and define the clipping box...
void SetEdgeAtt(Color_t color=1, Style_t style=1, Width_t width=1, Int_t n=0)
void Spectrum(Int_t nl, Double_t fmin, Double_t fmax, Int_t ic, Int_t idc, Int_t &irep)
Set Spectrum.
void MarchingCube(Double_t fiso, Double_t p[8][3], Double_t f[8], Double_t g[8][3], Int_t &nnod, Int_t &ntria, Double_t xyz[][3], Double_t grad[][3], Int_t itria[][3])
Topological decider for "Marching Cubes" algorithm Find set of triangles aproximating the isosurface ...
void BackBox(Double_t ang)
Draw back surfaces of surrounding box.
Mother of all ROOT objects.
Definition: TObject.h:58
typedef void((*Func_t)())
void FindVisibleLine(Double_t *p1, Double_t *p2, Int_t ntmax, Int_t &nt, Double_t *t)
Find visible part of a line ("RASTER SCREEN")
void FindPartEdge(Double_t *p1, Double_t *p2, Double_t f1, Double_t f2, Double_t fmin, Double_t fmax, Int_t &kpp, Double_t *pp)
Find part of edge where function defined on this edge has value from fmin to fmax ...
void MarchingCubeSurfacePenetration(Double_t a00, Double_t a10, Double_t a11, Double_t a01, Double_t b00, Double_t b10, Double_t b11, Double_t b01, Int_t &irep)
Check for surface penetration ("bottle neck")
static void SetF3ClippingBoxOff()
Static function Set the implicit function clipping box "off".
void DrawFaceRaster2(Int_t *icodes, Double_t *xyz, Int_t np, Int_t *iface, Double_t *tt)
Draw face - 2nd variant for "RASTER SCREEN" algorithm (draw face for stacked lego plot) ...
TF1 * f1
Definition: legend1.C:11
const Int_t kRAPIDITY
void ModifyScreen(Double_t *r1, Double_t *r2)
Modify SCREEN.
void SurfacePolar(Int_t iordr, Int_t na, Int_t nb, const char *chopt)
Draw surface in polar coordinates.
void SurfaceCartesian(Double_t ang, Int_t nx, Int_t ny, const char *chopt)
Draw surface in cartesian coordinate system.
void SetSurfaceFunction(SurfaceFunc_t pointer)
Store pointer to current surface function.
const Int_t kCARTESIAN
void GouraudFunction(Int_t ia, Int_t ib, Double_t *f, Double_t *t)
Find part of surface with luminosity in the corners.
const Int_t n
Definition: legend1.C:16
Line Attributes class.
Definition: TAttLine.h:32
void FillPolygonBorder(Int_t nn, Double_t *xy)
Fill a polygon including border ("RASTER SCREEN")
unsigned int r2[N_CITIES]
Definition: simanTSP.cxx:322