forums.silverfrost.com Forum Index forums.silverfrost.com
Welcome to the Silverfrost forums
 
 FAQFAQ   SearchSearch   MemberlistMemberlist   UsergroupsUsergroups   RegisterRegister 
 ProfileProfile   Log in to check your private messagesLog in to check your private messages   Log inLog in 

Request for extension
Goto page 1, 2  Next
 
Post new topic   Reply to topic    forums.silverfrost.com Forum Index -> ClearWin+
View previous topic :: View next topic  
Author Message
DanRRight



Joined: 10 Mar 2008
Posts: 2802
Location: South Pole, Antarctica

PostPosted: Wed Mar 15, 2023 3:07 pm    Post subject: Request for extension Reply with quote

Currently %PL allows to plot only 9 curves on the same graph.
As an example the plot looks like this:


This approximately covers 1/2 of evolution i intended to plot.

It is obvious that this limitation on number of curves is too restrictive. Would it be possible to triple or at least double that?
Back to top
View user's profile Send private message
Kenneth_Smith



Joined: 18 May 2012
Posts: 696
Location: Hamilton, Lanarkshire, Scotland.

PostPosted: Wed Mar 15, 2023 9:24 pm    Post subject: Reply with quote

Dan, take a look at the %pl[stacked] option.
Back to top
View user's profile Send private message Visit poster's website
DanRRight



Joined: 10 Mar 2008
Posts: 2802
Location: South Pole, Antarctica

PostPosted: Thu Mar 16, 2023 4:20 am    Post subject: Reply with quote

I know this option. It is extremely cumbersome way to deal with the multi-plot actually resembling a hack. Even those ways we discussed lately

http://forums.silverfrost.com/viewtopic.php?t=4547&highlight=multiplot

and which plot nicely

turn out to be inconvenient after some experience with them. I stopped using them. I am asking now about just the standard plotting option to extend it a bit. It also is not ideal, but it's straightforward, simple and stable like all plotting should be.

As a side note, i think that the ideal, shortest, simplest and ultimately intuitive and obvious plotting method of multiple curves should look exactly like this

Code:
I=winio@('%pl[file=settings.set, independent]', n,m,x,y)


That's all, nothing else, this is complete line. All adjustments should be inside the settings.set file.
Here n(m) is usual dimension of each of m plotted curves, and x and y are two-dimensional arrays x(n,m) and y(n,m) with n entries for each of m curves. The number of curves could be here infinite without any difficulties, conflicts or changes

How far is current standard method from ideal? Given its not punishing behavior it is pretty close to ideal. Currently i also use just one line to plot but it is long line like this

Code:
I=winio@('%pl[file=settings9.set, scale=linear, independent,n_graphs=9]', n, x(:,1),y(:,1),  x(:,2),y(:,2), x(:,3),y(:,3), x(:,4),y(:,4), x(:,5),y(:,5), x(:,6),y(:,6), x(:,7),y(:,7), x(:,8),y(:,8), x(:,9),y(:,9))


and even if i plot only say 3 curves all works OK. If extended to 20-30 i will just add more entries, it will be still universal without changes in number of plots in the source code
Back to top
View user's profile Send private message
PaulLaidler
Site Admin


Joined: 21 Feb 2005
Posts: 7876
Location: Salford, UK

PostPosted: Thu Mar 16, 2023 9:05 am    Post subject: Reply with quote

In the following sample code a callback function is used to make a single plot by caling draw_curvesD@. This call could be repeated for other sets of data.

Code:

      WINAPP
      MODULE mymod
      USE mswin
      INTEGER,PARAMETER::N=11
      DOUBLE PRECISION x(N),y(N)
      CONTAINS
      INTEGER FUNCTION cb()
      DOUBLE PRECISION xx(5),yy(5)
      CALL set_plot_mode@(1)
      CALL set_line_width@(2)
!     CALL set_line_style@(PS_DOT)
!     CALL draw_polylineD@(x,y,N,RGB@(0,0,255))   
      CALL draw_curvesD@(x,y,N,RGB@(0,0,255))   
!     CALL draw_beziersD@(x,y,N,RGB@(0,0,255))   
      CALL draw_line_betweenD@(0d0, 0d0, x(N), y(N), 255)
      CALL draw_line_betweenD@(0d0, 0d0, x(1), y(1), 255)
      CALL draw_rectangleD@(0.1d0, 0.1d0, 0.3d0, 0.15d0, RGB@(0,255,0))
      CALL draw_filled_ellipseD@(-0.1d0, 0.15d0, 0.04d0, 0.02d0, RGB@(0,255,0))
      CALL draw_charactersD@("Caption", 0.1d0, 0.2d0, 0)
      xx = (/0d0, 0.05d0, 0.05d0, 0d0, 0d0/)
      yy = (/0d0, 0d0, 0.035d0, 0.035d0, 0d0/)
      CALL draw_filled_polygonD@(xx,yy,5,RGB@(0,0,255))
      CALL draw_symbolsD@(x,y,N,11,5,255)
      CALL set_plot_mode@(0)
      cb = 2
      END FUNCTION
      END MODULE mymod
     
      PROGRAM main
      USE mymod
      INTEGER i
      DO i=1,N
        x(i)=0.1d0*(i-6)
        y(i)=x(i)*x(i)
      ENDDO
      i=winio@('%ca[Quadratic]%pv&')
      CALL winop@("%pl[title=Graph]") !graph title
      CALL winop@("%pl[y_max=0.25]")  !maximum y value on axis
      CALL winop@("%pl[x_array]")     !the data provides other max and min values
      CALL winop@("%pl[smoothing=5]")
      CALL winop@("%pl[link=user]")   !all lines and symbols drawn via set_plot_mode@
      CALL winop@("%pl[frame]")
      i=winio@('%^pl&',400,250,N,x,y,cb)
      i=winio@('%sf%ff%nl%cn%tt[OK]')
      END
Back to top
View user's profile Send private message AIM Address
DanRRight



Joined: 10 Mar 2008
Posts: 2802
Location: South Pole, Antarctica

PostPosted: Thu Mar 16, 2023 10:41 am    Post subject: Reply with quote

Paul,

This is exactly how i made (with the help of Ken) this nice colorful plot above. It took me several days by the way. It is serving as one of my general purpose plots in the program. But what i found when tried to duplicate it and create similar plot for other part of the same program - this required me more efforts than it worth to spend time on it because to make it work flawlessly definitely will need some debug time and such code is tens and tens of lines long.

So i ended with the standard plotting way (first my plot above) which if needed could be made similarly colorful via Settings (excluding these nice Legend numbering at the right hand side of the plot but that would take even more time)

This is why i ask to just to extend the standard %PL plotting to handle 20-30 graphs instead of 9 currently. I hope this would not require much efforts. I do plots like that for different tasks happens sometimes few times per week.
Back to top
View user's profile Send private message
PaulLaidler
Site Admin


Joined: 21 Feb 2005
Posts: 7876
Location: Salford, UK

PostPosted: Thu Mar 16, 2023 12:23 pm    Post subject: Reply with quote

Here is another sample. In this case there are 20 plots and the y-data is stored in a 2-dimensional array.

Code:
      WINAPP
      USE clrwin
      INTEGER M,N,i
      PARAMETER(N=11,M=20)
      REAL*8 xm,x(N),y(N,M)
      CHARACTER(len=24) str
      DO i=1,N
        xx = 0.1d0*(i-1)
        x(i) = xx
        xm = 0.1d0
! Replace this DO loop with data for each plot...
        DO j=0,M-1
         y(i,j+1) = xm*xx
         xm = xm + 0.1d0
        ENDDO
      ENDDO
      i=winio@("%ca[Multi-graph]&")
      CALL winop@("%pl[n_graphs=20]") !20=M
      CALL winop@("%pl[stacked]")
      CALL winop@("%pl[x_array]")
      CALL winop@("%pl[x_max=1.0]")
      CALL winop@("%pl[dx=0.1]")
      CALL winop@("%pl[smoothing=4]")
      !Various line colours going from blue to red...
      DO j=1,M
        i = 255*j/M
        write(str,"(a,z6.6,a)") "%pl[colour=#", RGB@(255-i,0,i), "]"
        CALL winop@(trim(str))
      ENDDO 
      i=winio@('%pl&',400,250,N,x,y)
      i=winio@('%ff%nl%cn%tt[OK]')
      END
Back to top
View user's profile Send private message AIM Address
DanRRight



Joined: 10 Mar 2008
Posts: 2802
Location: South Pole, Antarctica

PostPosted: Thu Mar 16, 2023 2:24 pm    Post subject: Reply with quote

Paul,
Thanks for an example. I tried before "stacked" option though some complexity stopped me. May be this was on my usual mischievous devilry saturdays Smile and i may look again at that but first my questions:

1) if stacked can do 20 graphs why regular %pl can not? What the difference?
2) does stacked handle "independent" X and Y similarly easy?
Back to top
View user's profile Send private message
PaulLaidler
Site Admin


Joined: 21 Feb 2005
Posts: 7876
Location: Salford, UK

PostPosted: Thu Mar 16, 2023 3:12 pm    Post subject: Reply with quote

1) Suffice to say that to provide the extension that you are asking for would not be quick and easy.

2) As I recall, only the Y values are stacked so you only get one set of X values.
Back to top
View user's profile Send private message AIM Address
PaulLaidler
Site Admin


Joined: 21 Feb 2005
Posts: 7876
Location: Salford, UK

PostPosted: Thu Mar 16, 2023 7:37 pm    Post subject: Reply with quote

I think that you could leave out [stacked] in the last sample above and present the data in a way that is not "stacked". The %pl line will have many arguments but you could, for example, use multiple Fortran continuation lines.
Back to top
View user's profile Send private message AIM Address
DanRRight



Joined: 10 Mar 2008
Posts: 2802
Location: South Pole, Antarctica

PostPosted: Fri Mar 17, 2023 6:15 pm    Post subject: Reply with quote

Do i correctly understand that something like this should work ?

Code:
I=winio@('%pl[file=settings9.set, scale=linear, stacked,independent,n_graphs=9]', n, &
& x(:,1),y(:,1),  x(:,2),y(:,2), x(:,3),y(:,3), x(:,4),y(:,4), &
& x(:,5),y(:,5), x(:,6),y(:,6), x(:,7),y(:,7), x(:,8),y(:,8), &
& x(:,9),y(:,9), x(:,10),y(:,10), x(:,11),y(:,11), x(:,12),y(:,12), &
& x(:,13),y(:,13), x(:,14),y(:,14), x(:,15),y(:,15), x(:,16),y(:,16), &
& x(:,17),y(:,17), x(:,18),y(:,18)......)
Back to top
View user's profile Send private message
PaulLaidler
Site Admin


Joined: 21 Feb 2005
Posts: 7876
Location: Salford, UK

PostPosted: Fri Mar 17, 2023 7:51 pm    Post subject: Reply with quote

Yes. I am suggesting that it may well work but you still need

Code:
CALL winop@("%pl[n_graphs=20]")


with 20 replaced by your value at the start.
Back to top
View user's profile Send private message AIM Address
DanRRight



Joined: 10 Mar 2008
Posts: 2802
Location: South Pole, Antarctica

PostPosted: Fri Mar 17, 2023 9:48 pm    Post subject: Reply with quote

It only plots in my case 1 graph.

So i returned to your demo and made some changes.
Here is the same code above i modified (see arrows in the text), this unfortunately does not plot anything. For simplicity I made it to plot just two graphs. And also for simplicity i made number of points equal in both plots

Did i made a mistake or three? Clearly i did, as at least this method (stacked+independent) plots something in my larger code
Code:
      WINAPP
      USE clrwin
      INTEGER M,N,i
      PARAMETER(N=11,M=2)          ! <-- changed M to 2
      REAL*8 xm,x(N,M),y(N,M)          !<-- changed X
      INTEGER :: Np(M)                ! <-- added Np
      CHARACTER(len=24) str

! Replace this DO loop with data for each plot...
      DO j=0,M-1
        DO i=1,N
          xx = 0.1d0*(i-1)
          x(i,j+1) = xx
          xm = 0.1d0
          y(i,j+1) = xm*xx
          xm = xm + 0.1d0
        ENDDO
      ENDDO
       
      Np(:)=N                             ! <-- added

      i=winio@("%ca[Multi-graph]&")
      CALL winop@("%pl[n_graphs=20]") !20=M
      CALL winop@("%pl[stacked]")
      CALL winop@("%pl[x_array]")
      CALL winop@("%pl[x_max=1.0]")
      CALL winop@("%pl[dx=0.1]")
      CALL winop@("%pl[smoothing=4]")
      CALL winop@("%pl[independent]")             !        <-- added
      !Various line colours going from blue to red...
      DO j=1,M
        i = 255*j/M
        write(str,"(a,z6.6,a)") "%pl[colour=#", RGB@(255-i,0,i), "]"
        CALL winop@(trim(str))
      ENDDO
      i=winio@('%pl&',400,250,Np,x(:,1),y(:,1),x(:,2),y(:,2))   ! <-- changed
      i=winio@('%ff%nl%cn%tt[OK]')
      END
Back to top
View user's profile Send private message
PaulLaidler
Site Admin


Joined: 21 Feb 2005
Posts: 7876
Location: Salford, UK

PostPosted: Fri Mar 17, 2023 10:28 pm    Post subject: Reply with quote

The number of graphs n_graphs is not 20 in your case. Also I suggested that you do not use [stacked]. The data presentation is probably wrong as well. That's all that I will contribute on this subject.
Back to top
View user's profile Send private message AIM Address
DanRRight



Joined: 10 Mar 2008
Posts: 2802
Location: South Pole, Antarctica

PostPosted: Fri Mar 17, 2023 11:55 pm    Post subject: Reply with quote

Thanks for finding my bugs in these demos and discussion which may shed the light on limiting amount of plots %pl can do

Then here is my last question: does %pl have 9 graphs limit on standard way of plotting or it was lifted lately (we discussed to take off this limit couple years ago)?

May be somebody here will look at that and find something else i missed.

Conclusion is pretty strange so far:
1) Looks like with Stacked it complains and swearing but plots 12 graphs ( 6 of which upside down ? May be still some other mistake is still here)

2) But surprisingly without staked it somehow manages to plot 12 ...which indicates that my original program of this post made with the standard %pl method exactly like in the last demo (and which refused to plot more than 9 ) has some error (??).
I use latest compiler 8.95

Here is the final modified Paul's demo code which supposed not to work for more than 9 graphs (and it indeed does not for more than 9 in my larger code why i started this thread and i do there exactly like is shown here) but it somehow unexplainable works in this demo

Code:
      WINAPP
      USE clrwin
      INTEGER M,N,i
      PARAMETER(N=11,M=12)
      REAL*8 xm,x(N,M),y(N,M)         !<-- changed X
      INTEGER :: Np(M)                ! <-- added Np
      CHARACTER(len=24) str

! Replace this DO loop with data for each plot...
      DO j=0,M-1                                       <-- changed a bit
          xm = 0.1d0 * (j+1)
        DO i=1,N
          xx = 0.1d0*(i-1)
          x(i,j+1) = xx

          y(i,j+1) = xm*xx
          xm = xm + 0.1d0
        ENDDO
      ENDDO
       
      Np(:)=N                             ! <-- added

      i=winio@("%ca[Multi-graph]&")
      CALL winop@("%pl[n_graphs=12]") 
!      CALL winop@("%pl[stacked]")
      CALL winop@("%pl[x_array]")
      CALL winop@("%pl[x_max=1.0]")
      CALL winop@("%pl[dx=0.1]")
      CALL winop@("%pl[smoothing=4]")
      CALL winop@("%pl[independent]")             !        <-- added
      !Various line colours going from blue to red...
      DO j=1,M
        i = 255*j/M
        write(str,"(a,z6.6,a)") "%pl[colour=#", RGB@(255-i,0,i), "]"
        CALL winop@(trim(str))
      ENDDO
      i=winio@('%pl&',400,250,Np, x(:,1),y(:,1), x(:,2),y(:,2), x(:,3),y(:,3), x(:,4),y(:,4),&
        & x(:,5),y(:,5), x(:,6),y(:,6), x(:,7),y(:,7), x(:,8),y(:,8), x(:,9),y(:,9), x(:,10),y(:,10),&
        & x(:,11),y(:,11), x(:,12),y(:,12) )


      i=winio@('%ff%nl%cn%tt[OK]')
      END
Back to top
View user's profile Send private message
PaulLaidler
Site Admin


Joined: 21 Feb 2005
Posts: 7876
Location: Salford, UK

PostPosted: Sat Mar 18, 2023 9:00 am    Post subject: Reply with quote

Yes you can have more than 10 plots (with or wihout [stacked]). Simply set n_graphs to the required number of graphs.

I have looked at the code for %pl and [stacked] with [independent] means both the x and the y arrays are stacked. But Dan, if you find [stacked] difficult to use why not go back to the old way of providing the data?
Back to top
View user's profile Send private message AIM Address
Display posts from previous:   
Post new topic   Reply to topic    forums.silverfrost.com Forum Index -> ClearWin+ All times are GMT + 1 Hour
Goto page 1, 2  Next
Page 1 of 2

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum


Powered by phpBB © 2001, 2005 phpBB Group