## Copyright (C) 1996, 1997 John W. Eaton
##
## This file is part of Octave.
##
## Octave is free software; you can redistribute it and/or modify it
## under the terms of the GNU General Public License as published by
## the Free Software Foundation; either version 2, or (at your option)
## any later version.
##
## Octave is distributed in the hope that it will be useful, but
## WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
## General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with Octave; see the file COPYING.  If not, write to the Free
## Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
## 02110-1301, USA.

## -*- texinfo -*-
## @deftypefn {Function File} {} axis (@var{limits})
## Set axis limits for plots.
##
## The argument @var{limits} should be a 2, 4, or 6 element vector.  The
## first and second elements specify the lower and upper limits for the x
## axis.  The third and fourth specify the limits for the y axis, and the
## fifth and sixth specify the limits for the z axis.
##
## Without any arguments, @code{axis} turns autoscaling on.  
##
## With one output argument, @code{x=axis} returns the current axes 
## (this is not yet implemented for automatic axes).
##
## The vector argument specifying limits is optional, and additional
## string arguments may be used to specify various axis properties.  For
## example,
##
## @example
## axis ([1, 2, 3, 4], "square");
## @end example
##
## @noindent
## forces a square aspect ratio, and
##
## @example
## axis ("labely", "tic");
## @end example
##
## @noindent
## turns tic marks on for all axes and tic mark labels on for the y-axis
## only.
##
## @noindent
## The following options control the aspect ratio of the axes.
##
## @table @code
## @item "square"
## Force a square aspect ratio.
## @item "equal"
## Force x distance to equal y-distance.
## @item "normal"
## Restore the balance.
## @end table
##
## @noindent
## The following options control the way axis limits are interpreted.
##
## @table @code
## @item "auto" 
## Set the specified axes to have nice limits around the data
## or all if no axes are specified.
## @item "manual" 
## Fix the current axes limits.
## @item "tight"
## Fix axes to the limits of the data (not implemented).
## @end table
##
## @noindent
## The option @code{"image"} is equivalent to @code{"tight"} and
## @code{"equal"}.
##
## @noindent
## The following options affect the appearance of tic marks.
##
## @table @code
## @item "on" 
## Turn tic marks and labels on for all axes.
## @item "off"
## Turn tic marks off for all axes.
## @item "tic[xyz]"
## Turn tic marks on for all axes, or turn them on for the
## specified axes and off for the remainder.
## @item "label[xyz]"
## Turn tic labels on for all axes, or turn them on for the 
## specified axes and off for the remainder.
## @item "nolabel"
## Turn tic labels off for all axes.
## @end table
## Note, if there are no tic marks for an axis, there can be no labels.
##
## @noindent
## The following options affect the direction of increasing values on
## the axes.
##
## @table @code
## @item "ij"
## Reverse y-axis, so lower values are nearer the top.
## @item "xy" 
## Restore y-axis, so higher values are nearer the top. 
## @end table
## 
## @end deftypefn

## Author: jwe

function curr_axis = axis (ax, varargin)

  ## This may not be correct if someone has used the gnuplot interface
  ## directly...

  global __current_axis__ = [-10, 10, -10, 10];

  ## To return curr_axis properly, octave needs to take control of scaling.
  ## It isn't hard to compute good axis limits:
  ##   scale = 10 ^ floor (log10 (max - min) - 1);
  ##   r = scale * [floor (min / scale), ceil (max / scale)];
  ## However, with axis("manual") there is little need to know the current
  ## limits.

  if (nargin == 0)
    if (nargout == 0)
      __gnuplot_raw__ ("set autoscale;\n");
    else
      curr_axis = __current_axis__;
    endif

  elseif (isstr (ax))
    ax = tolower (ax);
    len = length (ax);

    ## 'matrix mode' to reverse the y-axis
    if (strcmp (ax, "ij"))
      __gnuplot_raw__ ("set yrange [] reverse;\n"); 
    elseif (strcmp (ax, "xy"))
      __gnuplot_raw__ ("set yrange [] noreverse;\n");

      ## aspect ratio
    elseif (strcmp (ax, "image"))
      __gnuplot_raw__ ("set size ratio -1;\n"); 
      __gnuplot_raw__ ("set autoscale;\n"); ## XXX FIXME XXX should be the same as "tight"
    elseif (strcmp (ax, "equal"))
      __gnuplot_raw__ ("set size ratio -1;\n");
    elseif (strcmp (ax, "square"))
      __gnuplot_raw__ ("set size ratio 1;\n");
    elseif (strcmp (ax, "normal"))
      __gnuplot_raw__ ("set size noratio;\n");


      ## axis limits
    elseif (len >= 4 && strcmp (ax(1:4), "auto"))
      if (len > 4)
      	__gnuplot_raw__ (sprintf ("set autoscale %s;\n", ax(5:len)));
      else
	__gnuplot_raw__ ("set autoscale;\n");
      endif
    elseif (strcmp (ax, "manual"))
      ## fixes the axis limits, like axis(axis) should;
      __gnuplot_raw__ ("set xrange [] writeback;\n");
      __gnuplot_raw__ ("set yrange [] writeback;\n");
      __gnuplot_raw__ ("set zrange [] writeback;\n");
      ## XXX FIXME XXX if writeback were set in plot, no need to replot here.
      ## No semicolon (see replot.m).
      __gnuplot_replot__
      __gnuplot_raw__ ("set noautoscale x;\n");
      __gnuplot_raw__ ("set noautoscale y;\n");
      __gnuplot_raw__ ("set noautoscale z;\n");
    elseif (strcmp (ax, "tight"))
      ## XXX FIXME XXX if tight, plot must set ranges to limits of the
      ## all the data on the current plot, even if from a previous call.
      ## Instead, just let gnuplot do as it likes.
      __gnuplot_raw__ ("set autoscale;\n");


      ## tic marks
    elseif (strcmp (ax, "on"))
      __gnuplot_raw__ ("set xtics;\n");
      __gnuplot_raw__ ("set ytics;\n");
      __gnuplot_raw__ ("set ztics;\n");
      __gnuplot_raw__ ("set format;\n");
    elseif (strcmp (ax, "off"))
      __gnuplot_raw__ ("set noxtics;\n");
      __gnuplot_raw__ ("set noytics;\n");
      __gnuplot_raw__ ("set noztics;\n");
    elseif (strcmp (ax, "tic"))
      __gnuplot_raw__ ("set xtics;\n");
      __gnuplot_raw__ ("set ytics;\n");
      __gnuplot_raw__ ("set ztics;\n");
    elseif (len > 3 && strcmp (ax(1:3), "tic"))
      if (any (ax == "x"))
	__gnuplot_raw__ ("set xtics;\n");
      else
	__gnuplot_raw__ ("set noxtics;\n");
      endif
      if (any (ax == "y"))
	__gnuplot_raw__ ("set ytics;\n");
      else
	__gnuplot_raw__ ("set noytics;\n");
      endif
      if (any (ax == "z"))
	__gnuplot_raw__ ("set ztics;\n");
      else
	__gnuplot_raw__ ("set noztics;\n");
      endif
    elseif (strcmp (ax, "label"))
      __gnuplot_raw__ ("set format;\n");
    elseif (strcmp (ax, "nolabel"))
      __gnuplot_raw__ ("set format \"\\0\";\n");
    elseif (len > 5 && strcmp (ax(1:5), "label"))
      if (any (ax == "x"))
	__gnuplot_raw__ ("set format x;\n");
      else
	__gnuplot_raw__ ("set format x \"\\0\";\n");
      endif
      if (any (ax == "y"))
	__gnuplot_raw__ ("set format y;\n");
      else
	__gnuplot_raw__ ("set format y \"\\0\";\n");
      endif
      if (any (ax == "z"))
	__gnuplot_raw__ ("set format z;\n");
      else
	__gnuplot_raw__ ("set format z \"\\0\";\n");
      endif

    else
      warning ("unknown axis option '%s'", ax);
    endif

  elseif (isvector (ax))

    len = length (ax);

    if (len != 2 && len != 4 && len != 6)
      error ("axis: expecting vector with 2, 4, or 6 elements");
    endif

    __current_axis__ = reshape (ax, 1, len);

    if (len > 1)
      __gnuplot_raw__ (sprintf ("set xrange [%g:%g];\n", ax(1), ax(2)));
    endif

    if (len > 3)
      __gnuplot_raw__ (sprintf ("set yrange [%g:%g];\n", ax(3), ax(4)));
    endif

    if (len > 5)
      __gnuplot_raw__ (sprintf ("set zrange [%g:%g];\n", ax(5), ax(6)));
    endif

  else
    error ("axis: expecting no args, or a vector with 2, 4, or 6 elements");
  endif

  if (nargin > 1)
    axis (varargin{:});
  elseif (automatic_replot)
    ## No semicolon (see replot.m).
    __gnuplot_replot__
  endif
endfunction

%!demo
%! t=0:0.01:2*pi; x=sin(t);
%!
%! subplot(221);    title("normal plot");
%! plot(t, x, ";;");
%!
%! subplot(222);    title("square plot");
%! axis("square");  plot(t, x, ";;");
%!
%! subplot(223);    title("equal plot");
%! axis("equal");   plot(t, x, ";;");
%! 
%! subplot(224);    title("normal plot again");
%! axis("normal");  plot(t, x, ";;");

%!demo
%! t=0:0.01:2*pi; x=sin(t);
%!
%! subplot(121);   title("ij plot");
%! axis("ij");     plot(t, x, ";;");
%!
%! subplot(122);   title("xy plot");
%! axis("xy");     plot(t, x, ";;");

%!demo
%! t=0:0.01:2*pi; x=sin(t);
%!
%! subplot(331);   title("x tics & labels");
%! axis("ticx");   plot(t, x, ";;");
%!
%! subplot(332);   title("y tics & labels");
%! axis("ticy");   plot(t, x, ";;");
%!
%! subplot(334);     title("x & y tics, x labels");
%! axis("labelx","tic");   plot(t, x, ";;");
%!
%! subplot(335);     title("x & y tics, y labels");
%! axis("labely","tic");   plot(t, x, ";;");
%!
%! subplot(337);     title("x tics, no labels");
%! axis("nolabel","ticx");   plot(t, x, ";;");
%!
%! subplot(338);     title("y tics, no labels");
%! axis("nolabel","ticy");   plot(t, x, ";;");
%!
%! subplot(333);     title("no tics or labels");
%! axis("off");    plot(t, x, ";;");
%!
%! subplot(336);     title("all tics but no labels");
%! axis("nolabel","tic");    plot(t, x, ";;");
%!
%! subplot(339);     title("all tics & labels");
%! axis("on");       plot(t, x, ";;");

%!demo
%! t=0:0.01:2*pi; x=sin(t);
%!
%! subplot(321);    title("axes at [0 3 0 1]")
%! axis([0,3,0,1]); plot(t, x, ";;");
%!
%! subplot(322);    title("auto");
%! axis("auto");    plot(t, x, ";;");
%!
%! subplot(323);    title("manual");
%! plot(t, x, ";sine [0:2pi];"); hold on;
%! axis("manual");
%! plot(-3:3,-3:3, ";line (-3,-3)->(3,3);"); hold off;
%!
%! subplot(324);    title("axes at [0 3 0 1], then autox");
%! axis([0,3,0,1]); axis("autox");
%! plot(t, x, ";sine [0:2pi];");
%!
%! subplot(325);    title("axes at [3 6 0 1], then autoy");
%! axis([3,6,0,1]); axis("autoy");
%! plot(t, x, ";sine [0:2p];");
%!
%! subplot(326);    title("tight");
%! axis("tight");   plot(t, x, ";;");
%! % The last plot should not have any whitespace outside the data
%! % limits, but "tight" isn't implemented yet.