source: Sophya/trunk/SophyaPI/PIext/pihisto2d.cc@ 3572

Last change on this file since 3572 was 3559, checked in by ansari, 17 years ago

Amelioration gestion taille de fonte ds PIHisto, PINTuple ... en utilisant le flag AutoFontSize mis ds la classe PIDrawer (attributs autofontsize=, fixedfontsize), Reza 24/12/2008
File size: 34.0 KB
Line 
1#include <stdio.h>
2#include <stdlib.h>
3#include <pisysdep.h>
4#include PIAPP_H
5
6#include <iostream>
7#include <math.h>
8
9#include "sopnamsp.h"
10#include "pihisto2d.h"
11#include "srandgen.h"
12
13static int dbg = 0;
14
15
16//------ Implementation classe P2DHistoWrapper
17P2DHistoWrapper::P2DHistoWrapper(int_4 asx, int_4 asy)
18 : P2DArrayAdapter(asx, asy) ,
19 mScale(1.) , mOff(0.) , mRetFg(0)
20{
21}
22
23P2DHistoWrapper::~P2DHistoWrapper()
24{
25}
26
27
28double P2DHistoWrapper::MeanVal(int ix1, int ix2, int jy1, int jy2)
29{
30 int ec;
31 if (ix1>ix2) { ec=ix1; ix1=ix2; ix2=ec; }
32 if (jy1>jy2) { ec=jy1; jy1=jy2; jy2=ec; }
33 double ss = 0.;
34 for(int j=jy1; j<=jy2; j++)
35 for(int i=ix1; i<=ix2; i++) ss += (double)((*this)(i, j));
36 ss /= (double)((jy2-jy1+1)*(ix2-ix1+1));
37 return ss;
38}
39
40int P2DHistoWrapper::GetStatInfoAsText(vector<string> & /* text */ )
41{
42 return 0;
43}
44
45int P2DHistoWrapper::DecodeOptionString(vector<string> & opt, bool rmdecopt)
46{
47 if(opt.size() < 1) return(0);
48
49 vector<string> udopt; // On gardera ici les options non decodees
50 unsigned int k = 0;
51 int ndec = opt.size();
52 for( k=0; k<opt.size(); k++ ) {
53 string & opts = opt[k];
54 if(opts.substr(0,7) == "hs1") {
55 mScale = 1.; mOff = 0.;
56 }
57 else if(opts.substr(0,7) == "hscale=") {
58 mScale = atof(opts.substr(7).c_str());
59 }
60 else if(opts.substr(0,8) == "hoffset=") {
61 mOff = atof(opts.substr(8).c_str());
62 }
63 else if(opts.substr(0,8) == "hbincont") {
64 mRetFg = 0;
65 }
66 else if(opts.substr(0,7) == "hbinerr") {
67 mRetFg = 1;
68 }
69 else if(opts.substr(0,7) == "hbinent") {
70 mRetFg = 2;
71 }
72 else {
73 ndec--;
74 // S'il faut supprimer les options decodees
75 if (rmdecopt) udopt.push_back(opts);
76 }
77 }
78 // S'il faut supprimer les options decodees, on remplace l'argument opt
79 // par le vecteur des options non decodees.
80 if (rmdecopt) opt = udopt;
81 return(ndec);
82}
83
84int
85P2DHistoWrapper::OptionToString(vector<string> & opt) const
86{
87 char buff[64];
88 sprintf(buff, "hscale=%g", mScale); opt.push_back(buff);
89 sprintf(buff, "hoffset=%g", mOff); opt.push_back(buff);
90 if (mRetFg == 2) opt.push_back("hbinent");
91 else if (mRetFg == 1) opt.push_back("hbinerr");
92 else opt.push_back("hbincont");
93
94 return 1;
95}
96
97//++
98// Class PIHisto2D
99// Lib PIext
100// include pihisto2d.h
101//
102// Classes de dessin des histogrammes a 2 dimensions pour
103// objets *Histo2D*
104//--
105
106//++
107// Links Parents
108// PIDrawer
109//--
110
111//++
112// Titre Constructeur, méthodes
113//--
114
115//++
116PIHisto2D::PIHisto2D(P2DHistoWrapper* histowp, bool ad)
117//
118// Createur d'une classe de dessin pour l'histogramme 2D histo.
119//--
120: PIDrawer(), mHistoWp(histowp), mAdDO(ad), mLogScale(10.), mFPoints(0.5)
121{
122// mAdDO : Flag pour suppression automatique de mHistoWp
123// Attention: mFPoints n'est initialise que si on display par nuages de points
124// mLogScale n'est initialise que si on utilise une echelle log
125UseScale();
126UseColors();
127UseDisplay();
128UseDyn();
129UseFrac();
130SetStats();
131SetStatPosOffset();
132SetName("Histo2DDrw");
133// PIHisto2D has specific control tools
134mFgSpecContWind = true;
135}
136
137//++
138PIHisto2D::~PIHisto2D()
139//
140// Destructeur.
141//--
142{
143// La fonction virtuelle DeactivateControlWindow() doit etre appele
144// a ce niveau - En effet au niveau du destructeur de base, il
145// semble pointer sur la fonction de la classe de base
146// Reza - Octobre 2002
147// Desactivation totale de la fenetre de controle specialise
148// ---> parametre d'appel PIBaseWdgGen* wdg=NULL
149DeactivateControlWindow(NULL);
150if(mAdDO && mHistoWp!=NULL) delete mHistoWp;
151}
152
153//++
154void PIHisto2D::UseColors(bool fg,CMapId cmap,bool revcmap)
155//
156// Choix de la couleur si fg=true avec la color map cmap.
157// (pour la couleur cmap cf picmap.h).
158// Independamment du choix du display, la dynamique est
159// codee sur la color map donnant ainsi une double
160// information. Par exemple, carres de tailles variables
161// en couleur. Cette option est incontournable dans le cas
162// d'un display par des carres de taille fixe.
163// revcmap doit etre mis a "true" si on veut avoir une color map
164// inversee.
165//| -**- gestion dans H2WinArg par menu deroulant Black&White etc...
166//--
167{
168mFgCol = fg; mCmap = cmap; mRevCmap = revcmap;
169}
170
171//++
172void PIHisto2D::UseScale(unsigned short type,double logscale)
173//
174// Pour changer les echelles (lineaire ou logarithmique)
175//| Type = 0 : echelle lineaire
176//| = 1 : echelle log10
177//| -**- Explication du codage en type=0 (lineaire) :
178//| 1. [hmin,hmax] -> [0,1]
179//| h -> f = (h-hmin)/(hmax-hmin)
180//| 2. codage de f=[0,1] sur la dynamique du display choisi
181//| -**- Explication du codage en type=1 (logarithmique base 10) :
182//| 1. map lineaire entre 0 et 1:
183//| [hmin,hmax] -> [0,1]
184//| h -> f = (h-hmin)/(hmax-hmin)
185//| 2. transformation logarithmique de base 10 :
186//| [0,1] -> [0,1]
187//| f -> lf = log10(1.+f*(logscale-1))/log10(logscale)
188//| 3. codage de lf=[0,1] sur la dynamique du display choisi
189//| -**- gestion dans H2WinArg par menu deroulant Lineaire/Log10
190//| et "logscale" par saisie de valeur dans champ LogScal
191//--
192{
193if(type==0) mTypScal=0;
194else if(type==1) {mTypScal=1; if(logscale>1.) mLogScale=logscale;}
195else mTypScal=0;
196}
197
198//++
199void PIHisto2D::UseDisplay(unsigned short type,double fnpt)
200//
201// Type de Display
202//| Type = 0 : carres de tailles variables
203//| Type = 1 : nuages de points
204//| Le nombre de points a utiliser est fnpt*N
205//| ou N est le nombre de pixels ecran contenu
206//| dans un bin de l'histogramme.
207//| Type = 2 : code a la "hbook2" " .+123...9AB...YZ*"
208//| (cf detail PIHisto2D::HPrint2)
209//| Type = 3 : carres de taille fixe (couleur).
210//| -**- gestion dans H2WinArg par menu deroulant Carres_Var etc...
211//| et "fnpt" par saisie de valeur dans champ PerPt
212//--
213{
214if(type==0) mTypDisp = 0;
215else if(type==1) {
216 mTypDisp = 1;
217 if(fnpt<0.) mFPoints = 0.;
218 else if(fnpt>1.) mFPoints = 1.;
219 else mFPoints = fnpt;
220}
221else if(type==2) mTypDisp = 2;
222else if(type==3) mTypDisp = 3;
223else mTypDisp = 1;
224}
225
226//++
227void PIHisto2D::UseDyn(double fmin,double fmax)
228//
229// Gestion de la dynamique a representer:
230//| La dynamique va etre transformee de [hmin,hmax] vers [0,1] selon
231//| [hmin,hmax] -> [0,1]
232//| h -> f = (h-hmin)/(hmax-hmin)
233//| Par la suite, selon ce qui est demande, f va coder le display
234//| ou etre transforme en une autre echelle [0,1] (ex: echelle log10).
235//| Si hmax<=hmin, ils sont forces a la dynamique totale de l'histo2D.
236//| -**- gestion dans H2WinArg par saisie de valeurs dans champ Dyn
237//--
238{
239if( (mHistoWp) && (mHistoWp->NBinX()>0) && (mHistoWp->NBinY()>0) && (fmin>=fmax) ) {
240 double v,hmin,hmax;
241 hmin = hmax = (*mHistoWp)(0,0);
242 for (int_4 i=0; i<mHistoWp->NBinX(); i++)
243 for (int_4 j=0; j<mHistoWp->NBinY(); j++) {
244 v = (*mHistoWp)(i, j);
245 if(v<hmin) hmin = v;
246 if(v>hmax) hmax = v;
247 }
248 fmin = hmin; fmax = hmax;
249}
250
251if(fmin>=fmax) fmax = fmin+1.;
252mHMin = fmin; mHMax = fmax;
253}
254
255//++
256void PIHisto2D::UseFrac(double frmin,double frmax)
257//
258// Pour definir la fraction de la dynamique a dessiner:
259//| Selon le type de display (f=[0,1] cf PIHisto2D::UseDyn),
260//| - on ne dessine rien si f <= frmin dans les cas de display avec
261//| des nuages de points ou des carres de tailles variables.
262//| Pour un display "a la hbook2" on force frmin = 0.
263//| - frmax n'est utilise que pour la representation avec
264//| des carres de tailles variables: c'est la taille
265//| maximum que peut avoir le carre exprimee en unite
266//| de la taille du bin (ex: si frmax=0.8 le carre
267//| le + grand qui pourra etre dessine dans un bin
268//| aura une taille egale a 0.8*(taille du bin)).
269//| -**- gestion dans H2WinArg par saisie de valeurs dans champ Frac
270//--
271{
272if(frmin<0. || frmin>=1.) frmin = 0.;
273if(frmax<=0. || frmax>1. ) frmax = 1.;
274if(frmin>=frmax) {frmin=0.1; frmax=0.9;}
275mFracMin = frmin; mFracMax = frmax;
276}
277
278//++
279void PIHisto2D::Print(int lp)
280//
281// Print de l'etat des options du display.
282//--
283{
284printf("PIHisto2D::Print FgCol=%d Cmap=%d (Rev=%d) TypScal=%d TypDisp=%d (FPoints=%g)\n"
285 ,(int)mFgCol,(int)mCmap,(int)mRevCmap,mTypScal,mTypDisp,mFPoints);
286printf(" Dyn=%g,%g Frac=%g,%g LogSc=%g H=%p\n"
287 ,mHMin,mHMax,mFracMin,mFracMax,mLogScale,mHistoWp);
288fflush(stdout);
289}
290
291//++
292void PIHisto2D::UpdateLimits()
293//
294// Definition des tailles graphiques en fonction
295// des caracteristiques de l'histogramme a dessiner.
296//--
297{
298 if(!mHistoWp) return;
299 SetLimits(mHistoWp->XMin(), mHistoWp->XMax(), mHistoWp->YMin() , mHistoWp->YMax());
300}
301
302//++
303void PIHisto2D::ShowControlWindow(PIBaseWdgGen* wdg)
304//
305// Affichage de la fenetre de controle H2WinArg
306//--
307{
308 H2WinArg::SetCurrentPIHisto2D(this);
309 H2WinArg::ShowPIHisto2DTools(wdg);
310}
311
312//++
313void PIHisto2D::DeactivateControlWindow(PIBaseWdgGen* wdg)
314//
315// Desactivation de la fenetre de controle specialisee
316//--
317{
318 if (H2WinArg::GetCurrentPIHisto2D() == this) {
319 // si wdg != NULL, c'est un Detach (Drawer detache du PIBaseWdg
320 // si wdg == NULL, c'est un delete du PIHisto2D (du PIDrawer)
321 if ((wdg == NULL) || (H2WinArg::GetCurrentBaseWdg() == wdg)) {
322 H2WinArg::SetCurrentBaseWdg(NULL);
323 H2WinArg::SetCurrentPIHisto2D(NULL);
324 H2WinArg::HidePIHisto2DTools();
325 }
326 }
327 PIDrawer::DeactivateControlWindow(wdg);
328 return;
329}
330
331//++
332void PIHisto2D::Draw(PIGraphicUC* g,double xmin,double ymin,double xmax,double ymax)
333//
334// Dessin de l'histogramme.
335//| -**- Code de dessin selon choix des options:
336//| (detail voir UseColors UseScale UseDisplay UseDyn UseFrac)
337//| - [hmin,hmax] -> f=[0,1]
338//| (Choix hmin,hmax champ Dyn de H2WinArg)
339//| - Eventuellement ech Log -> re-codage log10 entre f=[0,1]
340//| (Choix menu deroulant et champ LogScal de H2WinArg)
341//| - Restriction de f=[0,1] -> f=[Frac(min),Frac(max)]
342//| (Choix champ Frac de H2WinArg)
343//| -**- Puis selon display:
344//| 0 carres variables, menu "Carres Var." de H2WinArg:
345//| if(f>Frac(min)) taille carre = f * Frac(max) * taille_du_bin
346//| 1 nuage de points, menu "....." et champ PerPt de H2WinArg:
347//| if(f>Frac(min)) npoints = f * PerPt * npoints_ecran_dans_bin
348//| 2 code hbook2, menu ".12..Z*" de H2WinArg:
349//| if(f>0) map de f=]0,1] dans ".+...Z*"
350//| 3 carres pleins, menu "Carres Pleins" et couleurs de H2WinArg):
351//| couleur = lut[ f * nombre_d_entree_dans_la_lut ]
352//--
353{
354if (axesFlags != kAxesNone) DrawAxes(g);
355
356if(!mHistoWp) return;
357if(mHistoWp->NBinX()<=0 || mHistoWp->NBinY()<=0) return;
358
359// Caracteristiques histogramme
360double dx = mHistoWp->WBinX(),dy = mHistoWp->WBinY();
361double p1dx,p1dy;
362g->DGrC2UC(1.,1.,p1dx,p1dy);
363
364// Gamme a representer entre [0,1] mais >=fracmin et scale fracmax
365double fracmin=FMin(), fracmax=FMax();
366double llscale = log10(LogScale());
367
368// gestion Couleurs.
369PIColors fgcoul = GetGraphicAtt().GetFgColor();
370PIColors bgcoul = GetGraphicAtt().GetBgColor();
371PIColorMap* cmap=NULL;
372int ncol = 0;
373if (mFgCol) {
374 cmap = new PIColorMap(mCmap);
375 cmap->ReverseColorIndex(mRevCmap);
376 ncol = cmap->NCol();
377 if(mTypDisp==3) fracmin=-1.;
378}
379
380// gestion epaisseur de ligne
381PILineAtt LineAtt = GetGraphicAtt().GetLineAtt();
382if(LineAtt == PI_NotDefLineAtt) GetGraphicAtt().SetLineAtt(PI_ThinLine);
383
384// gestion Markers ou plot avec des points.
385PIMarker Mk = GetGraphicAtt().GetMarker();
386int MkSz = GetGraphicAtt().GetMarkerSize();
387int npt = 1;
388if(mTypDisp==1) {
389 g->SelMarker(1,PI_DotMarker);
390 npt = (int)(NPixBin(g)*FPoints()); if(npt<=0) npt = 2;
391}
392
393// gestion Font.
394PIFontAtt FontAtt = GetGraphicAtt().GetFontAtt();
395int FontSize = GetGraphicAtt().GetFontSzPt();
396if(mTypDisp==2) {
397 double dxg,dyg,dg;
398 g->DUC2GrC(dx,dy,dxg,dyg);
399 dg =(dxg<dyg) ? dxg : dyg;
400 int npix = (int) (dg*0.9); if(npix<8) npix = 8;
401 g->SelFontSzPt(npix,FontAtt);
402 fracmin = 0;
403}
404
405// Plot de l'histogramme
406for(int i=0; i<mHistoWp->NBinX(); i++)
407 for(int j=0; j<mHistoWp->NBinY(); j++) {
408
409 r_8 left0,bottom0;
410 mHistoWp->BinLowEdge(i,j,left0,bottom0);
411
412 // Gestion de la dynamique a dessiner
413 double frac = ((*mHistoWp)(i,j)-HMin())/(HMax()-HMin());
414 if(frac<0.) continue;
415 if(mTypScal==1) { // echelle log10
416 frac = log10(1.+frac*(LogScale()-1.))/llscale;
417 if(frac<0.) continue;
418 }
419 if(frac<=fracmin) continue;
420 if(frac>1.) frac = 1.;
421 double fracred = frac * fracmax;
422
423 // Gestion de la couleur
424 int icol = 0;
425 if (cmap) {
426 icol = int(ncol*frac);
427 if(icol>=ncol) icol = ncol-1; else if(icol<0) icol=0;
428 g->SelForeground(*cmap,icol);
429 }
430
431// Pour ne pas dessiner en dehors des axes
432 if ( (left0+dx/2. < xmin) || (left0+dx/2. > xmax) ||
433 (bottom0+dy/2. < ymin) || (bottom0+dy/2. > ymax) ) continue;
434
435 // Dessin proprement dit selon le choix graphique.
436 if(mTypDisp==0) {
437 //..... carres de tailles variables
438 double left = left0 + 0.5*(1.-fracred)*dx, width = fracred*dx;
439 double bottom = bottom0 + 0.5*(1.-fracred)*dy, height = fracred*dy;
440 if (cmap) g->DrawFBox(left,bottom,width,height);
441 else g->DrawBox(left,bottom,width,height);
442 } else if(mTypDisp==1) {
443 //..... nuage de points .....
444 int ipt = int(npt *frac);
445 for(int k=0;k<ipt;k++) {
446 double x = left0 + frand01()*dx;
447 double y = bottom0 + frand01()*dy;
448 g->DrawMarker(x,y);
449 }
450 } else if(mTypDisp==2) {
451 //..... type hbook2/hprint .+23-Z*
452 char c[2];
453 c[0] = HPrint2(frac); c[1]='\0';
454 double x = left0 + dx/2.;
455 double y = bottom0 + dy/2.;
456 g->DrawString(x,y,c,PI_HorizontalCenter|PI_VerticalCenter);
457 } else if(mTypDisp==3) {
458 //..... carres de tailles fixes (avec gestion de continuite)
459 if (cmap) g->DrawFBox(left0,bottom0,dx+p1dx,dy+p1dy);
460 else g->DrawBox(left0,bottom0,dx+p1dx,dy+p1dy);
461 }
462
463}
464
465// Remise dans les conditions ulterieures pour la suite du graphique.
466GetGraphicAtt().SetMarkerAtt(MkSz,Mk);
467GetGraphicAtt().SetColAtt(fgcoul,bgcoul);
468g->SelFontSzPt(FontSize,FontAtt);
469GetGraphicAtt().SetLineAtt(LineAtt);
470if (cmap) delete cmap;
471
472// Fin du dessin, ecriture de la statistique.
473if(stats) DrawStats(g);
474}
475
476//++
477void PIHisto2D::GetClickInfo(string& info,double x,double y,double x0,double y0,bool fgdiff)
478//
479// Info specifique du drawer pour la position x,y
480//--
481{
482P2DHistoWrapper* h = HistoWrapper();
483if(h == NULL) return;
484
485int_4 i,j;
486i = floor((x-h->XMin())/h->WBinX());
487j = floor((y-h->YMin())/h->WBinY());
488
489if(i>=0 && i<h->NBinX() && j>=0 && j<h->NBinY()) {
490 char str[64];
491 if(fgdiff) {
492 int i0,j0;
493 i0 = floor((x0-h->XMin())/h->WBinX());
494 j0 = floor((y0-h->YMin())/h->WBinY());
495 if(i0>=0 && i0<h->NBinX() && j0>=0 && j0<h->NBinY()) {
496 sprintf(str," DV=%g",(*h)(i,j)-(*h)(i0,j0));
497 info += str;
498 } else {
499 info += " DV=?";
500 }
501 }
502 sprintf(str," V=%g",(*h)(i,j));
503 info += str;
504} else {
505 info += " V=?";
506}
507
508return;
509}
510
511//++
512void PIHisto2D::DrawStats(PIGraphicUC* g)
513//
514// Dessin des informations statistiques de l'histogramme.
515//--
516{
517 if (!mHistoWp) return;
518 if (GetGraphicAtt().GetLineAtt() == PI_NotDefLineAtt) g->SelLine(PI_ThinLine);
519// Reza(12/2008): utilisation de flag et fraction pour selection automatique taille de fonte
520 if (mAFSz) {
521 double fsz = (YMax() - YMin())*mFontSzF;
522 g->SelFontSz(fsz); // au lieu de 1/30*(ymax-ymin)
523 }
524
525 // La hauteur de la cellule
526 PIGrCoord a, d;
527 double cH = (double)g->GetFontHeight(a,d);
528 vector<string> lines;
529 int nlig = mHistoWp->GetStatInfoAsText(lines);
530 if (nlig < 1) return;
531
532 double cellHeight = nlig*1.2 * cH;
533
534 int idxll = 0;
535
536 int kl;
537 // on recherche la ligne la plus longue
538 for(kl=1; kl<nlig; kl++)
539 if ( lines[kl].length() > lines[idxll].length() ) idxll = kl;
540
541 double cellWidth = 1.1 * (double)g->CalcStringWidth(lines[idxll].c_str());
542
543
544 double ofpx = spoX*(XMax()-XMin());
545 double ofpy = spoY*(YMax()-YMin());
546
547 double xu, yu, cw;
548 // Les limites du cadre
549 xu = g->DeltaUCX(XMax(), -cellWidth);
550 yu = g->DeltaUCY(YMax(), -cellHeight);
551 double recw = XMax()-xu;
552 double rech = YMax()-yu;
553 xu += ofpx; yu += ofpy;
554 g->DrawBox(xu, yu, recw, rech);
555
556 // L'ecriture des labels
557 cw = (g->isAxeXDirRtoL()) ? -0.05*cellWidth : -0.95*cellWidth;
558 xu = g->DeltaUCX(XMax(),cw);
559
560 cw = (g->isAxeYDirUpDown()) ? -0.15*cH : -1.15*cH;
561 yu = g->DeltaUCY(YMax(),cw);
562 xu += ofpx; yu += ofpy;
563
564 for(kl=0; kl<nlig; kl++) {
565 g->DrawString(xu, yu, lines[kl].c_str() );
566 cw += -1.15*cH;
567 yu = g->DeltaUCY(YMax(),cw); yu += ofpy;
568 }
569
570
571 // printf("H[%d,%d] Dynamique: [%g,%g] Frac [%g,%g]\n"
572 // ,mHistoWp->NBinX(),mHistoWp->NBinY(),HMin(),HMax(),FMin(),FMax());
573}
574
575//++
576char PIHisto2D::HPrint2(double f)
577//
578// Codage des valeurs en caracteres (fct privee).
579//| f entre [0,1] mappee entre valeur=[0,37]
580//| si <0 alors =0, si >1 alors 1
581//| Display 4 ==> 4<=valeur<5
582//| C ==> 12<=valeur<13
583//| ==> valeur<=0
584//| * ==> valeur>=1
585//| . ==> 0<valeur<1
586//|------------------------------------------
587//| 1 2 3
588//| 01234567890123456789012345678901234567
589//| .+23456789ABCDEFGHIJKLMNOPQRSTUVWXYZ*
590//|------------------------------------------
591//--
592{
593char str[39] = " .+23456789ABCDEFGHIJKLMNOPQRSTUVWXYZ*";
594int i;
595if(f<=0.) i = 0;
596else if(f>=1.) i = 37;
597else { i = (int) (f*36.); i++;}
598if(i<0) i=0; else if (i>=38) i = 37;
599return str[i];
600}
601
602//++
603int PIHisto2D::NPixBin(PIGraphicUC* g)
604//
605// Nombre de pixels ecran dans un bin d'histogramme
606// (fct privee).
607//--
608{
609double dx = mHistoWp->WBinX(),dy = mHistoWp->WBinY();
610double dxg,dyg;
611g->DUC2GrC(dx,dy,dxg,dyg);
612int np = (int) dxg * (int) dyg;
613//printf("PIHisto2D::NPixBin H dx=%g dy=%g, G dx=%g dy=%g, np = %d\n"
614// ,dx,dy,dxg,dyg,np);
615return np;
616}
617
618//++
619int PIHisto2D::DecodeOptionString(vector<string> & opt, bool rmdecopt)
620//
621// Decodage des options
622//--
623{
624 // Decodage des options generales pidrawer
625 int optsz1 = opt.size();
626 if(optsz1<1) return(0);
627 int ndec1 = PIDrawer::DecodeOptionString(opt, rmdecopt);
628 if(optsz1-ndec1<1) return(ndec1); // si tout a ete decode
629
630 // Options generales pidrawer interessant le display Histo2D
631 bool rev;
632 if(GetGraphicAtt().GetColMapId(rev) != CMAP_OTHER) {
633 UseColors(true,GetGraphicAtt().GetColMapId(),rev);
634 } else UseColors(false);
635
636 // Decodage des options propres au display Histo2D
637 vector<string> udopt; // On gardera ici les options non decodees
638 unsigned int k = 0;
639 int ndec = opt.size();
640 bool listopt=false;
641 for( k=0; k<opt.size(); k++ ) {
642 string opts = opt[k];
643 if(opts=="h2help") {
644 string info; GetOptionsHelpInfo(info);
645 size_t q = info.find("PIHisto2D");
646 if(q<info.length()-1) cout<<info.substr(q)<<endl;
647 } else if(opts=="h2list") {
648 listopt=true;
649 } else if(opts=="sta" || opts=="stat" || opts=="stats") {
650 SetStats(true);
651 } else if( opts=="nsta" || opts=="nstat"
652 || opts=="nostat" || opts=="nostats") {
653 SetStats(false);
654 } else if(opts.substr(0,11) == "statposoff=") {
655 double xo=0., yo=0.;
656 sscanf(opts.substr(11).c_str(),"%lf,%lf",&xo, &yo);
657 SetStatPosOffset(xo, yo);
658 } else if(opts.substr(0,8)=="h2scale=") {
659 unsigned short t=TypScale(); double ls=LogScale();
660 if(opts.substr(8,3)=="lin") t=0;
661 else if(opts.substr(8,3)=="log")
662 {t=1; sscanf(opts.c_str(),"h2scale=log,%lf",&ls);}
663 UseScale(t,ls);
664 } else if(opts.substr(0,7)=="h2disp=") {
665 unsigned short t=TypDisplay(); double fpts=FPoints();
666 if(opts.substr(7,3)=="var") t=0;
667 else if(opts.substr(7,3)=="hbk") t=2;
668 else if(opts.substr(7,3)=="img") t=3;
669 else if(opts.substr(7,3)=="pts")
670 {t=1; sscanf(opts.c_str(),"h2disp=pts,%lf",&fpts);}
671 UseDisplay(t,fpts);
672 } else if(opts.substr(0,6)=="h2dyn=") {
673 double hmin=HMin(),hmax=HMax(); size_t q = opts.find(',');
674 sscanf(opts.c_str(),"h2dyn=%lf",&hmin);
675 if(q<opts.length()-1) sscanf(opts.substr(q+1).c_str(),"%lf",&hmax);
676 UseDyn(hmin,hmax);
677 } else if(opts.substr(0,7)=="h2frac=") {
678 double fmin=FMin(),fmax=FMax(); size_t q = opts.find(',');
679 sscanf(opts.c_str(),"h2frac=%lf",&fmin);
680 if(q<opts.length()-1) sscanf(opts.substr(q+1).c_str(),"%lf",&fmax);
681 UseFrac(fmin,fmax);
682 } else {
683 ndec--;
684 // S'il faut supprimer les options decodees
685 if (rmdecopt) udopt.push_back(opts);
686 }
687 }
688
689 // S'il faut supprimer les options decodees, on remplace l'argument opt
690 // par le vecteur des options non decodees.
691 if (rmdecopt) opt = udopt;
692
693 // Decodage des options par le histo-wrapper
694 int ndec2 = 0;
695 if ( ( optsz1-ndec1-ndec > 0) && (mHistoWp) ) {
696 ndec2 = mHistoWp->DecodeOptionString(opt, rmdecopt);
697 }
698
699 // Liste des options si demande
700 if(listopt) Print();
701
702 return(ndec+ndec1+ndec2);
703}
704
705int PIHisto2D::OptionToString(vector<string> & opt) const
706{
707 PIDrawer::OptionToString(opt);
708 char str[256];
709
710 if(stats) opt.push_back("stat"); else opt.push_back("nstat");
711
712 sprintf(str,"statposoff=%f,%f",spoX,spoY);
713 opt.push_back(str);
714
715 if(mTypDisp==0) sprintf(str,"h2disp=var");
716 else if(mTypDisp==1) sprintf(str,"h2disp=pts,%g",mFPoints);
717 else if(mTypDisp==2) sprintf(str,"h2disp=hbk");
718 else if(mTypDisp==3) sprintf(str,"h2disp=img");
719 opt.push_back(str);
720
721 if(mTypScal==0) sprintf(str,"h2scale=lin");
722 else if(mTypScal==1) sprintf(str,"h2disp=log,%g",mLogScale);
723 opt.push_back(str);
724
725 sprintf(str,"h2dyn=%g,%g",mHMin,mHMax);
726 opt.push_back(str);
727
728 sprintf(str,"h2frac=%g,%g",mFracMin,mFracMax);
729 opt.push_back(str);
730
731 if(mRevCmap) opt.push_back("revcmap");
732
733 // Les options du Histo2DWrapper :
734 if(mHistoWp) mHistoWp->OptionToString(opt);
735 return 1;
736}
737
738
739//++
740void PIHisto2D::GetOptionsHelpInfo(string& info)
741//
742// Help relatif au options
743//--
744{
745info += " ---- PIHisto2D options help info (see also ALT-O): \n" ;
746info += "- h2help: get this help text\n";
747info += "- h2list: list choosen options\n";
748info += "- sta,stat,stats: activate statistic display\n";
749info += " nsta,nstat,nostat,nostats: deactivate statistic display\n";
750info += "- h2disp=typ[,fracpts]: choose display type\n";
751info += " typ=var: variable size boxes\n";
752info += " typ=hbk: \"a la hbook2\"\n";
753info += " typ=img: image like (use \"h2col\" for color map)\n";
754info += " typ=pts: point clouds (fracpts=max possible fraction\n";
755info += " of used pixels per bin [0,1])\n";
756info += "- h2scale=lin/log[,logscale]: choose linear or logarithmic scale\n";
757info += "- h2dyn=[hmin][,hmax]: choose histogramme range for display\n";
758info += "- use graphic att. to define color table (ex: grey32,midas_heat,...)\n";
759info += " (see general graphicatt description)\n";
760info += "- use \"revcmap\" to reverse color table\n";
761info += "- h2frac=[fmin][,fmax]: choose sub-range display [0,1]\n";
762info += " ---- HistoWrapper options : \n" ;
763info += " hbincont: select bin content as Y value for display (default) \n";
764info += " hbinerr: select bin error as Y value for display \n";
765info += " hbinent: select bin entries as Y value for display \n";
766info += " hscale=value : multiplicative factor for Y value \n" ;
767info += " hoffset=value : additive coefficient for Y value \n" ;
768info += " hs1: set hscale=1 hoffset=0 (default) \n" ;
769// On recupere ensuite la chaine info de la classe de base
770PIDrawer::GetOptionsHelpInfo(info);
771return;
772}
773
774/////////////////////////////////////////////////////////////////
775// Classe H2WinArg
776/////////////////////////////////////////////////////////////////
777//++
778// Class H2WinArg
779// Lib PIext
780// include pihisto2d.h
781//
782// Fenêtre de dialogue pour le choix des options de tracé pour "PIHisto2D"
783// Classe de fenêtre de dialogue permettant de modifier interactivement
784// Les différents attributs de visualisation pour les *PIImage* .
785//--
786//++
787// Links Parents
788// PIWindow
789//--
790//++
791// Links Voir aussi
792// PIHisto2D
793// PIH2DWdg
794//--
795
796//++
797// Titre Constructeur, méthodes
798//--
799
800PIBaseWdgGen* H2WinArg::mBWdg = NULL;
801PIHisto2D* H2WinArg::mH2DDrw = NULL;
802static H2WinArg* cur_h2winarg = NULL;
803
804void H2WinArg::ShowPIHisto2DTools()
805{
806 if (cur_h2winarg == NULL) cur_h2winarg = new H2WinArg(PIApplicationGetApp());
807 cur_h2winarg->Show();
808}
809
810void H2WinArg::ShowPIHisto2DTools(PIBaseWdgGen* cbw)
811{
812 if (cur_h2winarg == NULL) cur_h2winarg = new H2WinArg(PIApplicationGetApp());
813 mBWdg = cbw;
814 cur_h2winarg->Show();
815}
816
817void H2WinArg::HidePIHisto2DTools()
818{
819 if (cur_h2winarg != NULL) cur_h2winarg->Hide();
820}
821
822void H2WinArg::SetCurrentBaseWdg(PIBaseWdgGen* cbw)
823{
824 mBWdg = cbw;
825}
826
827void H2WinArg::SetCurrentPIHisto2D(PIHisto2D* h2ddrw)
828{
829 mH2DDrw = h2ddrw;
830}
831
832PIBaseWdgGen* H2WinArg::GetCurrentBaseWdg()
833{
834 return(mBWdg);
835}
836
837PIHisto2D* H2WinArg::GetCurrentPIHisto2D()
838{
839 return(mH2DDrw);
840}
841
842//++
843H2WinArg::H2WinArg(PIMsgHandler* par)
844//
845// Creation de la fenetre de gestion des parametres
846// des dessins des histogrammes 2D. Cette fenetre de
847// dialogue est partagee par tous les widget de dessin
848// des histogrammes 2D. Pour faire apparaitre la fenetre
849// tapez ALT-O.
850//--
851//++
852//| - Menu 1: Choix du type de display
853//| Carres variables, nuages de points, caracteres a la hbook2
854//| et carres de tailles fixe (couleur ou niveauz de gris).
855//| - Menu 2: Choix du type d'echelle
856//| Lineaire ou logarithmique
857//| - Menu 3: Choix de la couleur
858//| noir et blanc, niveau de gris et couleurs diverses.
859//| - Champ texte Dyn: Pour donner la dynamique, si min>=max
860//| alors prend le min et le max de l'histogramme
861//| (cf PIHisto2D::UseDyn)
862//| - Champ texte Frac: fraction mini et maxi
863//| (cf PIHisto2D::UseFrac)
864//| - Champ texte LogScal: niveau de scaling pour le choix d'une
865//| echelle logarithmique (cf PIHisto2D::UseScale)
866//--
867//++
868//| - Curseur interactif PerPt: pourcentage de points a dessiner
869//| dans chaque bin (cf PIHisto2D::UseDisplay)
870//| - Bouton Apply: dessiner avec les options affichees
871//| - Bouton Dismiss: fermeture de la fenetre de dialogue.
872//| - Bouton Get: re-prendre les valeurs de display stoquees
873//| pour un histogramme donne.
874//| - Bouton Print: Imprimer les caracteristiques du display
875//| et de l'histogramme.
876//--
877: PIWindow((PIMsgHandler *)par,"H2D-Options",PIWK_dialog,250,260,150,150)
878, mFgCol(false), mCmap(CMAP_GREYINV32), mRevCmap(false)
879, mTypScal(0) , mLogScale(10.)
880, mTypDisp(0) , mFPoints(0.5)
881, mHMin(1.) , mHMax(-1.)
882, mFracMin(0.1), mFracMax(0.9)
883{
884if(dbg) printf("H2WinArg::H2WinArg %p par=%p\n",this,par);
885
886// Taille automatique
887int bsx, bsy;
888PIApplicationPrefCompSize(bsx, bsy); // environ 6 lettres
889int spx = (bsx>=10) ? bsx/10 : 1; // intervalle entre lettres X
890int spy = (bsy>=5) ? bsy/5 : 1; // intervalle entre lettres Y
891int wszx = 5*spx+bsx+int(2.5*bsx); // Taille fenetre en X
892int wszy = 11*spy+int(8.5*bsy); // Taille fenetre en Y
893SetSize(wszx, wszy);
894
895// menu du style de display des bins
896 int cpx = 2*spx, cpy = 2*spy;
897mOPop[0] = new PIOptMenu(this,"optmen-h2d-1",2*bsx,bsy,cpx,cpy);
898mOPop[0]->AppendItem("Carres Var." ,6101);
899mOPop[0]->AppendItem("....." ,6102);
900mOPop[0]->AppendItem(".+12..Z*" ,6103);
901mOPop[0]->AppendItem("Carres Pleins",6104);
902mOPop[0]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
903
904// Menu du choix de la dynamique
905 cpy += bsy+spy;
906mOPop[1] = new PIOptMenu(this,"optmen-h2d-2",2*bsx,bsy,cpx,cpy);
907mOPop[1]->AppendItem("Lineaire",6201);
908mOPop[1]->AppendItem("Log10" ,6202);
909mOPop[1]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
910
911// Menu du choix de la table des couleurs
912 cpy += bsy+spy;
913mOPop[2] = new PIOptMenu(this,"optmen-h2d-3",2*bsx,bsy,cpx,cpy);
914mOPop[2]->AppendItem("Black&White",7000);
915mCasc[0] = new PIMenu(mOPop[2]->Menu(),"PIStd-128Col");
916mCasc[1] = new PIMenu(mOPop[2]->Menu(), "PIUniCol32");
917mCasc[2] = new PIMenu(mOPop[2]->Menu(),"MIDAS-CMap");
918int kcc,nsct1=5,nsct2=9,nsct3=17,nsct4=PIColorMap::NumberStandardColorMaps()-2;
919for(kcc=0; kcc<nsct1; kcc++)
920 mOPop[2]->AppendItem(PIColorMap::GetStandardColorMapName(kcc).c_str(),7001+kcc);
921for(kcc=nsct1; kcc<nsct2; kcc++)
922 mCasc[0]->AppendItem(PIColorMap::GetStandardColorMapName(kcc).c_str(),7001+kcc);
923mOPop[2]->AppendPDMenu(mCasc[0]);
924for(kcc=nsct2; kcc<nsct3; kcc++)
925 mCasc[1]->AppendItem(PIColorMap::GetStandardColorMapName(kcc).c_str(),7001+kcc);
926mOPop[2]->AppendPDMenu(mCasc[1]);
927for(kcc=nsct3; kcc<nsct4; kcc++)
928 mCasc[2]->AppendItem(PIColorMap::GetStandardColorMapName(kcc).c_str(),7001+kcc);
929mOPop[2]->AppendPDMenu(mCasc[2]);
930for(kcc=nsct4; kcc<PIColorMap::NumberStandardColorMaps(); kcc++)
931 mOPop[2]->AppendItem(PIColorMap::GetStandardColorMapName(kcc).c_str(),7001+kcc);
932//mOPop[2]->SetValue(7000);
933mOPop[2]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
934
935// Reverse color map
936 cpy += bsy+spy;
937mCkb = new PICheckBox(this,"Reverse CMap",8001,2*bsx,bsy,cpx,cpy);
938mCkb->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
939
940// Labels et zones de saisie texte
941 cpy += bsy+spy;
942mLab[0] = new PILabel(this," Dyn: ",bsx,bsy,cpx,cpy);
943mLab[0]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
944mText[0] = new PIText(this,"Dynamique",int(2.5*bsx),bsy,cpx+bsx+spx,cpy);
945mText[0]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
946 cpy += bsy+spy;
947mLab[1] = new PILabel(this," Frac: ",bsx,bsy,cpx,cpy);
948mLab[1]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
949mText[1] = new PIText(this,"Fraction",int(2.5*bsx),bsy,cpx+bsx+spx,cpy);
950mText[1]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
951 cpy += bsy+spy;
952mLab[2] = new PILabel(this," LogScal: ",bsx,bsy,cpx,cpy);
953mLab[2]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
954mText[2] = new PIText(this,"LogScale",int(2.5*bsx),bsy,cpx+bsx+spx,cpy);
955mText[2]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
956
957// Labels et curseur mobile
958cpy += bsy+spy;
959mLab[3] = new PILabel(this," PerPt: ",bsx,bsy,cpx,cpy+int(0.25*bsy));
960mLab[3]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
961mPScal = new PIScale(this,"FracPoints",6401,kSDirLtoR
962 ,int(2.5*bsx),int(1.25*bsy),cpx+bsx+spx,cpy);
963mPScal->SetMinMax(0,100);
964mPScal->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
965
966//
967SetText();
968//
969
970// Boutons
971 cpx = 2*bsx+5*spx, cpy = 2*spy;
972mBut[0] = new PIButton(this,"Apply",6001,bsx,bsy,cpx,cpy);
973mBut[0]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
974 cpy += bsy+spy;
975mBut[1] = new PIButton(this,"Dismiss",6002,bsx,bsy,cpx,cpy);
976mBut[1]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
977 cpy += bsy+spy;
978mBut[2] = new PIButton(this,"Get",6003,bsx,bsy,cpx,cpy);
979mBut[2]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
980 cpy += bsy+spy;
981mBut[3] = new PIButton(this,"Print",6004,bsx,bsy,cpx,cpy);
982mBut[3]->SetBinding(PIBK_elastic,PIBK_elastic, PIBK_elastic,PIBK_elastic);
983// FinishCreate();
984}
985
986//++
987H2WinArg::~H2WinArg()
988//
989// Destructeur.
990//--
991{
992int i;
993if(dbg) printf("H2WinArg::~H2WinArg %p\n",this);
994for(i=0;i<3;i++) delete mCasc[i];
995for(i=0;i<3;i++) delete mOPop[i];
996for(i=0;i<4;i++) delete mBut[i];
997for(i=0;i<4;i++) delete mLab[i];
998for(i=0;i<3;i++) delete mText[i];
999delete mCkb;
1000delete mPScal;
1001}
1002
1003//++
1004void H2WinArg::Show()
1005//
1006// Initialisation sur ouverture ALT-O
1007//--
1008{
1009 if(dbg) printf("H2WinArg::Show() mH2DDrw=%p\n",mH2DDrw);
1010 // Pour recuperer les valeurs du Drawer sur lequel on fait ALT-O
1011 PIWindow::Show();
1012 return;
1013}
1014
1015//++
1016void H2WinArg::SetText()
1017//
1018// Gestion des fenetres de saisie de texte et des pop-menus.
1019//--
1020{
1021if(dbg) printf("H2WinArg::SetText()\n");
1022string sdum;
1023char str[256];
1024
1025sprintf(str,"%g %g",mHMin,mHMax);
1026mText[0]->SetText(str);
1027
1028sprintf(str,"%g %g",mFracMin,mFracMax);
1029mText[1]->SetText(str);
1030
1031sprintf(str,"%g",mLogScale);
1032mText[2]->SetText(str);
1033if(mTypScal==1) mText[2]->SetSensitive();
1034 else mText[2]->SetUnSensitive();
1035
1036if(mTypDisp==0) {sdum="Carres Var."; mOPop[0]->SetValueStr(sdum);}
1037else if(mTypDisp==1) {sdum="....."; mOPop[0]->SetValueStr(sdum);}
1038else if(mTypDisp==2) {sdum=".+12..Z*"; mOPop[0]->SetValueStr(sdum);}
1039else if(mTypDisp==3) {sdum="Carres Pleins"; mOPop[0]->SetValueStr(sdum);}
1040
1041if(mTypScal==0) {sdum="Lineaire"; mOPop[1]->SetValueStr(sdum);}
1042else if(mTypScal==1) {sdum="Log10"; mOPop[1]->SetValueStr(sdum);}
1043
1044if(!mFgCol) {mOPop[2]->SetValue(7000);}
1045else {
1046 for(int kk=0;kk<PIColorMap::NumberStandardColorMaps();kk++)
1047 if(mCmap == PIColorMap::GetStandardColorMapId(kk))
1048 {mOPop[2]->SetValue(7001+kk); break;}
1049}
1050mCkb->SetState(mRevCmap);
1051
1052mPScal->SetValue(int(mFPoints*100.));
1053if(mTypDisp==1) mPScal->SetSensitive();
1054 else mPScal->SetUnSensitive();
1055}
1056
1057//++
1058void H2WinArg::Process(PIMessage msg, PIMsgHandler* sender, void*)
1059//
1060// Gestions des messages.
1061//--
1062{
1063if(dbg) printf("H2WinArg::Process(%d-%d , %p ...) \n"
1064 ,(int)UserMsg(msg),(int)ModMsg(msg),sender);
1065
1066// if(!mH2Wdg) return;
1067if(!mBWdg) return;
1068// PIHisto2D* mpih = mH2Wdg->GetPIHisto();
1069PIHisto2D* mpih = mH2DDrw;
1070if(!mpih) return;
1071
1072int opt = UserMsg(msg);
1073
1074 if(opt == 6101) {mTypDisp = 0;}
1075else if(opt == 6102) {mTypDisp = 1;}
1076else if(opt == 6103) {mTypDisp = 2;}
1077else if(opt == 6104) {mTypDisp = 3;}
1078
1079else if(opt == 6201) {mTypScal = 0;}
1080else if(opt == 6202) {mTypScal = 1;}
1081
1082else if(opt == 7000) {mFgCol = false;}
1083else if(opt >= 7001 && opt <8000) {
1084 int k = opt-7001;
1085 mFgCol = true;
1086 mCmap = PIColorMap::GetStandardColorMapId(k);
1087}
1088
1089else if(opt == 8001) mRevCmap = mCkb->GetState();
1090
1091else if(opt == 6401) mFPoints = mPScal->GetValue()/100.;
1092
1093else if(opt==6001) {
1094 sscanf(mText[0]->GetText().c_str(),"%lf %lf",&mHMin,&mHMax);
1095 sscanf(mText[1]->GetText().c_str(),"%lf %lf",&mFracMin,&mFracMax);
1096 sscanf(mText[2]->GetText().c_str(),"%lf",&mLogScale);
1097 mpih->UseColors(mFgCol,mCmap,mRevCmap);
1098 mpih->UseScale(mTypScal,mLogScale);
1099 mpih->UseDisplay(mTypDisp,mFPoints);
1100 mpih->UseDyn(mHMin,mHMax);
1101 mpih->UseFrac(mFracMin,mFracMax);
1102 mBWdg->Refresh(); // On rafraichit le dessin (tout le PIScDrawWdg)
1103}
1104else if(opt==6002) this->Hide();
1105else if (opt==6003) {
1106 mFgCol = mpih->Color(); mCmap = mpih->ColMap();
1107 mRevCmap = mpih->IsColMapRev();
1108 mTypScal = mpih->TypScale(); mLogScale = mpih->LogScale();
1109 mTypDisp = mpih->TypDisplay(); mFPoints = mpih->FPoints();
1110 mHMin = mpih->HMin(); mHMax = mpih->HMax();
1111 mFracMin = mpih->FMin(); mFracMax = mpih->FMax();
1112}
1113else if(opt==6004) mpih->Print(2);
1114
1115SetText();
1116
1117if(dbg) {
1118 printf("H2WinArg::Process opt=%d col=%d,%d,%d scal=%d disp=%d npt=%g\n"
1119 ,opt,(int)mFgCol,(int)mCmap,(int)mRevCmap,mTypScal,mTypDisp,mFPoints);
1120 printf(" min,max= %g,%g frac= %g,%g logsc= %g\n"
1121 ,mHMin,mHMax,mFracMin,mFracMax,mLogScale);
1122}
1123
1124}
Note: See TracBrowser for help on using the repository browser.