savegeomtoh5.m
46f197fa7c5fgeom_testmumps

savegeomtoh5.m
View Options

			function savegeomtoh5(filename,splines,dist_extent,overwrite)
			if nargin<4
			overwrite=false;
			end
			if nargin<4
			epsge=1e-6;
			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);
			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);
			h5writeatt(filename,grp,'order',splines(i).order);
			h5writeatt(filename,grp,'dim',splines(i).dim);
			h5writeatt(filename,grp,'epsge',splines(i).epsge);
			h5writeatt(filename,grp,'epsce',splines(i).epsce);
			h5writeatt(filename,grp,'name',splines(i).name);
			h5create(filename,[grp,'/pos'],size(splines(i).points));
			h5write(filename,[grp,'/pos'],splines(i).points);
			end

			H5F.close(fid);
			end