From 12904eaa99bad389f62410e4cb1394b6636f1994 Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Wed, 11 Jan 2012 13:26:36 -0500
Subject: map2dat.m: added to main matlab lib and everything else linked to it

---
 matlab_helper_files/map2dat.m | 33 +++++++++++++++++++++++++++++++++
 1 file changed, 33 insertions(+)
 create mode 100644 matlab_helper_files/map2dat.m

(limited to 'matlab_helper_files/map2dat.m')

diff --git a/matlab_helper_files/map2dat.m b/matlab_helper_files/map2dat.m
new file mode 100644
index 0000000..4f4ca52
--- /dev/null
+++ b/matlab_helper_files/map2dat.m
@@ -0,0 +1,33 @@
+function map2dat(outfile, x,y,z, xskip, yskip)
+% saves 3D data in suitable way to be drawn by gnuplot
+% x,y - vectors of x,y values
+% z map of z values as used by Octave/Matlab
+% xskip, yskip - skip paprameters
+%  only every, xskip, yskip point will be written
+
+
+
+Nx=length(x);   
+Ny=length(y);   
+Nxs=ceil(Nx/xskip);
+Nys=ceil(Ny/yskip);
+points=zeros(1,3*Nxs*Nys);
+%points=[];
+tic; 
+for i=1:Nxs
+	for k=1:Nys
+		%points=[points x(i*xskip) y(k*yskip) z(k*yskip,i*xskip)];
+		points((i-1)*(Nys)*3+3*(k-1)+1) =   x(1+(i-1)*xskip);
+		points((i-1)*(Nys)*3+3*(k-1)+2) = y(1+(k-1)*yskip);
+		points((i-1)*(Nys)*3+3*(k-1)+3) = z(1+(k-1)*yskip,1+(i-1)*xskip);
+	end
+end
+disp('=== points formation complete ===');
+toc
+tic;
+%points
+fd = fopen(outfile, "wt");
+fprintf (fd, "%g %g %g\n", points);
+fclose(fd);
+disp('=== points saving complete ===');
+toc;
-- 
cgit v1.2.3