Mercurial > octave
view scripts/plot/draw/tetramesh.m @ 27978:a4268efb7334
maint: Eliminate single "Author:" lines from code base.
* color-picker.cc, color-picker.h, octave-qscintilla.cc, octave-qscintilla.h,
octave-txt-lexer.cc, octave-txt-lexer.h, balance.cc, kron.cc, lookup.cc,
ls-hdf5.cc, ls-mat5.cc, ls-oct-text.cc, sylvester.cc, tsearch.cc, jit-ir.cc,
jit-ir.h, jit-typeinfo.cc, jit-typeinfo.h, jit-util.cc, jit-util.h, pt-jit.cc,
pt-jit.h, Array-jit.cc, bsxfun.h, url-transfer.cc, url-transfer.h, acosd.m,
acotd.m, acscd.m, asecd.m, asind.m, atand.m, cosd.m, cotd.m, cscd.m, secd.m,
sind.m, tand.m, cart2pol.m, cart2sph.m, curl.m, del2.m, divergence.m, flip.m,
pol2cart.m, randi.m, sortrows.m, sph2cart.m, griddata3.m, griddatan.m,
rectint.m, voronoi.m, dialog.m, guidata.m, guihandles.m, __file_filter__.m,
__fltk_file_filter__.m, __get_funcname__.m, __is_function__.m,
__uigetdir_fltk__.m, __uigetfile_fltk__.m, __uiobject_split_args__.m,
__uiputfile_fltk__.m, uibuttongroup.m, uicontextmenu.m, uicontrol.m,
uigetdir.m, uigetfile.m, uimenu.m, uipanel.m, uipushtool.m, uiputfile.m,
uiresume.m, uitoggletool.m, uitoolbar.m, uiwait.m, doc.m, autumn.m, bone.m,
cmpermute.m, cmunique.m, cool.m, copper.m, cubehelix.m, flag.m, frame2im.m,
hot.m, hsv.m, hsv2rgb.m, im2frame.m, imfinfo.m, imformats.m, iscolormap.m,
jet.m, pink.m, prism.m, __imfinfo__.m, rainbow.m, rgb2hsv.m, spinmap.m,
spring.m, viridis.m, white.m, winter.m, importdata.m, strmatch.m, bandwidth.m,
isbanded.m, isdiag.m, istril.m, istriu.m, krylov.m, linsolve.m, logm.m,
lscov.m, subspace.m, bunzip2.m, compare_versions.m, edit.m, fullfile.m,
getfield.m, gunzip.m, inputParser.m, license.m, nargchk.m, narginchk.m,
nthargout.m, python.m, setfield.m, tar.m, unpack.m, ver.m, zip.m, odeplot.m,
fminbnd.m, fminunc.m, fsolve.m, fzero.m, glpk.m, humps.m, pathdef.m,
savepath.m, diffuse.m, lighting.m, material.m, shading.m, specular.m,
contour.m, ellipsoid.m, isocaps.m, isocolors.m, isonormals.m, isosurface.m,
light.m, pcolor.m, __interp_cube__.m, __marching_cube__.m, __patch__.m,
__stem__.m, reducepatch.m, reducevolume.m, ribbon.m, shrinkfaces.m, slice.m,
stem.m, surf.m, surfl.m, tetramesh.m, waterfall.m, __actual_axis_position__.m,
__next_line_color__.m, allchild.m, figure.m, findall.m, findobj.m, gcf.m,
hggroup.m, isprop.m, ndgrid.m, __add_default_menu__.m, saveas.m,
__splinefit__.m, prefdir.m, preferences.m, profexplore.m, profexport.m,
profile.m, profshow.m, setdiff.m, freqz.m, freqz_plot.m, sinc.m, unwrap.m,
pcr.m, qmr.m, sprandn.m, betaincinv.m, cosint.m, ellipke.m, factor.m,
gammaincinv.m, legendre.m, nchoosek.m, primes.m, invhilb.m, toeplitz.m,
erase.m, mat2str.m, strtrim.m, validatestring.m, fail.m, compare_plot_demos.m,
dump_demos.m, html_compare_plot_demos.m, addtodate.m, datenum.m, jit.tst:
Eliminate single "Author:" lines from code base.
author | Rik <rik@octave.org> |
---|---|
date | Tue, 21 Jan 2020 12:53:13 -0800 |
parents | bd51beb6205e |
children | 9342688e86b4 |
line wrap: on
line source
######################################################################## ## ## Copyright (C) 2012-2020 The Octave Project Developers ## ## See the file COPYRIGHT.md in the top-level directory of this ## distribution or <https://octave.org/copyright/>. ## ## 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 3 of the License, 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, see ## <https://www.gnu.org/licenses/>. ## ######################################################################## ## -*- texinfo -*- ## @deftypefn {} {} tetramesh (@var{T}, @var{X}) ## @deftypefnx {} {} tetramesh (@var{T}, @var{X}, @var{C}) ## @deftypefnx {} {} tetramesh (@dots{}, @var{property}, @var{val}, @dots{}) ## @deftypefnx {} {@var{h} =} tetramesh (@dots{}) ## Display the tetrahedrons defined in the m-by-4 matrix @var{T} as 3-D ## patches. ## ## @var{T} is typically the output of a Delaunay triangulation of a 3-D set ## of points. Every row of @var{T} contains four indices into the n-by-3 ## matrix @var{X} of the vertices of a tetrahedron. Every row in @var{X} ## represents one point in 3-D space. ## ## The vector @var{C} specifies the color of each tetrahedron as an index ## into the current colormap. The default value is 1:m where m is the number ## of tetrahedrons; the indices are scaled to map to the full range of the ## colormap. If there are more tetrahedrons than colors in the colormap then ## the values in @var{C} are cyclically repeated. ## ## Calling @code{tetramesh (@dots{}, "property", "value", @dots{})} passes all ## property/value pairs directly to the patch function as additional arguments. ## ## The optional return value @var{h} is a vector of patch handles where each ## handle represents one tetrahedron in the order given by @var{T}. ## A typical use case for @var{h} is to turn the respective patch ## @qcode{"visible"} property @qcode{"on"} or @qcode{"off"}. ## ## Type @code{demo tetramesh} to see examples on using @code{tetramesh}. ## @seealso{trimesh, delaunay, delaunayn, patch} ## @end deftypefn function h = tetramesh (varargin) [reg, prop] = parseparams (varargin); if (length (reg) < 2 || length (reg) > 3) print_usage (); endif T = reg{1}; X = reg{2}; if (! ismatrix (T) || columns (T) != 4) error ("tetramesh: T must be an n-by-4 matrix"); elseif (! ismatrix (X) || columns (X) != 3) error ("tetramesh: X must be an n-by-3 matrix"); endif size_T = rows (T); cmap = colormap (); if (length (reg) < 3) size_cmap = rows (cmap); C = mod ((1:size_T)' - 1, size_cmap) + 1; if (size_T < size_cmap && size_T > 1) ## expand to the available range of colors C = floor ((C - 1) * (size_cmap - 1) / (size_T - 1)) + 1; endif else C = reg{3}; if (! isvector (C) || size_T != length (C)) error ("tetramesh: C must be a vector of the same length as T"); endif endif hax = newplot (); hvec = zeros (size_T, 1); if (strcmp (graphics_toolkit (), "gnuplot")) ## Tiny reduction of the tetrahedron size to help gnuplot by ## avoiding identical faces with different colors for i = 1:size_T [th, p] = __shrink__ ([1 2 3 4], X(T(i, :), :), 1 - 1e-7); hvec(i) = patch ("Faces", th, "Vertices", p, "FaceColor", cmap(C(i), :), "FaceAlpha", 0.9, prop{:}); endfor else ## FLTK does not support FaceAlpha. for i = 1:size_T th = [1 2 3; 2 3 4; 3 4 1; 4 1 2]; hvec(i) = patch ("Faces", th, "Vertices", X(T(i, :), :), "FaceColor", cmap(C(i), :), "FaceAlpha", 1.0, prop{:}); endfor endif if (! ishold ()) set (hax, "view", [-37.5, 30]); endif if (nargout > 0) h = hvec; endif endfunction ## shrink the tetrahedron relative to its center of gravity function [tri, p] = __shrink__ (T, X, sf) midpoint = repmat (sum (X(T, :), 1) / 4, 12, 1); p = [X([1 2 3], :); X([2 3 4], :); X([3 4 1], :); X([4 1 2], :)]; p = sf * (p - midpoint) + midpoint; tri = reshape (1:12, 3, 4)'; endfunction %!demo %! clf; %! d = [-1 1]; %! [x,y,z] = meshgrid (d, d, d); %! x = [x(:); 0]; %! y = [y(:); 0]; %! z = [z(:); 0]; %! tetra = delaunay (x, y, z); %! X = [x(:) y(:) z(:)]; %! colormap (jet (64)); %! h = tetramesh (tetra, X); %! set (h(1:2:end), "visible", "off"); %! axis equal; %! view (30, 20); %! title ({"tetramesh() plot", ... %! "colormap = jet (64), every other tetrahedron invisible"}); %!demo %! clf; %! d = [-1 1]; %! [x,y,z] = meshgrid (d, d, d); %! x = [x(:); 0]; %! y = [y(:); 0]; %! z = [z(:); 0]; %! tetra = delaunay (x, y, z); %! X = [x(:) y(:) z(:)]; %! colormap (gray (256)); %! tetramesh (tetra, X, 21:20:241, "EdgeColor", "w"); %! axis equal; %! view (30, 20); %! title ({"tetramesh() plot", ... %! "colormap = gray (256) with white edges"});