From: Joachim G. <joa...@ip...> - 2004-05-13 18:11:07
Attachments:
x09f.f
|
Hello, after some days of debugging (frustration included), I think that the contour plotting routine used for plcon2 in the fortran interface has a bug in it. The attached example shows the problem and is taken from the expample x09f.f of version 5.3.0 and is the following: when you plot a function z(r) on the space x(r,theta),y(r,theta), then it depends on the storage of the grid values for plcon2 whether you see the contour lines or not. Specifically, x(r,theta) is stored in xg(i,j) (y and z the same way); if i corresponds to r, then you will get nice nested contours, if i corresponds to theta (thus j corresponds to r) then no contours are seen. If there is a dependence of z on r and theta, you will not notice this asymmetry in the coordinates. Has anyone an idea how to fix it, so that it does not depend on the order of the coordinates that you might or might not see contour lines? By the way, I love plplot!!! Really great! Joachim -- -----------------------------------------------@\\ Joachim Geiger @\\ Institute for Plasma Physics, Euratom Association\\\ Department E3 \\\__ Wendelsteinstr. 1 ( \\ \ D-17491 Greifswald, Germany \_ \\(_ Tel.: 03834/88-2327 ) \\ \ e-mail: joa...@ip... \ - | --------------------------------------------------- \____/ |
From: Alan W. I. <ir...@be...> - 2004-05-14 05:44:20
|
On 2004-05-13 20:11+0200 Joachim Geiger wrote: > (for version 5.3.0) > when you plot a function z(r) on the space x(r,theta),y(r,theta), > then it depends on the storage of the grid values for plcon2 > whether you see the contour lines or not. Specifically, > x(r,theta) is stored in xg(i,j) (y and z the same way); if > i corresponds to r, then you will get nice nested contours, > if i corresponds to theta (thus j corresponds to r) then > no contours are seen. To solve other contouring problems that have been found recently, there has been a replacement of our contouring algorithm since version 5.3.0, and it seems to have fixed the problem you have reported. To be specific, for the latest PLplot (from cvs) your fortran code produced two identical looking plots with all contours plotted. Could you please try the latest cvs snapshot version of PLplot from http://plplot.sourceforge.net/cvs-tarball/ to see whether that fixes the problem for you? Alan __________________________ Alan W. Irwin email: ir...@be... phone: 250-727-2902 Astronomical research affiliation with Department of Physics and Astronomy, University of Victoria (astrowww.phys.uvic.ca). Programming affiliations with the PLplot scientific plotting software package (plplot.org), the Yorick front-end to PLplot (yplot.sf.net), the Loads of Linux Links project (loll.sf.net), and the Linux Brochure Project (lbproject.sf.net). __________________________ Linux-powered Science __________________________ |
From: Joachim G. <joa...@ip...> - 2004-05-14 10:56:23
|
Thanks, this indeed solves the problem. Joachim Alan W. Irwin wrote: >On 2004-05-13 20:11+0200 Joachim Geiger wrote: > > > >>(for version 5.3.0) >>when you plot a function z(r) on the space x(r,theta),y(r,theta), >>then it depends on the storage of the grid values for plcon2 >>whether you see the contour lines or not. Specifically, >>x(r,theta) is stored in xg(i,j) (y and z the same way); if >>i corresponds to r, then you will get nice nested contours, >>if i corresponds to theta (thus j corresponds to r) then >>no contours are seen. >> >> > >To solve other contouring problems that have been found recently, there has >been a replacement of our contouring algorithm since version 5.3.0, and it >seems to have fixed the problem you have reported. To be specific, for the >latest PLplot (from cvs) your fortran code produced two identical looking >plots with all contours plotted. > >Could you please try the latest cvs snapshot version of PLplot from >http://plplot.sourceforge.net/cvs-tarball/ to see whether that fixes >the problem for you? > >Alan >__________________________ >Alan W. Irwin >email: ir...@be... >phone: 250-727-2902 > >Astronomical research affiliation with Department of Physics and Astronomy, >University of Victoria (astrowww.phys.uvic.ca). > >Programming affiliations with the PLplot scientific plotting software >package (plplot.org), the Yorick front-end to PLplot (yplot.sf.net), the >Loads of Linux Links project (loll.sf.net), and the Linux Brochure Project >(lbproject.sf.net). >__________________________ > >Linux-powered Science >__________________________ > > >------------------------------------------------------- >This SF.Net email is sponsored by: SourceForge.net Broadband >Sign-up now for SourceForge Broadband and get the fastest >6.0/768 connection for only $19.95/mo for the first 3 months! >http://ads.osdn.com/?ad_id=2562&alloc_id=6184&op=click >_______________________________________________ >Plplot-general mailing list >Plp...@li... >https://lists.sourceforge.net/lists/listinfo/plplot-general > > > -- -----------------------------------------------@\\ Joachim Geiger @\\ Institute for Plasma Physics, Euratom Association\\\ Department E3 \\\__ Wendelsteinstr. 1 ( \\ \ D-17491 Greifswald, Germany \_ \\(_ Tel.: 03834/88-2327 ) \\ \ e-mail: joa...@ip... \ - | --------------------------------------------------- \____/ |
From: Alan W. I. <ir...@be...> - 2004-05-14 14:03:32
|
On 2004-05-14 12:56+0200 Joachim Geiger wrote: > Thanks, this indeed solves the problem. That's great that Andrew Ross's new, contouring algorithm solved your problem. Alan __________________________ Alan W. Irwin email: ir...@be... phone: 250-727-2902 Astronomical research affiliation with Department of Physics and Astronomy, University of Victoria (astrowww.phys.uvic.ca). Programming affiliations with the PLplot scientific plotting software package (plplot.org), the Yorick front-end to PLplot (yplot.sf.net), the Loads of Linux Links project (loll.sf.net), and the Linux Brochure Project (lbproject.sf.net). __________________________ Linux-powered Science __________________________ |
From: Andrew R. <and...@us...> - 2004-05-14 18:07:27
|
On Fri, May 14, 2004 at 07:03:29AM -0700, Alan Irwin wrote: > On 2004-05-14 12:56+0200 Joachim Geiger wrote: > > > Thanks, this indeed solves the problem. > > That's great that Andrew Ross's new, contouring algorithm solved your > problem. I'm pleased to hear it too. Please let me know if you find any other problems with the new scheme. It seems to be a definite improvement on the old scheme but there may still be some teething problems. Andrew |