zcy
/
openEMS
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

matlab: Dump2VTK using scientific float print

pull/1/head
Thorsten Liebig 2010-10-14 12:20:19 +02:00
parent 27d88334b5
commit 1dacbeeb1a
1 changed files with 8 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
matlab/Dump2VTK.m
View File

@ -27,23 +27,23 @@ if (mesh.type==0) %write cartesian mesh to vtk
fprintf(fid,'DIMENSIONS %d %d %d\n',numel(x),numel(y),numel(z)); fprintf(fid,'DIMENSIONS %d %d %d\n',numel(x),numel(y),numel(z));
fprintf(fid,'X_COORDINATES %d float\n',numel(x)); fprintf(fid,'X_COORDINATES %d float\n',numel(x));
fprintf(fid,'%f',x(1)); fprintf(fid,'%e',x(1));
for n=2:numel(x) for n=2:numel(x)
fprintf(fid,' %f',x(n)); fprintf(fid,' %e',x(n));
end end
fprintf(fid,'\n'); fprintf(fid,'\n');
fprintf(fid,'Y_COORDINATES %d float\n',numel(y)); fprintf(fid,'Y_COORDINATES %d float\n',numel(y));
fprintf(fid,'%f',y(1)); fprintf(fid,'%e',y(1));
for n=2:numel(y) for n=2:numel(y)
fprintf(fid,' %f',y(n)); fprintf(fid,' %e',y(n));
end end
fprintf(fid,'\n'); fprintf(fid,'\n');
fprintf(fid,'Z_COORDINATES %d float\n',numel(z)); fprintf(fid,'Z_COORDINATES %d float\n',numel(z));
fprintf(fid,'%f',z(1)); fprintf(fid,'%e',z(1));
for n=2:numel(z) for n=2:numel(z)
fprintf(fid,' %f',z(n)); fprintf(fid,' %e',z(n));
end end
elseif (mesh.type==1) %write cylindrical mesh to vtk elseif (mesh.type==1) %write cylindrical mesh to vtk
@ -59,7 +59,7 @@ elseif (mesh.type==1) %write cylindrical mesh to vtk
for nz=1:numel(z) for nz=1:numel(z)
for ny=1:numel(y) for ny=1:numel(y)
for nx=1:numel(x) for nx=1:numel(x)
fprintf(fid,'%f %f %f\n',x(nx)*cos(y(ny)),x(nx)*sin(y(ny)),z(nz)); fprintf(fid,'%e %e %e\n',x(nx)*cos(y(ny)),x(nx)*sin(y(ny)),z(nz));
end end
end end
end end
@ -86,7 +86,7 @@ end
for nz=1:numel(z) for nz=1:numel(z)
for ny=1:numel(y) for ny=1:numel(y)
for nx=1:numel(x) for nx=1:numel(x)
fprintf(fid,'%f %f %f\n',fields(nx,ny,nz,1),fields(nx,ny,nz,2),fields(nx,ny,nz,3)); fprintf(fid,'%e %e %e\n',fields(nx,ny,nz,1),fields(nx,ny,nz,2),fields(nx,ny,nz,3));
end end
end end
end end