jwe/octave: scripts/plot/draw/private/__ezplot_

annotate scripts/plot/draw/private/ezplot.m @ 27924:1891570abac8

update Octave Project Developers copyright for the new year In files that have the "Octave Project Developers" copyright notice, update for 2020.

author	John W. Eaton <jwe@octave.org>
date	Mon, 06 Jan 2020 22:29:51 -0500
parents	b442ec6dda5c
children	bd51beb6205e

rev	line source
27924 1891570abac8 update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 27923 diff changeset	1 ## Copyright (C) 2007-2020 The Octave Project Developers
27923 b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 26671 diff changeset	2 ##
b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 26671 diff changeset	3 ## See the file COPYRIGHT.md in the top-level directory of this distribution
b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 26671 diff changeset	4 ## or <https://octave.org/COPYRIGHT.html/>.
b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 26671 diff changeset	5 ##
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	6 ##
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	7 ## This file is part of Octave.
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	8 ##
24534 194eb4bd202b maint: Update punctuation for GPL v3 license text. Rik <rik@octave.org> parents: 23220 diff changeset	9 ## Octave is free software: you can redistribute it and/or modify it
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	10 ## under the terms of the GNU General Public License as published by
24534 194eb4bd202b maint: Update punctuation for GPL v3 license text. Rik <rik@octave.org> parents: 23220 diff changeset	11 ## the Free Software Foundation, either version 3 of the License, or
22755 3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22351 diff changeset	12 ## (at your option) any later version.
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	13 ##
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	14 ## Octave is distributed in the hope that it will be useful, but
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	15 ## WITHOUT ANY WARRANTY; without even the implied warranty of
22755 3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22351 diff changeset	16 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22351 diff changeset	17 ## GNU General Public License for more details.
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	18 ##
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	19 ## You should have received a copy of the GNU General Public License
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	20 ## along with Octave; see the file COPYING. If not, see
24534 194eb4bd202b maint: Update punctuation for GPL v3 license text. Rik <rik@octave.org> parents: 23220 diff changeset	21 ## <https://www.gnu.org/licenses/>.
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	22
8812 7d48766c21a5 use consistent format for doc strings of internal functions John W. Eaton <jwe@octave.org> parents: 8046 diff changeset	23 ## -- texinfo --
20852 516bb87ea72e 2015 Code Sprint: remove class of function from docstring for all m-files. Rik <rik@octave.org> parents: 20711 diff changeset	24 ## @deftypefn {} {[@var{h}, @var{needusage}] =} __ezplot__ (@var{pltfunc}, @var{varargin})
8812 7d48766c21a5 use consistent format for doc strings of internal functions John W. Eaton <jwe@octave.org> parents: 8046 diff changeset	25 ## Undocumented internal function.
7d48766c21a5 use consistent format for doc strings of internal functions John W. Eaton <jwe@octave.org> parents: 8046 diff changeset	26 ## @end deftypefn
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	27
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	28 ## Overview: This function is the back-end for the 9 ez* plot functions.
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	29 ## As such, most of the function is actually dedicated to sorting
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	30 ## out the inputs and verifying that the particular ez* function
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	31 ## called was called correctly. The actual plotting occurs near
19593 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	32 ## the end in an unwind_protect block.
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	33
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	34 function [h, needusage] = __ezplot__ (pltfunc, varargin)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	35
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	36 ezfunc = ["ez" pltfunc];
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	37
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	38 [hax, varargin, nargin] = __plt_get_axis_arg__ (ezfunc, varargin{:});
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	39
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	40 ## Define outputs early in case of shorting out of function with return;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	41 h = [];
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	42 needusage = false;
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	43 if (nargin < 1)
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	44 needusage = true;
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	45 return;
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	46 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	47
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	48 iscontour = strncmp (pltfunc, "contour", 7);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	49
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	50 ## Defaults for ezplot
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	51 isplot = true;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	52 isplot3 = false;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	53 ispolar = false;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	54 nargs = 1;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	55 switch (pltfunc)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	56 case "plot"
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	57 ## defaults already set
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	58
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	59 case "plot3"
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	60 isplot = false;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	61 isplot3 = true;
19593 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	62
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	63 case "polar"
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	64 isplot = false;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	65 ispolar = true;
19593 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	66
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	67 otherwise
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	68 ## contour, mesh, surf plots
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	69 isplot = false;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	70 nargs = 2;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	71
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	72 endswitch
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	73
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	74 parametric = false;
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	75 fun = varargin{1};
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	76 if (ischar (fun))
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	77 if (exist (fun, "file") \|\| exist (fun, "builtin"))
20499 c0566df8cde3 Generate correct ezplot for 2-input functions (bug #46004). Rik <rik@octave.org> parents: 19697 diff changeset	78 fun = str2func (fun); # convert to function handle
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	79 else
20499 c0566df8cde3 Generate correct ezplot for 2-input functions (bug #46004). Rik <rik@octave.org> parents: 19697 diff changeset	80 fun = vectorize (inline (fun)); # convert to inline function
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	81 endif
20499 c0566df8cde3 Generate correct ezplot for 2-input functions (bug #46004). Rik <rik@octave.org> parents: 19697 diff changeset	82 endif
c0566df8cde3 Generate correct ezplot for 2-input functions (bug #46004). Rik <rik@octave.org> parents: 19697 diff changeset	83
c0566df8cde3 Generate correct ezplot for 2-input functions (bug #46004). Rik <rik@octave.org> parents: 19697 diff changeset	84 if (strcmp (typeinfo (fun), "inline function"))
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	85 argids = argnames (fun);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	86 if (isplot && length (argids) == 2)
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	87 nargs = 2;
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	88 elseif (numel (argids) != nargs)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	89 error ("%s: expecting a function of %d arguments", ezfunc, nargs);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	90 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	91 fun = vectorize (fun);
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	92 fstr = formula (fun);
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	93 if (isplot)
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	94 xarg = argids{1};
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	95 if (nargs == 2)
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	96 yarg = argids{2};
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	97 else
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	98 yarg = "";
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	99 endif
c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	100 elseif (isplot3)
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	101 xarg = "x";
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	102 yarg = "y";
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	103 elseif (isplot \|\| ispolar)
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	104 xarg = "";
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	105 yarg = "";
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	106 else
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	107 xarg = argids{1};
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	108 yarg = argids{2};
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	109 endif
25806 23483673ba43 Use is_function_handle instead of isa (x, "function_handle"). Rik <rik@octave.org> parents: 25054 diff changeset	110 elseif (is_function_handle (fun))
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	111 fstr = func2str (fun);
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	112 idx = index (fstr, ')');
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	113 if (idx != 0)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	114 args = regexp (fstr(3:(idx-1)), '\w+', 'match');
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	115 fstr = fstr(idx+2:end); # remove '@(x) ' from string name
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	116 else
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	117 args = {"x"};
20499 c0566df8cde3 Generate correct ezplot for 2-input functions (bug #46004). Rik <rik@octave.org> parents: 19697 diff changeset	118 try
c0566df8cde3 Generate correct ezplot for 2-input functions (bug #46004). Rik <rik@octave.org> parents: 19697 diff changeset	119 if (builtin ("nargin", fun) == 2)
c0566df8cde3 Generate correct ezplot for 2-input functions (bug #46004). Rik <rik@octave.org> parents: 19697 diff changeset	120 args{2} = "y";
c0566df8cde3 Generate correct ezplot for 2-input functions (bug #46004). Rik <rik@octave.org> parents: 19697 diff changeset	121 endif
c0566df8cde3 Generate correct ezplot for 2-input functions (bug #46004). Rik <rik@octave.org> parents: 19697 diff changeset	122 end_try_catch
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	123 endif
c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	124 if (isplot && length (args) == 2)
c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	125 nargs = 2;
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	126 elseif (numel (args) != nargs)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	127 error ("%s: expecting a function of %d arguments", ezfunc, nargs);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	128 endif
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	129 if (isplot)
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	130 xarg = args{1};
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	131 if (nargs == 2)
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	132 yarg = args{2};
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	133 else
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	134 yarg = "";
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	135 endif
c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	136 elseif (isplot3)
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	137 xarg = "x";
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	138 yarg = "y";
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	139 elseif (ispolar)
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	140 xarg = "";
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	141 yarg = "";
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	142 else
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	143 xarg = args{1};
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	144 yarg = args{2};
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	145 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	146 else
20711 7b608fadc663 Make error messages more specific about the variable and problem encountered. Rik <rik@octave.org> parents: 20499 diff changeset	147 error ("%s: F must be string, inline function, or function handle", ezfunc);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	148 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	149
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	150 if (nargin > 2 \|\| (nargin == 2 && isplot))
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	151 funx = fun;
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	152 fstrx = fstr;
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	153 funy = varargin{2};
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	154 if (ischar (funy) && ! strcmp (funy, "circ") && ! strcmp (funy, "animate"))
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	155 parametric = true;
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	156 if (exist (funy, "file") \|\| exist (funy, "builtin"))
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	157 funy = inline ([funy "(t)"]);
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	158 else
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	159 funy = vectorize (inline (funy));
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	160 endif
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	161 if (numel (argnames (funy)) != nargs)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	162 error ("%s: expecting a function of %d arguments", ezfunc, nargs);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	163 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	164 fstry = formula (funy);
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	165 elseif (strcmp (typeinfo (funy), "inline function"))
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	166 parametric = true;
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	167 if (numel (argnames (funy)) != nargs)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	168 error ("%s: expecting a function of %d arguments", ezfunc, nargs);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	169 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	170 funy = vectorize (funy);
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	171 fstry = formula (funy);
25806 23483673ba43 Use is_function_handle instead of isa (x, "function_handle"). Rik <rik@octave.org> parents: 25054 diff changeset	172 elseif (is_function_handle (funy))
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	173 parametric = true;
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	174 fstry = func2str (funy);
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	175 idx = index (fstry, ')');
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	176 if (idx != 0)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	177 args = regexp (fstry(3:(idx-1)), '\w+', 'match');
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	178 fstry = fstry(idx+2:end); # remove '@(x) ' from string name
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	179 else
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	180 args = {"y"};
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	181 endif
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	182 if (numel (args) != nargs)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	183 error ("%s: expecting a function of %d arguments", ezfunc, nargs);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	184 endif
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	185 endif
c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	186
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	187 if (! parametric && isplot3)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	188 needusage = true; # Can't call non-parametric ezplot3
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	189 return;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	190 elseif (parametric && isplot)
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	191 if (nargs == 2)
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	192 error ("%s: can not define a parametric function in this manner", ezfunc);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	193 else
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	194 xarg = "x";
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	195 yarg = "y";
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	196 endif
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	197 elseif (parametric)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	198 funz = varargin{3};
11149 fe3c3dfc07eb style fix: break lines before && and \|\|, not after John W. Eaton <jwe@octave.org> parents: 11032 diff changeset	199 if (ischar (funz) && ! strcmp (funz, "circ")
fe3c3dfc07eb style fix: break lines before && and \|\|, not after John W. Eaton <jwe@octave.org> parents: 11032 diff changeset	200 && ! strcmp (funz, "animate"))
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	201 if (exist (funz, "file") \|\| exist (funz, "builtin"))
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	202 funz = inline ([funz "(t)"]);
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	203 else
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	204 funz = vectorize (inline (funz));
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	205 endif
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	206 if (numel (argnames (funz)) > nargs)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	207 error ("%s: expecting a function of %d arguments", ezfunc, nargs);
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	208 endif
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	209 fstrz = formula (funz);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	210 elseif (strcmp (typeinfo (funz), "inline function"))
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	211 if (numel (argnames (funz)) != nargs)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	212 error ("%s: expecting a function of %d arguments", ezfunc, nargs);
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	213 endif
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	214 funz = vectorize (funz);
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	215 fstrz = formula (funz);
25806 23483673ba43 Use is_function_handle instead of isa (x, "function_handle"). Rik <rik@octave.org> parents: 25054 diff changeset	216 elseif (is_function_handle (funz))
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	217 fstrz = func2str (funz);
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	218 idx = index (fstrz, ')');
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	219 if (idx != 0)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	220 args = regexp (fstrz(3:(idx-1)), '\w+', 'match');
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	221 fstrz = fstrz(idx+2:end); # remove '@(x) ' from string name
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	222 else
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	223 args = {"z"};
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	224 endif
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	225 if (numel (args) != nargs)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	226 error ("%s: expecting a function of %d arguments", ezfunc, nargs);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	227 endif
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	228 else
20711 7b608fadc663 Make error messages more specific about the variable and problem encountered. Rik <rik@octave.org> parents: 20499 diff changeset	229 error ("%s: parametric plots require 3 functions", ezfunc);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	230 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	231 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	232 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	233
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	234 if ((isplot && nargs != 2) \|\| isplot3 \|\| ispolar)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	235 n = 500; # default for point-style functions like plot
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	236 else
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	237 n = 60; # default for meshgrid style functions like contour, surf
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	238 endif
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	239 domain = [];
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	240 circ = false;
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	241 animate = false;
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	242 if (parametric)
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	243 if (isplot)
c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	244 iarg = 3;
c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	245 else
c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	246 iarg = 4;
c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	247 endif
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	248 else
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	249 iarg = 2;
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	250 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	251 while (iarg <= nargin)
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	252 arg = varargin{iarg++};
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	253 if (ischar (arg) && strcmp (arg, "circ"))
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	254 circ = true;
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	255 elseif (ischar (arg) && strcmp (arg, "animate"))
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	256 animate = true;
21560 96c768d898f0 Use stricter input validation for ezplot N, DOM inputs (bug #46952). Rik <rik@octave.org> parents: 21385 diff changeset	257 elseif (isscalar (arg) && (n == 60 \|\| n == 500))
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	258 n = arg;
21560 96c768d898f0 Use stricter input validation for ezplot N, DOM inputs (bug #46952). Rik <rik@octave.org> parents: 21385 diff changeset	259 elseif (numel (arg) == 2 && isempty (domain))
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	260 domain = [arg(1) arg(2) arg(1) arg(2)];
21560 96c768d898f0 Use stricter input validation for ezplot N, DOM inputs (bug #46952). Rik <rik@octave.org> parents: 21385 diff changeset	261 elseif (numel (arg) == 4 && isempty (domain))
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	262 domain = arg(:).';
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	263 else
21560 96c768d898f0 Use stricter input validation for ezplot N, DOM inputs (bug #46952). Rik <rik@octave.org> parents: 21385 diff changeset	264 error ("%s: expecting scalar N, or 2-/4-element vector DOM", ezfunc);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	265 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	266 endwhile
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	267
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	268 if (circ && (iscontour \|\| isplot3 \|\| isplot))
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	269 needusage = true;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	270 return;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	271 elseif (circ && parametric)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	272 error ("%s: can not have both circular domain and parametric function",
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	273 ezfunc);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	274 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	275
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	276 if (animate && ! isplot3)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	277 error ("%s: animate option only valid for ezplot3", ezfunc);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	278 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	279
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	280 if (parametric)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	281 ## Make the label strings pretty by removing extra spaces between base
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	282 ## and exponent, the '.' in vectorized code, and the '*' for multiply.
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	283 fstrx = regexprep (regexprep (regexprep (fstrx,
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	284 '\s\.?(?:\^\|\\)\s','^'), '\.([/+-])', '$1'), '\s\.?\\s*', ' ');
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	285 fstry = regexprep (regexprep (regexprep (fstry,
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	286 '\s\.?(?:\^\|\\)\s','^'), '\.([/+-])', '$1'), '\s\.?\\s*', ' ');
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	287 if (isplot)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	288 fstr = ["x = " fstrx ", y = " fstry];
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	289 else
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	290 fstrz = regexprep (regexprep (regexprep (fstrz,
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	291 '\s\.?(?:\^\|\\)\s','^'), '\.([/+-])', '$1'), '\s\.?\\s*', ' ');
22308 b756cea4fdea __ezplot__.m: Fix Matlab incompatibilities. Rik <rik@octave.org> parents: 21561 diff changeset	292 fstr = ["x = " fstrx ", y = " fstry ", z = " fstrz];
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	293 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	294 else
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	295 fstr = regexprep (regexprep (regexprep (fstr,
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	296 '\s\.?(?:\^\|\\)\s','^'), '\.([/+-])', '$1'), '\s\.?\\s*', ' ');
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	297 if (isplot && nargs == 2)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	298 fstr = [fstr " = 0"]; # make title string of implicit function
22309 46d39db0c571 __ezplot__.m: Add 'r = ' to the title string for ezpolar plots. Rik <rik@octave.org> parents: 22308 diff changeset	299 elseif (ispolar)
46d39db0c571 __ezplot__.m: Add 'r = ' to the title string for ezpolar plots. Rik <rik@octave.org> parents: 22308 diff changeset	300 fstr = ["r = " fstr];
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	301 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	302 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	303
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	304 if (isempty (domain))
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	305 auto_domain = true;
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	306 if (isplot3 \|\| ispolar)
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	307 domain = [0, 2pi, 0, 2pi];
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	308 else
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	309 domain = [-2pi, 2pi, -2pi, 2pi];
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	310 endif
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	311 else
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	312 auto_domain = false;
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	313 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	314
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	315 auto_domain_done = false;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	316 do
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	317 domain_ok = true;
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	318
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	319 if ((isplot && nargs == 1) \|\| isplot3 \|\| ispolar)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	320 X = linspace (domain(1), domain(2), n);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	321 elseif (isplot && numel (domain) == 2)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	322 x = linspace (domain(1), domain(2), n);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	323 [X, Y] = meshgrid (x, x);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	324 elseif (circ)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	325 ## To plot on circular domain develop grid in polar coordinates
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	326 ## and then switch these to Cartesian coordinates.
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	327 cent = [domain(1) + domain(2), domain(3) + domain(4)] / 2;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	328 rmax = sqrt ((domain(2) - cent(1))^2 + (domain(4) - cent(2))^2);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	329 r = linspace (0, rmax, n);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	330 t = linspace (0, 2*pi, n);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	331 [T, R] = meshgrid (t, r);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	332 X = R .* cos (T) + cent(1);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	333 Y = R .* sin (T) + cent(2);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	334 domain = [-rmax+cent(1), +rmax+cent(1), -rmax+cent(2), +rmax+cent(2)];
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	335 else # contour, mesh, surf plots
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	336 x = linspace (domain(1), domain(2), n);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	337 y = linspace (domain(3), domain(4), n);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	338 [X, Y] = meshgrid (x, y);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	339 endif
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	340
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	341 if (parametric)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	342 if (isplot)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	343 XX = feval (funx, X);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	344 Z = feval (funy, X);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	345 X = XX;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	346 elseif (isplot3)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	347 Z = feval (funz, X);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	348 XX = feval (funx, X);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	349 YY = feval (funy, X);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	350 X = XX;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	351 Y = YY;
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	352 else
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	353 Z = feval (funz, X, Y);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	354 XX = feval (funx, X, Y);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	355 YY = feval (funy, X, Y);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	356 X = XX;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	357 Y = YY;
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	358
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	359 ## Eliminate the singularities
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	360 X = __eliminate_sing__ (X);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	361 Y = __eliminate_sing__ (Y);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	362 Z = __eliminate_sing__ (Z);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	363 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	364 else ## non-parametric plots
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	365 if (isplot && nargs == 2)
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	366 Z = feval (fun, X, Y);
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	367
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	368 ## Matlab returns line objects for this case and so can't call
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	369 ## contour directly as it returns patch objects to allow colormaps
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	370 ## to work with contours. Therefore recreate the lines from the
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	371 ## output for contourc, and store in cell arrays.
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	372 [c, ~] = contourc (X, Y, Z, [0, 0]);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	373
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	374 i = 1;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	375 XX = YY = {};
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	376 while (i < length (c))
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	377 clev = c(1,i);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	378 clen = c(2,i);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	379 XX = [XX, {c(1, i+1:i+clen)}];
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	380 YY = [YY, {c(2, i+1:i+clen)}];
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	381 i += clen+1;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	382 endwhile
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	383 else
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	384 if (ispolar)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	385 Z = feval (fun, X);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	386 ## FIXME: Why aren't singularities eliminated for polar plots?
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	387 elseif (isplot)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	388 Z = feval (fun, X);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	389 ## Eliminate the singularities
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	390 Z = __eliminate_sing__ (Z);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	391 domain = find_valid_domain (X, [], Z);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	392 elseif (iscontour)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	393 Z = feval (fun, X, Y);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	394 Z = __eliminate_sing__ (Z);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	395 else # mesh, surf plots
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	396 Z = feval (fun, X, Y);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	397 Z = __eliminate_sing__ (Z);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	398 if (circ)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	399 ## Use domain calculated at the start.
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	400 ## The X, Y grids are non-monotonic after conversion from polar
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	401 ## coordinates and find_valid_domain fails.
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	402
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	403 elseif (auto_domain && ! auto_domain_done)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	404 valid_domain = find_valid_domain (X, Y, Z);
22882 9b096bffc10d Fix plotting of constant functions with ezplot (bug #49766). Rik <rik@octave.org> parents: 22767 diff changeset	405 domain_ok = all (domain == valid_domain);
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	406 domain = valid_domain;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	407 auto_domain_done = true; # ensures only 1 round of do loop done
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	408 else
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	409 if (! auto_domain_done)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	410 domain = find_valid_domain (X, Y, Z);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	411 endif
17306 09543e9c8f40 Use explicit form of end (endif, endfor, etc.) in core m-files. Rik <rik@octave.org> parents: 17301 diff changeset	412 endif
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	413 endif
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	414 endif
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	415 endif
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	416 until (domain_ok)
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	417
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	418 ## Now, actually call the correct plot function with valid data and domain.
17211 87ba70043bfc Don't use ifelse in plot fcns to avoid unnecessary fcn evaluations. Rik <rik@octave.org> parents: 17167 diff changeset	419 oldfig = [];
17301 68bcac3c043a Correct inversion accidentally introduced in cset 87ba70043bfc. Rik <rik@octave.org> parents: 17211 diff changeset	420 if (! isempty (hax))
17211 87ba70043bfc Don't use ifelse in plot fcns to avoid unnecessary fcn evaluations. Rik <rik@octave.org> parents: 17167 diff changeset	421 oldfig = get (0, "currentfigure");
87ba70043bfc Don't use ifelse in plot fcns to avoid unnecessary fcn evaluations. Rik <rik@octave.org> parents: 17167 diff changeset	422 endif
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	423 unwind_protect
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	424 hax = newplot (hax);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	425 if (iscontour)
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	426 [~, h] = feval (pltfunc, hax, X, Y, Z);
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	427 elseif (isplot && nargs == 2)
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	428 h = zeros (length (XX), 1);
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	429 hold_state = get (hax, "nextplot");
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	430 for i = 1 : length (XX)
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	431 if (i == 1)
22308 b756cea4fdea __ezplot__.m: Fix Matlab incompatibilities. Rik <rik@octave.org> parents: 21561 diff changeset	432 h(1) = plot (hax, XX{1}, YY{1});
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	433 set (hax, "nextplot", "add");
22308 b756cea4fdea __ezplot__.m: Fix Matlab incompatibilities. Rik <rik@octave.org> parents: 21561 diff changeset	434 color = get (h(1), "color");
b756cea4fdea __ezplot__.m: Fix Matlab incompatibilities. Rik <rik@octave.org> parents: 21561 diff changeset	435 else
b756cea4fdea __ezplot__.m: Fix Matlab incompatibilities. Rik <rik@octave.org> parents: 21561 diff changeset	436 h(i) = plot (hax, XX{i}, YY{i}, "color", color);
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 9896 diff changeset	437 endif
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	438 endfor
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	439 set (hax, "nextplot", hold_state);
27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	440 axis (hax, domain);
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	441 elseif (isplot \|\| ispolar)
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	442 h = feval (pltfunc, hax, X, Z);
18718 a50d0d85dc73 Set axis equal for 2D parametric plots done with ezplot. Rik <rik@octave.org> parents: 17744 diff changeset	443 if (isplot)
a50d0d85dc73 Set axis equal for 2D parametric plots done with ezplot. Rik <rik@octave.org> parents: 17744 diff changeset	444 if (! parametric)
a50d0d85dc73 Set axis equal for 2D parametric plots done with ezplot. Rik <rik@octave.org> parents: 17744 diff changeset	445 axis (hax, domain);
a50d0d85dc73 Set axis equal for 2D parametric plots done with ezplot. Rik <rik@octave.org> parents: 17744 diff changeset	446 else
a50d0d85dc73 Set axis equal for 2D parametric plots done with ezplot. Rik <rik@octave.org> parents: 17744 diff changeset	447 axis ("equal");
a50d0d85dc73 Set axis equal for 2D parametric plots done with ezplot. Rik <rik@octave.org> parents: 17744 diff changeset	448 endif
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	449 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	450 elseif (isplot3)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	451 if (animate)
26671 fec69a310833 __ezplot__.m: for ezplot3, draw "animate" plot just once. Rik <rik@octave.org> parents: 26379 diff changeset	452 comet3 (hax, X, Y, Z);
fec69a310833 __ezplot__.m: for ezplot3, draw "animate" plot just once. Rik <rik@octave.org> parents: 26379 diff changeset	453 else
fec69a310833 __ezplot__.m: for ezplot3, draw "animate" plot just once. Rik <rik@octave.org> parents: 26379 diff changeset	454 h = feval (pltfunc, hax, X, Y, Z);
8046 c4482fc30c7f Add the ezplot function David Bateman <dbateman@free.fr> parents: 7540 diff changeset	455 endif
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	456 grid (hax, "on");
27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	457 zlabel (hax, "z");
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	458 else # mesh and surf plots
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	459 h = feval (pltfunc, hax, X, Y, Z);
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	460 ## FIXME: surf, mesh should really do a better job of setting zlim
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	461 if (! parametric)
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	462 axis (hax, domain);
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	463 endif
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	464 endif
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	465 xlabel (hax, xarg);
27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	466 ylabel (hax, yarg);
27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	467 title (hax, fstr);
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	468 unwind_protect_cleanup
17054 27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	469 if (! isempty (oldfig))
27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	470 set (0, "currentfigure", oldfig);
27b3a675ea6b __ezplot__.m: Update to use new __plt_get_axis_arg__. Rik <rik@octave.org> parents: 16993 diff changeset	471 endif
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	472 end_unwind_protect
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	473
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	474 endfunction
e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	475
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	476 ## Eliminate bad data (complex values, infinities, singularities)
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	477 function x = __eliminate_sing__ (x)
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	478 if (iscomplex (x))
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	479 x(imag (x) != 0) = NaN;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	480 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	481 x(isinf (x)) = NaN;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	482 ## High rates of curvature are treated as singularities
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	483 threshold = 0.2 * (max (x(:)) - min (x(:)));
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	484 x(abs (del2 (x)) > threshold) = NaN;
7337 e70789e0cd92 [project @ 2007-12-28 21:04:44 by jwe] jwe parents: diff changeset	485 endfunction
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	486
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	487 ## Find: 1) range of function where there are not NaN values,
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	488 ## 2) function is changing (not just flat surface)
22767 212333a97d8d maint: Remove dangling ';' from m-file function declarations. Rik <rik@octave.org> parents: 22755 diff changeset	489 function domain = find_valid_domain (X, Y, Z)
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	490
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	491 if (isvector (Z))
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	492 ## 2-D data for isplot
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	493 domain = [X(1) X(end)];
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	494
19593 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	495 ## Guess a range which includes the "mass" of the data by using a
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	496 ## median-based approach. The center 3/4 of the data is used to
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	497 ## determine the range of the data.
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	498 ## This seems to be vaguely what Matlab does, but can't be sure.
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	499 XX = sort (Z(isfinite (Z)));
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	500 if (length (X) > 4)
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	501 irlo = XX(fix (1/8 * length (XX)));
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	502 irhi = XX(fix (7/8 * length (XX)));
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	503 d = irhi - irlo;
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	504 domain(3) = max (XX(1) - d/8, irlo - d);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	505 domain(4) = min (XX(end) + d/8, irhi + d);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	506 else
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	507 domain(3:4) = [XX(1), XX(end)];
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	508 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	509
22882 9b096bffc10d Fix plotting of constant functions with ezplot (bug #49766). Rik <rik@octave.org> parents: 22767 diff changeset	510 ## Handle exceptional case of constant function
22910 23847979b91e maint: Strip trailing whitespace from source files. John W. Eaton <jwe@octave.org> parents: 22882 diff changeset	511 if (domain(3) == domain(4))
22882 9b096bffc10d Fix plotting of constant functions with ezplot (bug #49766). Rik <rik@octave.org> parents: 22767 diff changeset	512 domain(3) -= 1;
9b096bffc10d Fix plotting of constant functions with ezplot (bug #49766). Rik <rik@octave.org> parents: 22767 diff changeset	513 domain(4) += 1;
9b096bffc10d Fix plotting of constant functions with ezplot (bug #49766). Rik <rik@octave.org> parents: 22767 diff changeset	514 endif
9b096bffc10d Fix plotting of constant functions with ezplot (bug #49766). Rik <rik@octave.org> parents: 22767 diff changeset	515
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	516 else
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	517 ## 3-D data such as mesh, surf
19593 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	518 Zfinite = ! isnan (Z);
446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	519 Zrows = any (Zfinite, 2);
446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	520 rmin = find (Zrows, 1, "first");
446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	521 rmax = find (Zrows, 1, "last");
446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	522 Zcols = any (Zfinite, 1);
446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	523 cmin = find (Zcols, 1, "first");
446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	524 cmax = find (Zcols, 1, "last");
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	525
19593 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	526 ## Handle nasty case of all NaNs
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	527 if (isempty (rmin))
21385 89fa0694aa2e Fix check-missing-semicolon QA target (bug #47277). Rik <rik@octave.org> parents: 20852 diff changeset	528 rmin = 1; rmax = rows (Z);
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	529 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	530 if (isempty (cmin))
21385 89fa0694aa2e Fix check-missing-semicolon QA target (bug #47277). Rik <rik@octave.org> parents: 20852 diff changeset	531 cmin = 1; cmax = columns (Z);
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	532 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	533
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	534 if ( ! any (isnan (Z([rmin, rmax],:)(:)))
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	535 && ! any (isnan (Z(:, [cmin, cmax])(:))))
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	536 ## Exclude surfaces along borders which are flat (gradient =~ 0).
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	537 ## Technically, this calculation might be better done with actual
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	538 ## deltaX, deltaY values. But, data is usually meshgridded
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	539 ## (constant spacing) so working with deltaROW#, deltaCOL# is fine.
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	540 [Zx, Zy] = gradient (Z(rmin:rmax, cmin:cmax));
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	541 Zgrad = sqrt (Zx.^2 + Zy.^2);
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	542 slope = ((max (Z(:)) - min (Z(:)))
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	543 / sqrt ((rmax - rmin)^2 + (cmax - cmin)^2));
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	544 slope /= 125; # threshold for discarding points.
19593 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	545 Zrows = any (Zgrad > slope, 2);
446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	546 rmin += find (Zrows, 1, "first") - 1;
446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	547 rmax += find (Zrows, 1, "last") - rows (Zrows);
446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	548 Zcols = any (Zgrad > slope, 1);
446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	549 cmin += find (Zcols, 1, "first") - 1;
446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	550 cmax += find (Zcols, 1, "last") - columns (Zcols);
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	551 endif
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	552
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	553 domain = [X(1,cmin) X(1,cmax) Y(rmin,1) Y(rmax,1)];
78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	554 endif
19593 446c46af4b42 strip trailing whitespace from most source files John W. Eaton <jwe@octave.org> parents: 17744 diff changeset	555
16993 78f57b14535c Overhaul ez* family of plot functions. Rik <rik@octave.org> parents: 14868 diff changeset	556 endfunction

Mercurial > jwe > octave

annotate scripts/plot/draw/private/__ezplot__.m @ 27924:1891570abac8

annotate scripts/plot/draw/private/ezplot.m @ 27924:1891570abac8