ftp.nice.ch/pub/next/unix/developer/plplot.3.0.s.tar.gz#/plplot/src/plaxes.c

This is plaxes.c in view mode; [Download] [Up]

/* This draws a box around the current viewport. XOPT and YOPT are */
/* character strings which define the box as follows: */

/* A: Draw axis (X is horizontal line Y=0, Y is vertical line X=0) */
/* B: Draw bottom (X) or left (Y) edge of frame */
/* C: Draw top (X) or right (Y) edge of frame */
/* G: Draws a grid at the major tick interval */
/* I: Inverts tick marks */
/* L: Logarithmic axes, major ticks at decades, minor ticks at units */
/* N: Write numeric label at conventional location */
/* M: Write numeric label at unconventional location */
/* T: Draw major tick marks */
/* S: Draw minor tick marks */
/* V: (for Y only) Label vertically */

/* xtick, ytick are the major tick intervals required, zero for  */
/*              automatic selection */
/* nxsub, nysub are the number of subtick intervals in a major tick */
/*              interval */

#include "plplot.h"
#include <stdio.h>
#include <math.h>
#ifdef PLSTDC
#include <string.h>
#else
extern char *strcat();
extern char *strcpy();
#endif

#define betw(c,a,b) ((a<=c && c<=b) || (b<=c && c<=a))

static double xlog[8] = {0.301030, 0.477121, 0.602060, 0.698970,
                        0.778151, 0.845098, 0.903090, 0.954243};

void plaxes(x0,y0,xopt,xtick,nxsub,yopt,ytick,nysub)
FLOAT x0,y0;
char *xopt, *yopt;
double xtick, ytick;
int nxsub, nysub;
{
      static char string[40];
      int lax,lbx,lcx,lgx,lix,llx,lmx,lnx,lsx,ltx;
      int lay,lby,lcy,lgy,liy,lly,lmy,lny,lsy,lty,lvy;
      int xmajor, xminor, ymajor, yminor, xmode, xprec, xscal;
      int xorigin, yorigin;
      int ymode, yprec, yscal;
      int i, i1x, i2x, i3x, i4x, i1y, i2y, i3y, i4y;
      int nxsub1, nysub1;
      int lxmin, lxmax, lymin, lymax;
      int pxmin, pxmax, pymin, pymax;
      int vppxmi, vppxma, vppymi, vppyma;
      int level;
      double xpmm, ypmm, defmaj, defmin, htmaj, htmin;
      double xtick1, ytick1, vpwxmi, vpwxma, vpwymi, vpwyma;
      double pos, tn, tp, temp;

      glev(&level);
      if (level<3) plexit("Please set up window before calling plaxes.");

      /* Open  the clip limits to the subpage limits */
      gclp(&lxmin,&lxmax,&lymin,&lymax);
      gphy(&pxmin,&pxmax,&pymin,&pymax);
      sclp(pxmin,pxmax,pymin,pymax);

      gvpp(&vppxmi,&vppxma,&vppymi,&vppyma);

      /* convert world coordinates to physical	*/
      xorigin = wcpcx(x0); yorigin = wcpcy(y0);

      /* Tick and subtick sizes in device coords */
      gpixmm(&xpmm,&ypmm);
      gmaj(&defmaj,&htmaj);
      gmin(&defmin,&htmin);

      xmajor=MAX(round(htmaj*ypmm),1);
      ymajor=MAX(round(htmaj*xpmm),1);
      xminor=MAX(round(htmin*ypmm),1);
      yminor=MAX(round(htmin*xpmm),1);

      xtick1=xtick;
      nxsub1=nxsub;
      ytick1=ytick;
      nysub1=nysub;

      lax=stsearch(xopt,'a');
      lbx=stsearch(xopt,'b');
      lcx=stsearch(xopt,'c');
      lgx=stsearch(xopt,'g');
      lix=stsearch(xopt,'i');
      llx=stsearch(xopt,'l');
      lmx=stsearch(xopt,'m');
      lnx=stsearch(xopt,'n');
      lsx=stsearch(xopt,'s');
      ltx=stsearch(xopt,'t');

      lay=stsearch(yopt,'a');
      lby=stsearch(yopt,'b');
      lcy=stsearch(yopt,'c');
      lgy=stsearch(yopt,'g');
      liy=stsearch(yopt,'i');
      lly=stsearch(yopt,'l');
      lmy=stsearch(yopt,'m');
      lny=stsearch(yopt,'n');
      lsy=stsearch(yopt,'s');
      lty=stsearch(yopt,'t');
      lvy=stsearch(yopt,'v');

      gvpw(&vpwxmi,&vpwxma,&vpwymi,&vpwyma);

      lax=lax && (vpwymi*vpwyma<0.0) && !llx;
      lay=lay && (vpwxmi*vpwxma<0.0) && !lly;
      if (llx) xtick1=1.0;
      if (lly) ytick1=1.0;
      if (ltx || lgx)
          pldtik(vpwxmi,vpwxma,&xtick1,&nxsub1,&xmode,&xprec,&xscal);
      if (lty || lgy)
          pldtik(vpwymi,vpwyma,&ytick1,&nysub1,&ymode,&yprec,&yscal);

      /* Set up tick variables */
       i1x=xminor;
       i2x=xminor;
       i3x=xmajor;
       i4x=xmajor;

       i1y=yminor;
       i2y=yminor;
       i3y=ymajor;
       i4y=ymajor;

      /* draw boxing */
	if (lbx) {
		movphy(vppxmi, vppymi);
		draphy(vppxma, vppymi);
	}
	if (lcx) {
		movphy(vppxmi, vppyma);
		draphy(vppxma, vppyma);
	}
	if (lby) {
		movphy(vppxmi, vppymi);
		draphy(vppxmi, vppyma);
	}
	if (lcy) {
		movphy(vppxma, vppymi);
		draphy(vppxma, vppyma);
	}

      /* Draw the horizontal axis */
        movphy(vppxmi,yorigin);
        if (ltx) {
          tp=xtick1*floor(vpwxmi/xtick1);
          for(;;) {
             tn=tp+xtick1;
             if (lsx) {
               if (llx) {
                 for(i=0; i<=7;i++){
                   temp=tp+xlog[i];
                   if (betw(temp,vpwxmi,vpwxma))
                        plxtik(wcpcx(temp),yorigin,i1x,i2x);
                 }
               }
               else {
                 for ( i=1;i<=nxsub1-1;i++) {
                   temp=tp+i*(tn-tp)/nxsub1;
                   if (betw(temp,vpwxmi,vpwxma))
                        plxtik(wcpcx(temp),yorigin,i1x,i2x);
                 }
               }
             }
             if (!betw(tn,vpwxmi,vpwxma))
               break;
             plxtik(wcpcx(tn),yorigin,i3x,i4x);
             tp=tn;
          }
        }
        draphy(vppxma,yorigin);
	
      /* Draw vertical axis */
      if (lby) {
        movphy(xorigin,vppyma);
        if (lty) {
          tp=ytick1*(floor(vpwyma/ytick1)+1);
          for(;;) {
             tn=tp-ytick1;
             if (lsy) {
               if (lly) {
                 for(i=7;i>=0;i--) {
                   temp=tn+xlog[i];
                   if (betw(temp,vpwymi,vpwyma))
                        plytik(xorigin,wcpcy(temp),i1y,i2y);
                 }
               }
               else {
                 for(i=nysub1-1;i>=1;i--) {
                   temp=tn+i*(tp-tn)/nysub1;
                   if (betw(temp,vpwymi,vpwyma))
                        plytik(xorigin,wcpcy(temp),i1y,i2y);
                 }
               }
             }
             if (!betw(tn,vpwymi,vpwyma))
               break;
             plytik(xorigin,wcpcy(tn),i3y,i4y);
             tp=tn;
          }
        }
        draphy(xorigin,vppymi);
      }

      /* Draw grid in x direction */
      if (lgx) {
        tp=xtick1*(1.+floor(vpwxmi/xtick1+.5));
        for(tn=tp; betw(tn+xtick1/2.,vpwxmi,vpwxma); tn+=xtick1)
           pljoin(tn,vpwymi,tn,vpwyma);
      }

      /* Draw grid in y direction */
      if (lgy) {
        tp=ytick1*(1.+floor(vpwymi/ytick1+.5));
        for(tn=tp; betw(tn+ytick1/2.,vpwymi,vpwyma); tn+=ytick1)
           pljoin(vpwxmi,tn,vpwxma,tn);
      }

      /* Write horizontal label(s) */
      if ((lmx || lnx) && ltx) {
        tp=xtick1*(1.+floor(vpwxmi/xtick1));
        for(tn=tp; betw(tn,vpwxmi,vpwxma); tn+=xtick1) {
          if (!llx) 
              plform(tn,xscal,xprec,string);
          else {
            sprintf(string,"10#u%d",round(tn));
          }
          pos=(tn-vpwxmi)/(vpwxma-vpwxmi);
          if (lnx) plmtex("b",1.5,pos,0.5,string);
          if (lmx) plmtex("t",1.5,pos,0.5,string);
        }
        if(!llx && xmode) {
          sprintf(string,"E%+d",xscal);
          if (lnx) plmtex("b",3.0,1.,0.5,string);
          if (lmx) plmtex("t",3.0,1.,0.5,string);
        }
      }

      /* Write vertical label(s) */
      if ((lmy || lny) && lty) {
        tp=ytick1*(1.+floor(vpwymi/ytick1));
        for(tn=tp; betw(tn,vpwymi,vpwyma); tn+=ytick1) {
          if (!lly)
              plform(tn,yscal,yprec,string);
          else {
            sprintf(string,"10#u%d",round(tn));
          }
          pos=(tn-vpwymi)/(vpwyma-vpwymi);
          if (lny) {
            if (lvy)
              plmtex("lv",0.5,pos,1.0,string);
            else
              plmtex("l",1.5,pos,0.5,string);
          }
          if (lmy) {
            if (lvy)
              plmtex("rv",0.5,pos,0.0,string);
            else
              plmtex("r",1.5,pos,0.5,string);
          }
        }
        if(!lly && ymode) {
          sprintf(string,"E%+d",yscal);
          if (lny) {
            if (lvy)
              plmtex("t",2.0,0.0,1.0,string);
            else
              plmtex("l",3.0,1.0,0.5,string);
          }
          if (lmy) {
            if (lvy)
              plmtex("t",2.0,1.0,0.0,string);
            else
              plmtex("r",3.0,1.0,0.5,string);
          }
        }
      }
      /* Restore the clip limits to viewport edge */
      sclp(lxmin,lxmax,lymin,lymax);
}

These are the contents of the former NiCE NeXT User Group NeXTSTEP/OpenStep software archive, currently hosted by Marcel Waldvogel and Netfuture.ch.