From: <prn...@us...> - 2010-11-28 20:16:20
|
Revision: 7959 http://octave.svn.sourceforge.net/octave/?rev=7959&view=rev Author: prnienhuis Date: 2010-11-28 20:16:13 +0000 (Sun, 28 Nov 2010) Log Message: ----------- Will be moved tp fpl package Removed Paths: ------------- trunk/octave-forge/main/io/inst/savevtkvector.m Deleted: trunk/octave-forge/main/io/inst/savevtkvector.m =================================================================== --- trunk/octave-forge/main/io/inst/savevtkvector.m 2010-11-28 20:15:22 UTC (rev 7958) +++ trunk/octave-forge/main/io/inst/savevtkvector.m 2010-11-28 20:16:13 UTC (rev 7959) @@ -1,64 +0,0 @@ -## Copyright (C) 2010 Kurnia Wano, Levente Torok -##. -## 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. -##. -## You should have received a copy of the GNU General Public License -## along with Octave; see the file COPYING. If not, see -## - -## -*- texinfo -*- -## -## @deftypefn {Function File} savevtkvector ( @var{X}, @var{Y}, @var{Z}, @var{filename} ) -## savevtkvector Save a 3-D vector array in VTK format -## savevtkvector(X,Y,Z,filename) saves a 3-D vector of any size to -## filename in VTK format. X, Y and Z should be arrays of the same -## size, each storing speeds in the a single Cartesian directions. -## -## @end deftypefn -## -## -## Author: Kurnia Wano, Levente Torok <Tor...@gm...> -## Created: 2010-08-02 -## Updates: 2010-11-03 -## - -function savevtkvector(X, Y, Z, filename) - - if ((size(X) ~= size(Y)) | (size(X) ~= size(Z))) - error('Error: velocity arrays of unequal size\n'); - return; - end - [ny,nx,nz]=size(datax); - xx=1:size(datax,2); - yy=1:size(datax,1); - zz=1:size(datax,3); - datax=datax(:)’; - datay=datay(:)’; - dataz=dataz(:)’; - %% Header - fid=fopen(fullfile(outputd,filename),’w'); - fprintf(fid,’%s\n’,'# vtk DataFile Version 3.0′); - fprintf(fid,’%s\n’,’3D LFF extrapolation’); - fprintf(fid,’%s\n’,'ASCII’); - fprintf(fid,’%s\n’,'DATASET RECTILINEAR_GRID’); - fprintf(fid,’%s %1.0i %1.0i %1.0i\n’,'DIMENSIONS’,nx,ny,nz); - fprintf(fid,’%s %1.0i %s\n’,'X_COORDINATES’,nx,’float’); - fprintf(fid,’%1.0i ‘,xx); - fprintf(fid,’\n%s %1.0i %s\n’,'Y_COORDINATES’,ny,’float’); - fprintf(fid,’%1.0i ‘,yy); - fprintf(fid,’\n%s %1.0i %s\n’,'Z_COORDINATES’,nz,’float’); - fprintf(fid,’%1.0i ‘,zz); - %% Data - fprintf(fid,’\n%s %1.0i’,'POINT_DATA’,nx*ny*nz); - fprintf(fid,’\n%s\n’,'VECTORS BFIELD float’); - fprintf(fid,’%6.2f %6.2f %6.2f\n’,[datax;datay;dataz]); - fclose(fid); -return This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site. |