You can subscribe to this list here.
| 2001 |
Jan
|
Feb
(1) |
Mar
|
Apr
|
May
|
Jun
|
Jul
(1) |
Aug
|
Sep
|
Oct
|
Nov
|
Dec
|
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 2002 |
Jan
(1) |
Feb
|
Mar
|
Apr
|
May
|
Jun
|
Jul
(1) |
Aug
(1) |
Sep
|
Oct
|
Nov
(1) |
Dec
|
| 2003 |
Jan
|
Feb
|
Mar
|
Apr
|
May
|
Jun
|
Jul
(1) |
Aug
(1) |
Sep
|
Oct
(83) |
Nov
(57) |
Dec
(111) |
| 2004 |
Jan
(38) |
Feb
(121) |
Mar
(107) |
Apr
(241) |
May
(102) |
Jun
(190) |
Jul
(239) |
Aug
(158) |
Sep
(184) |
Oct
(193) |
Nov
(47) |
Dec
(68) |
| 2005 |
Jan
(190) |
Feb
(105) |
Mar
(99) |
Apr
(65) |
May
(92) |
Jun
(250) |
Jul
(197) |
Aug
(128) |
Sep
(101) |
Oct
(183) |
Nov
(186) |
Dec
(42) |
| 2006 |
Jan
(102) |
Feb
(122) |
Mar
(154) |
Apr
(196) |
May
(181) |
Jun
(281) |
Jul
(310) |
Aug
(198) |
Sep
(145) |
Oct
(188) |
Nov
(134) |
Dec
(90) |
| 2007 |
Jan
(134) |
Feb
(181) |
Mar
(157) |
Apr
(57) |
May
(81) |
Jun
(204) |
Jul
(60) |
Aug
(37) |
Sep
(17) |
Oct
(90) |
Nov
(122) |
Dec
(72) |
| 2008 |
Jan
(130) |
Feb
(108) |
Mar
(160) |
Apr
(38) |
May
(83) |
Jun
(42) |
Jul
(75) |
Aug
(16) |
Sep
(71) |
Oct
(57) |
Nov
(59) |
Dec
(152) |
| 2009 |
Jan
(73) |
Feb
(213) |
Mar
(67) |
Apr
(40) |
May
(46) |
Jun
(82) |
Jul
(73) |
Aug
(57) |
Sep
(108) |
Oct
(36) |
Nov
(153) |
Dec
(77) |
| 2010 |
Jan
(42) |
Feb
(171) |
Mar
(150) |
Apr
(6) |
May
(22) |
Jun
(34) |
Jul
(31) |
Aug
(38) |
Sep
(32) |
Oct
(59) |
Nov
(13) |
Dec
(62) |
| 2011 |
Jan
(114) |
Feb
(139) |
Mar
(126) |
Apr
(51) |
May
(53) |
Jun
(29) |
Jul
(41) |
Aug
(29) |
Sep
(35) |
Oct
(87) |
Nov
(42) |
Dec
(20) |
| 2012 |
Jan
(111) |
Feb
(66) |
Mar
(35) |
Apr
(59) |
May
(71) |
Jun
(32) |
Jul
(11) |
Aug
(48) |
Sep
(60) |
Oct
(87) |
Nov
(16) |
Dec
(38) |
| 2013 |
Jan
(5) |
Feb
(19) |
Mar
(41) |
Apr
(47) |
May
(14) |
Jun
(32) |
Jul
(18) |
Aug
(68) |
Sep
(9) |
Oct
(42) |
Nov
(12) |
Dec
(10) |
| 2014 |
Jan
(14) |
Feb
(139) |
Mar
(137) |
Apr
(66) |
May
(72) |
Jun
(142) |
Jul
(70) |
Aug
(31) |
Sep
(39) |
Oct
(98) |
Nov
(133) |
Dec
(44) |
| 2015 |
Jan
(70) |
Feb
(27) |
Mar
(36) |
Apr
(11) |
May
(15) |
Jun
(70) |
Jul
(30) |
Aug
(63) |
Sep
(18) |
Oct
(15) |
Nov
(42) |
Dec
(29) |
| 2016 |
Jan
(37) |
Feb
(48) |
Mar
(59) |
Apr
(28) |
May
(30) |
Jun
(43) |
Jul
(47) |
Aug
(14) |
Sep
(21) |
Oct
(26) |
Nov
(10) |
Dec
(2) |
| 2017 |
Jan
(26) |
Feb
(27) |
Mar
(44) |
Apr
(11) |
May
(32) |
Jun
(28) |
Jul
(75) |
Aug
(45) |
Sep
(35) |
Oct
(285) |
Nov
(99) |
Dec
(16) |
| 2018 |
Jan
(8) |
Feb
(8) |
Mar
(42) |
Apr
(35) |
May
(23) |
Jun
(12) |
Jul
(16) |
Aug
(11) |
Sep
(8) |
Oct
(16) |
Nov
(5) |
Dec
(8) |
| 2019 |
Jan
(9) |
Feb
(28) |
Mar
(4) |
Apr
(10) |
May
(7) |
Jun
(4) |
Jul
(4) |
Aug
|
Sep
(4) |
Oct
|
Nov
(23) |
Dec
(3) |
| 2020 |
Jan
(19) |
Feb
(3) |
Mar
(22) |
Apr
(17) |
May
(10) |
Jun
(69) |
Jul
(18) |
Aug
(23) |
Sep
(25) |
Oct
(11) |
Nov
(20) |
Dec
(9) |
| 2021 |
Jan
(1) |
Feb
(7) |
Mar
(9) |
Apr
|
May
(1) |
Jun
(8) |
Jul
(6) |
Aug
(8) |
Sep
(7) |
Oct
|
Nov
(2) |
Dec
(23) |
| 2022 |
Jan
(23) |
Feb
(9) |
Mar
(9) |
Apr
|
May
(8) |
Jun
(1) |
Jul
(6) |
Aug
(8) |
Sep
(30) |
Oct
(5) |
Nov
(4) |
Dec
(6) |
| 2023 |
Jan
(2) |
Feb
(5) |
Mar
(7) |
Apr
(3) |
May
(8) |
Jun
(45) |
Jul
(8) |
Aug
|
Sep
(2) |
Oct
(14) |
Nov
(7) |
Dec
(2) |
| 2024 |
Jan
(4) |
Feb
(4) |
Mar
|
Apr
(7) |
May
(2) |
Jun
(1) |
Jul
|
Aug
(5) |
Sep
|
Oct
|
Nov
(4) |
Dec
(14) |
| 2025 |
Jan
(22) |
Feb
(6) |
Mar
(5) |
Apr
(14) |
May
(6) |
Jun
(11) |
Jul
(19) |
Aug
|
Sep
(17) |
Oct
(1) |
Nov
(2) |
Dec
(18) |
| 2026 |
Jan
|
Feb
|
Mar
(5) |
Apr
|
May
(2) |
Jun
(1) |
Jul
(6) |
Aug
(1) |
Sep
|
Oct
|
Nov
|
Dec
|
|
From: Hans-Bernhard B. <HBB...@t-...> - 2008-06-01 22:06:32
|
Ethan A Merritt wrote: > One might plausibly want to toggle log scaling on the colorbar but you can't, > because the hotkey tries to toggle Z also (or instead?) I suspect that's "instead". If memory serves, when mouse interaction was added, pm3d still assumed z and cb axes to be the same thing. 4-column pm3d came later. The colour axis may not even have _had_ a log setting to be toggled, at that time. > I propose to restrict the effect of 'l' and 'L' to the X,Y axes in 2D > and the Z axis in 3D. We can add another hotkey, probably 'c', to > explicitly toggle log scaling on the colorbar in both 2D and 3D plots. I don't think we need a separate hotkey. 'l' only acts on the axis the mouse is close to. So it should act on the cb axis only if the colorbox is visible, and the mouse is inside it. |
|
From: James R. V. Z. <jr...@co...> - 2008-05-31 17:38:19
|
>> > Ethan A Merritt writes:
...
>> > Furthermore, this would present an excuse^H^H^Hopportunity to
>> > re-work how data values are stored internally. Right now setting
>> > log-scaling on an axis causes the corresponding data coordinate to
>> > be transformed and stored as the log. This creates great headaches
>> > if you want to toggle the log-scaling and redraw the plot from the
>> > data previously read in. I think it would be much cleaner and more
>> > flexible to store the original coordinate value on input, and only
>> > apply the log- or other scaling during plot generation.
>>
>> Yes! One reason I didn't take my patch any further was that the "last
>> minute scaling" revision ought to be done first.
>
>I had come to the opposite conclusion.
>I figured to leave the existing log-scale code in place, as you
>suggested above, while adding a new more general mecanism that worked
>on the original stored coordinates. Once the general mechanism was
>in place and working, the old log-scale code could be removed without
>ever having to modify it to work with "last-minute" scaling.
I think I see how that would work:
Currently:
Upon "set log x", all stored x values are transformed and on "unset
log x" they are inverse transformed.
Phase 1:
Add a scaling "newlog".
Add to AXIS a pointer to a function that transforms the data. The
function takes two parameters: the datum to be transformed, and a
double. For scaling "linear" or "log" the function is a no-op. For
"newlog", the extra parameter is the base of the logs. Otherwise the
extra parameter is ignored.
Always do "last-minute" scaling: When plotting a value, call the
function via the supplied pointer to transform it.
On "set log x" or "unset prob x" etc., update the function pointer.
On "(un)set log x", continue to (un)transform the stored data.
Phase 2 (after initial checkout):
Rename scaling "log" to "oldlog", and rename "newlog" to "log".
Phase 3 (after a major release):
Delete the "oldlog" commands and the code to (un)transform the stored
data.
The above only implements standard scalings (linear, log, probability,
and maybe weibull). To accommodate a user-defined scaling function,
the scaling step needs to be able to point to an action table. We
could use that mechanism to implement the standard scalings too -
building the action tables from static strings like "_linear(x,b)=x"
or "_log(x,b)=log(x)/log(b)". However, I'm worried that would be too
slow. The faster method would be to give the scaling function three
parameters (datum, extra, and action table), and let the standard
scaling functions ignore the third parameter.
- Jim Van Zandt
|
|
From: Ethan A M. <merritt@u.washington.edu> - 2008-05-31 15:10:37
|
On Saturday 31 May 2008 07:26, Allin Cottrell wrote: > I just noticed that if a gnuplot file contains a detailed > specification of an area to be filled with the filledcurve style, > the metapost terminal writes out a very long "fill" line, which > can exceed the length that mpost will process. > > The attached patchette arranges for such a fill line to be broken, > in the same way that a long "draw" line is broken. Applied to CVS -- Ethan A Merritt |
|
From: Allin C. <cot...@wf...> - 2008-05-31 14:27:03
|
I just noticed that if a gnuplot file contains a detailed specification of an area to be filled with the filledcurve style, the metapost terminal writes out a very long "fill" line, which can exceed the length that mpost will process. The attached patchette arranges for such a fill line to be broken, in the same way that a long "draw" line is broken. -- Allin Cottrell Department of Economics Wake Forest University, NC |
|
From: Ethan A M. <merritt@u.washington.edu> - 2008-05-31 04:38:31
|
Since I never use the logscale toggle keys ('l' and 'L') I hadn't realized
they are broken. The problem is that coloring options have proliferated, and
many kinds of plots now use the palette for things other than Z-coordinate.
For example, consider the last plot in heatmaps.dem.
One might plausibly want to toggle log scaling on the colorbar but you can't,
because the hotkey tries to toggle Z also (or instead?)
Many other examples of confusion are possible, also in 2D now that
palette coloring can be used in 2D also.
I propose to restrict the effect of 'l' and 'L' to the X,Y axes in 2D
and the Z axis in 3D. We can add another hotkey, probably 'c', to
explicitly toggle log scaling on the colorbar in both 2D and 3D plots.
Yes, this is a change.
Is that OK with everyone?
--
Ethan A Merritt
Biomolecular Structure Center
University of Washington, Seattle 98195-7742
|
|
From: Tim H. <hof...@hi...> - 2008-05-28 17:24:42
|
> A log file might be a good things to have. But I wonder if it would be worth aiming > for a portable solution - although we do not have the problem on other platforms. > (Some people would prefer to drop the windows console completly, I guess. IMHO > we should just make it just a bit more user friendly: readline support, more drag'n drop, > descent menu layout and a new button bar etc.) It definitly could need it a bit polishing. :) I don't think that one needs a portable solution. All other systems can redirect the error stream which is more natural there (would also be my favorite solution in windows, but this does not work). So I would leave it in the windows part, except someone really wants it for all systems. > The implementation for Windows via "MyPrintf" etc. would just be fine. Maybe putting the > code in TextPutStr() (wtext.c) would be even simpler - one has to check. In the My* functions there is still the file handle available. So I can check if the text is addressed to stderr and log only errors. Put apparently everything is coming via stderr, even the information at startup. I checked via (f==stderr) like it's done in the isterm() macro. Does someone know if this is right? Maybe it's a windows problem because the streams do not really exist. However since I call with a plotfile as parameter all the terminal parts like startup information and prompt aren't written anyway. > Btw. I would suggest -l <filename> instead of -log <filename> that's fine by me. Tim -- Tim Hoffmann Raum 052 HISKP Tel: +49-228-73-2942 Universitaet Bonn Fax: +49-228-73-2505 Nussallee 14-16 hof...@hi... 53115 Bonn http://www.hiskp.uni-bonn.de/gruppen/hikari |
|
From: Tim H. <tim...@un...> - 2008-05-28 17:05:07
|
Allin Cottrell wrote: > On Wed, 28 May 2008, Tim Hoffmann wrote: > wgnuplot 2>error.log > > That's supposed to work for win32 console programs at least. Right, but since wgnuplot is not a console program it does not. I tried but it creates just empty log files. As far as I could figure out, a WIN32 GUI application silently throws away everything written to stderr! However one can create a new console and redirect output to it. But it appears to be impossible to redirect it to the calling console. On the other hand, you could make it a console application and start the window yourself. But then you'll have a console poping up when you start the program from the menu. I think it's the same reason why there has to be pgnuplot.exe for piped input. So I suppose, a log file is the best choice in this case. Tim |
|
From: Allin C. <cot...@wf...> - 2008-05-28 15:26:39
|
On Wed, 28 May 2008, Tim Hoffmann wrote: > I am calling wgnuplot from an external program. Since there is no error > stream attached one cannot receive error messages. The exit code just > indicates that there is an error but not which. Therefore it would help, > if there was an option to write errors to a file, e.g.: > > wgnuplot -log error.log What happens if you do wgnuplot 2>error.log That's supposed to work for win32 console programs at least. Allin Cottrell |
|
From: Ethan M. <merritt@u.washington.edu> - 2008-05-28 13:26:50
|
On Wednesday 28 May 2008, Tim Hoffmann wrote: > Hello, > > I am calling wgnuplot from an external program. Since there is no error > stream attached one cannot receive error messages. Errors are sent to stderr, the usual error stream. In a normal operating environment (which may exclude MSWindows) you can receive the messages just fine. > The exit code just > indicates that there is an error but not which. The current CVS version also stores the most recent error message in the variable GPVAL_ERRMSG. You can print it to a file if you want to. > Therefore it would help, > if there was an option to write errors to a file, e.g.: > > wgnuplot -log error.log > > I think it should be quite simple. One could intercept the error > messages in MyFPrintF and similar functions. > > Before starting to implement, I would like to have some comments on the > idea form more experienced gnuplot developers. I am not sure there are any such people, with regard to Windows development. I gather that the Windows user interface (by which I mean Windows itself, not the gnuplot GUI) is fundamentally different from the stdin/stout/stderr model used by gnuplot and pretty much all other unix or C programs. -- Ethan Merritt (on the road) |
|
From: Tim H. <tim...@un...> - 2008-05-28 12:52:25
|
Hello, I am calling wgnuplot from an external program. Since there is no error stream attached one cannot receive error messages. The exit code just indicates that there is an error but not which. Therefore it would help, if there was an option to write errors to a file, e.g.: wgnuplot -log error.log I think it should be quite simple. One could intercept the error messages in MyFPrintF and similar functions. Before starting to implement, I would like to have some comments on the idea form more experienced gnuplot developers. Thanks, Tim |
|
From: Ethan M. <merritt@u.washington.edu> - 2008-05-27 22:02:30
|
On Sunday 25 May 2008, Philipp K. Janert wrote: > There seems to be a problem with the "line sample" > that is placed in the key for all styles that involve > boxes (w boxes, w boxerrorbars, w candlesticks, ...). > > For such styles, no line sample (actually: box sample) > is visible in the key. This problem is new with 4.3 and > seems to exist for all terminals (well, I tried wxt and eps). Fixed, thanks. Ethan > Best, > > Ph. > > ------------------------------------------------------------------------- > This SF.net email is sponsored by: Microsoft > Defy all challenges. Microsoft(R) Visual Studio 2008. > http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/ > _______________________________________________ > gnuplot-beta mailing list > gnu...@li... > https://lists.sourceforge.net/lists/listinfo/gnuplot-beta -- Ethan Merritt (on the road) |
|
From: <pl...@pi...> - 2008-05-26 01:51:27
|
On Mon, 26 May 2008 01:41:31 +0200, Philipp K. Janert <ja...@ie...> wrote: > The code does what all current gnuplot > smoothing algos do: they stop at the min > and max data point in the sample. I think > this is reasonable. > Best, > Ph. What I would suggest is that it only produces a plot line over the range where the data is complete. If the sample is large enough for this to be negligable , it won't notice anyway. If it is significant in relation to the data sample, the plot line will stop at the point where it is no longer mathematically valid. That would seem to be the correct thing to do. I see little justification for extending it beyond that range. It must be a trivial change to make if you accept the principal. best regards, Peter. |
|
From: <pl...@pi...> - 2008-05-26 00:52:05
|
On Mon, 26 May 2008 01:41:31 +0200, Philipp K. Janert <ja...@ie...> wrote: > On Saturday 24 May 2008 00:11, you wrote: >> On Sat, 24 May 2008 01:15:37 +0200, Philipp K. Janert <ja...@ie...> >> >> wrote: >> > I just submitted a patch (1970923) which >> > draws a smooth histogram-like curve >> > for a random collection of points, using >> > a Gaussian kernel density estimation >> > algorithm. >> > Demos are found here: >> > www.philipp-janert.com/kdensity >> >> very interesting. >> >> My initial impression on looking at your top left example is that there >> is >> a phase shift of +half a box in x most visible in the 0.01 and 0.05 >> plots. >> > The edge effect is actually in the histogram, > not in the kernel density (yet another advantage > of k-densities over histograms: the annoying > bin-placement problem goes away). hmm, never been much of a fan of bins and histograms, that's probably why. More for sociologists and economists. > > The code does what all current gnuplot > smoothing algos do: they stop at the min > and max data point in the sample. I think > this is reasonable. > Well I'm not sure that is comparable. IRRC all the "smoothing" algos (appart from unique) are splines , these are calculated over 4 data. In fact they would require just one point outside the data range at each end. I have not looked how they are dealt with but it is unlikely to be important for one point. However, techniques using a kernel require half the kernel width outside each end of the data range. I would guess by looking at your examples that the missing data are initialised as zero. Is that correct? Sorry to be a stickler for detail, it must be my rigourous physics training coming out. As people become less and less aware of what all these software tools are actually doing for them, it becomes more and more important that they do not introduce distortions. Don't think I knocking your efforts, I'm pretty impressed overall. best regards, Peter. > Best, > > Ph. > |
|
From: Philipp K. J. <ja...@ie...> - 2008-05-25 23:47:27
|
There seems to be a problem with the "line sample" that is placed in the key for all styles that involve boxes (w boxes, w boxerrorbars, w candlesticks, ...). For such styles, no line sample (actually: box sample) is visible in the key. This problem is new with 4.3 and seems to exist for all terminals (well, I tried wxt and eps). Best, Ph. |
|
From: Philipp K. J. <ja...@ie...> - 2008-05-25 23:41:33
|
On Saturday 24 May 2008 00:11, you wrote: > On Sat, 24 May 2008 01:15:37 +0200, Philipp K. Janert <ja...@ie...> > > wrote: > > I just submitted a patch (1970923) which > > draws a smooth histogram-like curve > > for a random collection of points, using > > a Gaussian kernel density estimation > > algorithm. > > Demos are found here: > > www.philipp-janert.com/kdensity > > very interesting. > > My initial impression on looking at your top left example is that there is > a phase shift of +half a box in x most visible in the 0.01 and 0.05 plots. > The edge effect is actually in the histogram, not in the kernel density (yet another advantage of k-densities over histograms: the annoying bin-placement problem goes away). The code does what all current gnuplot smoothing algos do: they stop at the min and max data point in the sample. I think this is reasonable. Best, Ph. |
|
From: Hans-Bernhard B. <HBB...@t-...> - 2008-05-25 09:28:57
|
[Cc: to list reinstated...] Philipp K. Janert wrote: > How does everybody else feel about this? I think it's quite bad. With the '=' operator in 'using' expressions, gnuplot shouldn't even have to be changed to compute that value. And yes, the proposed syntax is abominable. Yes, that's me speaking, who accepted the MAXITER=0 special case eleven years ago:-) |
|
From: <pl...@pi...> - 2008-05-24 07:16:22
|
On Sat, 24 May 2008 01:03:49 +0200, Philipp K. Janert <ja...@ie...> wrote: > > Ethan has asked me to take a look at this > patch. > > I guess the purpose of this patch is to allow > the user to enter a negative value for > MAX_ITER, whereupon the fit command will > calculate the residual using the starting > values of all fit parameters and exits. (That is > really zero iteration steps, but MAX_ITER==0 > has already been used to indicate unlimited > iteration). > > I am not very happy with this patch: > - I don't understand the purpose. Why do > you care about the unfitted residual? > - It leads to partially garbage output (the > errors in the fit parameters are not > defined for no iteration, etc) > - We are (ab-)using the MAX_ITER parameter > by introducing a magic value with > non-trivial semantics (which is never > a good idea in my mind). > > How does everybody else feel about this? > Who would find this patch useful, and why? > > (If it is useful, can we find a better way of > achieving the same goal, rather than shoe-horning > semantics into a negative iteration count?) > > Best, > > Ph. > Hi, I agree with your general comment , I find there is too much of this "magic number" and syntax tweeking approach in gnuplot. It makes it very difficult to remember and even harder to discover some very useful functions. It can get very aesoteric at times. best regards, Peter. |
|
From: <pl...@pi...> - 2008-05-24 07:10:59
|
On Sat, 24 May 2008 01:15:37 +0200, Philipp K. Janert <ja...@ie...> wrote: > I just submitted a patch (1970923) which > draws a smooth histogram-like curve > for a random collection of points, using > a Gaussian kernel density estimation > algorithm. > Demos are found here: > www.philipp-janert.com/kdensity very interesting. My initial impression on looking at your top left example is that there is a phase shift of +half a box in x most visible in the 0.01 and 0.05 plots. Considering that x=0 is in the middle of the first box it appears that the fits are responding in a way that aligns with the right of each box. It's a bit subjective due to the nature of the data but this is my impression for the peaks at 0.2 0.4 and 0.5 Maybe you could test this effect with a rapid change in the data. I also think there is an egde effect at the begining and end of the data. This is a common problem when applying this sort of technique to image data. How to deal with edges when the kernel goes outside the data. There are several "solutions" which involve falsely extening the data but applying a kernel to an incompete sample range is effectively filling it with zeros and is equally false. It's like a running mean cannot be meaningful upto the edges of the sample range since there are not enough samples to take the mean over. This also gives an artificial drop off at the edges. This is also rather marked near the origin in your lognormal example. In image processing it's just a case of prettying up the edges but in a scientific context this is clearly not appropriate. I think the only rigourous way to deal with this is not to plot the part where the data is incomplete. I hope the comments are useful. best regards, Peter. |
|
From: <so...@pi...> - 2008-05-24 06:36:43
|
On Sat, 24 May 2008 01:03:49 +0200, Philipp K. Janert <ja...@ie...> wrote: > > Ethan has asked me to take a look at this > patch. > > I guess the purpose of this patch is to allow > the user to enter a negative value for > MAX_ITER, whereupon the fit command will > calculate the residual using the starting > values of all fit parameters and exits. (That is > really zero iteration steps, but MAX_ITER==0 > has already been used to indicate unlimited > iteration). > > I am not very happy with this patch: > - I don't understand the purpose. Why do > you care about the unfitted residual? > - It leads to partially garbage output (the > errors in the fit parameters are not > defined for no iteration, etc) > - We are (ab-)using the MAX_ITER parameter > by introducing a magic value with > non-trivial semantics (which is never > a good idea in my mind). > > How does everybody else feel about this? > Who would find this patch useful, and why? > > (If it is useful, can we find a better way of > achieving the same goal, rather than shoe-horning > semantics into a negative iteration count?) > > Best, > > Ph. > Hi, I agree with your general comment , I find there is too much of this "magic number" and syntax tweeking approach in gnuplot. It makes it very difficult to remember and even harder to discover some very useful functions. It can get very aesoteric at times. best regards, Peter. |
|
From: Philipp K. J. <ja...@ie...> - 2008-05-23 23:15:36
|
I just submitted a patch (1970923) which draws a smooth histogram-like curve for a random collection of points, using a Gaussian kernel density estimation algorithm. Demos are found here: www.philipp-janert.com/kdensity The new method has the following advantages over the classic way of generating histograms using "smooth frequency": - the resulting histogram is a smooth curve, making the effect of binning less severe - it handles intermediate "bins" with no points in them gracefully. (smooth freq does so only if used "with boxes", but if you use "with lines" for example, the line will not drop to zero if an intermediate bin is empty) The method is invoked like a weighted smoothing algorithm: plot "data" u 1:(1):(1) smooth kdensity where the 2nd parameter is the weight of each point and the 3rd parameter is the bandwidth to be used. This patch complements the "smooth cumulative" algorithm as another way to visualize the distribution of a collection of random points. Comments and suggestions are welcome. Best, Ph. |
|
From: Philipp K. J. <ja...@ie...> - 2008-05-23 23:03:50
|
Ethan has asked me to take a look at this patch. I guess the purpose of this patch is to allow the user to enter a negative value for MAX_ITER, whereupon the fit command will calculate the residual using the starting values of all fit parameters and exits. (That is really zero iteration steps, but MAX_ITER==0 has already been used to indicate unlimited iteration). I am not very happy with this patch: - I don't understand the purpose. Why do you care about the unfitted residual? - It leads to partially garbage output (the errors in the fit parameters are not defined for no iteration, etc) - We are (ab-)using the MAX_ITER parameter by introducing a magic value with non-trivial semantics (which is never a good idea in my mind). How does everybody else feel about this? Who would find this patch useful, and why? (If it is useful, can we find a better way of achieving the same goal, rather than shoe-horning semantics into a negative iteration count?) Best, Ph. |
|
From: Ethan M. <merritt@u.washington.edu> - 2008-05-23 21:23:46
|
On Friday 23 May 2008 14:13, Philipp K. Janert wrote: > > I have applied and closed the patches: > 1962130 (smooth cumulative) > 1827826 (additional dgrid3d kernels) > > I suggest we also close patch > 632289 (dgrid3d: scaling and range) > because it is by now over 5 years out > of date and from what I can see has > been superseded by the just applied > patch 1827826. Any comments or > objections? OK by me. I'm sure it's not the only patch that is out of date. Ethan -- Ethan A Merritt |
|
From: Philipp K. J. <ja...@ie...> - 2008-05-23 21:13:46
|
I have applied and closed the patches: 1962130 (smooth cumulative) 1827826 (additional dgrid3d kernels) I suggest we also close patch 632289 (dgrid3d: scaling and range) because it is by now over 5 years out of date and from what I can see has been superseded by the just applied patch 1827826. Any comments or objections? Best, Ph. |
|
From: Thomas S. <t.s...@fz...> - 2008-05-23 10:10:50
|
> That's a bit strange, because apparently it segfaulted in the
> hidden-surface bookkeeping routines rather than anything having to do
> with bitmap management. I suppose one of those overflowed coordinates
> causes trouble in the depth-sort algorithm.
at least i have found what happens:
access to non-existing elements of 'quadtree' in procedure 'in_front' in
'hidden3d.c' resulting in useless values for 'listhead'
> for (grid_x = grid_x_low; grid_x <= grid_x_high; grid_x ++)
> for (grid_y = grid_y_low; grid_y <= grid_y_high; grid_y ++)
> for (listhead = quadtree[grid_x][grid_y]; <<====THERE====
but why it happens - don't know.
the following patch prevents segfaults and gives warnings:
------------------------------------------------------------------------------
--- hidden3d.c.orig 2008-03-29 10:28:01.000000000 +0100
+++ hidden3d.c 2008-05-23 11:56:46.000000000 +0200
@@ -1801,5 +1801,21 @@
grid_x_low = COORD_TO_TREECELL(xmin);
+ if (grid_x_low < 0) {
+ grid_x_low = 0;
+ fprintf(stderr,"in_front: grid_x_low set to 0\n");
+ }
grid_x_high = COORD_TO_TREECELL(xmax);
+ if (grid_x_high >= QUADTREE_GRANULARITY) {
+ grid_x_high = QUADTREE_GRANULARITY-1;
+ fprintf(stderr,"in_front: grid_x_high set to
QUADTREE_GRANULARITY-1\n");
+ }
grid_y_low = COORD_TO_TREECELL(ymin);
+ if (grid_y_low<0) {
+ grid_y_low = 0;
+ fprintf(stderr,"in_front: grid_y_low set to 0\n");
+ }
grid_y_high = COORD_TO_TREECELL(ymax);
+ if (grid_y_high >= QUADTREE_GRANULARITY) {
+ grid_y_high = QUADTREE_GRANULARITY-1;
+ fprintf(stderr,"in_front: grid_y_high set to
QUADTREE_GRANULARITY-1\n");
+ }
------------------------------------------------------------------------------
gnuplot> set term pbm large size 100,100
gnuplot> set out 'test.pbm'
gnuplot> load 'singulr.dem'
in_front: grid_x_high set to QUADTREE_GRANULARITY-1
in_front: grid_y_high set to QUADTREE_GRANULARITY-1
in_front: grid_x_high set to QUADTREE_GRANULARITY-1
in_front: grid_y_high set to QUADTREE_GRANULARITY-1
...
Hit return to continue (1)
in_front: grid_x_high set to QUADTREE_GRANULARITY-1
in_front: grid_y_high set to QUADTREE_GRANULARITY-1
...
it seems that the 'xmax' and 'ymax' values, set in the macro 'setup_edge',
are wrong sometimes.
--
View this message in context: http://www.nabble.com/DPU-414-terminal-tp17371497p17422742.html
Sent from the Gnuplot - Dev mailing list archive at Nabble.com.
|
|
From: Jonathan T. <J.T...@so...> - 2008-05-22 21:03:25
|
Hi,
On Thu, 22 May 2008, Ethan Merritt wrote:
[[about the PBM driver]]
> > -- I still use it fairly often
> > to generate individual frames for encoding into a movie with ppmtompeg
> > (which only groks a limited set of native input formats: PPM, PNM,
> > YUV, JPEG, and JMOVIE, but *not* PNG).
>
> That's not much of a reason.
> You could use the far more featureful png driver and interpose
> a png->pnm filter in the output:
> set term png truecolor enhanced font "verdana,11"
> set output '| convert png:- mypic.pnm'
Good idea, I'll play around the next time I make a movie. The
only problem is that crunching 3000 frames is already a cpu- and
disk-intensive process, and firing up ImageMagick for each frame
would make it considerably worse. What might be better would be
to do a pnmtopnm (part of netpbm) conversion on each frame, since
I already have a pnmchange (also part of netpbm) conversion there.
The netpbm programs are a lot lighter-weight than ImageMagick.
ciao,
--
-- Jonathan Thornburg (remove -animal to reply) <J.T...@so...>
School of Mathematics, U of Southampton, England
"Washing one's hands of the conflict between the powerful and the
powerless means to side with the powerful, not to be neutral."
-- quote by Freire / poster by Oxfam
|