Homec4science
Diffusion FENNECS lite (3Dmpimode) Browse

fennecslite/matlab/savegeomtoh5.m
2e0aa43721aa3Dmpimode

Last Change
Raw File

savegeomtoh5.m
View Options

function savegeomtoh5(filename,splines,dist_extent,overwrite)
% saves the geometry to be used with Fennecs
% splines is a structure array containing each individual boundary
% Overwrite defines if the .h5 file must be overwritten
% dist_extent is the distance in m on which the geometric weight goes from 0 to 1
% dist_extent can also be set in the job input file
% To look how to define each spline read the function
if nargin<4
overwrite=false;
end
group='/geometry_spl/';
try
fid=H5F.create(filename);
H5F.flush(fid,'H5F_SCOPE_LOCAL')
catch
if overwrite
delete(filename);
fid=H5F.create(filename);
H5F.flush(fid,'H5F_SCOPE_LOCAL')
else
error('cFile "%s" already exists\n',filename) ;
end
end
nbsplines=length(splines);
gid=H5G.create(fid,group,64);
h5writeatt(filename,group,'nbsplines',nbsplines); % total number of spline curves
h5writeatt(filename,group,'dist_extent',dist_extent);
if iscell(splines)
splines=cell2mat(splines);
end
for i=1:nbsplines
grp=sprintf('/geometry_spl/%2.2d',i);
gid=H5G.create(fid,grp,64);
h5writeatt(filename,grp,'Dirichlet_val',splines(i).Dval); % value of the electic potential in V at the boundary surface
h5writeatt(filename,grp,'order',splines(i).order); % order of the spline used to define the boundary
h5writeatt(filename,grp,'dim',splines(i).dim); % dimension of the spline curve
% epsge and epsce have no clear effect default values can be used
if isfield(splines(i),'epsge')
h5writeatt(filename,grp,'epsge',splines(i).epsge); % geometric presicion used by SISL for calculating distance to curve
else
h5writeatt(filename,grp,'epsge',1e-6); % geometric presicion used by SISL for calculating distance to curve
end
if isfield(splines(i),'epsce')
h5writeatt(filename,grp,'epsce',splines(i).epsce); % geometric presicion used by SISL for calculating distance to curve
else
h5writeatt(filename,grp,'epsce',1e-6); % geometric presicion used by SISL for calculating distance to curve
end
h5writeatt(filename,grp,'name',splines(i).name); % name of the boundary for readability of h5(not used in fennecs)
if isfield(splines(i),'type')
% boundary condition type (0 dirichlet fixed value)
% 1 Dirichlet with value depending on position
% 2 Natural
h5writeatt(filename,grp,'type',splines(i).type);
else
h5writeatt(filename,grp,'type',0);
end
if isfield(splines(i),'periodic')
% periodicity of the curve 1 closed periodic 0 open
h5writeatt(filename,grp,'periodic',splines(i).periodic);
else
h5writeatt(filename,grp,'periodic',0);
end
points=splines(i).points;
if size(splines(i).points,2)>splines(i).dim
points=points';
end
h5create(filename,[grp,'/pos'],size(points));
h5write(filename,[grp,'/pos'],points); % control points defining the curve
end
H5F.close(fid);
end
c4science · Help