Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]

Groups > comp.graphics.apps.gnuplot > #2553 > unrolled thread

difficult for loop synax :-|

Back to article view | Back to comp.graphics.apps.gnuplot

Contents

  difficult for loop synax :-| Kalin Kozhuharov <me.kalin@gmail.com> - 2014-08-19 11:57 -0700
    Re: difficult for loop synax :-| Ethan A Merritt <EAMerritt@gmail.com> - 2014-08-19 18:03 -0700
      Re: difficult for loop synax :-| Kalin Kozhuharov <me.kalin@gmail.com> - 2014-08-19 23:15 -0700
        Re: difficult for loop synax :-| Ethan A Merritt <sfeam@users.sourceforge.net> - 2014-08-20 10:15 -0700
          Re: difficult for loop synax :-| Kalin Kozhuharov <me.kalin@gmail.com> - 2014-08-20 11:50 -0700
            Re: difficult for loop synax :-| Ethan A Merritt <EAMerritt@gmail.com> - 2014-08-20 17:50 -0700
              Re: difficult for loop synax :-| Kalin Kozhuharov <me.kalin@gmail.com> - 2014-08-20 23:02 -0700

#2553 — difficult for loop synax :-|

I have been trying to plot two dependent functions per single for loop, but cannot get the syntax right...

This works but is a lot of typing and if I don't know the number of arguments cannot be expressed:

	plot  C=1, sin(C*x) w l, C=2, sin(C*x) w l, C=3, sin(C*x) w l;

Any ideas?

This works:

 plot for [n in "1 2 3"]  n*x, for [n in "1 2 3"] sin(n*x) w l;

but those are independent functions, so 3 times n*x is plotted, then 3 times sin(n*x) in that order

I need something like

	plot for [n in "1 2 3"]  {C=n, sin(C*x) w l};
i.e. set some parameter by a function, then plot something, then evaluate the parameter again, then plot second function...

Fiddling with multiplot is hard it seems, need to hardcode linestyles, etc.

Kalin.

[toc] | [next] | [standalone]

#2554

Kalin Kozhuharov wrote:

> I have been trying to plot two dependent functions per single for
> loop, but cannot get the syntax right...
> 
> This works but is a lot of typing and if I don't know the number of
> arguments cannot be expressed:
> 
> plot  C=1, sin(C*x) w l, C=2, sin(C*x) w l, C=3, sin(C*x) w l;
> 
> Any ideas?

Please clarify.   I see only a single function there, evaluated 
three times.   Where is the second dependent function?

> 
> This works:
> 
>  plot for [n in "1 2 3"]  n*x, for [n in "1 2 3"] sin(n*x) w l;
> 
> but those are independent functions, so 3 times n*x is plotted, then 3
> times sin(n*x) in that order
> 
> I need something like
> 
> plot for [n in "1 2 3"]  {C=n, sin(C*x) w l};
> i.e. set some parameter by a function, then plot something, then
> evaluate the parameter again, then plot second function...

I still see only a single function.
That command appears to be equivalent to 

    plot for [C=1:3] sin(C*x)

You must have some additional step in mind, I guess?

> Fiddling with multiplot is hard it seems, need to hardcode linestyles,
> etc.
> 
> Kalin.

[toc] | [prev] | [next] | [standalone]

#2555

On Wednesday, August 20, 2014 10:03:27 AM UTC+9, Ethan A Merritt wrote:
> Kalin Kozhuharov wrote:
> > I have been trying to plot two dependent functions per single for
> > loop, but cannot get the syntax right...
> > This works but is a lot of typing and if I don't know the number of
> > arguments cannot be expressed:
> 
> > plot  C=1, sin(C*x) w l, C=2, sin(C*x) w l, C=3, sin(C*x) w l;
> 
> Please clarify.   I see only a single function there, evaluated 
> three times.   Where is the second dependent function?
> 
Well, may be I oversimplified it, but the real thing is too complex anyway.
What about this:

reset; f1(x) = (C=1.134324*x,0); f2(x) = sin(C*x); set xrange [0:3.14];
plot a=f1(1), f2(x) title "n=1", a=f1(2), f2(x) title "n=2" , a=f1(3), f2(x) title "n=3";

How can I turn it into something like:
N="1 2 3"; plot for [n in N] a=f1(n), f2(x) title "n=" . n;

so that it produces the same result (all plots on same graph, automatically linestyled) when I don't know N. Currently this produces only one graph when n=3, because the for loop is only evaluated till the ",", then f2(x) is plotted once. I need some kind of syntactic "brackets".

I also tried:
N="1 2 3"; do for [n in N] {plot  a=f1(n), f2(x) title "n=" . n; pause .2; };
and saw how the plot is being generated, then overwritten by the next call of plot.

Cheers,
Kalin.

[toc] | [prev] | [next] | [standalone]

#2556

Kalin Kozhuharov wrote:

> On Wednesday, August 20, 2014 10:03:27 AM UTC+9, Ethan A Merritt wrote:
>> Kalin Kozhuharov wrote:
>> > I have been trying to plot two dependent functions per single for
>> > loop, but cannot get the syntax right...
>> > This works but is a lot of typing and if I don't know the number of
>> > arguments cannot be expressed:
>> 
>> > plot  C=1, sin(C*x) w l, C=2, sin(C*x) w l, C=3, sin(C*x) w l;
>> 
>> Please clarify.   I see only a single function there, evaluated
>> three times.   Where is the second dependent function?
>> 
> Well, may be I oversimplified it, but the real thing is too complex
> anyway. What about this:
> 
> reset; f1(x) = (C=1.134324*x,0); f2(x) = sin(C*x); set xrange [0:3.14];
> plot a=f1(1), f2(x) title "n=1", a=f1(2), f2(x) title "n=2" , a=f1(3),
> f2(x) title "n=3";

It is still not clear what you are aiming for.
The sequence of commands above does not make sense.
The function f1(x) is never plotted, so it plays no role except for
its side effect of setting C to a constant.
A simpler way to write it would be
   
   f1(x) = 1.134324*x
   plot for [n=1:3] sin(f1(n)*x) title sprintf("n=%d",n)

> How can I turn it into something like:
> N="1 2 3"; plot for [n in N] a=f1(n), f2(x) title "n=" . n;
> when I don't know N

Same as above, replacing 3 with N

   plot for [n=1:N] ...

Of course N still has to come from somewhere, I assume from some
earlier part of the script.

> so that it produces the same result (all plots on same graph,
> automatically linestyled) when I don't know N. Currently this produces
> only one graph when n=3, because the for loop is only evaluated till the
> ",", then f2(x) is plotted once. I need some kind of syntactic "brackets".
> 
> I also tried:
> N="1 2 3"; do for [n in N] {plot  a=f1(n), f2(x) title "n=" . n; pause .2;
> }; and saw how the plot is being generated, then overwritten by the next
> call of plot.
> 
> Cheers,
> Kalin.
-- 
Ethan A Merritt

[toc] | [prev] | [next] | [standalone]

#2557

Hello Ethan,

Thank you for your time! I was trying to make self-contained example, but because it was simple enough you "optimized" it ;-) "Premature optimization is the root of all evil", said D. Knuth :-D

On Thursday, August 21, 2014 2:15:51 AM UTC+9, Ethan A Merritt wrote:
>> Well, may be I oversimplified it, but the real thing is too complex
>> anyway. What about this:
>> 
>> reset; f1(x) = (C=1.134324*x,0); f2(x) = sin(C*x); set xrange [0:3.14];
>> plot a=f1(1), f2(x) title "n=1", a=f1(2), f2(x) title "n=2" , a=f1(3),
>> f2(x) title "n=3";
> 
> It is still not clear what you are aiming for.
> The sequence of commands above does not make sense.
> 
> The function f1(x) is never plotted, so it plays no role except for
> its side effect of setting C to a constant.
>
Well, yes I'd have to admit I am using side-effects here.
And my code is turning to a small mess because I am trying to have gnuplot do statistics in-line, instead of preprocessing the data.

My current code is here: https://github.com/thinrope/fixed_sensor_visualization/blob/master/timeplot_all.gpl#L26

plot for [ID in IDs] "tmp/" . ID . ".data" u 1:2 axes x1y2 w l lw 2 title graph_title(ID);

What I wanted to do is add another SMA filter on each datafile (or 2 with different window) and still have it plot similar graph.

so having defined SMA as in https://github.com/thinrope/fixed_sensor_visualization/blob/master/timeplot.gpl#L46 I was planning to write something along:

plot for [ID in IDs] a=sma_24_init(0), "tmp/" . ID . ".data" u 1:(sma_24($2)) axes x1y2 w l lw 2 title graph_title(ID);
 
IDs are list of numbers coming form the environment/Makefile, but can be strings. 

May be possible with `set table` and `do for [ID in IDs]{init; plot;}` but I'd rather avoid it.


Alternatively, something along

init_OK=0; f3(x,n) = (init_OK ? a=f1(n): 0, f2(x));
plot for [n=1:3] a=f3(x,n);
         undefined variable: x
:-|

Kalin.

[toc] | [prev] | [next] | [standalone]

#2558

Kalin Kozhuharov wrote:

> Hello Ethan,
> 
> Thank you for your time! I was trying to make self-contained example,
> but because it was simple enough you "optimized" it ;-) "Premature
> optimization is the root of all evil", said D. Knuth :-D
> 
> On Thursday, August 21, 2014 2:15:51 AM UTC+9, Ethan A Merritt wrote:
>>> Well, may be I oversimplified it, but the real thing is too complex
>>> anyway. What about this:
>>> 
>>> reset; f1(x) = (C=1.134324*x,0); f2(x) = sin(C*x); set xrange
>>> [0:3.14]; plot a=f1(1), f2(x) title "n=1", a=f1(2), f2(x) title
>>> "n=2" , a=f1(3), f2(x) title "n=3";
>> 
>> It is still not clear what you are aiming for.
>> The sequence of commands above does not make sense.
>> 
>> The function f1(x) is never plotted, so it plays no role except for
>> its side effect of setting C to a constant.
>>
> Well, yes I'd have to admit I am using side-effects here.
> And my code is turning to a small mess because I am trying to have
> gnuplot do statistics in-line, instead of preprocessing the data.
> 
> My current code is here:
> https://github.com/thinrope/fixed_sensor_visualization/blob/master/timeplot_all.gpl#L26
> 
> plot for [ID in IDs] "tmp/" . ID . ".data" u 1:2 axes x1y2 w l lw 2
> title graph_title(ID);
> 
> What I wanted to do is add another SMA filter on each datafile (or 2
> with different window) and still have it plot similar graph.
> 
> so having defined SMA as in
> https://github.com/thinrope/fixed_sensor_visualization/blob/master/timeplot.gpl#L46
> I was planning to write something along:
> 
> plot for [ID in IDs] a=sma_24_init(0), "tmp/" . ID . ".data" u
> 1:(sma_24($2)) axes x1y2 w l lw 2 title graph_title(ID);

I think that statement would work if you just remove the unnecessary comma
after the definition a=<something>.
Doesn't it?

	Ethan




> IDs are list of numbers coming form the environment/Makefile, but can
> be strings.
> 
> May be possible with `set table` and `do for [ID in IDs]{init; plot;}`
> but I'd rather avoid it.
> 
> 
> Alternatively, something along
> 
> init_OK=0; f3(x,n) = (init_OK ? a=f1(n): 0, f2(x));
> plot for [n=1:3] a=f3(x,n);
>          undefined variable: x
> :-|
> 
> Kalin.

[toc] | [prev] | [next] | [standalone]

#2559

On Thursday, August 21, 2014 9:50:52 AM UTC+9, Ethan A Merritt wrote:
>> plot for [ID in IDs] a=sma_24_init(0), "tmp/" . ID . ".data" u
>> 1:(sma_24($2)) axes x1y2 w l lw 2 title graph_title(ID);
> I think that statement would work if you just remove the unnecessary comma
> after the definition a=<something>.
> 
> Doesn't it?
>
LoL, LoL, LoL! Yes it works now. It was _only_ the extra coma that broke it...
Thank you very much.

The data is now better visualized in http://gamma.tar.bz/nGeigies/

Kalin.

[toc] | [prev] | [standalone]

Back to top | Article view | comp.graphics.apps.gnuplot

csiph-web