/usr/share/plplot_octave/support/__pl_mesh.m is in octave-plplot 5.10.0+dfsg-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 | ## Copyright (C) 1998-2003 Joao Cardoso.
##
## This program 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 of the License, or (at your
## option) any later version.
##
## This program 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.
##
## This file is part of plplot_octave.
function __pl_mesh(caller, varargin)
global __pl
global DRAW_LINEXY MAG_COLOR BASE_CONT
global TOP_CONT SURF_CONT DRAW_SIDES FACETED MESH
strm = __pl_init;
switch (caller)
case ("mesh")
__pl.type(strm) = 100;
case ("meshc")
__pl.type(strm) = 101;
case ("meshz")
__pl.type(strm) = 102;
case ("surf")
__pl.type(strm) = 103;
case ("surfc")
__pl.type(strm) = 104;
case ("surfl")
__pl.type(strm) = 105;
otherwise
error("__pl_mesh: FIXME")
endswitch
if (nargin == 2)
z = varargin{1};
if (ismatrix (z))
[rz, cz] = size(z);
x = (1:cz)';
y = (1:rz)';
else
error ("mesh: argument must be a matrix.\n");
endif
elseif (nargin == 4 )
x = varargin{1};
y = varargin{2};
z = varargin{3};
if (isvector (x) && isvector (y) && ismatrix (z))
xlen = length (x);
ylen = length (y);
if (xlen == columns (z) && ylen == rows (z))
if (rows (y) == 1)
y = y';
endif
if (rows (x) == 1)
x = x';
endif
else
msg = "mesh: rows (z) must be the same as length (x) and";
msg = sprintf ("%s\ncolumns (z) must be the same as length (y).\n", msg);
error (msg);
endif
elseif (ismatrix (x) && ismatrix (y) && ismatrix (z))
error("x,y and z all matrices not yet implemented.\n")
else
error ("mesh: x and y must be vectors and z must be a matrix.\n");
endif
else
help (caller);
return
endif
if (!ishold)
__pl.items(strm) = 1;
__pl.lxm(strm) = __pl.lym(strm) = __pl.lzm(strm) = realmax;
__pl.lxM(strm) = __pl.lyM(strm) = __pl.lzM(strm) = -realmax;
endif
items = __pl.items(strm);
## find the max/min x/y values. Currently reset at ??
__pl.lxm(strm) = min([__pl.lxm(strm), min(min(x))]);
__pl.lxM(strm) = max([__pl.lxM(strm), max(max(x))]);
__pl.lym(strm) = min([__pl.lym(strm), min(min(y))]);
__pl.lyM(strm) = max([__pl.lyM(strm), max(max(y))]);
__pl.lzm(strm) = min([__pl.lzm(strm), min(min(z))]);
__pl.lzM(strm) = max([__pl.lzM(strm), max(max(z))]);
## kludge, use "fmt" as plot type. But __pl_plotit still uses __pl.type
__pl.x{items, strm} = x;
__pl.y{items, strm} = y;
__pl.z{items, strm} = z;
__pl.fmt{items, strm} = __pl.type(strm);
__pl.items(strm) = __pl.items(strm) + 1;
__pl_meshplotit;
endfunction
|