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

matlab: field interpolation on given lines

pull/1/head
Thorsten Liebig 2010-10-12 16:32:47 +02:00
parent 5f033659ea
commit 58d8fc1ec5
1 changed files with 12 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
matlab/GetField_Interpolation.m
View File

@ -1,5 +1,5 @@
function [field_i mesh_i] = GetField_Interpolation(field, mesh, numLines, varargin) function [field_i mesh_i] = GetField_Interpolation(field, mesh, lines, varargin)
% [field_i mesh_i] = GetField_Interpolation(field, mesh, numLines, varargin) % [field_i mesh_i] = GetField_Interpolation(field, mesh, lines, varargin)
% %
% example: % example:
% [field mesh] = ReadHDF5Dump('Et.h5'); % [field mesh] = ReadHDF5Dump('Et.h5');
@ -15,7 +15,7 @@ function [field_i mesh_i] = GetField_Interpolation(field, mesh, numLines, vararg
% %
% See also ReadHDF5Dump ReadHDF5FieldData ReadHDF5Mesh % See also ReadHDF5Dump ReadHDF5FieldData ReadHDF5Mesh
if (~isnumeric(numLines) || numel(numLines)~=3) if ((~iscell(lines) && ~isnumeric(lines)) || numel(lines)~=3)
error('openEMS:GetField_Interpolation: numLines for interpolation must be a vector...'); error('openEMS:GetField_Interpolation: numLines for interpolation must be a vector...');
end end
@ -23,9 +23,15 @@ x = mesh.lines{1};
y = mesh.lines{2}; y = mesh.lines{2};
z = mesh.lines{3}; z = mesh.lines{3};
x_i = linspace(x(1),x(end),numLines(1)); if (isnumeric(lines))
y_i = linspace(y(1),y(end),numLines(2)); x_i = linspace(x(1),x(end),lines(1));
z_i = linspace(z(1),z(end),numLines(3)); y_i = linspace(y(1),y(end),lines(2));
z_i = linspace(z(1),z(end),lines(3));
else
x_i = lines{1};
y_i = lines{2};
z_i = lines{3};
end
field_i = field; field_i = field;
mesh_i = mesh; mesh_i = mesh;